/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   D E S   R A C I N E S   N - I E M E S   D E   L ' U N I T E                                                  */
/*        D A N S   L E   P L A N   C O M P L E X E  :                                                                               */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrc/racN.01$K' :                                                                                               */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 20010602212015).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

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

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        P A R A M E T R E S  :                                                                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   EXPOSANT                                                                                                                      \
                    TROIS                                                                                                               \
                                        /* Exposant.                                                                                 */
#define   X_DEPART_GAUCHE                                                                                                               \
                    PARE(-1.0)
#define   X_DEPART_DROITE                                                                                                               \
                    PARE(1.0)
#define   Y_DEPART_BAS                                                                                                                  \
                    PARE(-1.0)
#define   Y_DEPART_HAUT                                                                                                                 \
                    PARE(1.0)
                                        /* Definition de la fenetre de depart (vue d'ensemble).                                      */
#define   X_ARRIVEE_GAUCHE                                                                                                              \
                    PARE(-1.0)
#define   X_ARRIVEE_DROITE                                                                                                              \
                    PARE(1.0)
#define   Y_ARRIVEE_BAS                                                                                                                 \
                    PARE(-1.0)
#define   Y_ARRIVEE_HAUT                                                                                                                \
                    PARE(1.0)
                                        /* Definition de la fenetre d'arrivee (dans les filaments).                                  */

#define   NOMBRE_D_ITERATIONS                                                                                                           \
                    SOUS(BLANC,NOIR)                                                                                                    \
                                        /* Nombre d'iterations demandees.                                                            */
#define   PRECISION                                                                                                                     \
                    PARE(0.00001)                                                                                                       \
                                        /* Pour calculer le rapport de reduction.                                                    */
#include  xci/sequence.01.I"

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   D E S   R A C I N E S   N - I E M E S   D E   L ' U N I T 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(Positive,INIT(exposant,EXPOSANT));
                                        /* Exposant.                                                                                 */

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

     DEFV(Logical,INIT(visualisation_arguments_de_sortie,VRAI));
                                        /* Indicateur de choix entre la visualisation du nombre d'iterations ('FAUX'),               */
                                        /* et l'argument de sortie ('VRAI').                                                         */
                                        /*                                                                                           */
                                        /* On notera le 20061205165402 que c'est le fait que l'on visualise a priori l'argument      */
                                        /* de sortie qui permet de visualiser la racine (parmi les N) vers laquelle on converge,     */
                                        /* puisque qu'elles sont sur le cercle trigonometrique. Ainsi, les N arguments possibles     */
                                        /* sont tres differents les uns des autres et se distinguent par les couleurs. C'est la      */
                                        /* raison pour laquelle a disparu :                                                          */
                                        /*                                                                                           */
                                        /*                  GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie);          */
                                        /*                                                                                           */
                                        /* a cette date...                                                                           */
     /*..............................................................................................................................*/
     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.                                                       */

     GET_ARGUMENTSi(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("exposant=""e=",exposant);
                         GET_ARGUMENT_I("iterations=",nombre_d_iterations);
                         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("marquer_racines=""mr="
                                       ,Iracines_n_iemes_de_l_unite_dans_C_____marquer_les_racines_de_l_unite
                                        );
                         GET_ARGUMENT_P("niveau_racines=""nr=",Iracines_n_iemes_de_l_unite_dans_C_____niveau_des_racines_de_l_unite);
                                        /* Arguments introduits le 20200623121508...                                                 */

                         GET_ARGUMENT_L("gF_ITERATION_DANS_C_____compatibilite_20211015=""compatibilite_20211015="
                                       ,gF_ITERATION_DANS_C_____compatibilite_20211015
                                        );
                                        /* Introduit le 20211015104719.                                                              */

                         GET_ARGUMENT_L("editer_trajectoires=""et=",gF_ITERATION_DANS_C_____editer_les_trajectoires);
                         GET_ARGUMENT_L("normaliser_trajectoires=""nt=",gF_ITERATION_DANS_C_____normaliser_les_trajectoires);

                         GET_ARGUMENT_I("translation_X_trajectoires=""tXt="
                                       ,gF_ITERATION_DANS_C_____translation_X_d_editions_des_trajectoires
                                        );
                         GET_ARGUMENT_I("Xmin_trajectoires=""Xmint=",gF_ITERATION_DANS_C_____Xmin__d_editions_des_trajectoires);
                         GET_ARGUMENT_I("Xmax_trajectoires=""Xmaxt=",gF_ITERATION_DANS_C_____Xmax__d_editions_des_trajectoires);
                         GET_ARGUMENT_F("maille_X_trajectoires=""mXt=",gF_ITERATION_DANS_C_____maille_X_d_editions_des_trajectoires);

                         GET_ARGUMENT_I("translation_Y_trajectoires=""tYt="
                                       ,gF_ITERATION_DANS_C_____translation_Y_d_editions_des_trajectoires
                                        );
                         GET_ARGUMENT_I("Ymin_trajectoires=""Ymint=",gF_ITERATION_DANS_C_____Ymin__d_editions_des_trajectoires);
                         GET_ARGUMENT_I("Ymax_trajectoires=""Ymaxt=",gF_ITERATION_DANS_C_____Ymax__d_editions_des_trajectoires);
                         GET_ARGUMENT_F("maille_Y_trajectoires=""mYt=",gF_ITERATION_DANS_C_____maille_Y_d_editions_des_trajectoires);
                                        /* Introduits le 20211010113059 et completes le 20211012130202 par la possibilite de         */
                                        /* "resserer" le maillage des points dont on visualise les trajectoires...                   */

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DES_ITERATIONS_1;
                                        /* Introduit sous cette forme le 20070108142834...                                           */
                         )
                    );

     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...                                                */
               CALi(Inoir(ImageR));
                                        /* Initialisation de l'image Resultat a priori (introduit le 20061205151138)...              */

               CALS(Iracines_n_iemes_de_l_unite_dans_C(ImageR
                                                      ,ADRESSE(coin_bas_gauche_de_depart_puis_courant)
                                                      ,ADRESSE(coin_haut_droite_de_depart_puis_courant)
                                                      ,exposant
                                                      ,nombre_d_iterations
                                                      ,visualisation_arguments_de_sortie
                                                       )
                    );

               EGAL(nom_image,chain_Aconcaten2_sauf_nom_pipe(nom_imageR,chain_numero(numero_d_image,nombre_de_chiffres)));
               CALi(gIupdate_image(NOTL(iterations_dans_C_HC_HHC_____etendre_la_precision)
                                  ,nom_imageR
                                  ,ImageR,iterations_dans_C_HC_HHC_____IFmageR
                                   )
                    );
               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.