/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N E V A R I E T E D E C A L A B I - Y A U 8 - D I M E N S I O N N E L L E : */
/* */
/* */
/* Author of '$xrs/CalabiYau.42$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20230210121926). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* I N T E R F A C E ' listG ' : */
/* */
/* */
/* :Debut_listG: */
/* :Fin_listG: */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D I R E C T I V E S S P E C I F I Q U E S D E C O M P I L A T I O N : */
/* */
/*************************************************************************************************************************************/
@define PRAGMA_CPP_____MODULE_NON_DERIVABLE
@define PRAGMA_CL_____MODULE_NON_OPTIMISABLE
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* F I C H I E R S D ' I N C L U D E S : */
/* */
/*************************************************************************************************************************************/
#include INCLUDES_BASE
#include image_image_IMAGESF_EXT
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N S D E B A S E E T U N I V E R S E L L E S : */
/* */
/*************************************************************************************************************************************/
#include xrv/dimensionnement.01.vv.I"
#include xrk/attractor.11.I"
#include xrv/particule.41.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 3 */
/* D E F I N I T I O N D E L ' E S P A C E P H Y S I Q U E D A N S R ( D E B U T ) : */
/* */
/* */
/* Nota : */
/* */
/* Les extrema des coordonnees {x,y,z} */
/* ainsi que ceux de leurs differentielles */
/* {dx,dy,dz} sont fixees un peu arbitrairement */
/* et sans etre parametrees. */
/* */
/* */
/*************************************************************************************************************************************/
#define hXmin_ESPACE \
PARE(-1.0)
#define hYmin_ESPACE \
PARE(-1.0)
#define hZmin_ESPACE \
PARE(-1.0)
/* Definition du "coin" inferieur-gauche-arriere de l'espace physique. */
#define hXmax_ESPACE \
PARE(+1.0)
#define hYmax_ESPACE \
PARE(+1.0)
#define hZmax_ESPACE \
PARE(+1.0)
/* Definition du "coin" superieur-droite-avant de l'espace physique. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 3 */
/* D E F I N I T I O N D E L ' E S P A C E P H Y S I Q U E D A N S R ( D E B U T ) : */
/* */
/*************************************************************************************************************************************/
#include xrk/attractor.12.I"
#define dXmin_ESPACE \
PARE(-1.00)
#define dYmin_ESPACE \
PARE(-1.00)
#define dZmin_ESPACE \
PARE(-1.00)
/* Definition des minima des differentielles {dx,dy,dz}. */
#define dXmax_ESPACE \
PARE(+1.00)
#define dYmax_ESPACE \
PARE(+1.00)
#define dZmax_ESPACE \
PARE(+1.00)
/* Definition des maxima des differentielles {dx,dy,dz}. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N E V A R I E T E D E C A L A B I - Y A U 8 - D I M E N S I O N N E L L E : */
/* */
/* */
/* Definition de l'hyper-volume ('v $xrs/CalabiYau.42$K') : */
/* */
/* Une variete hyper-complexe de Calabi-Yau */
/* est definie implicitement par : */
/* */
/* n1 n2 */
/* Q + Q = 1 */
/* 1 2 */
/* */
/* ou 'Q1' et 'Q2' sont deux nombres hyper-complexes. */
/* Cette variete est donc a 8 dimensions, et la */
/* surface a 3 dimensions que l'on visualisera */
/* sera definie par : */
/* */
/* */
/* X = (A .R(Q )) + (A .I(Q )) + (A .J(Q )) + (A .K(Q )) + (A .R(Q )) + (A .I(Q )) + (A .J(Q )) + (A .K(Q )) */
/* X1 1 X2 1 X3 1 X4 1 X5 2 X6 2 X7 2 X8 2 */
/* */
/* */
/* Y = (A .R(Q )) + (A .I(Q )) + (A .J(Q )) + (A .K(Q )) + (A .R(Q )) + (A .I(Q )) + (A .J(Q )) + (A .K(Q )) */
/* Y1 1 Y2 1 Y3 1 Y4 1 Y5 2 Y6 2 Y7 2 Y8 2 */
/* */
/* */
/* Z = (A .R(Q )) + (A .I(Q )) + (A .J(Q )) + (A .K(Q )) + (A .R(Q )) + (A .I(Q )) + (A .J(Q )) + (A .K(Q )) */
/* Z1 1 Z2 1 Z3 1 Z4 1 Z5 2 Z6 2 Z7 2 Z8 2 */
/* */
/* */
/* soit, par defaut (qui fait que 'v $xiirs/CAYA.A1.0129' et 'v $xiirs/CAYA.O3' sont semblables bien qu'ayant ete */
/* generees par '$xrs/CalabiYau.12$X' et '$xrs/CalabiYau.42$X' respectivement) : */
/* */
/* */
/* X = R(Q ) + I(Q ) */
/* 1 2 */
/* */
/* */
/* Y = I(Q ) + I(Q ) */
/* 1 2 */
/* */
/* */
/* Z = R(Q ) + I(Q ) */
/* 2 2 */
/* */
/* */
/* (ou 'R(...)', 'I(...)', J(...)' et 'K(...)' designent respectivement */
/* les parties Reelle, Imaginaire, Jmaginaire et Kmaginaire) avec : */
/* */
/* u E [ 0 , p/2 ] */
/* */
/* v E [ -1 , +1 ] */
/* */
/* w E [ -1 , +1 ] */
/* */
/* t E [ -1 , +1 ] */
/* */
/* (ou 'p' designe 'pi'). */
/* */
/* */
/*************************************************************************************************************************************/
#include xrs/CalabiYau.42.I"
#include xrs/CalabiYau.44.I"
#include xrs/CalabiYau.45.I"
#include xrs/CalabiYau.47.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* F O N C T I O N S D E V I S U A L I S A T I O N E T D ' I N T E R P O L A T I O N : */
/* */
/*************************************************************************************************************************************/
#include xrs/hyper_volumes.21.I"
#define PONDERATION_ROUGE_VERTE_BLEUE_POUR_dcx_dcy_dcz \
FU
#define PRENDRE_LA_VALEUR_ABSOLUE_DES_DIFFERENTIELLES \
VRAI
#include xrs/CalabiYau.46.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S T R O I S F O N C T I O N S ' F ' : */
/* */
/*************************************************************************************************************************************/
#define EDITER_LES_PROJECTIONS \
FAUX
DEFV(Local,DEFV(Logical,INIT(editer_les_projections,EDITER_LES_PROJECTIONS)));
/* Faut-il editer les projections ('VRAI') ou pas ('FAUX') ? Ceci a ete introduit le */
/* 20230317105152, la valeur par defaut garantissant la compatibilite anterieure... */
#define FAIRE_UNE_PROJECTION_PARALLELE \
VRAI
DEFV(Local,DEFV(Logical,INIT(faire_une_projection_parallele,FAIRE_UNE_PROJECTION_PARALLELE)));
/* Faut-il faire une projection parallele ('VRAI') ou perspective ('FAUX') ? Ceci a ete */
/* introduit le 20230316163011, la valeur par defaut garantissant la compatibilite */
/* anterieure... */
#define VALIDER_LA_QUALITE_DU_CALCUL_DE_Q1_ET_DE_Q2 \
FAUX
#define EPSILON_DE_VALIDATION_DU_CALCUL_DE_Q1_ET_DE_Q2 \
tgEPSILON
DEFV(Local,DEFV(Logical,INIT(valider_la_qualite_du_calcul_de_q1_et_de_q2,VALIDER_LA_QUALITE_DU_CALCUL_DE_Q1_ET_DE_Q2)));
DEFV(Local,DEFV(Float,INIT(epsilon_de_validation_du_calcul_de_q1_et_de_q2,EPSILON_DE_VALIDATION_DU_CALCUL_DE_Q1_ET_DE_Q2)));
/* Faut-il valider le calcul de {q1,q2} ('VRAI') ou pas ('FAUX'). */
/* */
/* L'epsilon de validation a ete introduit le 20230321103102... */
#include xrs/CalabiYau.41.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* P O U R D E S R A I S O N S D E C O M P A T I B I L I T E : */
/* */
/*************************************************************************************************************************************/
#include xrk/integr.1B.vv.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S I N I T I A L I S A T I O N S : */
/* */
/*************************************************************************************************************************************/
#include xrk/attractor.18.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N E V A R I E T E D E C A L A B I - Y A U 8 - D I M E N S I O N N E L L E : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/hyper_volumes.41.I"
#include xrs/hyper_volumes.22.I"
/*..............................................................................................................................*/
INITIALISATIONS_GENERALES;
/* Initialisations generales faites au tout debut... */
iTRANSFORMAT_31(liste_PARAMETRE_N1,PARAMETRE_N1_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_N2,PARAMETRE_N2_IMPLICITE);
/* Initialisation des parametres 'n1' et 'n2'. */
#include xrv/champs_5.1A.I"
#include xrs/hyper_volumes.13.I"
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(CONTROLE_DES_COMPATIBILITES_DES_FONCTIONS_DES_NOMBRES_HYPER_COMPLEXES;
/* Mis sous cette forme le 20230322100547... */
PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_FICHIER_PARAMETRE("PARAMETRE_N1=",_N1);
GET_FICHIER_PARAMETRE("PARAMETRE_N2=",_N2);
GET_ARGUMENT_F("HRA=",HReelle(parametre_qA));
GET_ARGUMENT_F("HIA=",HImaginaire(parametre_qA));
GET_ARGUMENT_F("HJA=",HJmaginaire(parametre_qA));
GET_ARGUMENT_F("HKA=",HKmaginaire(parametre_qA));
GET_ARGUMENT_F("HRB=",HReelle(parametre_qB));
GET_ARGUMENT_F("HIB=",HImaginaire(parametre_qB));
GET_ARGUMENT_F("HJB=",HJmaginaire(parametre_qB));
GET_ARGUMENT_F("HKB=",HKmaginaire(parametre_qB));
GET_ARGUMENT_L("editer_projections=""ep=",editer_les_projections);
/* Argument introduit le 20230317105152... */
/* */
/* On notera le 20230317114657 que l'on a interet a definir alors les arguments suivants : */
/* */
/* ModeFlot=f */
/* ChiffresFlot=6 */
/* */
/* pour ameliorer la lisibilite des editions... */
GET_ARGUMENT_F("xR1=",coefficient_de_projection_parallele_01_8D_3D_x1_X);
GET_ARGUMENT_F("xI1=",coefficient_de_projection_parallele_01_8D_3D_x2_X);
GET_ARGUMENT_F("xJ1=",coefficient_de_projection_parallele_01_8D_3D_x3_X);
GET_ARGUMENT_F("xK1=",coefficient_de_projection_parallele_01_8D_3D_x4_X);
GET_ARGUMENT_F("xR2=",coefficient_de_projection_parallele_01_8D_3D_x5_X);
GET_ARGUMENT_F("xI2=",coefficient_de_projection_parallele_01_8D_3D_x6_X);
GET_ARGUMENT_F("xJ2=",coefficient_de_projection_parallele_01_8D_3D_x7_X);
GET_ARGUMENT_F("xK2=",coefficient_de_projection_parallele_01_8D_3D_x8_X);
GET_ARGUMENT_F("yR1=",coefficient_de_projection_parallele_01_8D_3D_x1_Y);
GET_ARGUMENT_F("yI1=",coefficient_de_projection_parallele_01_8D_3D_x2_Y);
GET_ARGUMENT_F("yJ1=",coefficient_de_projection_parallele_01_8D_3D_x3_Y);
GET_ARGUMENT_F("yK1=",coefficient_de_projection_parallele_01_8D_3D_x4_Y);
GET_ARGUMENT_F("yR2=",coefficient_de_projection_parallele_01_8D_3D_x5_Y);
GET_ARGUMENT_F("yI2=",coefficient_de_projection_parallele_01_8D_3D_x6_Y);
GET_ARGUMENT_F("yJ2=",coefficient_de_projection_parallele_01_8D_3D_x7_Y);
GET_ARGUMENT_F("yK2=",coefficient_de_projection_parallele_01_8D_3D_x8_Y);
GET_ARGUMENT_F("zR1=",coefficient_de_projection_parallele_01_8D_3D_x1_Z);
GET_ARGUMENT_F("zI1=",coefficient_de_projection_parallele_01_8D_3D_x2_Z);
GET_ARGUMENT_F("zJ1=",coefficient_de_projection_parallele_01_8D_3D_x3_Z);
GET_ARGUMENT_F("zK1=",coefficient_de_projection_parallele_01_8D_3D_x4_Z);
GET_ARGUMENT_F("zR2=",coefficient_de_projection_parallele_01_8D_3D_x5_Z);
GET_ARGUMENT_F("zI2=",coefficient_de_projection_parallele_01_8D_3D_x6_Z);
GET_ARGUMENT_F("zJ2=",coefficient_de_projection_parallele_01_8D_3D_x7_Z);
GET_ARGUMENT_F("zK2=",coefficient_de_projection_parallele_01_8D_3D_x8_Z);
GET_ARGUMENT_L("ProjectionParallele=""pparallele=""ppar=",faire_une_projection_parallele);
GET_ARGUMENT_N("ProjectionPerspective=""pperspective=""pper=",faire_une_projection_parallele);
GET_ARGUMENT_F("O1=",observateur_de_projection_perspective_01_8D_X1);
GET_ARGUMENT_F("O2=",observateur_de_projection_perspective_01_8D_X2);
GET_ARGUMENT_F("O3=",observateur_de_projection_perspective_01_8D_X3);
GET_ARGUMENT_F("O4=",observateur_de_projection_perspective_01_8D_X4);
GET_ARGUMENT_F("O5=",observateur_de_projection_perspective_01_8D_X5);
GET_ARGUMENT_F("O6=",observateur_de_projection_perspective_01_8D_X6);
GET_ARGUMENT_F("O7=",observateur_de_projection_perspective_01_8D_X7);
GET_ARGUMENT_F("O8=",observateur_de_projection_perspective_01_8D_X8);
GET_ARGUMENT_F("PA=",plan_de_projection_perspective_01_8D_A);
GET_ARGUMENT_F("PB=",plan_de_projection_perspective_01_8D_B);
GET_ARGUMENT_F("PC=",plan_de_projection_perspective_01_8D_C);
GET_ARGUMENT_F("PD=",plan_de_projection_perspective_01_8D_D);
GET_ARGUMENT_F("PE=",plan_de_projection_perspective_01_8D_E);
GET_ARGUMENT_F("PF=",plan_de_projection_perspective_01_8D_F);
GET_ARGUMENT_F("PG=",plan_de_projection_perspective_01_8D_G);
GET_ARGUMENT_F("PH=",plan_de_projection_perspective_01_8D_H);
GET_ARGUMENT_F("PI=",plan_de_projection_perspective_01_8D_I);
/* Les arguments relatifs au choix de la projection ont ete introduit le 20230316163011 */
/* et le 20230316163223... */
GET_ARGUMENT_L("puissance_editer=",FHCpuissance_____editer_Argument_et_resultat);
/* Argument introduit le 20230321150520... */
GET_ARGUMENT_L("valider=""qualite=",valider_la_qualite_du_calcul_de_q1_et_de_q2);
GET_ARGUMENT_F("espilon_validation=""epsval=",epsilon_de_validation_du_calcul_de_q1_et_de_q2);
/* Argument introduit le 20230321103252... */
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v_w_t;
/* Argument introduit le 20230213121431 car oublie anterieurement... */
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;
PROCESS_ARGUMENTS_DE_GESTION_DE_LA_RELATIVITE_RESTREINTE;
)
);
#include xrs/hyper_volumes.12.I"
/* Visualisation et evolution de l'hyper-volume au cours du temps... */
Test(IL_FAUT(valider_la_qualite_du_calcul_de_q1_et_de_q2))
/* Test introduit le 20230403171036 suite a 'v $ximcd/operator$FON 2023032410053'... */
Bblock
PRINT_ATTENTION("la mutiplication des hyper-complexes n'etant pas commutative, cette validation n'a pas de sens");
Eblock
ATes
Bblock
Eblock
ETes
VISUALISATION_GENERALE_DE_L_HYPER_VOLUME(BLOC(Bblock
EGAL(parametre_n1,sPARAMETRE_N1(numero_de_la_periode_courante));
EGAL(parametre_n2,sPARAMETRE_N2(numero_de_la_periode_courante));
/* Calcul de la valeur des parametres 'n1' et 'n2'. */
Test(N_EST_PAS_ENTIER(parametre_n1))
Bblock
PRINT_ATTENTION("le parametre 'N1' n'est pas entier");
Eblock
ATes
Bblock
Eblock
ETes
Test(N_EST_PAS_ENTIER(parametre_n2))
Bblock
PRINT_ATTENTION("le parametre 'N2' n'est pas entier");
Eblock
ATes
Bblock
Eblock
ETes
Eblock
)
,BLOC(Bblock
DEFV(Float,INIT(minimum_k1
,FLOT(RESZ(NEUT(INTE(parametre_n1)),INTE(parametre_n1)))
)
);
DEFV(Float,INIT(maximum_k1
,FLOT(RESZ(PRED(ARRX(parametre_n1)),ARRX(parametre_n1)))
)
);
DEFV(Float,INIT(minimum_k2
,FLOT(RESZ(NEUT(INTE(parametre_n2)),INTE(parametre_n2)))
)
);
DEFV(Float,INIT(maximum_k2
,FLOT(RESZ(PRED(ARRX(parametre_n2)),ARRX(parametre_n2)))
)
);
/* Afin de faire evoluer 'k1' et 'k2'. */
DEFV(Float,INIT(facteur_1_de_la_translation_de_uv,FLOT__UNDEF));
DEFV(Float,INIT(facteur_2_de_la_translation_de_uv,FLOT__UNDEF));
/* Pour tester le changement des associations {u,v} <--> {1,2}. */
EGAL(facteur_1_de_la_translation_de_uv,FZERO);
DoIn(parametre_k1,minimum_k1,maximum_k1,FU)
Bblock
EGAL(facteur_2_de_la_translation_de_uv,FZERO);
DoIn(parametre_k2,minimum_k2,maximum_k2,FU)
/* Boucles sur les N1xN2 "patches" composant la variete de Calabi-Yau. */
Bblock
GENER_HYPER_VOLUME_AVEC_DIFFERENCES____(Fx_CalabiYau_4
,Fy_CalabiYau_4
,Fz_CalabiYau_4
,BLOC(VIDE;)
);
/* Il n'y a rien de specifique a calculer pour le point courant... */
Eblock
EDoI
Eblock
EDoI
Eblock
)
,Fx_CalabiYau_4,Fy_CalabiYau_4,Fz_CalabiYau_4
/* Fonctions definissant parametriquement la surface... */
,Pxyz_CalabiYau_4
);
RETU_Commande;
Eblock
ECommande