/*************************************************************************************************************************************/
/* */
/* C A L C U L D ' U N Z O O M S U R U N E N S E M B L E D E J U L I A */
/* D A N S L E P L A N C O M P L E X E : */
/* */
/* */
/* Author of '$xrc/julia.01$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 1987??????????). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 20070129112426... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 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(0.0)
#define X_ARRIVEE_DROITE \
PARE(0.0)
#define Y_ARRIVEE_BAS \
PARE(0.0)
#define Y_ARRIVEE_HAUT \
PARE(0.0)
/* Definition de la fenetre d'arrivee (dans les filaments). */
#define GENERATEUR_REEL \
PARE(0.0)
#define GENERATEUR_IMAGINAIRE \
PARE(1.0)
/* Generateur de l'ensemble : on prend "i". */
#define NOMBRE_D_ITERATIONS \
SOUS(BLANC,NOIR) \
/* Nombre d'iterations demandees. */
#define PRECISION \
PARE(0.00001) \
/* Pour calculer le rapport de reduction. */
#define VISUALISATION_ARGUMENTS_DE_SORTIE \
FAUX \
/* Faut-il visualiser l'argument de sortie ('VRAI') ou les equipotentielles ('FAUX'). */
#include xci/sequence.01.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* M A C R O S U T I L E S : */
/* */
/*************************************************************************************************************************************/
#include xrc/ITERATION.11.I"
/* Introduit le 20070108142434... */
#define CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_JULIA \
Bblock \
Test(IL_FAUT(iterations_dans_C_HC_HHC_____editer_la_fenetre_courante_avant_et_apres_homothetie)) \
Bblock \
CAL3(Prme1("fenetre carree numero %d : ",numero_d_image)); \
Eblock \
ATes \
Bblock \
Eblock \
ETes \
\
CALS(Ijulia_dans_C(ImageR \
,ADRESSE(coin_bas_gauche_de_depart_puis_courant) \
,ADRESSE(coin_haut_droite_de_depart_puis_courant) \
,ADRESSE(point_argument) \
,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_image \
,ImageR,iterations_dans_C_HC_HHC_____IFmageR \
) \
); \
CALZ_FreCC(nom_image); \
Eblock
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* C A L C U L D ' U N Z O O M S U R U N E N S E M B L E D E J U L I A */
/* D A N S L E P L A N C O M P L E X E : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
DEFV(complexe,coin_bas_gauche_de_depart_puis_courant);
/* Definition du point situe en bas a gauche de la fenetre de depart, puis de la fenetre */
/* courante au cours du zoom. */
DEFV(complexe,coin_haut_droite_de_depart_puis_courant);
/* Definition du point situe en haut a droite de la fenetre de depart, puis de la fenetre */
/* courante au cours du zoom. */
DEFV(complexe,coin_bas_gauche_d_arrivee);
/* Definition du point situe en bas a gauche de la fenetre d'arrivee, */
DEFV(complexe,coin_haut_droite_d_arrivee);
/* Definition du point situe en haut a droite de la fenetre d'arrivee. */
DEFV(complexe,point_argument);
/* Nombre complexe generateur. */
DEFV(Positive,INIT(nombre_d_iterations,NOMBRE_D_ITERATIONS));
/* Nombre maximal d'iterations a effectuer. */
DEFV(Float,INIT(precision,PRECISION));
/* Pour calculer le rapport de reduction. */
DEFV(Float,INIT(rapport_de_reduction,FLOT__UNDEF));
/* Rapport de passage d'une fenetre a l'autre. */
DEFV(Logical,INIT(visualisation_arguments_de_sortie,VISUALISATION_ARGUMENTS_DE_SORTIE));
/* Indicateur de choix entre la visualisation du nombre d'iterations ('FAUX'), */
/* et l'argument de sortie ('VRAI'). */
DEFV(Int,INIT(premiere_image,PREMIERE_IMAGE_ABSOLUE));
/* Numero de la premiere image (ce parametre a ete introduit le 20030330103500). */
DEFV(Int,INIT(nombre_d_images,NOMBRE_D_IMAGES));
/* Nombre d'images a generer. */
DEFV(Int,INIT(numero_d_image,UNDEF));
/* Numero de l'image courante (celle-ci n'est pas necessairement generee : voir */
/* le parametre 'PREMIERE_IMAGE'). */
DEFV(Int,INIT(multiplicite_du_numero_d_image,MULTIPLICITE_DU_NUMERO_D_IMAGE));
/* "Multiplicite" du numero de l'image courante. Ce nouveau parametre a ete introduit */
/* le 20030330102112 afin de pouvoir choisir iterativement les coordonnees des coins */
/* d'arrivee et ce le plus vite possible et donc en calculant le moins d'images possible. */
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 20150228111740...). */
/*..............................................................................................................................*/
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. */
Cinitialisation(point_argument,GENERATEUR_REEL,GENERATEUR_IMAGINAIRE);
/* Initialisation du nombre complexe generateur. */
GET_ARGUMENTSi(nombre_d_arguments
,BLOC(GET_ARGUMENT_C("imageR=""R=",nom_imageR);
GET_ARGUMENT_I("chiffres=",nombre_de_chiffres);
GET_ARGUMENT_I("premiere=",premiere_image);
GET_ARGUMENT_I("n=""images=""N=",nombre_d_images);
GET_ARGUMENT_C("imageI=""I=",nom_imageI);
GET_ARGUMENT_C("imageF=""F=",nom_imageF);
/* Arguments introduits le 20150228111740... */
GET_ARGUMENT_I("multiplicite=""divisible=",multiplicite_du_numero_d_image);
GET_ARGUMENT_I("iterations=",nombre_d_iterations);
GET_ARGUMENT_F("exposantJ=""exposant=",Ijulia_dans_C_____exposant);
/* Arguments introduits le 20091118121636... */
/* */
/* Le 20100608131913, passage de 'GET_ARGUMENT_I(...)' a 'GET_ARGUMENT_F(...)'. */
GET_ARGUMENT_F("seuilJ=""seuil=",Ijulia_dans_C_____seuil_du_carre_du_module);
/* Arguments introduits le 20091118183717... */
GET_ARGUMENT_L("editer=",iterations_dans_C_HC_HHC_____editer_la_fenetre_courante_avant_et_apres_homothetie);
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 20100610111603... */
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("r=""gr=",Reelle(point_argument));
GET_ARGUMENT_F("i=""gi=",Imaginaire(point_argument));
GET_ARGUMENT_F("precision=",precision);
GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie);
GET_ARGUMENT_F("ponderation_point_argument_A=""ppA=",Ijulia_dans_C_____ponderation_point_argument_A);
GET_ARGUMENT_F("ponderation_point_courant_C=""ppC=",Ijulia_dans_C_____ponderation__point_courant_C);
/* Introduits le 20130825093449... */
GET_ARGUMENT_L("gF_ITERATION_DANS_C_____compatibilite_20211015=""compatibilite_20211015="
,gF_ITERATION_DANS_C_____compatibilite_20211015
);
/* Introduit le 20211015104913. */
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 20190314184506 et complete le 20190315100802 avec la normalisation... */
/* Le 20211014082427 furent ajoutees les parametres d'encadrement des points du maillage... */
PROCESS_ARGUMENTS_DE_PARAMETRAGE_DES_ITERATIONS_1;
/* Introduit sous cette forme le 20070108142434... */
CONTROLE_DE_L_ARITHMETIQUE_FOLKLORIQUE_DES_NOMBRES_COMPLEXES;
/* Cette extension a ete introduite le 20181211112231... */
CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_COMPLEXES;
/* Cette extension a ete introduite le 20070129112426... */
CONTROLE_DU_REPLIEMENT_DE_L_ESPACE_DES_NOMBRES_COMPLEXES;
/* Cette extension a ete introduite le 20121104075753... */
)
);
EGAL(numero_d_image,PREMIERE_IMAGE_ABSOLUE);
/* Numero de l'image courante (celle-ci n'est pas necessairement generee : voir */
/* le parametre 'PREMIERE_IMAGE'). */
CALCUL_DU_RAPPORT_DE_REDUCTION;
Test(IFOU(IFNE_chaine(nom_imageI,NOM_PIPE),IFNE_chaine(nom_imageF,NOM_PIPE)))
/* Test introduit le 20150228111740 et qui permet de sauvergarder l'etat Final et de */
/* pouvoir le reutiliser ulterieurement comme etat Initial... */
Bblock
BMalloc(Ijulia_dans_C_____etat_initial,imageJ);
BMalloc(Ijulia_dans_C_____etat_final__,imageJ);
/* 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(Ijulia_dans_C_____charger_____l_etat_initial,VRAI);
Test(PAS_D_ERREUR(CODE_ERROR(IloadJ_image(Ijulia_dans_C_____etat_initial,nom_imageI))))
Bblock
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
CALi(IJinitialisation(Ijulia_dans_C_____etat_initial,ADRESSE(C_____nombre_complexe__0__0)));
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
Test(IFNE_chaine(nom_imageF,NOM_PIPE))
Bblock
EGAL(Ijulia_dans_C_____sauvegarder_l_etat_final__,VRAI);
Eblock
ATes
Bblock
Eblock
ETes
CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_JULIA;
Test(IL_FAUT(Ijulia_dans_C_____sauvegarder_l_etat_final__))
Bblock
CALi(IupdateJ_image(nom_imageF,Ijulia_dans_C_____etat_final__));
Eblock
ATes
Bblock
Eblock
ETes
EMalloc(Ijulia_dans_C_____etat_final__,imageJ);
EMalloc(Ijulia_dans_C_____etat_initial,imageJ);
Eblock
ATes
Bblock
DoIn(numero_d_image,premiere_image,LSTX(premiere_image,nombre_d_images),PAS_DES_IMAGES)
Bblock
Test(IFGE(numero_d_image,PREMIERE_IMAGE))
Bblock
/* Afin de se positionner dans la sequence... */
Test(DIVISIBLE(numero_d_image,multiplicite_du_numero_d_image))
Bblock
/* Afin de faire le moins de calculs possible... */
CALi(Inoir(ImageR));
/* Initialisation de l'image Resultat. */
CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_JULIA;
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
REDUCTION_DE_LA_FENETRE_COURANTE;
Eblock
EDoI
Eblock
ETes
RETU_Commande;
Eblock
ECommande