/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D ' U N E   H Y P E R - S P H E R E  :                                                               */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrs/hyper_sphere.11$K' :                                                                                       */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 20191125111450).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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"

#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_HYPER_SPHERE_11:                                                                                             */
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D ' U N E   H Y P E R - S P H E R E  :                                                               */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Definition de la variete ('v $xrs/hyper_sphere.11$K') :                                                                    */
/*                                                                                                                                   */
/*                    L'Hyper-Sphere est definie dans un espace reel                                                                 */
/*                  a 4 dimensions. Elle est parametree a l'aide de                                                                  */
/*                  3 coordonnees {u,v,w}. Le point courant {x1,x2,x3,x4}                                                            */
/*                  est alors defini par :                                                                                           */
/*                                                                                                                                   */
/*                                      x  = A.Xcartesienne_4D(R,u,v,w)                                                              */
/*                                       1                                                                                           */
/*                                                                                                                                   */
/*                                      x  = B.Ycartesienne_4D(R,u,v,w)                                                              */
/*                                       2                                                                                           */
/*                                                                                                                                   */
/*                                      x  = C.Zcartesienne_4D(R,u,v,w)                                                              */
/*                                       3                                                                                           */
/*                                                                                                                                   */
/*                                      x  = D.Tcartesienne_4D(R,u,v,w)                                                              */
/*                                       4                                                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*                  avec pour les parametres {u,v,w} :                                                                               */
/*                                                                                                                                   */
/*                                      u E [ 0 , p ]                                                                                */
/*                                                                                                                                   */
/*                                      v E [ 0 , 2.p ]                                                                              */
/*                                                                                                                                   */
/*                                      w E [ 0 , p ]                                                                                */
/*                                                                                                                                   */
/*                  (ou 'p' designe 'pi').                                                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
/*        :Fin_listMN_HYPER_SPHERE_11:                                                                                               */

#include  xrs/hyper_sphere.12.I"
#include  xrs/hyper_sphere.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   EDITER_LES_MESSAGES_DE_MAUVAIS_ENCADREMENT_DES_NIVEAUX_INTERPOLES                                                             \
                    FAUX                                                                                                                \
                                        /* Introduit le 20191125134930...                                                            */
#include  xrs/hyper_sphere.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/hyper_sphere.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 E   H Y P E R - S P H E R E  :                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
#include  xrs/volumes.41.I"

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

     iTRANSFORMAT_31(liste_PARAMETRE_R,PARAMETRE_R_IMPLICITE);
                                        /* Initialisation du parametre 'R'.                                                          */
     iTRANSFORMAT_31(liste_PARAMETRE_A,PARAMETRE_A_IMPLICITE);
     iTRANSFORMAT_31(liste_PARAMETRE_B,PARAMETRE_B_IMPLICITE);
     iTRANSFORMAT_31(liste_PARAMETRE_C,PARAMETRE_C_IMPLICITE);
     iTRANSFORMAT_31(liste_PARAMETRE_D,PARAMETRE_D_IMPLICITE);
                                        /* Initialisation des parametres {A,B,C,D} introduits le 20191126123726...                   */

#include  xrv/champs_5.1A.I"
#include  xrs/volumes.13.I"

     GET_ARGUMENTSv(nombre_d_arguments
                   ,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES;

                         GET_FICHIER_PARAMETRE("PARAMETRE_R=",_R);

                         GET_FICHIER_PARAMETRE("PARAMETRE_A=",_A);
                         GET_FICHIER_PARAMETRE("PARAMETRE_B=",_B);
                         GET_FICHIER_PARAMETRE("PARAMETRE_C=",_C);
                         GET_FICHIER_PARAMETRE("PARAMETRE_D=",_D);
                                        /* Introduits le 20191126123726...                                                           */

                         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;

                         PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v_w;
                                        /* Introduit le 20221221185642...                                                            */

                         PROCESS_ARGUMENTS_DE_PROJECTION_DU_VOLUME;

                         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_R,sPARAMETRE_R(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'R'.                                                     */
                                                       EGAL(parametre_A,sPARAMETRE_A(numero_de_la_periode_courante));
                                                       EGAL(parametre_B,sPARAMETRE_B(numero_de_la_periode_courante));
                                                       EGAL(parametre_C,sPARAMETRE_C(numero_de_la_periode_courante));
                                                       EGAL(parametre_D,sPARAMETRE_D(numero_de_la_periode_courante));
                                        /* Calcul de la valeur des parametres {A,B,C,D}.                                             */
                                                       Eblock
                                                       )
                                                 ,Fx_HyperSphere_1,Fy_HyperSphere_1,Fz_HyperSphere_1
                                        /* Fonctions definissant parametriquement le volume...                                       */
                                                 ,Pxyz_HyperSphere_1
                                                 ,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.