/*************************************************************************************************************************************/ /* */ /* C A L C U L D E L ' E N S E M B L E D E M A N D E L B R O T */ /* D A N S L ' E S P A C E H Y P E R - C O M P L E X E : */ /* */ /* */ /* Author of '$xrc/mandel.41$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 20070129112331... */ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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 (introduit le 20091121115929) ? */ \ /* La valeur par defaut ('VRAI') est logique pour un calcul 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(-2.0) #define X_DEPART_DROITE \ PARE(0.5) #define Y_DEPART_BAS \ PARE(-1.25) #define Y_DEPART_HAUT \ PARE(1.25) /* 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 /* Definition des composantes 'J' et 'K' des quaternions. */ #define NOMBRE_D_ITERATIONS \ HUIT \ /* Nombre d'iterations demandees. */ #define VISUALISATION_ARGUMENTS_DE_SORTIE \ FAUX #define PONDERATION_ARGUMENTT \ FU #define PONDERATION_ARGUMENTP \ FZERO #define PONDERATION_ARGUMENTA \ FZERO /* Faut-il visualiser les arguments de sortie ('VRAI'), auquel cas les trois parametres */ /* suivants sont les ponderations des 3 angles {T,P,A}, 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" /* Introduit le 20091120123336... */ #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. */ #define CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_MANDELBROT \ Bblock \ 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); \ /* Calcul de la composante 'J' des quaternions lorsqu'elle est variable. */ \ vTRANSFORMAT_31(composante_K,sCOMPOSANTE_K,numero_de_la_periode_courante,fichier_COMPOSANTE_K); \ /* Calcul de la composante 'K' des quaternions lorsqu'elle est variable. */ \ \ CALS(Imandelbrot_dans_HC(ImageR \ ,ADRESSE(coin_bas_gauche_de_depart_puis_courant) \ ,ADRESSE(coin_haut_droite_de_depart_puis_courant) \ ,composante_J \ ,composante_K \ ,nombre_d_iterations \ ,visualisation_arguments_de_sortie \ ,ponderation_argumentT,ponderation_argumentP,ponderation_argumentA \ ) \ ); \ /* Calcul de l'ensemble de Mandelbrot dans l'espace hyper-Complexe. */ \ \ GENERATION_D_UN_TYPE_D_IMAGE(VRAI \ ,nom_imageR \ ,ImageR \ ); \ /* Generation des images demandees. */ \ Eblock /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* D E F I N I T I O N S D E S R O T A T I O N S T R I D I M E N S I O N N E L L E S : */ /* */ /*************************************************************************************************************************************/ #include xrc/julia.41.I" /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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 : */ /* */ /*************************************************************************************************************************************/ 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) /* Definition du fichier des composantes 'J' et 'K' des quaternions. */ #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) /* Formule generale definissant les variations des composantes 'J' et 'K' des quaternions. */ #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 M A N D E L B R O T */ /* D A N S L ' E S P A C E 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 (introduit le 20091121115929) ? */ /* La valeur par defaut ('VRAI') est logique pour un calcul 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)); /* Definition des composantes 'J' et 'K' des quaternions. */ 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_argumentT,PONDERATION_ARGUMENTT)); DEFV(Float,INIT(ponderation_argumentP,PONDERATION_ARGUMENTP)); DEFV(Float,INIT(ponderation_argumentA,PONDERATION_ARGUMENTA)); /* Faut-il visualiser les arguments de sortie ('VRAI'), auquel cas les trois parametres */ /* suivants sont les ponderations des 3 angles {T,P,A}, 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. */ 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 20150228115203...). */ /*..............................................................................................................................*/ 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. */ DEBUT_DE_L_INITIALISATION_DES_INTERPOLATIONS; FIN_DE_L_INITIALISATION_DES_INTERPOLATIONS; /* Initialisations generales pour l'interpolation des rotations... */ 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); /* Initialisation des composantes 'J' et 'K' des quaternions. */ 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_ARGUMENTS_GEOMETRIQUES_DE_ROTATION; PROCESS_ARGUMENTS_GEOMETRIQUES_DE_ROTATION_2_HC; /* Introduit le 20091222140116... */ 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 ); GET_ARGUMENT_C("imageR=""R=",nom_imageR); GET_ARGUMENT_I("chiffres=",nombre_de_chiffres); GET_ARGUMENT_C("imageI=""I=",nom_imageI); GET_ARGUMENT_C("imageF=""F=",nom_imageF); /* Arguments introduits le 20150228115203... */ GET_ARGUMENT_I("iterations=",nombre_d_iterations); GET_ARGUMENT_F("exposantM=""exposant=",Imandelbrot_dans_HC_____exposant); /* Arguments introduits le 20091118121214... */ /* */ /* Le 20100608132532, passage de 'GET_ARGUMENT_I(...)' a 'GET_ARGUMENT_F(...)'. */ GET_ARGUMENT_F("seuilM=""seuil=",Imandelbrot_dans_HC_____seuil_du_carre_du_module); /* Arguments introduits le 20091118183519... */ GET_ARGUMENT_L("fenetres_depart_arrivee_identiques=""fdai=" ,les_fenetres_de_depart_et_d_arrivee_sont_identiques ); /* Arguments introduits le 20091121115929... */ 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)); 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)); /* Introduction des arguments {"xbg=","ybg=","xhd=","yhd="} le 20091122120041... */ /* */ /* Les arguments suivants de 'v $xrc/mandel.41$K' : */ /* */ /* OX=... OY=... EX=... EY=... */ /* fenetres_depart_arrivee_identiques=VRAI */ /* editer_fenetre_courante=VRAI */ /* */ /* s'utiliseront conjointement avec les arguments suivants de 'v $xci/display$K' : */ /* */ /* ModeFlot='f' */ /* ChiffresFlot=... */ /* */ /* editer_XY_curseur=VRAI */ /* normaliser_XY_curseur=VRAI */ /* DX_curseur=... DY_curseur=... */ /* OX=... OY=... EX=... EY=... */ /* */ /* Supposons que '$xrc/mandel.41$X' en soit a : */ /* */ /* xbg=+0.0014 xhd=+0.0016 ybg=+0.0045 yhd=+0.0047 */ /* --- --- --- --- */ /* */ /* Supposons qu'alors '$xci/display$X' renvoie : */ /* */ /* curseur=(+0.00153333,+0.00467777) */ /* ----- ----- */ /* */ /* on aura alors interet a choisir pour l'etape suivante de '$xrc/mandel.41$X' : */ /* */ /* xbg=+0.00152 xhd=+0.00154 ybg=+0.00466 yhd=+0.00468 */ /* ---- ---- ---- ---- */ /* */ /* et ainsi {+0.00152,+0.00154} sont de part et d'autre de +0.00153, de meme que */ /* ---- ---- ---- */ /* {+0.00466,+0.00468} sont de part et d'autre de +0.00467 (dans les deux cas on fait */ /* ---- ---- ---- */ /* donc moins et plus 0.00001). */ GET_ARGUMENT_F("precision=",precision); GET_ARGUMENT_F("J=",composante_J); GET_ARGUMENT_F("K=",composante_K); GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie); GET_ARGUMENT_F("ponderationT=",ponderation_argumentT); GET_ARGUMENT_F("ponderationP=",ponderation_argumentP); GET_ARGUMENT_F("ponderationA=",ponderation_argumentA); PROCESS_ARGUMENTS_DE_PARAMETRAGE_DES_ITERATIONS_1; /* Introduit le 20091120122821... */ GET_ARGUMENT_L("F_ITERATION_DANS_HC_____compatibilite_20110825=""compatibilite_20110825=" ,F_ITERATION_DANS_HC_____compatibilite_20110825 ); /* Introduit le 20110826123121. */ CONTROLE_DE_L_ARITHMETIQUE_FOLKLORIQUE_DES_NOMBRES_COMPLEXES; /* Cette extension a ete introduite le 20181211112338... */ CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_HYPER_COMPLEXES; /* Cette extension a ete introduite le 20070129112331... */ CONTROLE_DE_LA_TRANFORMATION_DE_LA_SUITE_DE_POINTS_DANS_HC; /* Cette extension a ete introduite le 20100513181745... */ CONTROLE_DU_REPLIEMENT_DE_L_ESPACE_DES_NOMBRES_HYPER_COMPLEXES; /* Cette extension a ete introduite le 20121104080422... */ ) ); Test(EST_VRAI(les_fenetres_de_depart_et_d_arrivee_sont_identiques)) /* Test introduit le 20091121115929... */ 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; Test(IFOU(IFNE_chaine(nom_imageI,NOM_PIPE),IFNE_chaine(nom_imageF,NOM_PIPE))) /* Test introduit le 20150228120012 et qui permet de sauvergarder l'etat Final et de */ /* pouvoir le reutiliser ulterieurement comme etat Initial... */ Bblock BMalloc(Imandelbrot_dans_HC_____etat_initial,imageHJ); BMalloc(Imandelbrot_dans_HC_____etat_final__,imageHJ); /* 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(Imandelbrot_dans_HC_____charger_____l_etat_initial,VRAI); Test(PAS_D_ERREUR(CODE_ERROR(IloadHJ_image(Imandelbrot_dans_HC_____etat_initial,nom_imageI)))) Bblock Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; CALi(IHJinitialisation(Imandelbrot_dans_HC_____etat_initial,ADRESSE(HC_____quaternion__0__0__0__0))); Eblock ETes Eblock ATes Bblock Eblock ETes Test(IFNE_chaine(nom_imageF,NOM_PIPE)) Bblock EGAL(Imandelbrot_dans_HC_____sauvegarder_l_etat_final__,VRAI); Eblock ATes Bblock Eblock ETes CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_MANDELBROT; Test(IL_FAUT(Imandelbrot_dans_HC_____sauvegarder_l_etat_final__)) Bblock CALi(IupdateHJ_image(nom_imageF,Imandelbrot_dans_HC_____etat_final__)); Eblock ATes Bblock Eblock ETes EMalloc(Imandelbrot_dans_HC_____etat_final__,imageJ); EMalloc(Imandelbrot_dans_HC_____etat_initial,imageJ); Eblock ATes Bblock 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 CHANGEMENT_DU_POINT_DE_VUE_GEOMETRIQUE; /* Mise en place des rotations tridimensionnelles dans les Quaternions. */ 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... */ CALCUL_ET_ENREGISTREMENT_D_UNE_IMAGE_DE_L_ENSEMBLE_DE_MANDELBROT; Eblock ATes Bblock Eblock ETes REDUCTION_DE_LA_FENETRE_COURANTE; Eblock EDoI Eblock ETes RETU_Commande; Eblock ECommande