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

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        F I C H I E R S   D ' I N C L U D E S  :                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  INCLUDES_BASE
#include  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 ? La valeur par defaut ('VRAI') */ \
                                        /* est logique pour un calcul hyper-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(0.0)
#define   X_DEPART_DROITE                                                                                                               \
                    PARE(1.0)
#define   Y_DEPART_BAS                                                                                                                  \
                    PARE(-0.5)
#define   Y_DEPART_HAUT                                                                                                                 \
                    PARE(0.5)
                                        /* 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
#define   COMPOSANTE_S                                                                                                                  \
                    FZERO
#define   COMPOSANTE_T                                                                                                                  \
                    FZERO
#define   COMPOSANTE_U                                                                                                                  \
                    FZERO
#define   COMPOSANTE_V                                                                                                                  \
                    FZERO
                                        /* Definition des composantes {J,K,S,T,U,V} des Octonions.                                   */

#define   PARTIE_REELLE_DU_TAUX_DE_CROISSANCE                                                                                           \
                    PARE(3.1)
#define   PARTIE_IMAGINAIRE_DU_TAUX_DE_CROISSANCE                                                                                       \
                    PARE(0.0)
#define   PARTIE_JMAGINAIRE_DU_TAUX_DE_CROISSANCE                                                                                       \
                    PARE(0.0)
#define   PARTIE_KMAGINAIRE_DU_TAUX_DE_CROISSANCE                                                                                       \
                    PARE(0.0)
#define   PARTIE_SMAGINAIRE_DU_TAUX_DE_CROISSANCE                                                                                       \
                    PARE(0.0)
#define   PARTIE_TMAGINAIRE_DU_TAUX_DE_CROISSANCE                                                                                       \
                    PARE(0.0)
#define   PARTIE_UMAGINAIRE_DU_TAUX_DE_CROISSANCE                                                                                       \
                    PARE(0.0)
#define   PARTIE_VMAGINAIRE_DU_TAUX_DE_CROISSANCE                                                                                       \
                    PARE(0.0)
                                        /* Taux de croissance de l'ensemble.                                                         */

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

#define   VISUALISATION_ARGUMENTS_DE_SORTIE                                                                                             \
                    FAUX
#define   PONDERATION_ARGUMENTA1                                                                                                        \
                    FU
#define   PONDERATION_ARGUMENTA2                                                                                                        \
                    FZERO
#define   PONDERATION_ARGUMENTA3                                                                                                        \
                    FZERO
#define   PONDERATION_ARGUMENTA4                                                                                                        \
                    FZERO
#define   PONDERATION_ARGUMENTA5                                                                                                        \
                    FZERO
#define   PONDERATION_ARGUMENTA6                                                                                                        \
                    FZERO
#define   PONDERATION_ARGUMENTA7                                                                                                        \
                    FZERO
                                        /* Faut-il visualiser les arguments de sortie ('VRAI'), auquel cas les sept parametres       */
                                        /* suivants sont les ponderations des 7 angles {A1,A2,A3,A4,A5,A6,A7}, 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"

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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  :                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrv/champs_5.41.I"

#include  xrq/nucleon.Lf.2.I"

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)
dfTRANSFORMAT_31(liste_COMPOSANTE_S,fichier_COMPOSANTE_S,COMPOSANTE_S_IMPLICITE,COMPOSANTE_S)
dfTRANSFORMAT_31(liste_COMPOSANTE_T,fichier_COMPOSANTE_T,COMPOSANTE_T_IMPLICITE,COMPOSANTE_T)
dfTRANSFORMAT_31(liste_COMPOSANTE_U,fichier_COMPOSANTE_U,COMPOSANTE_U_IMPLICITE,COMPOSANTE_U)
dfTRANSFORMAT_31(liste_COMPOSANTE_V,fichier_COMPOSANTE_V,COMPOSANTE_V_IMPLICITE,COMPOSANTE_V)
                                        /* Definition du fichier des composantes {J,K,S,T,U,V} des Octonions.                        */
#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)
#define   sCOMPOSANTE_S(numero_de_la_periode)                                                                                           \
                    sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_S)
#define   sCOMPOSANTE_T(numero_de_la_periode)                                                                                           \
                    sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_T)
#define   sCOMPOSANTE_U(numero_de_la_periode)                                                                                           \
                    sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_U)
#define   sCOMPOSANTE_V(numero_de_la_periode)                                                                                           \
                    sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_V)
                                        /* Formule generale definissant les variations des composantes {J,K,S,T,U,V} des Octonions.  */

#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   V E R H U L S T                                                            */
/*        D A N S   L ' E S P A C E   H Y P E R - 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 ? La valeur par defaut ('VRAI') */
                                        /* est logique pour un calcul hyper-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));
     DEFV(Float,INIT(composante_S,COMPOSANTE_S));
     DEFV(Float,INIT(composante_T,COMPOSANTE_T));
     DEFV(Float,INIT(composante_U,COMPOSANTE_U));
     DEFV(Float,INIT(composante_V,COMPOSANTE_V));
                                        /* Definition des composantes {J,K,S,T,U,V} des Octonions.                                   */
     DEFV(hyper_hyper_complexe,taux_de_croissance);
                                        /* Taux de croissance de l'ensemble.                                                         */
     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_argumentA1,PONDERATION_ARGUMENTA1));
     DEFV(Float,INIT(ponderation_argumentA2,PONDERATION_ARGUMENTA2));
     DEFV(Float,INIT(ponderation_argumentA3,PONDERATION_ARGUMENTA3));
     DEFV(Float,INIT(ponderation_argumentA4,PONDERATION_ARGUMENTA4));
     DEFV(Float,INIT(ponderation_argumentA5,PONDERATION_ARGUMENTA5));
     DEFV(Float,INIT(ponderation_argumentA6,PONDERATION_ARGUMENTA6));
     DEFV(Float,INIT(ponderation_argumentA7,PONDERATION_ARGUMENTA7));
                                        /* Faut-il visualiser les arguments de sortie ('VRAI'), auquel cas les sept parametres       */
                                        /* suivants sont les ponderations des 7 angles {A1,A2,A3,A4,A5,A6,A7}, 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.                                                                   */
     /*..............................................................................................................................*/
     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.                                                       */

     HHCinitialisation(taux_de_croissance
                      ,PARTIE_REELLE_DU_TAUX_DE_CROISSANCE
                      ,PARTIE_IMAGINAIRE_DU_TAUX_DE_CROISSANCE
                      ,PARTIE_JMAGINAIRE_DU_TAUX_DE_CROISSANCE
                      ,PARTIE_KMAGINAIRE_DU_TAUX_DE_CROISSANCE
                      ,PARTIE_SMAGINAIRE_DU_TAUX_DE_CROISSANCE
                      ,PARTIE_TMAGINAIRE_DU_TAUX_DE_CROISSANCE
                      ,PARTIE_UMAGINAIRE_DU_TAUX_DE_CROISSANCE
                      ,PARTIE_VMAGINAIRE_DU_TAUX_DE_CROISSANCE
                       );
                                        /* Definition du taux de croissance implicite.                                               */

     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);
     iTRANSFORMAT_31(liste_COMPOSANTE_S,COMPOSANTE_S_IMPLICITE);
     iTRANSFORMAT_31(liste_COMPOSANTE_T,COMPOSANTE_T_IMPLICITE);
     iTRANSFORMAT_31(liste_COMPOSANTE_U,COMPOSANTE_U_IMPLICITE);
     iTRANSFORMAT_31(liste_COMPOSANTE_V,COMPOSANTE_V_IMPLICITE);
                                        /* Initialisation des composantes {J,K,S,T,U,V} des Octonions.                               */

     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_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
                                                  );
                         PROCESS_ARGUMENT_FICHIER("COMPOSANTE_S="
                                                 ,fichier_COMPOSANTE_S
                                                 ,liste_COMPOSANTE_S
                                                 ,COMPOSANTE_S_IMPLICITE
                                                 ,gTRANSFORMAT_31
                                                  );
                         PROCESS_ARGUMENT_FICHIER("COMPOSANTE_T="
                                                 ,fichier_COMPOSANTE_T
                                                 ,liste_COMPOSANTE_T
                                                 ,COMPOSANTE_T_IMPLICITE
                                                 ,gTRANSFORMAT_31
                                                  );
                         PROCESS_ARGUMENT_FICHIER("COMPOSANTE_U="
                                                 ,fichier_COMPOSANTE_U
                                                 ,liste_COMPOSANTE_U
                                                 ,COMPOSANTE_U_IMPLICITE
                                                 ,gTRANSFORMAT_31
                                                  );
                         PROCESS_ARGUMENT_FICHIER("COMPOSANTE_V="
                                                 ,fichier_COMPOSANTE_V
                                                 ,liste_COMPOSANTE_V
                                                 ,COMPOSANTE_V_IMPLICITE
                                                 ,gTRANSFORMAT_31
                                                  );

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

                         GET_ARGUMENT_I("iterations=",nombre_d_iterations);

                         GET_ARGUMENT_L("fenetres_depart_arrivee_identiques=""fdai="
                                       ,les_fenetres_de_depart_et_d_arrivee_sont_identiques
                                        );
                         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));
                                        /* Introduction des arguments {"xbg=","ybg=","xhd=","yhd="} le 20091122115932...             */
                         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_F("J=",composante_J);
                         GET_ARGUMENT_F("K=",composante_K);
                         GET_ARGUMENT_F("S=",composante_S);
                         GET_ARGUMENT_F("T=",composante_T);
                         GET_ARGUMENT_F("U=",composante_U);
                         GET_ARGUMENT_F("V=",composante_V);

                         GET_ARGUMENT_F("Rtaux=""Rr=""RR=",HHReelle(taux_de_croissance));
                         GET_ARGUMENT_F("Itaux=""Ir=""IR=",HHImaginaire(taux_de_croissance));
                         GET_ARGUMENT_F("Jtaux=""Jr=""JR=",HHJmaginaire(taux_de_croissance));
                         GET_ARGUMENT_F("Ktaux=""Kr=""KR=",HHKmaginaire(taux_de_croissance));
                         GET_ARGUMENT_F("Staux=""Sr=""SR=",HHSmaginaire(taux_de_croissance));
                         GET_ARGUMENT_F("Ttaux=""Tr=""TR=",HHTmaginaire(taux_de_croissance));
                         GET_ARGUMENT_F("Utaux=""Ur=""UR=",HHUmaginaire(taux_de_croissance));
                         GET_ARGUMENT_F("Vtaux=""Vr=""VR=",HHVmaginaire(taux_de_croissance));

                         GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie);
                         GET_ARGUMENT_F("ponderationA1=",ponderation_argumentA1);
                         GET_ARGUMENT_F("ponderationA2=",ponderation_argumentA2);
                         GET_ARGUMENT_F("ponderationA3=",ponderation_argumentA3);
                         GET_ARGUMENT_F("ponderationA4=",ponderation_argumentA4);
                         GET_ARGUMENT_F("ponderationA5=",ponderation_argumentA5);
                         GET_ARGUMENT_F("ponderationA6=",ponderation_argumentA6);
                         GET_ARGUMENT_F("ponderationA7=",ponderation_argumentA7);
                         GET_ARGUMENT_L("simple=",Iverhulst_dans_HHC_____utiliser_la_forme_la_plus_simple);

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DES_ITERATIONS_1;
                         GET_ARGUMENT_L("F_ITERATION_DANS_HHC_____compatibilite_20110825=""compatibilite_20110825="
                                       ,F_ITERATION_DANS_HHC_____compatibilite_20110825
                                        );

                         CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_HYPER_COMPLEXES;

                                        /* Le 20121103105322, 'CONTROLE_DE_LA_TRANFORMATION_DE_LA_SUITE_DE_POINTS_DANS_HC' a ete     */
                                        /* supprime puisqu'il concernait 'HC' et non pas 'HHC'...                                    */

                         CONTROLE_DU_REPLIEMENT_DE_L_ESPACE_DES_NOMBRES_HYPER_HYPER_COMPLEXES;
                                        /* Cette extension a ete introduite le 20121104081923...                                     */
                         )
                    );

     Test(EST_VRAI(les_fenetres_de_depart_et_d_arrivee_sont_identiques))
          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;

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

               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);
               vTRANSFORMAT_31(composante_K,sCOMPOSANTE_K,numero_de_la_periode_courante,fichier_COMPOSANTE_K);
               vTRANSFORMAT_31(composante_S,sCOMPOSANTE_S,numero_de_la_periode_courante,fichier_COMPOSANTE_S);
               vTRANSFORMAT_31(composante_T,sCOMPOSANTE_T,numero_de_la_periode_courante,fichier_COMPOSANTE_T);
               vTRANSFORMAT_31(composante_U,sCOMPOSANTE_U,numero_de_la_periode_courante,fichier_COMPOSANTE_U);
               vTRANSFORMAT_31(composante_V,sCOMPOSANTE_V,numero_de_la_periode_courante,fichier_COMPOSANTE_V);
                                        /* Calcul des composantes {J,K,S,T,U,V} des Octonions lorsqu'elles sont variables.           */

               CALS(Iverhulst_dans_HHC(ImageR
                                      ,ADRESSE(coin_bas_gauche_de_depart_puis_courant)
                                      ,ADRESSE(coin_haut_droite_de_depart_puis_courant)
                                      ,composante_J
                                      ,composante_K
                                      ,composante_S
                                      ,composante_T
                                      ,composante_U
                                      ,composante_V
                                      ,ADRESSE(taux_de_croissance)
                                      ,nombre_d_iterations
                                      ,visualisation_arguments_de_sortie
                                      ,ponderation_argumentA1
                                      ,ponderation_argumentA2
                                      ,ponderation_argumentA3
                                      ,ponderation_argumentA4
                                      ,ponderation_argumentA5
                                      ,ponderation_argumentA6
                                      ,ponderation_argumentA7
                                       )
                    );
                                        /* Calcul de l'ensemble de Verhulst dans l'espace hyper-hyper-Complexe.                      */

               GENERATION_D_UN_TYPE_D_IMAGE(VRAI
                                           ,nom_imageR
                                           ,ImageR
                                            );
                                        /* Generation des images demandees.                                                          */
               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.