/*************************************************************************************************************************************/
/* */
/* 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.32$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20230202182100). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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_CPP_____AUTORISER_LE_GooF
#define PRAGMA_CPP_____SI_LE_GooF_EST_ACTIVABLE_ALORS_FAIRE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_FLOTTANTS
@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 la surface ('v $xrs/CalabiYau.32$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.32$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 ] */
/* */
/* (ou 'p' designe 'pi'). */
/* */
/* */
/*************************************************************************************************************************************/
#include xrs/CalabiYau.32.I"
#include xrs/CalabiYau.34.I"
#include xrs/CalabiYau.35.I"
#include xrs/CalabiYau.37.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/surfaces.21.I"
#define PONDERATION_ROUGE_VERTE_BLEUE_POUR_dcx_dcy_dcz \
FU
#define PRENDRE_LA_VALEUR_ABSOLUE_DES_DIFFERENTIELLES \
VRAI
#include xrs/CalabiYau.36.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 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 20230321103719... */
#include xrs/CalabiYau.31.I"
#define COMPATIBILITE_2009110209 \
VRAI
DEFV(Local,DEFV(Logical,INIT(compatibilite_2009110209,COMPATIBILITE_2009110209)));
/* Introduit le 20230203170910 en notant bien que la date de cette compatibilite vient */
/* de 'v $xrs/CalabiYau.12$K COMPATIBILITE_2009110209' et qu'elle conservee ici en l'etat... */
/* Je ne suis pas sur que cela soit utile, mais on ne sait jamais... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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/surfaces.41.I"
#include xrs/surfaces.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/surfaces.13.I"
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(GET_ARGUMENT_L("compatibilite_2009110209=",compatibilite_2009110209);
/* Argument introduit le 20230203170910... */
CONTROLE_DES_COMPATIBILITES_DES_FONCTIONS_DES_NOMBRES_HYPER_COMPLEXES;
/* Mis sous cette forme le 20230322100423... */
PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_FICHIER_PARAMETRE("PARAMETRE_N1=",_N1);
GET_FICHIER_PARAMETRE("PARAMETRE_N2=",_N2);
GET_ARGUMENT_F("u_prime=""up=",coordonnee_u_prime);
GET_ARGUMENT_F("v_prime=""vp=",coordonnee_v_prime);
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_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("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 20230321103719... */
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_TEXTURE;
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_NORMALE;
PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;
PROCESS_ARGUMENTS_DE_DIFFERENTIATION_PAR_DIFFERENCES;
PROCESS_ARGUMENTS_D_ACTIVATION_DE_L_ARITHMETIQUE_ETENDUE_ET_NON_BASIQUE_POUR_Fx_Fy_Fz;
PROCESS_ARGUMENTS_DE_GESTION_DE_LA_RELATIVITE_RESTREINTE;
)
);
#include xrs/surfaces.12.I"
/* Visualisation et evolution de la surface au cours du temps... */
VISUALISATION_GENERALE_DE_LA_SURFACE(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
EGAL(facteur_de_la_translation_de_u_effectif_lors_d_un_lister
,facteur_1_de_la_translation_de_uv
);
EGAL(facteur_de_la_translation_de_v_effectif_lors_d_un_lister
,facteur_2_de_la_translation_de_uv
);
GENERATION_D_UNE_IMAGE_DE_LA_SURFACE_AVEC_DIFFERENCES____(Fx_CalabiYau_3
,Fy_CalabiYau_3
,Fz_CalabiYau_3
,BLOC(VIDE;)
);
/* Il n'y a rien de specifique a calculer pour le point courant... */
INCR(facteur_2_de_la_translation_de_uv
,COND(IL_FAUT(compatibilite_2009110209),FZERO,FU)
);
Eblock
EDoI
INCR(facteur_1_de_la_translation_de_uv
,COND(IL_FAUT(compatibilite_2009110209),FZERO,FU)
);
Eblock
EDoI
Eblock
)
,Fx_CalabiYau_3,Fy_CalabiYau_3,Fz_CalabiYau_3
/* Fonctions definissant parametriquement la surface... */
,Pxyz_CalabiYau_3
);
RETU_Commande;
Eblock
ECommande