/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D ' U N   T O R E   E L L I P T I Q U E   E T   O S C I L L A N T  :                                 */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrs/tore.21$K' :                                                                                               */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 20090605180604).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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
#define   PRAGMA_CPP_____SI_LE_GooF_EST_ACTIVABLE_ALORS_FAIRE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_FLOTTANTS
@define   PRAGMA_CL_____MODULE_NON_OPTIMISABLE

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        F I C H I E R S   D ' I N C L U D E S  :                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  INCLUDES_BASE
#include  image_image_IMAGESF_EXT

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

#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.                                                                                        */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrs/tore.18.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 )  :                           */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrk/attractor.12.I"

#include  xrs/tore.19.I"

/*===================================================================================================================================*/
/*        :Debut_listMN_TORE_21:                                                                                                     */
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D ' U N   T O R E  :                                                                                 */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Definition de la surface ('v $xrs/tore.21$K') :                                                                            */
/*                                                                                                                                   */
/*                    Le tore elliptique est defini parametriquement                                                                 */
/*                  en fonction des deux parametres 'u'                                                                              */
/*                  (appele aussi 'distance polaire' ou                                                                              */
/*                  'theta' ou encore 'latitude') et 'v'                                                                             */
/*                  (appele aussi 'longitude' ou 'phi') :                                                                            */
/*                                                                                                                                   */
/*                                      F (u,v) = (R11+R21.cos(u)).cos(v)                                                            */
/*                                       x                                                                                           */
/*                                                                                                                                   */
/*                                      F (u,v) = (R12+R22.cos(u)).sin(v)                                                            */
/*                                       y                                                                                           */
/*                                                                                                                                   */
/*                                      F (u,v) = R23.sin(u) + KO.cos(AO.v+BO)                                                       */
/*                                       z                                                                                           */
/*                                                                                                                                   */
/*                  avec :                                                                                                           */
/*                                                                                                                                   */
/*                                      u E [ 0 , 2.p ]                                                                              */
/*                                                                                                                                   */
/*                                      v E [ 0 , 2.p ]                                                                              */
/*                                                                                                                                   */
/*                  (ou 'p' designe 'pi').                                                                                           */
/*                                                                                                                                   */
/*                    Ainsi, la coordonnee 'u' decrit la                                                                             */
/*                  petite ellipse et 'v' la grande ellipse.                                                                         */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
/*        :Fin_listMN_TORE_21:                                                                                                       */

#include  xrs/tore.22.I"
#include  xrs/tore.14.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/surfaces.21.I"

#include  xrs/tore.23.I"

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D E S   T R O I S   F O N C T I O N S   ' F '  :                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrs/tore.21.I"
                                        /* ATTENTION, il est imperatif que les fonctions derivables formellement soient definies     */
                                        /* dans un fichier a part afin qu'elles ne soient pas l'objet d'une double definition. En    */
                                        /* effet, on trouve dans 'v $xcc/cpp$Z' :                                                    */
                                        /*                                                                                           */
                                        /*                  $CA  $module$w                                                           */
                                        /*                                 | $xcp/cpp$X          ... -c$PASSE_D -e$PASSE_5           */
                                        /*                                                                              > $fichierR  */
                                        /*                                                                                           */
                                        /* Ainsi, si par exemple la definition de 'RAYON(h,theta)' etait faite localement dans le    */
                                        /* fichier '$xrf/recursif.11$K', cette definition apparaitrait deux fois dans le fichier     */
                                        /* '$module$w' : une premiere fois en tete, provenant de '$PASSE_D' ou l'on ne trouve que    */
                                        /* les definitions du type '-define ...', puis une deuxieme fois a la suite, la ou se trouve */
                                        /* concatene au fichier '$module$W' des '-define ...' le fichier '$fichierA' a compiler...   */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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   T O R E  :                                                                                 */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
#include  xrs/surfaces.41.I"
                                        /* Introduit le 20150528140325...                                                            */

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

     iTRANSFORMAT_31(liste_PARAMETRE_R11,PARAMETRE_R11_IMPLICITE);
                                        /* Initialisation du parametre 'R11'.                                                        */
     iTRANSFORMAT_31(liste_PARAMETRE_R12,PARAMETRE_R12_IMPLICITE);
                                        /* Initialisation du parametre 'R12'.                                                        */

     iTRANSFORMAT_31(liste_PARAMETRE_R21,PARAMETRE_R21_IMPLICITE);
                                        /* Initialisation du parametre 'R21'.                                                        */
     iTRANSFORMAT_31(liste_PARAMETRE_R22,PARAMETRE_R22_IMPLICITE);
                                        /* Initialisation du parametre 'R22'.                                                        */
     iTRANSFORMAT_31(liste_PARAMETRE_R23,PARAMETRE_R23_IMPLICITE);
                                        /* Initialisation du parametre 'R23'.                                                        */

     iTRANSFORMAT_31(liste_PARAMETRE_KO,PARAMETRE_KO_IMPLICITE);
                                        /* Initialisation du parametre 'KO' (introduite le 20120605171629).                          */
     iTRANSFORMAT_31(liste_PARAMETRE_AO,PARAMETRE_AO_IMPLICITE);
                                        /* Initialisation du parametre 'AO' (introduite le 20120605171629).                          */
     iTRANSFORMAT_31(liste_PARAMETRE_BO,PARAMETRE_BO_IMPLICITE);
                                        /* Initialisation du parametre 'BO' (introduite le 20120605171629).                          */

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

     GET_ARGUMENTSv(nombre_d_arguments
                   ,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES;

                         GET_FICHIER_PARAMETRE("PARAMETRE_R11=",_R11);
                         GET_FICHIER_PARAMETRE("PARAMETRE_R12=",_R12);

                         GET_FICHIER_PARAMETRE("PARAMETRE_R21=",_R21);
                         GET_FICHIER_PARAMETRE("PARAMETRE_R22=",_R22);
                         GET_FICHIER_PARAMETRE("PARAMETRE_R23=",_R23);

                         GET_FICHIER_PARAMETRE("PARAMETRE_KO=",_KO);
                         GET_FICHIER_PARAMETRE("PARAMETRE_AO=",_AO);
                         GET_FICHIER_PARAMETRE("PARAMETRE_BO=",_BO);
                                        /* Arguments {KO,AO,BO} introduits le 20120605171629...                                      */

                         PROCESS_ARGUMENTS_DE_VISUALISATION;

                         PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;

                         PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;

                         PROCESS_ARGUMENTS_DE_PROJECTION_DE_LA_SURFACE;

                         PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_TEXTURE;

                         PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_NORMALE;

                         PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;

                         PROCESS_ARGUMENTS_DE_DEFINITION_D_UNE_DROITE_DU_PLAN_UV;

                         PROCESS_ARGUMENTS_D_ACTIVATION_DE_L_ARITHMETIQUE_ETENDUE_ET_NON_BASIQUE_POUR_Fx_Fy_Fz;

                         PROCESS_ARGUMENTS_DE_GESTION_DE_LA_RELATIVITE_RESTREINTE;
                         )
                    );

#include  xrs/surfaces.12.I"
                                        /* Visualisation et evolution de la surface au cours du temps...                             */

     VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENTIATION(BLOC(Bblock
                                                           EGAL(parametre_R11,sPARAMETRE_R11(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'R11'.                                                   */
                                                           EGAL(parametre_R12,sPARAMETRE_R12(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'R12'.                                                   */

                                                           EGAL(parametre_R21,sPARAMETRE_R21(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'R21'.                                                   */
                                                           EGAL(parametre_R22,sPARAMETRE_R22(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'R22'.                                                   */
                                                           EGAL(parametre_R23,sPARAMETRE_R23(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'R23'.                                                   */

                                                           EGAL(parametre_KO,sPARAMETRE_KO(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'KO'.                                                    */
                                                           EGAL(parametre_AO,sPARAMETRE_AO(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'AO'.                                                    */
                                                           EGAL(parametre_BO,sPARAMETRE_BO(numero_de_la_periode_courante));
                                        /* Calcul de la valeur du parametre 'BO'.                                                    */
                                                           Eblock
                                                           )
                                                     ,Fx_Tore_2,Fy_Tore_2,Fz_Tore_2
                                        /* Fonctions definissant parametriquement la surface...                                      */
                                                     ,Pxyz_Tore_2
                                                     ,BLOC(VIDE;)
                                        /* Il n'y a rien de specifique a calculer pour le point courant...                           */
                                                      );

     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.