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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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 20111209173051...                                                            */

#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_31:                                                                          */
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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.31$K') :                                                  */
/*                                                                                                                                   */
/*                    Le ruban de Mobius tridimensionnel est                                                                         */
/*                  defini dans un espace reel a 4 dimensions. 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 {x1,x2,x3,x4} est alors                                                            */
/*                  defini par :                                                                                                     */
/*                                                                                                                                   */
/*                                                      v            v                                                               */
/*                                      x  = R1.[u.cos(---) - w.sin(---)]                                                            */
/*                                       1              d            d                                                               */
/*                                                                                                                                   */
/*                                                      v            v                                                               */
/*                                      x  = R2.[u.sin(---) + w.cos(---)]                                                            */
/*                                       2              d            d                                                               */
/*                                                                                                                                   */
/*                                      x  = R3.sin(v)                                                                               */
/*                                       3                                                                                           */
/*                                                                                                                                   */
/*                                      x  = R4.cos(u)                                                                               */
/*                                       4                                                                                           */
/*                                                                                                                                   */
/*                  avec en general :                                                                                                */
/*                                                                                                                                   */
/*                                      R1 = +2                                                                                      */
/*                                      R2 = +2                                                                                      */
/*                                      R3 = +1                                                                                      */
/*                                      R4 = +1                                                                                      */
/*                                      d  = +4                                                                                      */
/*                                                                                                                                   */
/*                  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_31:                                                                            */

#include  xrs/Mobius3D.32.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.33.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.31.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 20150528140510...                                                            */

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

     iTRANSFORMAT_31(liste_PARAMETRE_R1,PARAMETRE_R1_IMPLICITE);
     iTRANSFORMAT_31(liste_PARAMETRE_R2,PARAMETRE_R2_IMPLICITE);
     iTRANSFORMAT_31(liste_PARAMETRE_R3,PARAMETRE_R3_IMPLICITE);
     iTRANSFORMAT_31(liste_PARAMETRE_R4,PARAMETRE_R4_IMPLICITE);
                                        /* Initialisation des parametres 'R1', 'R2', 'R3' et 'R4'.                                   */

     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_R1=",_R1);
                         GET_FICHIER_PARAMETRE("PARAMETRE_R2=",_R2);
                         GET_FICHIER_PARAMETRE("PARAMETRE_R3=",_R3);
                         GET_FICHIER_PARAMETRE("PARAMETRE_R4=",_R4);

                         GET_FICHIER_PARAMETRE("PARAMETRE_D=",_D);

                         GET_ARGUMENT_F("xx1=",coefficient_de_projection_parallele_01_4D_3D_x1_X);
                         GET_ARGUMENT_F("xx2=",coefficient_de_projection_parallele_01_4D_3D_x2_X);
                         GET_ARGUMENT_F("xx3=",coefficient_de_projection_parallele_01_4D_3D_x3_X);
                         GET_ARGUMENT_F("xx4=",coefficient_de_projection_parallele_01_4D_3D_x4_X);

                         GET_ARGUMENT_F("yx1=",coefficient_de_projection_parallele_01_4D_3D_x1_Y);
                         GET_ARGUMENT_F("yx2=",coefficient_de_projection_parallele_01_4D_3D_x2_Y);
                         GET_ARGUMENT_F("yx3=",coefficient_de_projection_parallele_01_4D_3D_x3_Y);
                         GET_ARGUMENT_F("yx4=",coefficient_de_projection_parallele_01_4D_3D_x4_Y);

                         GET_ARGUMENT_F("zx1=",coefficient_de_projection_parallele_01_4D_3D_x1_Z);
                         GET_ARGUMENT_F("zx2=",coefficient_de_projection_parallele_01_4D_3D_x2_Z);
                         GET_ARGUMENT_F("zx3=",coefficient_de_projection_parallele_01_4D_3D_x3_Z);
                         GET_ARGUMENT_F("zx4=",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 20050623090420.                                                              */

                         PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v_w;
                                        /* Introduit le 20221221185604...                                                            */

                         PROCESS_ARGUMENTS_DE_PROJECTION_DU_VOLUME;
                                        /* Introduit le 20050628090816.                                                              */

                         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_R1,sPARAMETRE_R1(numero_de_la_periode_courante));
                                                       EGAL(parametre_R2,sPARAMETRE_R2(numero_de_la_periode_courante));
                                                       EGAL(parametre_R3,sPARAMETRE_R3(numero_de_la_periode_courante));
                                                       EGAL(parametre_R4,sPARAMETRE_R4(numero_de_la_periode_courante));
                                        /* Calcul de la valeur des parametres 'R1', 'R2', 'R3' et 'R4'.                              */
                                                       EGAL(parametre_d,sPARAMETRE_D(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'd'.                                                     */
                                                       Eblock
                                                       )
                                                 ,Fx_Mobius3D_3,Fy_Mobius3D_3,Fz_Mobius3D_3
                                        /* Fonctions definissant parametriquement le volume...                                       */
                                                 ,Pxyz_Mobius3D_3
                                                 ,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



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.