/*************************************************************************************************************************************/
/* */
/* 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 4 - D I M E N S I O N N E L L E : */
/* */
/* */
/* Exemple : */
/* */
/* L'image 'v $xiirs/CAYA.61.0129' */
/* en est un exemple... */
/* */
/* */
/* Author of '$xrs/CalabiYau.12$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20010528145410). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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
/* Le 'GooF' fut introduit le 20080126104007... */
@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
/* Introduit le 20041229180623 pour 'v $xrs/surfaces.12$I IFmageA4'... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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"
/* Introduit le 20111209173154... */
#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 4 - D I M E N S I O N N E L L E : */
/* */
/* */
/* Definition de la surface ('v $xrs/CalabiYau.12$K') : */
/* */
/* Une variete complexe de Calabi-Yau */
/* est definie implicitement par : */
/* */
/* n1 n2 */
/* Z + Z = 1 */
/* 1 2 */
/* */
/* ou 'Z1' et 'Z2' sont deux nombres complexes. */
/* Cette variete est donc a 4 dimensions, et la */
/* surface a 3 dimensions que l'on visualisera */
/* sera definie par : */
/* */
/* */
/* F (u,v) = A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) */
/* x xR1 1 xI1 1 xR2 2 xI2 2 */
/* */
/* */
/* F (u,v) = A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) */
/* y yR1 1 yI1 1 yR2 2 yI2 2 */
/* */
/* */
/* F (u,v) = A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) */
/* z zR1 1 zI1 1 zR2 2 zI2 2 */
/* */
/* */
/* soit, par defaut : */
/* */
/* */
/* F (u,v) = R(Z ) + I(Z ) */
/* x 1 2 */
/* */
/* */
/* F (u,v) = I(Z ) + I(Z ) */
/* y 1 2 */
/* */
/* */
/* F (u,v) = R(Z ) + I(Z ) */
/* z 2 2 */
/* */
/* */
/* (ou 'R(...)' et 'I(...)' designent respectivement */
/* les parties Reelle et Imaginaire) avec : */
/* */
/* u E [ 0 , p/2 ] */
/* */
/* v E [ -1 , +1 ] */
/* */
/* (ou 'p' designe 'pi'). */
/* */
/* */
/*************************************************************************************************************************************/
#include xrs/CalabiYau.12.I"
#include xrs/CalabiYau.14.I"
#include xrs/CalabiYau.15.I"
#include xrs/CalabiYau.17.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 \
/* Afin d'assurer la compatibilite anterieure (introduit le 20060130113630)... */
#define PRENDRE_LA_VALEUR_ABSOLUE_DES_DIFFERENTIELLES \
VRAI \
/* Introduit le 20060131120842 pour le probleme 'v $xrs/surfaces.14$I 20060131113621'... */
#include xrs/CalabiYau.16.I"
#define COMPATIBILITE_2009110209 \
VRAI
DEFV(Local,DEFV(Logical,INIT(compatibilite_2009110209,COMPATIBILITE_2009110209)));
/* Permet de proceder a des 'lister_u_v=VRAI' compatibles a ceux qui furent effectues */
/* anterieurement au 20091102092405. */
/* */
/* Le 20091103090312, principalement grace a l'image 'v $xiirs/CAYA.A1.5x5.1', je note */
/* que les "patches" ne se raccordent pas comme je le croyais. En particulier, sur cette */
/* image, des symetries bord a bord entre couple de "patches" voisins sont visibles. Cela */
/* rend impossible la mise en place d'un systeme de coordonnees {u,v} globales par */
/* translation (par des multiples des periodes respectives de 'u' et de 'v') des coordonnees */
/* {u,v} locales a chaque "patch". Au passage, c'est cela qui explique (enfin !), les sortes */
/* de "creneaux" qui etaient observes sur le "bord" de 'v $xiirs/CAYA.A1.0129'. C'est donc */
/* pourquoi, la valeur par defaut passe de 'FAUX' a 'VRAI' a cette date, meme si le */
/* dispositif est conserve, on ne sait jamais. En conclusion, cela condamne l'utilisation */
/* de 'v $xrs/SurfParti.01$Z' avec '$xrs/CalabiYau.12$X' comme generateur ainsi que je */
/* l'esperais (provisoirement ?). L'image 'v $xiirs/CAYA.A1.5x5.6' montre elle-aussi tres */
/* clairement qu'en fait il semble que l'on passe d'un patch a son voisin par une rotation */
/* de pi/2 (mais est-ce aussi simple ?). Donc plus que de translater les coordonnees {u,v}, */
/* il conviendrait de les permuter, mais malheureusement, elles ne sont pas definies dans */
/* le meme segment : */
/* */
/* u E [0 , p/2] */
/* v E [-1 , +1]. */
/* */
/* d'ou un gros probleme... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 */
/* 20230317105202, 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 20230310122110, la valeur par defaut garantissant la compatibilite */
/* anterieure... */
#define VALIDER_LA_QUALITE_DU_CALCUL_DE_Z1_ET_DE_Z2 \
FAUX
#define EPSILON_DE_VALIDATION_DU_CALCUL_DE_Z1_ET_DE_Z2 \
tgEPSILON
DEFV(Local,DEFV(Logical,INIT(valider_la_qualite_du_calcul_de_z1_et_de_z2,VALIDER_LA_QUALITE_DU_CALCUL_DE_Z1_ET_DE_Z2)));
DEFV(Local,DEFV(Float,INIT(epsilon_de_validation_du_calcul_de_z1_et_de_z2,EPSILON_DE_VALIDATION_DU_CALCUL_DE_Z1_ET_DE_Z2)));
/* Faut-il valider le calcul de {z1,z2} ('VRAI') ou pas ('FAUX'). */
/* */
/* L'epsilon de validation a ete introduit le 20230321104411... */
#include xrs/CalabiYau.11.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 4 - 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"
/* Introduit le 20150528134956... */
#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 20091102092405... */
CONTROLE_DES_COMPATIBILITES_DES_FONCTIONS_DES_NOMBRES_COMPLEXES;
/* Mis sous cette forme le 20230322100156... */
PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_FICHIER_PARAMETRE("PARAMETRE_N1=",_N1);
GET_FICHIER_PARAMETRE("PARAMETRE_N2=",_N2);
GET_ARGUMENT_F("RA=",Reelle(parametre_zA));
GET_ARGUMENT_F("IA=",Imaginaire(parametre_zA));
GET_ARGUMENT_F("RB=",Reelle(parametre_zB));
GET_ARGUMENT_F("IB=",Imaginaire(parametre_zB));
GET_ARGUMENT_L("editer_projections=""ep=",editer_les_projections);
/* Argument introduit le 20230317105202... */
/* */
/* On notera le 20230317114731 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_4D_3D_x1_X);
GET_ARGUMENT_F("xI1=",coefficient_de_projection_parallele_01_4D_3D_x2_X);
GET_ARGUMENT_F("xR2=",coefficient_de_projection_parallele_01_4D_3D_x3_X);
GET_ARGUMENT_F("xI2=",coefficient_de_projection_parallele_01_4D_3D_x4_X);
GET_ARGUMENT_F("yR1=",coefficient_de_projection_parallele_01_4D_3D_x1_Y);
GET_ARGUMENT_F("yI1=",coefficient_de_projection_parallele_01_4D_3D_x2_Y);
GET_ARGUMENT_F("yR2=",coefficient_de_projection_parallele_01_4D_3D_x3_Y);
GET_ARGUMENT_F("yI2=",coefficient_de_projection_parallele_01_4D_3D_x4_Y);
GET_ARGUMENT_F("zR1=",coefficient_de_projection_parallele_01_4D_3D_x1_Z);
GET_ARGUMENT_F("zI1=",coefficient_de_projection_parallele_01_4D_3D_x2_Z);
GET_ARGUMENT_F("zR2=",coefficient_de_projection_parallele_01_4D_3D_x3_Z);
GET_ARGUMENT_F("zI2=",coefficient_de_projection_parallele_01_4D_3D_x4_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_4D_X1);
GET_ARGUMENT_F("O2=",observateur_de_projection_perspective_01_4D_X2);
GET_ARGUMENT_F("O3=",observateur_de_projection_perspective_01_4D_X3);
GET_ARGUMENT_F("O4=",observateur_de_projection_perspective_01_4D_X4);
GET_ARGUMENT_F("PA=",plan_de_projection_perspective_01_4D_A);
GET_ARGUMENT_F("PB=",plan_de_projection_perspective_01_4D_B);
GET_ARGUMENT_F("PC=",plan_de_projection_perspective_01_4D_C);
GET_ARGUMENT_F("PD=",plan_de_projection_perspective_01_4D_D);
GET_ARGUMENT_F("PE=",plan_de_projection_perspective_01_4D_E);
/* Les arguments relatifs au choix de la projection ont ete introduit le 20230310122110... */
GET_ARGUMENT_L("valider=""qualite=",valider_la_qualite_du_calcul_de_z1_et_de_z2);
GET_ARGUMENT_F("espilon_validation=""epsval=",epsilon_de_validation_du_calcul_de_z1_et_de_z2);
/* Argument introduit le 20230321103252... */
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;
/* Introduit le 20060130112542. */
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
/* Introduit le 20060130113630. */
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_TEXTURE;
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_NORMALE;
/* Introduit le 20060130113630. */
PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;
PROCESS_ARGUMENTS_DE_DIFFERENTIATION_PAR_DIFFERENCES;
/* Introduit le 20060131090342... */
PROCESS_ARGUMENTS_D_ACTIVATION_DE_L_ARITHMETIQUE_ETENDUE_ET_NON_BASIQUE_POUR_Fx_Fy_Fz;
/* Introduit le 20080126104007... */
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));
/* Introduits le 20091102141848 pour tester le changement des associations {u,v} <--> {1,2}. */
EGAL(facteur_1_de_la_translation_de_uv,FZERO);
/* A priori et quel que soit 'compatibilite_2009110209' (introduit le 20091102092405). */
DoIn(parametre_k1,minimum_k1,maximum_k1,FU)
Bblock
EGAL(facteur_2_de_la_translation_de_uv,FZERO);
/* A priori et quel que soit 'compatibilite_2009110209' (introduit le 20091102092405). */
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_1
,Fy_CalabiYau_1
,Fz_CalabiYau_1
,BLOC(VIDE;)
);
/* Il n'y a rien de specifique a calculer pour le point courant... */
/* */
/* Le 20060130113630, 'GENERATION_D_UNE_IMAGE_DE_LA_SURFACE_SANS_DIFFERENTIATION...)' est */
/* devenu 'GENERATION_D_UNE_IMAGE_DE_LA_SURFACE_AVEC_DIFFERENCES____(...)'. */
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_1,Fy_CalabiYau_1,Fz_CalabiYau_1
/* Fonctions definissant parametriquement la surface... */
,Pxyz_CalabiYau_1
);
RETU_Commande;
Eblock
ECommande