/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        R O T A T I O N   T R I D I M E N S I O N N E L L E   D E   T R O I S   F I C H I E R S   N U M E R I Q U E S              */
/*        E T   T R A N S F O R M A T I O N   A U T O - S I M I L A I R E  :                                                         */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrv/ROTATION.01$K' :                                                                                           */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 1998??????????).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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
                                        /* Ceci a ete introduit le 20070727142541.                                                   */
#define   PRAGMA_CPP_____SI_LE_GooF_EST_ACTIVABLE_ALORS_FAIRE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_FLOTTANTS
                                        /* Ceci a ete introduit le 20070727142541.                                                   */

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        V A L E U R S   I M P L I C I T E S   D E S   P A R A M E T R E S  :                                                       */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   GENERER_UNE_STRUCTURE_AUTOSIMILAIRE                                                                                           \
                    FAUX                                                                                                                \
                                        /* Faut-il generer une structure autosimilaire ('VRAI') ou bien une "simple" transformation  */ \
                                        /* via une certaine matrice ('FAUX'). Ceci fut introduit le 20070719100226 et la valeur par  */ \
                                        /* defaut garantit la compatibilite anterieure...                                            */

#define   ALPHA_X                                                                                                                       \
                    FU
#define   ALPHA_Y                                                                                                                       \
                    FU
#define   ALPHA_Z                                                                                                                       \
                    FU
                                        /* Ponderations des {Xi,Yi,Zi} si 'IL_FAUT(generer_une_structure_autosimilaire)'. Ceci fut   */
                                        /* introduit le 20070720210623...                                                            */

#define   FAIRE_UNE_ROTATION_TRIDIMENSIONNELLE                                                                                          \
                    VRAI                                                                                                                \
                                        /* Faut-il faire une rotation tridimensionnelle ('VRAI') ou bien une transformation          */ \
                                        /* arbitraire ('FAUX') definie par les 3x3 elements de sa matrice, ce qui ouvre la porte,    */ \
                                        /* par exemple, a des 'IFS' tridimensionnelles...                                            */

#define   PONDERATION_OX_IMPLICITE                                                                                                      \
                    FZERO
#define   PONDERATION_OY_IMPLICITE                                                                                                      \
                    FZERO
#define   PONDERATION_OZ_IMPLICITE                                                                                                      \
                    FZERO
                                        /* Ponderations de selections des coordonnees apres rotation.                                */
#define   PONDERATION_NIVEAU_IMPLICITE                                                                                                  \
                    FZERO
                                        /* Ponderation de selection des niveaux (introduit le 20231102174138).                       */

#define   ROTATION_OX_IMPLICITE                                                                                                         \
                    FZERO
#define   ROTATION_OY_IMPLICITE                                                                                                         \
                    FZERO
#define   ROTATION_OZ_IMPLICITE                                                                                                         \
                    FZERO
                                        /* Angles de rotation.                                                                       */
#define   TRANSLATION_OX_IMPLICITE                                                                                                      \
                    FZERO
#define   TRANSLATION_OY_IMPLICITE                                                                                                      \
                    FZERO
#define   TRANSLATION_OZ_IMPLICITE                                                                                                      \
                    FZERO
                                        /* Translation de l'origine.                                                                 */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D E S   T R O I S   F I C H I E R S  :                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrv/ARITHMET.1d.I"
                                        /* Passage a l'allocation dynamique le 20060214181321...                                     */
#include  xrv/ARITHMET.21.I"
#include  xrv/champs_5.41.I"

#define   X_IMPLICITE                                                                                                                   \
                    FZERO
#define   Y_IMPLICITE                                                                                                                   \
                    FZERO
#define   Z_IMPLICITE                                                                                                                   \
                    FZERO

#define   NIVEAU_IMPLICITE                                                                                                              \
                    FLOT__BLANC

#define   SELECTION_DES_POINTS_D_ANCRAGE_IMPLICITE                                                                                      \
                    FLOT(VRAI)
#define   SELECTION_DES_POINTS_A_ANCRER_IMPLICITE                                                                                       \
                    FLOT(VRAI)
#define   DUPLIQUER_LA_VALEUR_IMPLICITE                                                                                                 \
                    FLOT(FAUX)

gGENERATION_D_UN_FICHIER(fichier_LISTE_X,liste_initiale_des_X);
gGENERATION_D_UN_FICHIER(fichier_LISTE_Y,liste_initiale_des_Y);
gGENERATION_D_UN_FICHIER(fichier_LISTE_Z,liste_initiale_des_Z);

gGENERATION_D_UN_FICHIER(fichier_LISTE_NIVEAU,liste_initiale_des_NIVEAU);

gGENERATION_D_UN_FICHIER(fichier_des_selections_des_points_d_ancrage,liste_des_selections_des_points_d_ancrage);
gGENERATION_D_UN_FICHIER(fichier_des_selections_des_points_a_ancrer,liste_des_selections_des_points_a_ancrer);
gGENERATION_D_UN_FICHIER(fichier_des_duplications_des_valeurs,liste_des_duplications_des_valeurs);
                                        /* Definition en memoire des fichiers de coordonnees, ainsi que du fichier de                */
                                        /* selection des points d'ancrage introduit le 20070719100226. Le fichier de                 */
                                        /* selection des points a ancrer fut introduit le 20070721104407. Le fichier de              */
                                        /* duplication des valeurs fut introduit le 20070721181021...                                */

#define   ELEMENT_DU_FICHIER_LISTE_X(index)                                                                                             \
                    gELEMENT_DU_FICHIER(liste_initiale_des_X,index)
#define   ELEMENT_DU_FICHIER_LISTE_Y(index)                                                                                             \
                    gELEMENT_DU_FICHIER(liste_initiale_des_Y,index)
#define   ELEMENT_DU_FICHIER_LISTE_Z(index)                                                                                             \
                    gELEMENT_DU_FICHIER(liste_initiale_des_Z,index)

#define   ELEMENT_DU_FICHIER_LISTE_NIVEAU(index)                                                                                        \
                    gELEMENT_DU_FICHIER(liste_initiale_des_NIVEAU,index)

#define   ELEMENT_DE_SELECTION_DES_POINTS_D_ANCRAGE(index)                                                                              \
                    LOGI(gELEMENT_DU_FICHIER(liste_des_selections_des_points_d_ancrage,index))
#define   ELEMENT_DE_SELECTION_DES_POINTS_A_ANCRER(index)                                                                               \
                    LOGI(gELEMENT_DU_FICHIER(liste_des_selections_des_points_a_ancrer,index))
#define   ELEMENT_DE_SELECTION_DE_DUPLICATION_DES_VALEURS(index)                                                                        \
                    LOGI(gELEMENT_DU_FICHIER(liste_des_duplications_des_valeurs,index))
                                        /* Acces a un element courant des fichiers de coordonnees, ainsi que du fichier de           */
                                        /* selection des points d'ancrage introduit le 20070719100226. Le fichier de                 */
                                        /* selection des points a ancrer fut introduit le 20070721104407. Le fichier de              */
                                        /* duplication des valeurs fut introduit le 20070721181021...                                */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        R O T A T I O N   T R I D I M E N S I O N N E L L E   D E   T R O I S   F I C H I E R S   N U M E R I Q U E S              */
/*        E T   T R A N S F O R M A T I O N   A U T O - S I M I L A I R E  :                                                         */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
#include  xrv/ARITHMET.22.I"
#include  xci/valeurs.03.I"

     DEFV(Logical,INIT(generer_une_structure_autosimilaire,GENERER_UNE_STRUCTURE_AUTOSIMILAIRE));
                                        /* Faut-il generer une structure autosimilaire ('VRAI') ou bien une "simple" transformation  */
                                        /* via une certaine matrice ('FAUX'). Ceci fut introduit le 20070719100226 et la valeur par  */
                                        /* defaut garantit la compatibilite anterieure...                                            */

     DEFV(Float,INIT(alpha_X,ALPHA_X));
     DEFV(Float,INIT(alpha_Y,ALPHA_Y));
     DEFV(Float,INIT(alpha_Z,ALPHA_Z));
                                        /* Ponderations des {Xi,Yi,Zi} si 'IL_FAUT(generer_une_structure_autosimilaire)'. Ceci fut   */
                                        /* introduit le 20070720210623...                                                            */

     DEFV(Logical,INIT(faire_une_rotation_tridimensionnelle,FAIRE_UNE_ROTATION_TRIDIMENSIONNELLE));
                                        /* Faut-il faire une rotation tridimensionnelle ('VRAI') ou bien une transformation          */
                                        /* arbitraire ('FAUX') definie par les 3x3 elements de sa matrice, ce qui ouvre la porte,    */
                                        /* par exemple, a des 'IFS' tridimensionnelles...                                            */

     DEFV(Float,INIT(ponderation_OX,PONDERATION_OX_IMPLICITE));
     DEFV(Float,INIT(ponderation_OY,PONDERATION_OY_IMPLICITE));
     DEFV(Float,INIT(ponderation_OZ,PONDERATION_OZ_IMPLICITE));
                                        /* Ponderations de selections des coordonnees apres rotation.                                */
     DEFV(Float,INIT(ponderation_NIVEAU,PONDERATION_NIVEAU_IMPLICITE));
                                        /* Ponderation de selection des niveaux (introduit le 20231102174138).                       */

     DEFV(Float,INIT(angle_de_ROTATION_OX,ROTATION_OX_IMPLICITE));
     DEFV(Float,INIT(angle_de_ROTATION_OY,ROTATION_OY_IMPLICITE));
     DEFV(Float,INIT(angle_de_ROTATION_OZ,ROTATION_OZ_IMPLICITE));
                                        /* Angles de rotation.                                                                       */
     DEFV(Float,INIT(translation_OX,TRANSLATION_OX_IMPLICITE));
     DEFV(Float,INIT(translation_OY,TRANSLATION_OY_IMPLICITE));
     DEFV(Float,INIT(translation_OZ,TRANSLATION_OZ_IMPLICITE));
                                        /* Translation de l'origine.                                                                 */
     /*..............................................................................................................................*/
#include  xrv/champs_5.1A.I"
                                        /* Ceci fut introduit le 20070103172621...                                                   */

     INITIALISATION_TRANSFORMATION;
                                        /* Mise en place de la transformation neutre...                                              */

     GET_ARGUMENTS_(nombre_d_arguments
                   ,BLOC(PROCESS_ARGUMENT_I("nombre_elements=""ne=",nombre_d_elements
                                           ,BLOC(VIDE;)
                                           ,BLOC(Bblock
                                                 PRINT_AVERTISSEMENT("'ne=' doit etre defini avant toute entree de fichiers");
                                                 Eblock
                                                 )
                                            );

                         PROCESS_ARGUMENTS_DE_DEFINITION_DES_FICHIERS_01;

                         PROKESF_ARGUMENT_FICHIER("LISTE_X="
                                                 ,fichier_LISTE_X
                                                 ,liste_initiale_des_X
                                                 ,X_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                         PROKESF_ARGUMENT_FICHIER("LISTE_Y="
                                                 ,fichier_LISTE_Y
                                                 ,liste_initiale_des_Y
                                                 ,Y_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                         PROKESF_ARGUMENT_FICHIER("LISTE_Z="
                                                 ,fichier_LISTE_Z
                                                 ,liste_initiale_des_Z
                                                 ,Z_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );

                         PROKESF_ARGUMENT_FICHIER("LISTE_NIVEAU="
                                                 ,fichier_LISTE_NIVEAU
                                                 ,liste_initiale_des_NIVEAU
                                                 ,NIVEAU_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                                        /* Parametre introduit le 20231102174138...                                                  */

                         PROKESF_ARGUMENT_FICHIER("SELECTION_ANCRAGE=""ANCRAGE="
                                                 ,fichier_des_selections_des_points_d_ancrage
                                                 ,liste_des_selections_des_points_d_ancrage
                                                 ,SELECTION_DES_POINTS_D_ANCRAGE_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                                        /* Parametre introduit le 20070719100226...                                                  */
                         PROKESF_ARGUMENT_FICHIER("SELECTION_A_ANCRER=""A_ANCRER="
                                                 ,fichier_des_selections_des_points_a_ancrer
                                                 ,liste_des_selections_des_points_a_ancrer
                                                 ,SELECTION_DES_POINTS_A_ANCRER_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                                        /* Parametre introduit le 20070719100226...                                                  */

                         PROKESF_ARGUMENT_FICHIER("DUPLIQUER_VALEUR=""DUPLIQUER="
                                                 ,fichier_des_duplications_des_valeurs
                                                 ,liste_des_duplications_des_valeurs
                                                 ,DUPLIQUER_LA_VALEUR_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                                        /* Parametre introduit le 20070721181021...                                                  */

                         GET_ARGUMENT_F("px=""pX=""Pond1=",ponderation_OX);
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_pasX'...                     */
                         GET_ARGUMENT_F("py=""pY=""Pond2=",ponderation_OY);
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_pasY'...                     */
                         GET_ARGUMENT_F("pz=""pZ=""Pond3=",ponderation_OZ);
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_pasZ'...                     */
                         GET_ARGUMENT_F("pniveau=""pn=""pNIVEAU=""pN=""Pond4=",ponderation_NIVEAU);

                         GET_ARGUMENT_L("transformation_autosimilaire=""structure_autosimilaire=""self=""auto="
                                       ,generer_une_structure_autosimilaire
                                        );
                                        /* Parametre introduit le 20070719100226...                                                  */
                         GET_ARGUMENT_N("transformation_simple=""structure_simple=""simple=",generer_une_structure_autosimilaire);
                                        /* Parametre introduit le 20070722092531...                                                  */

                         GET_ARGUMENT_F("alphaX=""ax=""aX=",alpha_X);
                         GET_ARGUMENT_F("alphaY=""ay=""aY=",alpha_Y);
                         GET_ARGUMENT_F("alphaZ=""az=""aZ=",alpha_Z);
                                        /* Parametres introduits le 20070720210623...                                                */

                         GET_ARGUMENT_L("rotation=",faire_une_rotation_tridimensionnelle);
                         GET_ARGUMENT_N("transformation=""transformation_quelconque=",faire_une_rotation_tridimensionnelle);
                                        /* Parametre introduit le 20070719164210...                                                  */

                         PROCESS_ARGUMENT_F("ROTATION_OX=",angle_de_ROTATION_OX
                                           ,BLOC(VIDE;)
                                           ,BLOC(EGAL(vecteurs_____ordre_de_la_ROTATION_0X,nombre_d_arguments_recuperes);
                                                 )
                                            );
                         PROCESS_ARGUMENT_F("ROTATION_OY=",angle_de_ROTATION_OY
                                           ,BLOC(VIDE;)
                                           ,BLOC(EGAL(vecteurs_____ordre_de_la_ROTATION_0Y,nombre_d_arguments_recuperes);
                                                 )
                                            );
                         PROCESS_ARGUMENT_F("ROTATION_OZ=",angle_de_ROTATION_OZ
                                           ,BLOC(VIDE;)
                                           ,BLOC(EGAL(vecteurs_____ordre_de_la_ROTATION_0Z,nombre_d_arguments_recuperes);
                                                 )
                                            );

                         GET_ARGUMENT_F("a11=",ASD2(vecteurs_____matrix_3D,cx,cx));
                         GET_ARGUMENT_F("a12=",ASD2(vecteurs_____matrix_3D,cx,cy));
                         GET_ARGUMENT_F("a13=",ASD2(vecteurs_____matrix_3D,cx,cz));

                         GET_ARGUMENT_F("a21=",ASD2(vecteurs_____matrix_3D,cy,cx));
                         GET_ARGUMENT_F("a22=",ASD2(vecteurs_____matrix_3D,cy,cy));
                         GET_ARGUMENT_F("a23=",ASD2(vecteurs_____matrix_3D,cy,cz));

                         GET_ARGUMENT_F("a31=",ASD2(vecteurs_____matrix_3D,cz,cx));
                         GET_ARGUMENT_F("a32=",ASD2(vecteurs_____matrix_3D,cz,cy));
                         GET_ARGUMENT_F("a33=",ASD2(vecteurs_____matrix_3D,cz,cz));

                         GET_ARGUMENT_F("tx=""tX=""TRANSLATION_OX=",translation_OX);
                                        /* Introduit le 20060516104923...                                                            */
                                        /*                                                                                           */
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_translationX'...             */
                         GET_ARGUMENT_F("ty=""tY=""TRANSLATION_OY=",translation_OY);
                                        /* Introduit le 20060516104923...                                                            */
                                        /*                                                                                           */
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_translationY'...             */
                         GET_ARGUMENT_F("tz=""tZ=""TRANSLATION_OZ=",translation_OZ);
                                        /* Introduit le 20060516104923...                                                            */
                                        /*                                                                                           */
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_translationZ'...             */

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_3;
                                        /* Cette procedure fut introduite le 20070103172621...                                       */

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_1;
                                        /* Cette procedure fut introduite le 20061226192152...                                       */

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_5;
                                        /* Cette procedure fut introduite le 20211005104829...                                       */
                         )
                    );

     Test(IL_FAUT(faire_une_rotation_tridimensionnelle))
          Bblock
          GENERATION_DE_LA_MATRICE_DE_ROTATION(vecteurs_____ordre_de_la_ROTATION_0X
                                              ,T_ROTATION_X(angle_de_ROTATION_OX);
                                              ,vecteurs_____ordre_de_la_ROTATION_0Y
                                              ,T_ROTATION_Y(angle_de_ROTATION_OY);
                                              ,vecteurs_____ordre_de_la_ROTATION_0Z
                                              ,T_ROTATION_Z(angle_de_ROTATION_OZ);
                                               );
                                        /* Mise en place du produit des trois rotations relatives aux trois axes.                    */
          Eblock
     ATes
          Bblock
                                        /* Dans le cas ou la transformation est arbitraire, les coefficients de la matrice sont      */
                                        /* supposes corrects apres leur entree (totale ou partielle...).                             */
          Eblock
     ETes

     Test(IL_NE_FAUT_PAS(generer_une_structure_autosimilaire))
          Bblock
          Test(IFNE_chaine(fichier_des_selections_des_points_d_ancrage,NOM_PIPE_Local))
                                        /* Test introduit le 20070719132807...                                                       */
               Bblock
               PRINT_ATTENTION("la selection des points d'ancrage n'est utilisee que pour la generation de structures autosimilaires");
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          Test(IFNE_chaine(fichier_des_selections_des_points_a_ancrer,NOM_PIPE_Local))
                                        /* Test introduit le 20070721104407...                                                       */
               Bblock
               PRINT_ATTENTION("la selection des points a ancrer n'est utilisee que pour la generation de structures autosimilaires");
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          Test(IFNE_chaine(fichier_des_duplications_des_valeurs,NOM_PIPE_Local))
                                        /* Test introduit le 20070721181021...                                                       */
               Bblock
               PRINT_ATTENTION("la duplication des valeurs n'est utilisee que pour la generation de structures autosimilaires");
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          gOPERATION_SUR_LES_FICHIERS(BLOC(
                                           DEFV(Float,INIT(Xf,ELEMENT_DU_FICHIER_LISTE_X(index)));
                                           DEFV(Float,INIT(Yf,ELEMENT_DU_FICHIER_LISTE_Y(index)));
                                           DEFV(Float,INIT(Zf,ELEMENT_DU_FICHIER_LISTE_Z(index)));
                                        /* Recuperation des coordonnees {X,Y,Z} courantes dans les fichiers.                         */
                                           DEFV(Float,INIT(niveau,ELEMENT_DU_FICHIER_LISTE_NIVEAU(index)));
                                        /* Recuperation du niveau (introduit le 20231102174138).                                     */
                                           )
                                     ,USe_GooF(LIZ4(ponderation_OX,TRANSFORMATION_GEOMETRIQUE_3D_Fx(Xf,Yf,Zf,translation_OX)
                                                   ,ponderation_OY,TRANSFORMATION_GEOMETRIQUE_3D_Fy(Xf,Yf,Zf,translation_OY)
                                                   ,ponderation_OZ,TRANSFORMATION_GEOMETRIQUE_3D_Fz(Xf,Yf,Zf,translation_OZ)
                                                   ,ponderation_NIVEAU
                                                   ,niveau
                                                    )
                                               )
                                     ,EDITER_LA_VALEUR_RESULTANTE_DANS_gOPERATION_SUR_LES_FICHIERS
                                     ,nombre_d_exemplaires_du_resultat_de_l_operation_sur_les_valeurs_courantes
                                      );
                                        /* Rotation tridimensionnelle.                                                               */
          Eblock
     ATes
          Bblock
          DEFV(Int,INIT(indexB,UNDEF));
                                        /* Index de la boucle "de Base"...                                                           */

          DoIn(indexB
              ,PREMIER_ELEMENT_D_UN_FICHIER
              ,DERNIER_ELEMENT_D_UN_FICHIER
              ,I
               )
               Bblock
               Test(EST_VRAI(ELEMENT_DE_SELECTION_DES_POINTS_D_ANCRAGE(indexB)))
                    Bblock
                    gOPERATION_SUR_LES_FICHIERS(BLOC(
                                                     DEFV(Float,INIT(Xf,ELEMENT_DU_FICHIER_LISTE_X(index)));
                                                     DEFV(Float,INIT(Yf,ELEMENT_DU_FICHIER_LISTE_Y(index)));
                                                     DEFV(Float,INIT(Zf,ELEMENT_DU_FICHIER_LISTE_Z(index)));
                                        /* Recuperation des coordonnees {X,Y,Z} courantes dans les fichiers.                         */
                                                     DEFV(Float,INIT(niveau,ELEMENT_DU_FICHIER_LISTE_NIVEAU(index)));
                                        /* Recuperation du niveau (introduit le 20231102174138) et ainsi, les points "fils"          */
                                        /* auront le meme niveau que le point "pere"...                                              */

                                                     DEFV(Logical,INIT(le_point_courant_est_a_ancrer
                                                                      ,ELEMENT_DE_SELECTION_DES_POINTS_A_ANCRER(index)
                                                                       )
                                                          );
                                        /* Le point courant est-il a ancrer ?                                                        */
                                                     DEFV(Logical,INIT(le_point_courant_est_a_dupliquer
                                                                      ,ELEMENT_DE_SELECTION_DE_DUPLICATION_DES_VALEURS(index)
                                                                       )
                                                          );
                                        /* Le point courant est-il a dupliquer ?                                                     */
                                                     )
                                               ,USe_GooF(LIZ4(ponderation_OX
                                                             ,TRANSFORMATION_GEOMETRIQUE_3D_Fx(Xf,Yf,Zf
                                                                                              ,AXPB(alpha_X
                                                                                                   ,ELEMENT_DU_FICHIER_LISTE_X(indexB)
                                                                                                   ,translation_OX
                                                                                                    )
                                                                                               )
                                                             ,ponderation_OY
                                                             ,TRANSFORMATION_GEOMETRIQUE_3D_Fy(Xf,Yf,Zf
                                                                                              ,AXPB(alpha_Y
                                                                                                   ,ELEMENT_DU_FICHIER_LISTE_Y(indexB)
                                                                                                   ,translation_OY
                                                                                                    )
                                                                                               )
                                                             ,ponderation_OZ
                                                             ,TRANSFORMATION_GEOMETRIQUE_3D_Fz(Xf,Yf,Zf
                                                                                              ,AXPB(alpha_Z
                                                                                                   ,ELEMENT_DU_FICHIER_LISTE_Z(indexB)
                                                                                                   ,translation_OZ
                                                                                                    )
                                                                                               )
                                                             ,ponderation_NIVEAU
                                                             ,niveau
                                                              )
                                                         )
                                               ,EST_VRAI(le_point_courant_est_a_ancrer)
                                               ,ADD2(nombre_d_exemplaires_du_resultat_de_l_operation_sur_les_valeurs_courantes
                                                    ,COND(EST_VRAI(le_point_courant_est_a_dupliquer)
                                                         ,UN
                                                         ,ZERO
                                                          )
                                                     )
                                                );
                                        /* Transformation autosimilaire (introduite le 20070719100226...).                           */
                                        /*                                                                                           */
                                        /* On notera donc que par defaut le nombre d'exemplaires edites par defaut est de 1 aux      */
                                        /* extremites du segment de 'index' et de 2 a l'interieur. Ceci est en prevision d'une       */
                                        /* exploitation de ces sorties par 'v $xrv/particule.10$K' avec les parametres :             */
                                        /*                                                                                           */
                                        /*                  isoles=FAUX                                                              */
                                        /*                  chainer=FAUX                                                             */
                                        /*                                                                                           */
                                        /* qui permettent de visualiser des "segments" formes de couples de points...                */
                                        /*                                                                                           */
                                        /* Le principe est le suivant : soit la liste de points {Pi}. Pour chaque point 'Pi', on     */
                                        /* genere ici la liste {Pij} (possedant le meme nombre 'N' d'elements que la liste {Pi})     */
                                        /* ou points 'Pij' sont definis par :                                                        */
                                        /*                                                                                           */
                                        /*                  X   = a  .X  + F (X ,Y ,Z ) + T                                          */
                                        /*                   ij    x   i    x  j  j  j     x                                         */
                                        /*                                                                                           */
                                        /*                  Y   = a  .Y  + F (X ,Y ,Z ) + T                                          */
                                        /*                   ij    y   i    y  j  j  j     y                                         */
                                        /*                                                                                           */
                                        /*                  Z   = a  .Z  + F (X ,Y ,Z ) + T                                          */
                                        /*                   ij    z   i    z  j  j  j     z                                         */
                                        /*                                                                                           */
                                        /* pour :                                                                                    */
                                        /*                                                                                           */
                                        /*                  j E [1,n]                                                                */
                                        /*                                                                                           */
                                        /* ou 'N' designe le nombre d'elements de la liste {Pi} et ou {ax,ay,az} sont les            */
                                        /* {alpha_X,alpha_Y,alpha_Z} introduits le 20070720210623 et qui permettent de traiter       */
                                        /* ainsi autre chose que des coordonnees et par exemple des derivees ou encore des rayons... */
                                        /*                                                                                           */
                                        /* La liste de "selection des points d'ancrage" permet d'ignorer certains points 'Pi'...     */
                                        /* La liste de "selection des points a ancrer" permet d'ignorer certains points 'Pj'...      */
                                        /*                                                                                           */
                                        /* Dans le cas ou {Fx,Fy,Fz} definit une application contractante, par exemple :             */
                                        /*                                                                                           */
                                        /*                  rotation=FAUX                                                            */
                                        /*                                                                                           */
                                        /*                  a11=1/2                                                                  */
                                        /*                  a22=1/2                                                                  */
                                        /*                  a33=0                                                                    */
                                        /*                                                                                           */
                                        /*                  Tx=-X2/(2x2)                                                             */
                                        /*                  Ty=0                                                                     */
                                        /*                  Tz=0                                                                     */
                                        /*                                                                                           */
                                        /* la liste {Pij} pour 'i' donne est une version reduite (d'un rapport 1/2) de la liste      */
                                        /* {Pi} qui est "accrochee" au point 'Pi'. Pour que les points {Pij} soient, par exemple,    */
                                        /* "centres" sur chaque point 'Pi', il faut utiliser a bon escient les translations          */
                                        /* {Tx,Ty,Tz}. Soit par exemple la liste {Pi} suivante :                                     */
                                        /*                                                                                           */
                                        /*                  P1 _______________ P3                                                    */
                                        /*                     \      |      /                                                       */
                                        /*                       \    |    /                                                         */
                                        /*                         \  |  /                                                           */
                                        /*                           \|/                                                             */
                                        /*                                                                                           */
                                        /*                            P2                                                             */
                                        /*                                                                                           */
                                        /* en inhibant (via la liste de "selection des points d'ancrage") le point 'P2', cela        */
                                        /* donnera avec la matrice '{Aij}' ci dessus :                                               */
                                        /*                                                                                           */
                                        /*                  P11 _______ P13    P31 _______ P33                                       */
                                        /*                      \  |  /            \  |  /                                           */
                                        /*                        \|/                \|/                                             */
                                        /*                                                                                           */
                                        /*                        P12                P32                                             */
                                        /*                                                                                           */
                                        /*                        P1 _______________ P3                                              */
                                        /*                           \      |      /                                                 */
                                        /*                             \    |    /                                                   */
                                        /*                               \  |  /                                                     */
                                        /*                                 \|/                                                       */
                                        /*                                                                                           */
                                        /*                                  P2                                                       */
                                        /*                                                                                           */
                                        /* cette liste {Pij} (pour tout 'i' et tout 'j') etant autosimilaire. Elle sera editee       */
                                        /* selon :                                                                                   */
                                        /*                                                                                           */
                                        /*                  P11,P12,P12,P13                                                          */
                                        /*                  P31,P32,P32,P33                                                          */
                                        /*                                                                                           */
                                        /* avec la demande de duplication de 'P12' et de 'P32', et inhibition du point 'P2' en       */
                                        /* ne figurant pas dans la liste de "selection des points d'ancrage"...                      */
                    Eblock
               ATes
                    Bblock
                    Eblock
               ETes
               Eblock
          EDoI
          Eblock
     ETes

     lGENERATION_D_UN_FICHIER(liste_des_duplications_des_valeurs,DUPLIQUER_LA_VALEUR_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_des_selections_des_points_a_ancrer,SELECTION_DES_POINTS_A_ANCRER_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_des_selections_des_points_d_ancrage,SELECTION_DES_POINTS_D_ANCRAGE_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_initiale_des_Z,Z_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_initiale_des_Y,Y_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_initiale_des_X,X_IMPLICITE);

     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.