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