/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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.11'                                                                                     */
/*                  en est un exemple...                                                                                             */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrs/Mobius3D.11$K' :                                                                                           */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 20040326090514).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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_____AUTORISER_LE_GooF
                                        /* Introduit le 20040424172633 afin de pouvoir generer une image de type 'v $xiirs/MOB3.11'  */
                                        /* sans l'anomalie decrite dans 'v $xiirs/$Fnota 20040330115702'...                          */
@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 20111209172246...                                                            */

#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_11:                                                                          */
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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.11$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). Posons :                                                                                            */
/*                                                                                                                                   */
/*                                                        v            v                                                             */
/*                                      M(u,v,w) = u.cos(---) - w.sin(---) + tM                                                      */
/*                                                        d            d                                                             */
/*                                                                                                                                   */
/*                                                        v            v                                                             */
/*                                      N(u,v,w) = u.sin(---) + w.cos(---) + tN                                                      */
/*                                                        d            d                                                             */
/*                                                                                                                                   */
/*                  avec en general :                                                                                                */
/*                                                                                                                                   */
/*                                      d  = +2                                                                                      */
/*                                      tM = -1                                                                                      */
/*                                      tN = 0                                                                                       */
/*                                                                                                                                   */
/*                  Le point courant {Z1,Z2} est alors defini par :                                                                  */
/*                                                                                                                                   */
/*                                                     i.v                                                                           */
/*                                      Z  = M(u,v,w).e                                                                              */
/*                                       1                                                                                           */
/*                                                                                                                                   */
/*                                                        v                                                                          */
/*                                                     i.---                                                                         */
/*                                                        d                                                                          */
/*                                      Z  = N(u,v,w).e                                                                              */
/*                                       2                                                                                           */
/*                                                                                                                                   */
/*                  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_11:                                                                            */

#include  xrs/Mobius3D.12.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...                                                 */

#define   ACTIVER__S_IL_EST_ACTIVABLE__ET__SI_UTILE__GooF_____COMPATIBILITE_20080124                                                    \
                    VRAI                                                                                                                \
                                        /* Pour activer le 'GooF' (introduit le 20080124183238...).                                  */
#include  xrs/Mobius3D.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 '  :                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrs/Mobius3D.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   " 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 20150528140422...                                                            */

#include  xrs/volumes.22.I"
     /*..............................................................................................................................*/
     INITIALISATIONS_GENERALES;
                                        /* Initialisations generales faites au tout debut...                                         */

     iTRANSFORMAT_31(liste_PARAMETRE_TM,PARAMETRE_TM_IMPLICITE);
     iTRANSFORMAT_31(liste_PARAMETRE_TN,PARAMETRE_TN_IMPLICITE);
                                        /* Initialisation des parametres 'tM' et 'tN'.                                               */

     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_TM=",_TM);
                         GET_FICHIER_PARAMETRE("PARAMETRE_TN=",_TN);

                         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 20050623090323.                                                              */

                         PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v_w;
                                        /* Introduit le 20221221183547...                                                            */

                         PROCESS_ARGUMENTS_DE_VISUALISATION_LIES_A_L_USAGE_DU_GooF;
                                        /* Introduit le 20080125084357.                                                              */

                         PROCESS_ARGUMENTS_DE_PROJECTION_DU_VOLUME;
                                        /* Introduit le 20050628090707.                                                              */

                         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_tM,sPARAMETRE_TM(numero_de_la_periode_courante));
                                                       EGAL(parametre_tN,sPARAMETRE_TN(numero_de_la_periode_courante));
                                                       EGAL(parametre_d,sPARAMETRE_D(numero_de_la_periode_courante));
                                        /* Calcul de la valeur des parametres 'tM', 'tN' et 'd'.                                     */
                                                       Eblock
                                                       )
                                                 ,Fx_Mobius3D_1,Fy_Mobius3D_1,Fz_Mobius3D_1
                                        /* Fonctions definissant parametriquement le volume...                                       */
                                                 ,Pxyz_Mobius3D_1
                                                 ,BLOC(VIDE;)
                                        /* Il n'y a rien de specifique a calculer pour le point courant...                           */
                                                  );
                                        /* Passe le 20040327165609 de 'VISUALISATION_DU_VOLUME_SANS_DIFFERENTIATION(...)' a          */
                                        /* VISUALISATION_DU_VOLUME_AVEC_DIFFERENCES____(...)'.                                       */
                                        /*                                                                                           */
                                        /* 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



Copyright © Jean-François COLONNA, 2019-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2019-2024.