/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   G E N E R A L   D E   L ' E N S E M B L E   D E   M A N D E L B R O T                                        */
/*        A V E C   M I S E   E N   M O N T A G N E S   D E   C H A Q U E   I M A G E                                                */
/*        D A N S   L E   P L A N   C O M P L E X E  :                                                                               */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrc/mandel.02$K' :                                                                                             */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 1988??????????).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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   GENERER_LES_GET_ARGUMENTS_DE_CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_FLOTTANTS
                                        /* Introduit le 20070129111928...                                                            */

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        P A R A M E T R E S  :                                                                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   X_DEPART_GAUCHE                                                                                                               \
                    PARE(-2.0)
#define   X_DEPART_DROITE                                                                                                               \
                    PARE(0.5)
#define   Y_DEPART_BAS                                                                                                                  \
                    PARE(-1.25)
#define   Y_DEPART_HAUT                                                                                                                 \
                    PARE(1.25)
                                        /* Definition de la fenetre de depart (vue d'ensemble).                                      */
#define   X_ARRIVEE_GAUCHE                                                                                                              \
                    PARE(-0.234254)
#define   X_ARRIVEE_DROITE                                                                                                              \
                    PARE(-0.234254)
#define   Y_ARRIVEE_BAS                                                                                                                 \
                    PARE(-0.828110)
#define   Y_ARRIVEE_HAUT                                                                                                                \
                    PARE(-0.828110)
                                        /* Definition de la fenetre d'arrivee (dans les filaments).                                  */
#define   OMBRES_PORTEES                                                                                                                \
                    VRAI
#define   LARGEUR_ZONE_PENOMBRE                                                                                                         \
                    PARE(12.0)
                                        /* Pour                                                                                      */
#define   DEPTH_CUEING                                                                                                                  \
                    VRAI
#define   MIN_DEPTH_CUEING                                                                                                              \
                    FZERO
                                        /*      generer                                                                              */
#define   AVION     FAUX                                                                                                                \
                                        /*              la surface                                                                   */
#define   ALIAS     VRAI                                                                                                                \
                                        /*                         montagneuse.                                                      */
#define   TRANS_OX  0.0
#define   TRANS_OY  0.80000
                                        /* Translation de la montagne.                                                               */
#define   X_SOURCE  6.4453
#define   Y_SOURCE  2.7343                                                                                                              \
                                        /* Position de la source lumineuse.                                                          */
#define   MODULATION_DU_CHAMP                                                                                                           \
                    FU                                                                                                                  \
                                        /* Modulation du champ de Mandelbrot.                                                        */
#define   NOMBRE_D_ITERATIONS                                                                                                           \
                    SOUS(BLANC,NOIR)                                                                                                    \
                                        /* Nombre d'iterations demandees.                                                            */
#define   PRECISION                                                                                                                     \
                    PARE(0.00001)                                                                                                       \
                                        /* Pour calculer le rapport de reduction.                                                    */
#define   VISUALISATION_ARGUMENTS_DE_SORTIE                                                                                             \
                    FAUX                                                                                                                \
                                        /* Faut-il visualiser l'argument de sortie ('VRAI') ou les equipotentielles ('FAUX').        */

#include  xci/sequence.01.I"

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        M A C R O S   U T I L E S  :                                                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrc/ITERATION.11.I"
                                        /* Introduit le 20091122182507...                                                            */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   G E N E R A L   D E   L ' E N S E M B L E   D E   M A N D E L B R O T                                        */
/*        A V E C   M I S E   E N   M O N T A G N E S   D E   C H A Q U E   I M A G E                                                */
/*        D A N S   L E   P L A N   C O M P L E X E  :                                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
     DEFV(complexe,coin_bas_gauche_de_depart_puis_courant);
                                        /* Definition du point situe en bas a gauche de la fenetre de depart, puis de la fenetre     */
                                        /* courante au cours du zoom.                                                                */
     DEFV(complexe,coin_haut_droite_de_depart_puis_courant);
                                        /* Definition du point situe en haut a droite de la fenetre de depart, puis de la fenetre    */
                                        /* courante au cours du zoom.                                                                */
     DEFV(complexe,coin_bas_gauche_d_arrivee);
                                        /* Definition du point situe en bas a gauche de la fenetre d'arrivee,                        */
     DEFV(complexe,coin_haut_droite_d_arrivee);
                                        /* Definition du point situe en haut a droite de la fenetre d'arrivee.                       */
     DEFV(Positive,INIT(nombre_d_iterations,NOMBRE_D_ITERATIONS));
                                        /* Nombre maximal d'iterations a effectuer.                                                  */
     DEFV(Float,INIT(precision,PRECISION));
                                        /* Pour calculer le rapport de reduction.                                                    */
     DEFV(Float,INIT(rapport_de_reduction,FLOT__UNDEF));
                                        /* Rapport de passage d'une fenetre a l'autre.                                               */
     DEFV(Logical,INIT(visualisation_arguments_de_sortie,VISUALISATION_ARGUMENTS_DE_SORTIE));
                                        /* Indicateur de choix entre la visualisation du nombre d'iterations ('FAUX'),               */
                                        /* et l'argument de sortie ('VRAI').                                                         */
     DEFV(Logical,INIT(vue_d_avion,AVION));
                                        /* Pour savoir si l'on genere une vue de cote, ou bien une vue d'avion...                    */
     DEFV(Logical,INIT(ombres_portees,OMBRES_PORTEES));
     DEFV(Float,INIT(largeur_zone_penombre,LARGEUR_ZONE_PENOMBRE));
                                        /* Parametres de l'ombrage.                                                                  */
     DEFV(Logical,INIT(depth_cueing,DEPTH_CUEING));
     DEFV(Float,INIT(min_depth_cueing,MIN_DEPTH_CUEING));
                                        /* Parametres du depth-cueing.                                                               */
     DEFV(deltaF_2D,translation_montagne);
                                        /* Translation de la surface montagneuse.                                                    */
     DEFV(Float,INIT(modulation_du_champ,MODULATION_DU_CHAMP));
                                        /* Modulation du champ de Mandelbrot.                                                        */
     DEFV(pointF_2D,source_lumineuse);
                                        /* Position de la source lumineuse.                                                          */
     DEFV(Int,INIT(nombre_d_images,NOMBRE_D_IMAGES));
                                        /* Nombre d'images a generer.                                                                */
     DEFV(Int,INIT(numero_d_image,PREMIERE_IMAGE_ABSOLUE));
                                        /* Numero de l'image courante (celle-ci n'est pas necessairement generee : voir              */
                                        /* le parametre 'PREMIERE_IMAGE').                                                           */
     DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE));
                                        /* Nom de la sequence d'images a generer.                                                    */
     DEFV(Int,INIT(nombre_de_chiffres,NOMBRE_DE_CHIFFRES));
                                        /* Nombre de chiffres codant le numero des images de la sequence a generer.                  */
     DEFV(CHAR,INIT(POINTERc(nom_image),NOM_UNDEF));
                                        /* Nom courant des images.                                                                   */
     /*..............................................................................................................................*/
     EGAL(Z_Buffer_____valeur_initiale,VALEUR_INITIALE_DU_Z_BUFFER_POUR_LES_MONTAGNES);
                                        /* Introduit le 20060414171417...                                                            */
     EGAL(Imontagnes_precises_____visualiser_la_falaise_avant_de_la_montagne,VRAI);
                                        /* Introduit le 20060414173128...                                                            */

     Cinitialisation(coin_bas_gauche_de_depart_puis_courant,X_DEPART_GAUCHE,Y_DEPART_BAS);
     Cinitialisation(coin_haut_droite_de_depart_puis_courant,X_DEPART_DROITE,Y_DEPART_HAUT);
                                        /* Definition de la fenetre de depart.                                                       */
     Cinitialisation(coin_bas_gauche_d_arrivee,X_ARRIVEE_GAUCHE,Y_ARRIVEE_BAS);
     Cinitialisation(coin_haut_droite_d_arrivee,X_ARRIVEE_DROITE,Y_ARRIVEE_HAUT);
                                        /* Definition de la fenetre d'arrivee.                                                       */
     INITIALISATION_ACCROISSEMENT_2D(translation_montagne,TRANS_OX,TRANS_OY);
     INITIALISATION_POINT_2D(source_lumineuse,X_SOURCE,Y_SOURCE);
                                        /* Definition de la montagne...                                                              */

     GET_ARGUMENTSv(nombre_d_arguments
                   ,BLOC(GET_ARGUMENT_C("imageR=""R=",nom_imageR);
                         GET_ARGUMENT_I("chiffres=",nombre_de_chiffres);
                         GET_ARGUMENT_I("n=""images=""N=",nombre_d_images);
                         GET_ARGUMENT_I("iterations=",nombre_d_iterations);
                         GET_ARGUMENT_F("exposantM=""exposant=",Imandelbrot_dans_C_____exposant);
                                        /* Arguments introduits le 20091118122305...                                                 */
                                        /*                                                                                           */
                                        /* Le 20100608132420, passage de 'GET_ARGUMENT_I(...)' a 'GET_ARGUMENT_F(...)'.              */
                         GET_ARGUMENT_F("seuilM=""seuil=",Imandelbrot_dans_C_____seuil_du_carre_du_module);
                                        /* Arguments introduits le 20091118184202...                                                 */
                         GET_ARGUMENT_F("xbgd=""xbgD=",Reelle(coin_bas_gauche_de_depart_puis_courant));
                         GET_ARGUMENT_F("ybgd=""ybgD=",Imaginaire(coin_bas_gauche_de_depart_puis_courant));
                         GET_ARGUMENT_F("xhdd=""xhdD=",Reelle(coin_haut_droite_de_depart_puis_courant));
                         GET_ARGUMENT_F("yhdd=""yhdD=",Imaginaire(coin_haut_droite_de_depart_puis_courant));
                         GET_ARGUMENT_F("xbga=""xbgA=",Reelle(coin_bas_gauche_d_arrivee));
                         GET_ARGUMENT_F("ybga=""ybgA=",Imaginaire(coin_bas_gauche_d_arrivee));
                         GET_ARGUMENT_F("xhda=""xhdA=",Reelle(coin_haut_droite_d_arrivee));
                         GET_ARGUMENT_F("yhda=""yhdA=",Imaginaire(coin_haut_droite_d_arrivee));
                         GET_ARGUMENT_F("precision=",precision);
                         GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie);
                         GET_ARGUMENT_F("modulation=",modulation_du_champ);
                         GET_ARGUMENT_L("depth-cueing=",depth_cueing);
                         GET_ARGUMENT_L("avion=",vue_d_avion);
                         GET_ARGUMENT_L("Imontagnes_precises_____compatibilite_20061220=""compatibilite_20061220="
                                       ,Imontagnes_precises_____compatibilite_20061220
                                        );
                                        /* Introduit le 20061220110516...                                                            */
                         GET_ARGUMENT_F("fc1=""facteur_cretes_1="
                                       ,Imontagnes_precises_____facteur_1_interpolation_toutes_lignes_de_crete
                                        );
                         GET_ARGUMENT_F("fc2=""facteur_cretes_2="
                                       ,Imontagnes_precises_____facteur_2_interpolation_toutes_lignes_de_crete
                                        );
                                        /* Introduit le 20061220120334 et etendu le 20061220135323 a deux facteurs...                */

                         CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_COMPLEXES;
                                        /* Cette extension a ete introduite le 20070129111928...                                     */
                         )
                    );

     CALi(Inoir(ImageR));
                                        /* Initialisation de l'image Resultat.                                                       */

     CALS(Iblanc(ImageA2));
                                        /* A2 : texture a appliquer sur la montagne.                                                 */

     CALCUL_DU_RAPPORT_DE_REDUCTION;

     DoIn(numero_d_image,PREMIERE_IMAGE_ABSOLUE,LSTX(PREMIERE_IMAGE_ABSOLUE,nombre_d_images),PAS_DES_IMAGES)
          Bblock
          Test(IFGE(numero_d_image,PREMIERE_IMAGE))
               Bblock
                                        /* Afin de se positionner dans la sequence...                                                */
               CALS(Imandelbrot_dans_C(ImageA1
                                      ,ADRESSE(coin_bas_gauche_de_depart_puis_courant)
                                      ,ADRESSE(coin_haut_droite_de_depart_puis_courant)
                                      ,nombre_d_iterations
                                      ,visualisation_arguments_de_sortie
                                       )
                    );
                                        /* A1 : partie courante de l'ensemble de Mandelbrot.                                         */
               CALS(Iinit_Z_Buffer());
               CALS(Imontagnes(ImageR
                              ,modulation_du_champ,ImageA1
                              ,ImageA2
                              ,VRAI
                              ,ADRESSE(translation_montagne)
                              ,ombres_portees,largeur_zone_penombre
                              ,ADRESSE(source_lumineuse)
                              ,depth_cueing,min_depth_cueing
                              ,vue_d_avion
                              ,ALIAS
                               )
                    );
                                        /* R  : montagnes correspondant a la partie courante de l'ensemble de Mandelbrot.            */

               EGAL(nom_image,chain_Aconcaten2_sauf_nom_pipe(nom_imageR,chain_numero(numero_d_image,nombre_de_chiffres)));
               CALi(Iupdate_image(nom_image,ImageR));
               CALZ_FreCC(nom_image);
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          REDUCTION_DE_LA_FENETRE_COURANTE;
          Eblock
     EDoI

     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.