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

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

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N S   D E   B A S E   E T   U N I V E R S E L L E S  :                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrk/attractor.11.I"

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        P A R A M E T R E S  :                                                                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   LES_FENETRES_DE_DEPART_ET_D_ARRIVEE_SONT_IDENTIQUES                                                                           \
                    VRAI                                                                                                                \
                                        /* Les fenetres de Depart et d'Arrivee sont-elles identiques (introduit le 20091121115929) ? */ \
                                        /* La valeur par defaut ('VRAI') est logique pour un calcul hyper-complexe puisque dans ce   */ \
                                        /* cas, en general, les differentes fenetres correspondent a differentes coupes dans l'objet */ \
                                        /* et ces fenetres doivent donc etres alors identiques...                                    */

#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.                                                       */
#define   X_ARRIVEE_GAUCHE                                                                                                              \
                    X_DEPART_GAUCHE
#define   X_ARRIVEE_DROITE                                                                                                              \
                    X_DEPART_DROITE
#define   Y_ARRIVEE_BAS                                                                                                                 \
                    Y_DEPART_BAS
#define   Y_ARRIVEE_HAUT                                                                                                                \
                    Y_DEPART_HAUT
                                        /* Definition de la fenetre d'arrivee.                                                       */

#define   COMPOSANTE_J                                                                                                                  \
                    FZERO
#define   COMPOSANTE_K                                                                                                                  \
                    FZERO
                                        /* Definition des composantes 'J' et 'K' des quaternions.                                    */

#define   NOMBRE_D_ITERATIONS                                                                                                           \
                    HUIT                                                                                                                \
                                        /* Nombre d'iterations demandees.                                                            */

#define   VISUALISATION_ARGUMENTS_DE_SORTIE                                                                                             \
                    FAUX
#define   PONDERATION_ARGUMENTT                                                                                                         \
                    FU
#define   PONDERATION_ARGUMENTP                                                                                                         \
                    FZERO
#define   PONDERATION_ARGUMENTA                                                                                                         \
                    FZERO
                                        /* Faut-il visualiser les arguments de sortie ('VRAI'), auquel cas les trois parametres      */
                                        /* suivants sont les ponderations des 3 angles {T,P,A}, ou les equipotentielles ('FAUX').    */

#define   PRECISION                                                                                                                     \
                    PARE(0.00001)                                                                                                       \
                                        /* Pour calculer le rapport de reduction.                                                    */

#define   NOMBRE_D_IMAGES                                                                                                               \
                    UN                                                                                                                  \
                                        /* Nombre d'images a generer.                                                                */
#define   PREMIERE_IMAGE                                                                                                                \
                    NUMERO_DE_LA_PREMIERE_PERIODE_DE_LA_SIMULATION                                                                      \
                                        /* Numero de la premiere image a generer : ce parametre permet de faire une reprise sur une  */ \
                                        /* sequence interrompue. On notera qu'habituellement sa valeur est 'PREMIERE_IMAGE_ABSOLUE', */ \
                                        /* mais que pour des raisons de compatibilite avec les differents modules inclus dans ce     */ \
                                        /* programme, on lui substitue 'NUMERO_DE_LA_PREMIERE_PERIODE_DE_LA_SIMULATION' partout...   */

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

#define   GENERATION_D_UN_TYPE_D_IMAGE(generer_ce_type_d_image,nom_du_type_d_image,image)                                               \
                    Bblock                                                                                                              \
                    Test(IL_FAUT(generer_ce_type_d_image))                                                                              \
                         Bblock                                                                                                         \
                         EGAL(nom_image                                                                                                 \
                             ,chain_Aconcaten2_sauf_nom_pipe(nom_du_type_d_image                                                        \
                                                            ,chain_numero(numero_d_image,nombre_de_chiffres)                            \
                                                             )                                                                          \
                              );                                                                                                        \
                         CALi(Iupdate_image(nom_image,image));                                                                          \
                         CALZ_FreCC(nom_image);                                                                                         \
                         Eblock                                                                                                         \
                    ATes                                                                                                                \
                         Bblock                                                                                                         \
                         Eblock                                                                                                         \
                    ETes                                                                                                                \
                    Eblock                                                                                                              \
                                        /* Procedure generant si necessaire un certain type d'image.                                 */

#define   CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_MANDELBROT                                                              \
                    Bblock                                                                                                              \
                    vTRANSFORMAT_31(nombre_d_iterations                                                                                 \
                                   ,sNOMBRE_D_ITERATIONS,numero_de_la_periode_courante,fichier_NOMBRE_D_ITERATIONS                      \
                                    );                                                                                                  \
                                        /* Calcul du nombre d'iterations lorsqu'il est variable.                                     */ \
                                                                                                                                        \
                    vTRANSFORMAT_31(composante_J,sCOMPOSANTE_J,numero_de_la_periode_courante,fichier_COMPOSANTE_J);                     \
                                        /* Calcul de la composante 'J' des quaternions lorsqu'elle est variable.                     */ \
                    vTRANSFORMAT_31(composante_K,sCOMPOSANTE_K,numero_de_la_periode_courante,fichier_COMPOSANTE_K);                     \
                                        /* Calcul de la composante 'K' des quaternions lorsqu'elle est variable.                     */ \
                                                                                                                                        \
                    CALS(Imandelbrot_dans_HC(ImageR                                                                                     \
                                            ,ADRESSE(coin_bas_gauche_de_depart_puis_courant)                                            \
                                            ,ADRESSE(coin_haut_droite_de_depart_puis_courant)                                           \
                                            ,composante_J                                                                               \
                                            ,composante_K                                                                               \
                                            ,nombre_d_iterations                                                                        \
                                            ,visualisation_arguments_de_sortie                                                          \
                                            ,ponderation_argumentT,ponderation_argumentP,ponderation_argumentA                          \
                                             )                                                                                          \
                         );                                                                                                             \
                                        /* Calcul de l'ensemble de Mandelbrot dans l'espace hyper-Complexe.                          */ \
                                                                                                                                        \
                    GENERATION_D_UN_TYPE_D_IMAGE(VRAI                                                                                   \
                                                ,nom_imageR                                                                             \
                                                ,ImageR                                                                                 \
                                                 );                                                                                     \
                                        /* Generation des images demandees.                                                          */ \
                    Eblock

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N S   D E S   R O T A T I O N S   T R I D I M E N S I O N N E L L E S  :                                 */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrc/julia.41.I"

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N S   D E   L ' I N T E R P O L A T I O N   D E S   P A R A M E T R E S  :                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
dfTRANSFORMAT_31(liste_NOMBRE_D_ITERATIONS,fichier_NOMBRE_D_ITERATIONS,NOMBRE_D_ITERATIONS_IMPLICITE,NOMBRE_D_ITERATIONS)
                                        /* Definition du fichier des nombres d'iterations.                                           */
#define   sNOMBRE_D_ITERATIONS(numero_de_la_periode)                                                                                    \
                    sTRANSFORMAT_31(numero_de_la_periode,liste_NOMBRE_D_ITERATIONS)
                                        /* Formule generale definissant les variations du nombre d'iterations.                       */

dfTRANSFORMAT_31(liste_COMPOSANTE_J,fichier_COMPOSANTE_J,COMPOSANTE_J_IMPLICITE,COMPOSANTE_J)
dfTRANSFORMAT_31(liste_COMPOSANTE_K,fichier_COMPOSANTE_K,COMPOSANTE_K_IMPLICITE,COMPOSANTE_K)
                                        /* Definition du fichier des composantes 'J' et 'K' des quaternions.                         */
#define   sCOMPOSANTE_J(numero_de_la_periode)                                                                                           \
                    sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_J)
#define   sCOMPOSANTE_K(numero_de_la_periode)                                                                                           \
                    sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_K)
                                        /* Formule generale definissant les variations des composantes 'J' et 'K' des quaternions.   */

#define   numero_de_la_periode_courante                                                                                                 \
                    numero_d_image
#define   nombre_de_periodes_de_la_simulation                                                                                           \
                    nombre_d_images
                                        /* Pour assurer la compatibilite avec '$xrq/nucleon.Lf.2.I'...                               */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   D E   L ' E N S E M B L E   D E   M A N D E L B R O T                                                        */
/*        D A N S   L ' E S P A C E   H Y P E R - C O M P L E X E   :                                                                */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
     DEFV(Logical,INIT(les_fenetres_de_depart_et_d_arrivee_sont_identiques,LES_FENETRES_DE_DEPART_ET_D_ARRIVEE_SONT_IDENTIQUES));
                                        /* Les fenetres de Depart et d'Arrivee sont-elles identiques (introduit le 20091121115929) ? */
                                        /* La valeur par defaut ('VRAI') est logique pour un calcul hyper-complexe puisque dans ce   */
                                        /* cas, en general, les differentes fenetres correspondent a differentes coupes dans l'objet */
                                        /* et ces fenetres doivent donc etres alors identiques...                                    */
     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(Float,INIT(composante_J,COMPOSANTE_J));
     DEFV(Float,INIT(composante_K,COMPOSANTE_K));
                                        /* Definition des composantes 'J' et 'K' des quaternions.                                    */

     DEFV(Positive,INIT(nombre_d_iterations,NOMBRE_D_ITERATIONS));
                                        /* Nombre maximal d'iterations a effectuer.                                                  */

     DEFV(Logical,INIT(visualisation_arguments_de_sortie,VISUALISATION_ARGUMENTS_DE_SORTIE));
     DEFV(Float,INIT(ponderation_argumentT,PONDERATION_ARGUMENTT));
     DEFV(Float,INIT(ponderation_argumentP,PONDERATION_ARGUMENTP));
     DEFV(Float,INIT(ponderation_argumentA,PONDERATION_ARGUMENTA));
                                        /* Faut-il visualiser les arguments de sortie ('VRAI'), auquel cas les trois parametres      */
                                        /* suivants sont les ponderations des 3 angles {T,P,A}, ou les equipotentielles ('FAUX').    */

     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(Int,INIT(nombre_d_images,NOMBRE_D_IMAGES));
                                        /* Nombre d'images a generer.                                                                */
     DEFV(Int,INIT(numero_d_image,NUMERO_DE_LA_PREMIERE_PERIODE_DE_LA_SIMULATION));
                                        /* 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.                                                                   */

     DEFV(CHAR,INIC(POINTERc(nom_imageI),NOM_PIPE));
     DEFV(CHAR,INIC(POINTERc(nom_imageF),NOM_PIPE));
                                        /* Noms des eventuelles 'imageJ' definissant respectivement les etats Initial et Final       */
                                        /* (ceci a ete introduit le 20150228115203...).                                              */
     /*..............................................................................................................................*/
     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.                                                       */

     DEBUT_DE_L_INITIALISATION_DES_INTERPOLATIONS;
     FIN_DE_L_INITIALISATION_DES_INTERPOLATIONS;
                                        /* Initialisations generales pour l'interpolation des rotations...                           */

     iTRANSFORMAT_31(liste_NOMBRE_D_ITERATIONS,NOMBRE_D_ITERATIONS_IMPLICITE);
                                        /* Initialisation du nombre d'iterations.                                                    */

     iTRANSFORMAT_31(liste_COMPOSANTE_J,COMPOSANTE_J_IMPLICITE);
     iTRANSFORMAT_31(liste_COMPOSANTE_K,COMPOSANTE_K_IMPLICITE);
                                        /* Initialisation des composantes 'J' et 'K' des quaternions.                                */

     GET_ARGUMENTSi(nombre_d_arguments
                   ,BLOC(PROCESS_ARGUMENT_I("n=""images=""N=",nombre_d_images
                                           ,BLOC(VIDE;)
                                           ,BLOC(PRINT_AVERTISSEMENT("'n=''images=''N=' doit etre avant toute entree de fichiers");)
                                            );

                         PROCESS_ARGUMENTS_GEOMETRIQUES_DE_ROTATION;

                         PROCESS_ARGUMENTS_GEOMETRIQUES_DE_ROTATION_2_HC;
                                        /* Introduit le 20091222140116...                                                            */

                         PROCESS_ARGUMENT_FICHIER("NOMBRE_D_ITERATIONS="
                                                 ,fichier_NOMBRE_D_ITERATIONS
                                                 ,liste_NOMBRE_D_ITERATIONS
                                                 ,NOMBRE_D_ITERATIONS_IMPLICITE
                                                 ,gTRANSFORMAT_31
                                                  );

                         PROCESS_ARGUMENT_FICHIER("COMPOSANTE_J="
                                                 ,fichier_COMPOSANTE_J
                                                 ,liste_COMPOSANTE_J
                                                 ,COMPOSANTE_J_IMPLICITE
                                                 ,gTRANSFORMAT_31
                                                  );
                         PROCESS_ARGUMENT_FICHIER("COMPOSANTE_K="
                                                 ,fichier_COMPOSANTE_K
                                                 ,liste_COMPOSANTE_K
                                                 ,COMPOSANTE_K_IMPLICITE
                                                 ,gTRANSFORMAT_31
                                                  );

                         GET_ARGUMENT_C("imageR=""R=",nom_imageR);
                         GET_ARGUMENT_I("chiffres=",nombre_de_chiffres);

                         GET_ARGUMENT_C("imageI=""I=",nom_imageI);
                         GET_ARGUMENT_C("imageF=""F=",nom_imageF);
                                        /* Arguments introduits le 20150228115203...                                                 */

                         GET_ARGUMENT_I("iterations=",nombre_d_iterations);

                         GET_ARGUMENT_F("exposantM=""exposant=",Imandelbrot_dans_HC_____exposant);
                                        /* Arguments introduits le 20091118121214...                                                 */
                                        /*                                                                                           */
                                        /* Le 20100608132532, passage de 'GET_ARGUMENT_I(...)' a 'GET_ARGUMENT_F(...)'.              */
                         GET_ARGUMENT_F("seuilM=""seuil=",Imandelbrot_dans_HC_____seuil_du_carre_du_module);
                                        /* Arguments introduits le 20091118183519...                                                 */

                         GET_ARGUMENT_L("fenetres_depart_arrivee_identiques=""fdai="
                                       ,les_fenetres_de_depart_et_d_arrivee_sont_identiques
                                        );
                                        /* Arguments introduits le 20091121115929...                                                 */
                         GET_ARGUMENT_F("xbg=""xbgd=""xbgD=",Reelle(coin_bas_gauche_de_depart_puis_courant));
                         GET_ARGUMENT_F("ybg=""ybgd=""ybgD=",Imaginaire(coin_bas_gauche_de_depart_puis_courant));
                         GET_ARGUMENT_F("xhd=""xhdd=""xhdD=",Reelle(coin_haut_droite_de_depart_puis_courant));
                         GET_ARGUMENT_F("yhd=""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));
                                        /* Introduction des arguments {"xbg=","ybg=","xhd=","yhd="} le 20091122120041...             */
                                        /*                                                                                           */
                                        /* Les arguments suivants de 'v $xrc/mandel.41$K' :                                          */
                                        /*                                                                                           */
                                        /*                  OX=... OY=... EX=... EY=...                                              */
                                        /*                  fenetres_depart_arrivee_identiques=VRAI                                  */
                                        /*                  editer_fenetre_courante=VRAI                                             */
                                        /*                                                                                           */
                                        /* s'utiliseront conjointement avec les arguments suivants de 'v $xci/display$K' :           */
                                        /*                                                                                           */
                                        /*                  ModeFlot='f'                                                             */
                                        /*                  ChiffresFlot=...                                                         */
                                        /*                                                                                           */
                                        /*                  editer_XY_curseur=VRAI                                                   */
                                        /*                  normaliser_XY_curseur=VRAI                                               */
                                        /*                  DX_curseur=... DY_curseur=...                                            */
                                        /*                  OX=... OY=... EX=... EY=...                                              */
                                        /*                                                                                           */
                                        /* Supposons que '$xrc/mandel.41$X' en soit a :                                              */
                                        /*                                                                                           */
                                        /*                  xbg=+0.0014  xhd=+0.0016  ybg=+0.0045  yhd=+0.0047                       */
                                        /*                         ---          ---          ---          ---                        */
                                        /*                                                                                           */
                                        /* Supposons qu'alors '$xci/display$X' renvoie :                                             */
                                        /*                                                                                           */
                                        /*                  curseur=(+0.00153333,+0.00467777)                                        */
                                        /*                              -----       -----                                            */
                                        /*                                                                                           */
                                        /* on aura alors interet a choisir pour l'etape suivante de '$xrc/mandel.41$X' :             */
                                        /*                                                                                           */
                                        /*                  xbg=+0.00152 xhd=+0.00154 ybg=+0.00466 yhd=+0.00468                      */
                                        /*                         ----         ----         ----         ----                       */
                                        /*                                                                                           */
                                        /* et ainsi {+0.00152,+0.00154} sont de part et d'autre de +0.00153, de meme que             */
                                        /*              ----     ----                                 ----                           */
                                        /* {+0.00466,+0.00468} sont de part et d'autre de +0.00467 (dans les deux cas on fait        */
                                        /*     ----     ----                                 ----                                    */
                                        /* donc moins et plus 0.00001).                                                              */

                         GET_ARGUMENT_F("precision=",precision);

                         GET_ARGUMENT_F("J=",composante_J);
                         GET_ARGUMENT_F("K=",composante_K);

                         GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie);
                         GET_ARGUMENT_F("ponderationT=",ponderation_argumentT);
                         GET_ARGUMENT_F("ponderationP=",ponderation_argumentP);
                         GET_ARGUMENT_F("ponderationA=",ponderation_argumentA);

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DES_ITERATIONS_1;
                                        /* Introduit le 20091120122821...                                                            */
                         GET_ARGUMENT_L("F_ITERATION_DANS_HC_____compatibilite_20110825=""compatibilite_20110825="
                                       ,F_ITERATION_DANS_HC_____compatibilite_20110825
                                        );
                                        /* Introduit le 20110826123121.                                                              */

                         CONTROLE_DE_L_ARITHMETIQUE_FOLKLORIQUE_DES_NOMBRES_COMPLEXES;
                                        /* Cette extension a ete introduite le 20181211112338...                                     */
                         CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_HYPER_COMPLEXES;
                                        /* Cette extension a ete introduite le 20070129112331...                                     */

                         CONTROLE_DE_LA_TRANFORMATION_DE_LA_SUITE_DE_POINTS_DANS_HC;
                                        /* Cette extension a ete introduite le 20100513181745...                                     */

                         CONTROLE_DU_REPLIEMENT_DE_L_ESPACE_DES_NOMBRES_HYPER_COMPLEXES;
                                        /* Cette extension a ete introduite le 20121104080422...                                     */
                         )
                    );

     Test(EST_VRAI(les_fenetres_de_depart_et_d_arrivee_sont_identiques))
                                        /* Test introduit le 20091121115929...                                                       */
          Bblock
          Cegal(coin_bas_gauche_d_arrivee,coin_bas_gauche_de_depart_puis_courant);
          Cegal(coin_haut_droite_d_arrivee,coin_haut_droite_de_depart_puis_courant);
                                        /* Definition de la fenetre d'Arrivee identique a la fenetre de Depart.                      */
          Eblock
     ATes
          Bblock
          Eblock
     ETes

     CALCUL_DU_RAPPORT_DE_REDUCTION;

     Test(IFOU(IFNE_chaine(nom_imageI,NOM_PIPE),IFNE_chaine(nom_imageF,NOM_PIPE)))
                                        /* Test introduit le 20150228120012 et qui permet de sauvergarder l'etat Final et de         */
                                        /* pouvoir le reutiliser ulterieurement comme etat Initial...                                */
          Bblock
          BMalloc(Imandelbrot_dans_HC_____etat_initial,imageHJ);
          BMalloc(Imandelbrot_dans_HC_____etat_final__,imageHJ);
                                        /* Il faut declarer ces images a priori (meme en cas de non utilisation et ceci a cause      */
                                        /* des validations {BblockV,EblockV}...).                                                    */

          Test(IFGT(nombre_d_images,UN))
               Bblock
               PRINT_ATTENTION("dans le mode 'Initial/Final' une seule image sera generee");
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          CALi(Inoir(ImageR));

          Test(IFNE_chaine(nom_imageI,NOM_PIPE))
               Bblock
               EGAL(Imandelbrot_dans_HC_____charger_____l_etat_initial,VRAI);

               Test(PAS_D_ERREUR(CODE_ERROR(IloadHJ_image(Imandelbrot_dans_HC_____etat_initial,nom_imageI))))
                    Bblock
                    Eblock
               ATes
                    Bblock
                    Test__CODE_ERREUR__ERREUR07;

                    CALi(IHJinitialisation(Imandelbrot_dans_HC_____etat_initial,ADRESSE(HC_____quaternion__0__0__0__0)));
                    Eblock
               ETes
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          Test(IFNE_chaine(nom_imageF,NOM_PIPE))
               Bblock
               EGAL(Imandelbrot_dans_HC_____sauvegarder_l_etat_final__,VRAI);
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_MANDELBROT;

          Test(IL_FAUT(Imandelbrot_dans_HC_____sauvegarder_l_etat_final__))
               Bblock
               CALi(IupdateHJ_image(nom_imageF,Imandelbrot_dans_HC_____etat_final__));
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          EMalloc(Imandelbrot_dans_HC_____etat_final__,imageJ);
          EMalloc(Imandelbrot_dans_HC_____etat_initial,imageJ);
          Eblock
     ATes
          Bblock
          DoIn(numero_d_image
              ,NUMERO_DE_LA_PREMIERE_PERIODE_DE_LA_SIMULATION
              ,LSTX(NUMERO_DE_LA_PREMIERE_PERIODE_DE_LA_SIMULATION,nombre_d_images)
              ,I
               )
               Bblock
               CHANGEMENT_DU_POINT_DE_VUE_GEOMETRIQUE;
                                        /* Mise en place des rotations tridimensionnelles dans les Quaternions.                      */

               Test(IFGE(numero_d_image,PREMIERE_IMAGE))
                    Bblock
                                        /* Afin de se positionner dans la sequence...                                                */
                    CALi(Inoir(ImageR));
                                        /* Initialisation de l'image Resultat a priori...                                            */

                    CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_MANDELBROT;
                    Eblock
               ATes
                    Bblock
                    Eblock
               ETes

               REDUCTION_DE_LA_FENETRE_COURANTE;
               Eblock
          EDoI
          Eblock
     ETes

     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.