/*************************************************************************************************************************************/ /* */ /* D E F I N I T I O N D ' U N " R U B A N D E M O B I U S " T R I D I M E N S I O N N E L : */ /* */ /* */ /* Exemple : */ /* */ /* L'image 'v $xiirs/MOB3.21' */ /* en est un exemple... */ /* */ /* */ /* Author of '$xrs/Mobius3D.21$K' : */ /* */ /* Jean-Francois COLONNA (LACTAMME, 20040402110222). */ /* */ /*************************************************************************************************************************************/ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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 /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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 20111209172801... */ #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(-2.0) #define hYmin_ESPACE \ PARE(-2.0) #define hZmin_ESPACE \ PARE(-2.0) /* Definition du "coin" inferieur-gauche-arriere de l'espace physique. */ #define hXmax_ESPACE \ PARE(+2.0) #define hYmax_ESPACE \ PARE(+2.0) #define hZmax_ESPACE \ PARE(+2.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(-2.0) #define dYmin_ESPACE \ PARE(-2.0) #define dZmin_ESPACE \ PARE(-2.0) /* Definition des minima des differentielles {dx,dy,dz}. */ #define dXmax_ESPACE \ PARE(+2.0) #define dYmax_ESPACE \ PARE(+2.0) #define dZmax_ESPACE \ PARE(+2.0) /* Definition des maxima des differentielles {dx,dy,dz}. */ /*===================================================================================================================================*/ /* :Debut_listMN_RUBAN_DE_MOBIUS_TRIDIMENSIONNEL_21: */ /*************************************************************************************************************************************/ /* */ /* D E F I N I T I O N D ' U N " R U B A N D E M O B I U S " T R I D I M E N S I O N N E L : */ /* */ /* */ /* Definition de la variete (due a Patrice Jeener, 'v $xrs/Mobius3D.21$K') : */ /* */ /* Le ruban de Mobius tridimensionnel est */ /* defini dans un espace reel a 4 dimensions */ /* (ou complexe a 2 dimensions {Z1,Z2}). Il est */ /* parametre a l'aide de 3 coordonnees {u,v,w} */ /* (en assimilant le ruban de Mobius tridimensionnel */ /* a un tore, la coordonnee 'v' decrit un grand cercle */ /* de ce tore). Le point courant {Z1,Z2} est alors */ /* defini par : */ /* */ /* i.v */ /* Z = (u+t).e */ /* 1 */ /* */ /* v */ /* i.--- */ /* d */ /* Z = w.e */ /* 2 */ /* */ /* avec en general : */ /* */ /* d = +2 */ /* t = +1 */ /* */ /* avec pour les parametres {u,v,w} : */ /* */ /* u E [ -p/8 , +p/8 ] */ /* */ /* v E [ -p , +p ] */ /* */ /* w E [ -p/8 , +p/8 ] */ /* */ /* (ou 'p' designe 'pi'). */ /* */ /* */ /*************************************************************************************************************************************/ /* :Fin_listMN_RUBAN_DE_MOBIUS_TRIDIMENSIONNEL_21: */ #include xrs/Mobius3D.22.I" #include xrs/Mobius3D.14.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/volumes.21.I" #define RENORMALISER_ARBITRAIREMENT_LES_DIFFERENTIELLES \ FAUX \ /* Pour ne pas renormaliser arbitrairement les differentielles, et ce afin de prendre en */ \ /* compte leur extrema veritables dans 'RENORMALISATION_D_UNE_LISTE_DE_POINTS'... */ #include xci/sequence.01.I" /* ATTENTION, on definit ainsi le symbole 'DERNIERE_IMAGE' qui ne sert a rien ici, puisque */ /* c'est en effet 'Zmax' qui joue ce role... */ #include xrs/Mobius3D.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 ' : */ /* */ /*************************************************************************************************************************************/ #include xrs/Mobius3D.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 N " R U B A N D E M O B I U S " T R I D I M E N S I O N N E L : */ /* */ /*************************************************************************************************************************************/ BCommande(nombre_d_arguments,arguments) /*-----------------------------------------------------------------------------------------------------------------------------------*/ Bblock #include xrs/volumes.41.I" /* Introduit le 20150528140446... */ #include xrs/volumes.22.I" /*..............................................................................................................................*/ INITIALISATIONS_GENERALES; /* Initialisations generales faites au tout debut... */ iTRANSFORMAT_31(liste_PARAMETRE_T,PARAMETRE_T_IMPLICITE); /* Initialisation du parametre 't'. */ iTRANSFORMAT_31(liste_PARAMETRE_D,PARAMETRE_D_IMPLICITE); /* Initialisation du parametre 'd'. */ #include xrv/champs_5.1A.I" #include xrs/volumes.13.I" GET_ARGUMENTSv(nombre_d_arguments ,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES; GET_FICHIER_PARAMETRE("PARAMETRE_T=",_T); GET_FICHIER_PARAMETRE("PARAMETRE_D=",_D); 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); PROCESS_ARGUMENTS_DE_DERIVATION_PARTIELLE_NUMERIQUE; PROCESS_ARGUMENTS_DE_COLORIAGE; PROCESS_ARGUMENTS_DE_VISUALISATION; PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES; /* Introduit le 20050623090354. */ PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v_w; /* Introduit le 20221221185529... */ PROCESS_ARGUMENTS_DE_PROJECTION_DU_VOLUME; /* Introduit le 20050628090749. */ PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES; PROCESS_ARGUMENTS_DE_GESTION_DE_LA_RELATIVITE_RESTREINTE; ) ); #include xrs/volumes.12.I" /* Visualisation et evolution du volume au cours du temps... */ VISUALISATION_DU_VOLUME_AVEC_DIFFERENCES____(BLOC(Bblock EGAL(parametre_t,sPARAMETRE_T(numero_de_la_periode_courante)); EGAL(parametre_d,sPARAMETRE_D(numero_de_la_periode_courante)); /* Calcul de la valeur des parametres 't' et 'd'. */ Eblock ) ,Fx_Mobius3D_2,Fy_Mobius3D_2,Fz_Mobius3D_2 /* Fonctions definissant parametriquement le volume... */ ,Pxyz_Mobius3D_2 ,BLOC(VIDE;) /* Il n'y a rien de specifique a calculer pour le point courant... */ ); /* On notera que 'VISUALISATION_DU_VOLUME_AVEC_DIFFERENCES_ET_AVEC_TEXTURAGE____(...)' */ /* n'est pas utilise ici car, en effet, 'dimZ' a une valeur importante (en general la */ /* valeur standard 'k___dimZ') qui fait que la capacite memoire est insuffisante pour */ /* accueillir les differents albums ('albumF') de texturage. Au passage, son utilisation */ /* demanderait aussi 'PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_TEXTURE' dans le */ /* 'GET_ARGUMENTSv(...)' ci-dessus (et rien d'autre...). */ RETU_Commande; Eblock ECommande