/*************************************************************************************************************************************/
/* */
/* C A L C U L D ' U N Z O O M S U R L ' E X P O S A N T D E L Y A P U N O V */
/* D E L ' 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.21$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 20070129112657... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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.3)
#define GENERATEUR_IMAGINAIRE \
PARE(0.3)
/* Generateur de l'ensemble : valeur arbitraire afin d'etre sur d'avoir un Julia connexe... */
#define NOMBRE_D_ITERATIONS \
CENT \
/* Nombre d'iterations demandees. */
#define RENORMALISATION_GLOBALE \
FAUX \
/* Faut-il faire une renormalisation locale a chaque image ('FAUX') ou globale ('VRAI') a */ \
/* l'ensemble des images... */
#define MODE_DE_CALCUL_DE_L_EXPOSANT_DE_LYAPUNOV \
CALCUL_DE_L_EXPOSANT_DE_LYAPUNOV_A_PARTIR_DE_LA_PARTIE_REELLE_DE_LA_DERIVEE \
/* A priori, on calculera l'exposant de Lyapunov a partir de la partie Reelle. */
#define VISUALISER_LES_EXPOSANTS_NEGATIFS \
VRAI \
/* A priori, on etudie les zones stables, */
#define VISUALISER_LES_EXPOSANTS_POSITIFS \
FAUX \
/* Et on ignore les zones chaotiques... */
#define FACTEUR_MULTIPLICATIF \
NEGA(FU) \
/* Facteur multiplicatif destine a faire une "inversion video"... */
#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 20091122182415... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* C A L C U L D ' U N Z O O M S U R L ' E X P O S A N T D E L Y A P U N O V */
/* D E L ' 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(Logical,INIT(renormalisation_globale,RENORMALISATION_GLOBALE));
/* Faut-il faire une renormalisation locale a chaque image ('FAUX') ou globale ('VRAI') a */
/* l'ensemble des images... */
DEFV(Int,INIT(mode_de_calcul_de_l_exposant_de_Lyapunov,MODE_DE_CALCUL_DE_L_EXPOSANT_DE_LYAPUNOV));
/* A priori, on calculera l'exposant de Lyapunov a partir de la partie Reelle. */
DEFV(Logical,INIT(visualiser_les_exposants_negatifs,VISUALISER_LES_EXPOSANTS_NEGATIFS));
/* A priori, on etudie les zones stables, */
DEFV(Logical,INIT(visualiser_les_exposants_positifs,VISUALISER_LES_EXPOSANTS_POSITIFS));
/* Et on ignore les zones chaotiques... */
DEFV(Float,INIT(facteur_multiplicatif,FACTEUR_MULTIPLICATIF));
/* Facteur multiplicatif destine a faire une "inversion video"... */
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. */
/*..............................................................................................................................*/
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("n=""images=""N=",nombre_d_images);
GET_ARGUMENT_I("iterations=",nombre_d_iterations);
GET_ARGUMENT_L("renormalisation_globale=""globale=",renormalisation_globale);
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_F("r=""gr=",Reelle(point_argument));
GET_ARGUMENT_F("i=""gi=",Imaginaire(point_argument));
GET_ARGUMENT_I("mode=",mode_de_calcul_de_l_exposant_de_Lyapunov);
GET_ARGUMENT_L("negatifs=",visualiser_les_exposants_negatifs);
GET_ARGUMENT_L("positifs=",visualiser_les_exposants_positifs);
GET_ARGUMENT_F("facteur=",facteur_multiplicatif);
CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_COMPLEXES;
/* Cette extension a ete introduite le 20070129112657... */
)
);
begin_nouveau_block
Bblock
DEFV(genere_Float,INIT(niveau_minimum_de_la_sequence,F_INFINI));
DEFV(genere_Float,INIT(niveau_maximum_de_la_sequence,F_MOINS_L_INFINI));
/* Afin de rechercher les niveaux extrema globaux de l'ensemble de la sequence... */
BDEFV(imageF,ensemble_de_Julia);
/* Image flottante dans laquelle on trouve la couche courante. */
CALi(Inoir(ImageR));
/* Initialisation de l'image Resultat. */
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... */
EGAL(nom_image
,chain_Aconcaten2_sauf_nom_pipe(nom_imageR,chain_numero(numero_d_image,nombre_de_chiffres))
);
Test(IL_FAUT(renormalisation_globale))
Bblock
DEFV(genere_Float,INIT(niveau_minimum_de_l_image,FLOT__NIVEAU_UNDEF));
DEFV(genere_Float,INIT(niveau_maximum_de_l_image,FLOT__NIVEAU_UNDEF));
/* Afin de rechercher les niveaux extrema de l'image courante. */
CALS(IFjulia_Lyapunov_dans_C(ensemble_de_Julia
,ADRESSE(coin_bas_gauche_de_depart_puis_courant)
,ADRESSE(coin_haut_droite_de_depart_puis_courant)
,ADRESSE(point_argument)
,nombre_d_iterations
,mode_de_calcul_de_l_exposant_de_Lyapunov
,visualiser_les_exposants_negatifs
,visualiser_les_exposants_positifs
,facteur_multiplicatif
)
);
CALS(IFnivo_extrema(ensemble_de_Julia
,ADRESSE(niveau_minimum_de_l_image)
,ADRESSE(niveau_maximum_de_l_image)
)
);
/* Recherche des extrema locaux a l'image courante. */
EGAL(niveau_minimum_de_la_sequence
,MIN2(niveau_minimum_de_la_sequence
,niveau_minimum_de_l_image
)
);
EGAL(niveau_maximum_de_la_sequence
,MAX2(niveau_maximum_de_la_sequence
,niveau_maximum_de_l_image
)
);
/* Recherche des extrema globaux a la sequence entiere. */
CALi(IupdateF_image(nom_image,ensemble_de_Julia));
Eblock
ATes
Bblock
CALS(Ijulia_Lyapunov_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
,mode_de_calcul_de_l_exposant_de_Lyapunov
,visualiser_les_exposants_negatifs
,visualiser_les_exposants_positifs
,facteur_multiplicatif
)
);
CALi(Iupdate_image(nom_image,ImageR));
Eblock
ETes
CALZ_FreCC(nom_image);
Eblock
ATes
Bblock
Eblock
ETes
REDUCTION_DE_LA_FENETRE_COURANTE;
Eblock
EDoI
Test(IL_FAUT(renormalisation_globale))
Bblock
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))
/* Afin de se positionner dans la sequence... */
Bblock
EGAL(nom_image
,chain_Aconcaten2_sauf_nom_pipe(nom_imageR,chain_numero(numero_d_image,nombre_de_chiffres))
);
Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(ensemble_de_Julia,nom_image))))
Bblock
CALS(Ifloat_std(ImageR
,ensemble_de_Julia
,niveau_minimum_de_la_sequence
,niveau_maximum_de_la_sequence
)
);
/* Renormalisation globale des images de la sequence... */
CALi(Iupdate_image(nom_image,ImageR));
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
CALZ_FreCC(nom_image);
Eblock
ATes
Bblock
Eblock
ETes
Eblock
EDoI
Eblock
ATes
Bblock
Eblock
ETes
EDEFV(imageF,ensemble_de_Julia);
/* Image flottante dans laquelle on trouve la couche courante. */
Eblock
end_nouveau_block
RETU_Commande;
Eblock
ECommande