/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D U " P A T C H " { K1 , K2 } */
/* D E L A 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 { N1 , N2 } : */
/* */
/* */
/* Author of '$xrs/CalabiYau.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20010525151422). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 20080126103913... */
@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 20041229180323 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 20111209173327... */
#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.11$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 20060130115940)... */
#define PRENDRE_LA_VALEUR_ABSOLUE_DES_DIFFERENTIELLES \
VRAI \
/* Introduit le 20060131120756 pour le probleme 'v $xrs/surfaces.14$I 20060131113621'... */
#include xrs/CalabiYau.13.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 */
/* 20230320110647, 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 20230310122105, 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 20230321104611... */
#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 " P A T C H " { K1 , K2 } */
/* D E L A 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 { N1 , N2 } : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/surfaces.41.I"
/* Introduit le 20150528134910... */
#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'. */
iTRANSFORMAT_31(liste_PARAMETRE_K1,PARAMETRE_K1_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_K2,PARAMETRE_K2_IMPLICITE);
/* Initialisation des parametres 'k1' et 'k2'. */
#include xrv/champs_5.1A.I"
#include xrs/surfaces.13.I"
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(CONTROLE_DES_COMPATIBILITES_DES_FONCTIONS_DES_NOMBRES_COMPLEXES;
/* Mis sous cette forme le 20230322100044... */
PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_FICHIER_PARAMETRE("PARAMETRE_N1=",_N1);
GET_FICHIER_PARAMETRE("PARAMETRE_N2=",_N2);
GET_FICHIER_PARAMETRE("PARAMETRE_K1=",_K1);
GET_FICHIER_PARAMETRE("PARAMETRE_K2=",_K2);
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 20230320110647... */
/* */
/* 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 20230310122105... */
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 20230321104611... */
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;
/* Introduit le 20060117145200. */
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
/* Introduit le 20050615085652. */
PROCESS_ARGUMENTS_DE_PROJECTION_DE_LA_SURFACE;
/* Introduit le 20041126210153. */
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_TEXTURE;
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_NORMALE;
PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;
PROCESS_ARGUMENTS_DE_DEFINITION_D_UNE_DROITE_DU_PLAN_UV;
/* Introduit le 20060831113024... */
PROCESS_ARGUMENTS_DE_DIFFERENTIATION_PAR_DIFFERENCES;
/* Introduit le 20060131090152... */
PROCESS_ARGUMENTS_D_ACTIVATION_DE_L_ARITHMETIQUE_ETENDUE_ET_NON_BASIQUE_POUR_Fx_Fy_Fz;
/* Introduit le 20080126103913... */
PROCESS_ARGUMENTS_DE_GESTION_DE_LA_RELATIVITE_RESTREINTE;
)
);
#include xrs/surfaces.12.I"
/* Visualisation et evolution de la surface au cours du temps... */
VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENCES____(BLOC(Bblock
EGAL(parametre_n1,sPARAMETRE_N1(numero_de_la_periode_courante));
EGAL(parametre_n2,sPARAMETRE_N2(numero_de_la_periode_courante));
EGAL(parametre_k1,sPARAMETRE_K1(numero_de_la_periode_courante));
EGAL(parametre_k2,sPARAMETRE_K2(numero_de_la_periode_courante));
/* Calcul de la valeur des parametres 'n1', 'n2', 'k1' et 'k2'. */
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
Test(N_EST_PAS_ENTIER(parametre_k1))
Bblock
PRINT_ATTENTION("le parametre 'K1' n'est pas entier");
Eblock
ATes
Bblock
Eblock
ETes
Test(N_EST_PAS_ENTIER(parametre_k2))
Bblock
PRINT_ATTENTION("le parametre 'K2' n'est pas entier");
Eblock
ATes
Bblock
Eblock
ETes
Test(IFET(IFET(EST_ENTIER(parametre_n1)
,EST_ENTIER(parametre_n2)
)
,IFET(EST_ENTIER(parametre_k1)
,EST_ENTIER(parametre_k2)
)
)
)
Bblock
DEFV(Int,INIT(minimum_k1
,RESZ(NEUT(INTE(parametre_n1)),INTE(parametre_n1))
)
);
DEFV(Int,INIT(maximum_k1
,RESZ(PRED(INTE(parametre_n1)),INTE(parametre_n1))
)
);
DEFV(Int,INIT(minimum_k2
,RESZ(NEUT(INTE(parametre_n2)),INTE(parametre_n2))
)
);
DEFV(Int,INIT(maximum_k2
,RESZ(PRED(INTE(parametre_n2)),INTE(parametre_n2))
)
);
/* Afin de valider 'k1' et 'k2'. */
Test(IFEXff(INTE(parametre_k1),minimum_k1,maximum_k1))
Bblock
PRINT_ERREUR("le parametre 'K1' est incorrect");
CAL1(Prer3("il vaut %d et est donc hors de [%d,%d]\n"
,INTE(parametre_k1)
,minimum_k1
,maximum_k1
)
);
Eblock
ATes
Bblock
Eblock
ETes
Test(IFEXff(INTE(parametre_k2),minimum_k2,maximum_k2))
Bblock
PRINT_ERREUR("le parametre 'K2' est incorrect");
CAL1(Prer3("il vaut %d et est donc hors de [%d,%d]\n"
,INTE(parametre_k2)
,minimum_k2
,maximum_k2
)
);
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
Eblock
)
,Fx_CalabiYau_1,Fy_CalabiYau_1,Fz_CalabiYau_1
/* Fonctions definissant parametriquement la surface... */
,Pxyz_CalabiYau_1
,BLOC(VIDE;)
/* Il n'y a rien de specifique a calculer pour le point courant... */
);
/* Le 20060130115940, 'VISUALISATION_DE_LA_SURFACE_SANS_DIFFERENTIATION...)' est */
/* devenu 'VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENCES____(...)'. */
RETU_Commande;
Eblock
ECommande