/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        G E N E R A T I O N   D E   F R A C T A L S   " C A R R E S "   ( V E R S I O N   2 )  :                                   */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrf/FractalCarre.01$K' :                                                                                       */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 20231103173722).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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   ALPHA_X                                                                                                                       \
                    FDU
#define   ALPHA_Y                                                                                                                       \
                    FDU
#define   ALPHA_Z                                                                                                                       \
                    FDU
                                        /* Ponderations des {Xi,Yi,Zi}.                                                              */

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

#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"
#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);
                                        /* Definition en memoire des fichiers de coordonnees et de niveaux.                          */

#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)
                                        /* Acces a un element courant des fichiers de coordonnees et de niveaux.                     */

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

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

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

     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"

     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
                                                  );

                         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_F("alphaX=""ax=""aX=",alpha_X);
                         GET_ARGUMENT_F("alphaY=""ay=""aY=",alpha_Y);
                         GET_ARGUMENT_F("alphaZ=""az=""aZ=",alpha_Z);

                         GET_ARGUMENT_L("rotation=",faire_une_rotation_tridimensionnelle);
                         GET_ARGUMENT_N("transformation=""transformation_quelconque=",faire_une_rotation_tridimensionnelle);

                         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);
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_translationX'...             */
                         GET_ARGUMENT_F("ty=""tY=""TRANSLATION_OY=",translation_OY);
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_translationY'...             */
                         GET_ARGUMENT_F("tz=""tZ=""TRANSLATION_OZ=",translation_OZ);
                                        /* 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;

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_1;

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_5;
                         )
                    );

     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

     begin_nouveau_block
          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
               DEFV(Float,INIT(XfB,ELEMENT_DU_FICHIER_LISTE_X(indexB)));
               DEFV(Float,INIT(YfB,ELEMENT_DU_FICHIER_LISTE_Y(indexB)));
               DEFV(Float,INIT(ZfB,ELEMENT_DU_FICHIER_LISTE_Z(indexB)));
                                        /* Recuperation des coordonnees {X,Y,Z} courantes dans les fichiers.                         */
               DEFV(Float,INIT(niveauB,ELEMENT_DU_FICHIER_LISTE_NIVEAU(indexB)));
                                        /* Recuperation du niveau et ainsi, les points "fils" auront le meme niveau que le point     */
                                        /* "pere"...                                                                                 */

               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.                         */
                                                )
                                          ,USe_GooF(LIZ4(ponderation_OX
                                                        ,TRANSFORMATION_GEOMETRIQUE_3D_Fx(XfB,YfB,ZfB
                                                                                         ,AXPB(alpha_X
                                                                                              ,Xf
                                                                                              ,translation_OX
                                                                                               )
                                                                                          )
                                                        ,ponderation_OY
                                                        ,TRANSFORMATION_GEOMETRIQUE_3D_Fy(XfB,YfB,ZfB
                                                                                         ,AXPB(alpha_Y
                                                                                              ,Yf
                                                                                              ,translation_OY
                                                                                               )
                                                                                          )
                                                        ,ponderation_OZ
                                                        ,TRANSFORMATION_GEOMETRIQUE_3D_Fz(XfB,YfB,ZfB
                                                                                         ,AXPB(alpha_Z
                                                                                              ,Zf
                                                                                              ,translation_OZ
                                                                                               )
                                                                                          )
                                                        ,ponderation_NIVEAU
                                                        ,niveauB
                                                         )
                                                    )
                                          ,EDITER_LA_VALEUR_RESULTANTE_DANS_gOPERATION_SUR_LES_FICHIERS
                                          ,nombre_d_exemplaires_du_resultat_de_l_operation_sur_les_valeurs_courantes
                                           );
               Eblock
          EDoI
          Eblock
     end_nouveau_block

     lGENERATION_D_UN_FICHIER(liste_initiale_des_NIVEAU,Z_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, 2023-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2023-2024.