/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D U " P A T C H " { K1 , K2 , K3 } */
/* D E L A V A R I E T E D E C A L A B I - Y A U 6 - D I M E N S I O N N E L L E { N1 , N2 , N3 } : */
/* */
/* */
/* Author of '$xrs/CalabiYau.21$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20010531085531). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 20080126104041... */
@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 20041229180645 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 20111209172434... */
#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 6 - D I M E N S I O N N E L L E : */
/* */
/* */
/* Definition de la surface ('v $xrs/CalabiYau.21$K') : */
/* */
/* Une variete complexe de Calabi-Yau */
/* est definie implicitement par : */
/* */
/* n1 n2 n3 */
/* Z + Z + Z = 1 */
/* 1 2 3 */
/* */
/* ou 'Z1', 'Z2' et 'Z3' sont trois nombres complexes. */
/* Cette variete est donc a 6 dimensions, et la */
/* surface a 3 dimensions que l'on visualisera */
/* sera definie par : */
/* */
/* */
/* F (u1,v1,u2,v2) = A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) */
/* x xR1 1 xI1 1 xR2 2 xI2 2 xR3 3 xI3 3 */
/* */
/* */
/* F (u1,v1,u2,v2) = A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) */
/* y yR1 1 yI1 1 yR2 2 yI2 2 yR3 3 yI3 3 */
/* */
/* */
/* F (u1,v1,u2,v2) = A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) + A R(Z ) + A I(Z ) */
/* z zR1 1 zI1 1 zR2 2 zI2 2 zR3 3 zI3 3 */
/* */
/* */
/* soit, par defaut : */
/* */
/* */
/* F (u1,v1,u2,v2) = R(Z ) + I(Z ) */
/* x 1 2 */
/* */
/* */
/* F (u1,v1,u2,v2) = I(Z ) + R(Z ) */
/* y 1 3 */
/* */
/* */
/* F (u1,v1,u2,v2) = R(Z ) + I(Z ) */
/* z 2 3 */
/* */
/* */
/* (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' et ou '?' represente */
/* 'phi' et 'theta'). */
/* */
/* */
/*************************************************************************************************************************************/
#include xrs/CalabiYau.22.I"
#include xrs/CalabiYau.14.I"
#include xrs/CalabiYau.25.I"
#include xrs/CalabiYau.27.I"
#define ALPHA_u1 \
FU
DEFV(Local,DEFV(Float,INIT(alpha_u1,ALPHA_u1)));
#define BETA_u1 \
FZERO
DEFV(Local,DEFV(Float,INIT(beta_u1,BETA_u1)));
#define ALPHA_v1 \
FU
DEFV(Local,DEFV(Float,INIT(alpha_v1,ALPHA_v1)));
#define BETA_v1 \
FZERO
DEFV(Local,DEFV(Float,INIT(beta_v1,BETA_v1)));
#define ALPHA_u2 \
FU
DEFV(Local,DEFV(Float,INIT(alpha_u2,ALPHA_u2)));
#define BETA_u2 \
FZERO
DEFV(Local,DEFV(Float,INIT(beta_u2,BETA_u2)));
#define ALPHA_v2 \
FU
DEFV(Local,DEFV(Float,INIT(alpha_v2,ALPHA_v2)));
#define BETA_v2 \
FZERO
DEFV(Local,DEFV(Float,INIT(beta_v2,BETA_v2)));
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 20060130120254)... */
#define PRENDRE_LA_VALEUR_ABSOLUE_DES_DIFFERENTIELLES \
VRAI \
/* Introduit le 20060131120905 pour le probleme 'v $xrs/surfaces.14$I 20060131113621'... */
#include xrs/CalabiYau.23.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_Z1_DE_Z2_ET_DE_Z3 \
FAUX
#define EPSILON_DE_VALIDATION_DU_CALCUL_DE_Z1_DE_Z2_ET_DE_Z3 \
tgEPSILON
DEFV(Local,DEFV(Logical,INIT(valider_la_qualite_du_calcul_de_z1_de_z2_et_de_z3,VALIDER_LA_QUALITE_DU_CALCUL_DE_Z1_DE_Z2_ET_DE_Z3)));
DEFV(Local,DEFV(Float,INIT(epsilon_de_validation_du_calcul_de_z1_de_z2_et_de_z3,EPSILON_DE_VALIDATION_DU_CALCUL_DE_Z1_DE_Z2_ET_DE_Z3)));
/* Faut-il valider le calcul de {z1,z2,z3} ('VRAI') ou pas ('FAUX'). */
/* */
/* L'epsilon de validation a ete introduit le 20230321104236... */
#include xrs/CalabiYau.21.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 , K3 } */
/* D E L A V A R I E T E D E C A L A B I - Y A U 6 - D I M E N S I O N N E L L E { N1 , N2 , N3 } : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/surfaces.41.I"
/* Introduit le 20150528135027... */
#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);
iTRANSFORMAT_31(liste_PARAMETRE_N3,PARAMETRE_N3_IMPLICITE);
/* Initialisation des parametres 'n1', 'n2' et 'n3'. */
iTRANSFORMAT_31(liste_PARAMETRE_K1,PARAMETRE_K1_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_K2,PARAMETRE_K2_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_K3,PARAMETRE_K3_IMPLICITE);
/* Initialisation des parametres 'k1', 'k2' et 'k3'. */
#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 20230322100242... */
PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_FICHIER_PARAMETRE("PARAMETRE_N1=",_N1);
GET_FICHIER_PARAMETRE("PARAMETRE_N2=",_N2);
GET_FICHIER_PARAMETRE("PARAMETRE_N3=",_N3);
GET_FICHIER_PARAMETRE("PARAMETRE_K1=",_K1);
GET_FICHIER_PARAMETRE("PARAMETRE_K2=",_K2);
GET_FICHIER_PARAMETRE("PARAMETRE_K3=",_K3);
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_F("RC=",Reelle(parametre_zC));
GET_ARGUMENT_F("IC=",Imaginaire(parametre_zC));
GET_ARGUMENT_F("au1=",alpha_u1);
GET_ARGUMENT_F("bu1=",beta_u1);
GET_ARGUMENT_F("av1=",alpha_v1);
GET_ARGUMENT_F("bv1=",beta_v1);
GET_ARGUMENT_F("au2=",alpha_u2);
GET_ARGUMENT_F("bu2=",beta_u2);
GET_ARGUMENT_F("av2=",alpha_v2);
GET_ARGUMENT_F("bv2=",beta_v2);
GET_ARGUMENT_F("xR1=",coefficient_de_projection_parallele_01_6D_3D_x1_X);
GET_ARGUMENT_F("xI1=",coefficient_de_projection_parallele_01_6D_3D_x2_X);
GET_ARGUMENT_F("xR2=",coefficient_de_projection_parallele_01_6D_3D_x3_X);
GET_ARGUMENT_F("xI2=",coefficient_de_projection_parallele_01_6D_3D_x4_X);
GET_ARGUMENT_F("xR3=",coefficient_de_projection_parallele_01_6D_3D_x5_X);
GET_ARGUMENT_F("xI3=",coefficient_de_projection_parallele_01_6D_3D_x6_X);
GET_ARGUMENT_F("yR1=",coefficient_de_projection_parallele_01_6D_3D_x1_Y);
GET_ARGUMENT_F("yI1=",coefficient_de_projection_parallele_01_6D_3D_x2_Y);
GET_ARGUMENT_F("yR2=",coefficient_de_projection_parallele_01_6D_3D_x3_Y);
GET_ARGUMENT_F("yI2=",coefficient_de_projection_parallele_01_6D_3D_x4_Y);
GET_ARGUMENT_F("yR3=",coefficient_de_projection_parallele_01_6D_3D_x5_Y);
GET_ARGUMENT_F("yI3=",coefficient_de_projection_parallele_01_6D_3D_x6_Y);
GET_ARGUMENT_F("zR1=",coefficient_de_projection_parallele_01_6D_3D_x1_Z);
GET_ARGUMENT_F("zI1=",coefficient_de_projection_parallele_01_6D_3D_x2_Z);
GET_ARGUMENT_F("zR2=",coefficient_de_projection_parallele_01_6D_3D_x3_Z);
GET_ARGUMENT_F("zI2=",coefficient_de_projection_parallele_01_6D_3D_x4_Z);
GET_ARGUMENT_F("zR3=",coefficient_de_projection_parallele_01_6D_3D_x5_Z);
GET_ARGUMENT_F("zI3=",coefficient_de_projection_parallele_01_6D_3D_x6_Z);
GET_ARGUMENT_L("valider=""qualite=",valider_la_qualite_du_calcul_de_z1_de_z2_et_de_z3);
GET_ARGUMENT_F("espilon_validation=""epsval=",epsilon_de_validation_du_calcul_de_z1_de_z2_et_de_z3);
/* Argument introduit le 20230321104236... */
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;
/* Introduit le 20060117145222. */
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
/* Introduit le 20050615085715. */
PROCESS_ARGUMENTS_DE_PROJECTION_DE_LA_SURFACE;
/* Introduit le 20041126211107. */
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 20060831113131... */
PROCESS_ARGUMENTS_DE_DIFFERENTIATION_PAR_DIFFERENCES;
/* Introduit le 20060131090229... */
PROCESS_ARGUMENTS_D_ACTIVATION_DE_L_ARITHMETIQUE_ETENDUE_ET_NON_BASIQUE_POUR_Fx_Fy_Fz;
/* Introduit le 20080126104041... */
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_n3,sPARAMETRE_N3(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));
EGAL(parametre_k3,sPARAMETRE_K3(numero_de_la_periode_courante));
/* Calcul de la valeur des parametres 'n1', 'n2', ,'n3, 'k1', k2' et 'k3'. */
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_n3))
Bblock
PRINT_ATTENTION("le parametre 'N3' 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(N_EST_PAS_ENTIER(parametre_k3))
Bblock
PRINT_ATTENTION("le parametre 'K3' n'est pas entier");
Eblock
ATes
Bblock
Eblock
ETes
Test(IFET(I3ET(EST_ENTIER(parametre_n1)
,EST_ENTIER(parametre_n2)
,EST_ENTIER(parametre_n3)
)
,I3ET(EST_ENTIER(parametre_k1)
,EST_ENTIER(parametre_k2)
,EST_ENTIER(parametre_k3)
)
)
)
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))
)
);
DEFV(Int,INIT(minimum_k3
,RESZ(NEUT(INTE(parametre_n3)),INTE(parametre_n3))
)
);
DEFV(Int,INIT(maximum_k3
,RESZ(PRED(INTE(parametre_n3)),INTE(parametre_n3))
)
);
/* Afin de valider 'k1', k2' et 'k3'. */
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
Test(IFEXff(INTE(parametre_k3),minimum_k3,maximum_k3))
Bblock
PRINT_ERREUR("le parametre 'K3' est incorrect");
CAL1(Prer3("il vaut %d et est donc hors de [%d,%d]\n"
,INTE(parametre_k3)
,minimum_k3
,maximum_k3
)
);
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
Eblock
)
,Fx_CalabiYau_2,Fy_CalabiYau_2,Fz_CalabiYau_2
/* Fonctions definissant parametriquement la surface... */
,Pxyz_CalabiYau_2
,BLOC(VIDE;)
/* Il n'y a rien de specifique a calculer pour le point courant... */
);
/* Le 20060131090039, 'VISUALISATION_DE_LA_SURFACE_SANS_DIFFERENTIATION...)' est */
/* devenu 'VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENCES____(...)'. */
RETU_Commande;
Eblock
ECommande