/*************************************************************************************************************************************/ /* */ /* C A L C U L D E L ' I T E R A T I O N D E P O L Y N O M E */ /* D A N S L ' E S P A C E H Y P E R - H Y P E R - C O M P L E X E : */ /* */ /* */ /* Author of '$xrc/polynome.81$K' : */ /* */ /* Jean-Francois COLONNA (LACTAMME, 20111118212316). */ /* */ /*************************************************************************************************************************************/ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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 PRAGMA_CL_____MODULE_NON_OPTIMISABLE /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* F I C H I E R S D ' I N C L U D E S : */ /* */ /*************************************************************************************************************************************/ #include INCLUDES_BASE #include image_image_IMAGESF_EXT #include maths_compl_fonct_ITERATIONS_EXT #include image_image_QUAD_IMAGE_EXT #include image_image_DI_ALBUM_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 : */ /* */ /*************************************************************************************************************************************/ #include xci/sequence.01.I" /* ATTENTION, on definit ainsi le symbole 'DERNIERE_IMAGE' qui ne sert a rien ici, puisque */ /* c'est en effet 'Zmax' qui joue ce role... */ =define PREMIERE_IMAGE_DE_L_ALBUM \ PREMIERE_IMAGE #undef PREMIERE_IMAGE /* Ceci vient du fait que 'v $xci/sequence.01$I PREMIERE_IMAGE' definit le parametre */ /* 'PREMIERE_IMAGE' (='PREMIERE_IMAGE_ABSOLUE'='ZERO'), alors qu'un peu plus loin, il */ /* convient de le redefinir (='NUMERO_DE_LA_PREMIERE_PERIODE_DE_LA_SIMULATION' */ /* ='PREMIERE_ITERATION_D_UN_Komp'='UN'). */ #define AVERTIR_S_IL_Y_A_PARAMETRAGE_DU_PRODUIT_NON_THEORIQUE_AVEC_PLUS_D_UNE_IMAGE_CALCULEE \ VRAI \ /* Indicateur introduit ne sachant pas comment faire autrement... */ #define LES_FENETRES_DE_DEPART_ET_D_ARRIVEE_SONT_IDENTIQUES \ FAUX \ /* Les fenetres de Depart et d'Arrivee sont-elles identiques ? */ \ /* */ \ /* La valeur par defaut ('VRAI') est logique pour un calcul hyper-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(-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. */ #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 #define COMPOSANTE_S \ FZERO #define COMPOSANTE_T \ FZERO #define COMPOSANTE_U \ FZERO #define COMPOSANTE_V \ FZERO /* Definition des composantes {J,K,S,T,U,V} des octonions. */ #define GENERATEUR_REEL \ PARE(0.0) #define GENERATEUR_IMAGINAIRE \ PARE(1.0) #define GENERATEUR_JMAGINAIRE \ PARE(0.0) #define GENERATEUR_KMAGINAIRE \ PARE(0.0) #define GENERATEUR_SMAGINAIRE \ PARE(0.0) #define GENERATEUR_TMAGINAIRE \ PARE(0.0) #define GENERATEUR_UMAGINAIRE \ PARE(0.0) #define GENERATEUR_VMAGINAIRE \ PARE(0.0) /* Generateur de l'ensemble : on prend "i". */ #define PONDERATION_DU_POINT_ARGUMENT \ FZERO #define PONDERATION_DU_POINT_COURANT \ FU /* Ponderations initialisees sur un ensemble de type "Mandelbrot"... */ #define SEUIL_DU_CARRE_DU_MODULE \ EXP2(DOUB(UNITE)) \ /* Seuil du carre du module definissant l'arret des iterations. */ #define NOMBRE_D_ITERATIONS \ HUIT \ /* Nombre d'iterations demandees. */ #define VISUALISATION_ARGUMENTS_DE_SORTIE \ FAUX #define PONDERATION_ARGUMENTA1 \ FU #define PONDERATION_ARGUMENTA2 \ FZERO #define PONDERATION_ARGUMENTA3 \ FZERO #define PONDERATION_ARGUMENTA4 \ FZERO #define PONDERATION_ARGUMENTA5 \ FZERO #define PONDERATION_ARGUMENTA6 \ FZERO #define PONDERATION_ARGUMENTA7 \ 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 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 : */ /* */ /*************************************************************************************************************************************/ #define definition_de_exposa1_Rho \ IFmageA1 #define definition_de_ponder1_A1_ \ IFmageA2 #define definition_de_ponder1_A2_ \ IFmageA3 #define definition_de_ponder1_A3_ \ IFmageA4 #define definition_de_ponder1_A4_ \ IFmageA5 #define definition_de_ponder1_A5_ \ IFmageA6 #define definition_de_ponder1_A6_ \ IFmageA7 #define definition_de_ponder1_A7_ \ IFmageA8 #include xrc/ITERATION.11.I" /* Introduit le 20070108143853... */ #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(gIupdate_image(NOTL(iterations_dans_C_HC_HHC_____etendre_la_precision) \ ,nom_image \ ,image,iterations_dans_C_HC_HHC_____IFmageR \ ) \ ); \ CALZ_FreCC(nom_image); \ Eblock \ ATes \ Bblock \ Eblock \ ETes \ Eblock \ /* Procedure generant si necessaire un certain type d'image. */ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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) dfTRANSFORMAT_31(liste_COMPOSANTE_S,fichier_COMPOSANTE_S,COMPOSANTE_S_IMPLICITE,COMPOSANTE_S) dfTRANSFORMAT_31(liste_COMPOSANTE_T,fichier_COMPOSANTE_T,COMPOSANTE_T_IMPLICITE,COMPOSANTE_T) dfTRANSFORMAT_31(liste_COMPOSANTE_U,fichier_COMPOSANTE_U,COMPOSANTE_U_IMPLICITE,COMPOSANTE_U) dfTRANSFORMAT_31(liste_COMPOSANTE_V,fichier_COMPOSANTE_V,COMPOSANTE_V_IMPLICITE,COMPOSANTE_V) /* Definition du fichier des composantes {J,K,S,T,U,V} des octonions. */ #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) #define sCOMPOSANTE_S(numero_de_la_periode) \ sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_S) #define sCOMPOSANTE_T(numero_de_la_periode) \ sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_T) #define sCOMPOSANTE_U(numero_de_la_periode) \ sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_U) #define sCOMPOSANTE_V(numero_de_la_periode) \ sTRANSFORMAT_31(numero_de_la_periode,liste_COMPOSANTE_V) /* Formule generale definissant les variations des composantes {J,K,S,T,U,V} des octonions. */ #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 ' I T E R A T I O N D E P O L Y N O M E */ /* D A N S L ' E S P A C E H Y P E R - H Y P E R - C O M P L E X E : */ /* */ /*************************************************************************************************************************************/ BCommande(nombre_d_arguments,arguments) /*-----------------------------------------------------------------------------------------------------------------------------------*/ Bblock DEFV(CHAR,INIC(POINTERc(nom_imageA),NOM_PIPE)); /* Nom de l'album definissant le polynome... */ DEFV(CHAR,INIC(POINTERc(nom_postfixe),NOM_UNDEF_VIDE)); /* Nom d'un eventuel postfixe a placer derriere <nom_imageA><numero> (par exemple '$ROUGE'). */ DEFV(Int,INIT(premiere_image_de_l_album,PREMIERE_IMAGE_DE_L_ALBUM)); /* Numero de la premiere image, */ DEFV(Int,INIT(pas_des_images_de_l_album,PAS_DES_IMAGES)); /* Pas de passage d'un numero d'image a une autre. */ DEFV(Int,INIT(nombre_de_chiffres_de_l_album,NOMBRE_DE_CHIFFRES)); /* Nombre de chiffres codant le numero des images de la serie... */ DEFV(Logical,INIT(avertir_s_il_y_a_parametrage_du_produit_non_theorique_avec_plus_d_une_image_calculee ,AVERTIR_S_IL_Y_A_PARAMETRAGE_DU_PRODUIT_NON_THEORIQUE_AVEC_PLUS_D_UNE_IMAGE_CALCULEE ) ); /* Indicateur introduit ne sachant pas comment faire autrement... */ DEFV(CHAR,INIC(POINTERc(nom_definition_de_exposa1_Rho),NOM_PIPE)); DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_A1_),NOM_PIPE)); DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_A2_),NOM_PIPE)); DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_A3_),NOM_PIPE)); DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_A4_),NOM_PIPE)); DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_A5_),NOM_PIPE)); DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_A6_),NOM_PIPE)); DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_A7_),NOM_PIPE)); /* Parametrage du "faux produit". */ 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 ? */ /* */ /* La valeur par defaut ('VRAI') est logique pour un calcul hyper-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)); DEFV(Float,INIT(composante_S,COMPOSANTE_S)); DEFV(Float,INIT(composante_T,COMPOSANTE_T)); DEFV(Float,INIT(composante_U,COMPOSANTE_U)); DEFV(Float,INIT(composante_V,COMPOSANTE_V)); /* Definition des composantes {J,K,S,T,U,V} des octonions. */ DEFV(hyper_hyper_complexe,point_argument); /* Nombre hyper-hyper-complexe generateur. */ DEFV(Float,INIT(ponderation_du_point_argument,PONDERATION_DU_POINT_ARGUMENT)); DEFV(Float,INIT(ponderation_du_point_courant,PONDERATION_DU_POINT_COURANT)); /* Ponderations initialisees sur un ensemble de typ "Mandelbrot"... */ DEFV(Positive,INIT(nombre_d_iterations,NOMBRE_D_ITERATIONS)); /* Nombre maximal d'iterations a effectuer. */ DEFV(Float,INIT(seuil_du_carre_du_module,SEUIL_DU_CARRE_DU_MODULE)); /* Seuil du carre du module definissant l'arret des iterations. */ DEFV(Logical,INIT(visualisation_arguments_de_sortie,VISUALISATION_ARGUMENTS_DE_SORTIE)); DEFV(Float,INIT(ponderation_argumentA1,PONDERATION_ARGUMENTA1)); DEFV(Float,INIT(ponderation_argumentA2,PONDERATION_ARGUMENTA2)); DEFV(Float,INIT(ponderation_argumentA3,PONDERATION_ARGUMENTA3)); DEFV(Float,INIT(ponderation_argumentA4,PONDERATION_ARGUMENTA4)); DEFV(Float,INIT(ponderation_argumentA5,PONDERATION_ARGUMENTA5)); DEFV(Float,INIT(ponderation_argumentA6,PONDERATION_ARGUMENTA6)); DEFV(Float,INIT(ponderation_argumentA7,PONDERATION_ARGUMENTA7)); /* Faut-il visualiser les arguments de sortie ('VRAI'), auquel cas les trois parametres */ /* suivants sont les ponderations des 7 angles {A1,A2,A3,A4,A5,A6,A7} 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. */ /*..............................................................................................................................*/ 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. */ HHCinitialisation(point_argument ,GENERATEUR_REEL ,GENERATEUR_IMAGINAIRE ,GENERATEUR_JMAGINAIRE ,GENERATEUR_KMAGINAIRE ,GENERATEUR_SMAGINAIRE ,GENERATEUR_TMAGINAIRE ,GENERATEUR_UMAGINAIRE ,GENERATEUR_VMAGINAIRE ); /* Initialisation du nombre hyper-hyper-complexe generateur. */ 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 octonions. */ 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_HHC; /* Introduit le 20091222140039... */ 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 ); PROCESS_ARGUMENT_FICHIER("COMPOSANTE_S=" ,fichier_COMPOSANTE_S ,liste_COMPOSANTE_S ,COMPOSANTE_S_IMPLICITE ,gTRANSFORMAT_31 ); PROCESS_ARGUMENT_FICHIER("COMPOSANTE_T=" ,fichier_COMPOSANTE_T ,liste_COMPOSANTE_T ,COMPOSANTE_T_IMPLICITE ,gTRANSFORMAT_31 ); PROCESS_ARGUMENT_FICHIER("COMPOSANTE_U=" ,fichier_COMPOSANTE_U ,liste_COMPOSANTE_U ,COMPOSANTE_U_IMPLICITE ,gTRANSFORMAT_31 ); PROCESS_ARGUMENT_FICHIER("COMPOSANTE_V=" ,fichier_COMPOSANTE_V ,liste_COMPOSANTE_V ,COMPOSANTE_V_IMPLICITE ,gTRANSFORMAT_31 ); GET_ARGUMENT_C("imageA=""A=",nom_imageA); GET_ARGUMENT_C("postfixeA=",nom_postfixe); GET_ARGUMENT_I("premiereA=",premiere_image_de_l_album); GET_ARGUMENT_I("pasA=",pas_des_images_de_l_album); GET_ARGUMENT_I("chiffresA=",nombre_de_chiffres_de_l_album); GET_ARGUMENT_I("premier_exposant_polynome=""pep=""pe=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____premier_exposant_du_polynome ); GET_ARGUMENT_L("avertir_parametrage_produit_plus_d_une_image=""avertir=""apppui=" ,avertir_s_il_y_a_parametrage_du_produit_non_theorique_avec_plus_d_une_image_calculee ); GET_ARGUMENT_C("imageAeRho=""AeRho=""AeR=",nom_definition_de_exposa1_Rho); GET_ARGUMENT_L("partie_entiere_ApA1=""peApA1=""iApA1=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____prendre_la_partie_entiere_de_ponder1_A1_ ); GET_ARGUMENT_C("imageApA1=""ApA1=",nom_definition_de_ponder1_A1_); GET_ARGUMENT_L("partie_entiere_ApA2=""peApA2=""iApA2=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____prendre_la_partie_entiere_de_ponder1_A2_ ); GET_ARGUMENT_C("imageApA2=""ApA2=",nom_definition_de_ponder1_A2_); GET_ARGUMENT_L("partie_entiere_ApA3=""peApA3=""iApA3=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____prendre_la_partie_entiere_de_ponder1_A3_ ); GET_ARGUMENT_C("imageApA3=""ApA3=",nom_definition_de_ponder1_A3_); GET_ARGUMENT_L("partie_entiere_ApA4=""peApA4=""iApA4=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____prendre_la_partie_entiere_de_ponder1_A4_ ); GET_ARGUMENT_C("imageApA4=""ApA4=",nom_definition_de_ponder1_A4_); GET_ARGUMENT_L("partie_entiere_ApA5=""peApA5=""iApA5=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____prendre_la_partie_entiere_de_ponder1_A5_ ); GET_ARGUMENT_C("imageApA5=""ApA5=",nom_definition_de_ponder1_A5_); GET_ARGUMENT_L("partie_entiere_ApA6=""peApA6=""iApA6=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____prendre_la_partie_entiere_de_ponder1_A6_ ); GET_ARGUMENT_C("imageApA6=""ApA6=",nom_definition_de_ponder1_A6_); GET_ARGUMENT_L("partie_entiere_ApA7=""peApA7=""iApA7=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____prendre_la_partie_entiere_de_ponder1_A7_ ); GET_ARGUMENT_C("imageApA7=""ApA7=",nom_definition_de_ponder1_A7_); GET_ARGUMENT_C("imageR=""R=",nom_imageR); GET_ARGUMENT_I("chiffresR=",nombre_de_chiffres); GET_ARGUMENT_I("iterations=",nombre_d_iterations); GET_ARGUMENT_L("fenetres_depart_arrivee_identiques=""fdai=" ,les_fenetres_de_depart_et_d_arrivee_sont_identiques ); 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 20091122120116... */ GET_ARGUMENT_F("precision=",precision); GET_ARGUMENT_F("J=",composante_J); GET_ARGUMENT_F("K=",composante_K); GET_ARGUMENT_F("S=",composante_S); GET_ARGUMENT_F("T=",composante_T); GET_ARGUMENT_F("U=",composante_U); GET_ARGUMENT_F("V=",composante_V); GET_ARGUMENT_F("r=""gr=",HHReelle(point_argument)); GET_ARGUMENT_F("i=""gi=",HHImaginaire(point_argument)); GET_ARGUMENT_F("j=""gj=",HHJmaginaire(point_argument)); GET_ARGUMENT_F("k=""gk=",HHKmaginaire(point_argument)); GET_ARGUMENT_F("s=""gs=",HHSmaginaire(point_argument)); GET_ARGUMENT_F("t=""gt=",HHTmaginaire(point_argument)); GET_ARGUMENT_F("u=""gu=",HHUmaginaire(point_argument)); GET_ARGUMENT_F("v=""gv=",HHVmaginaire(point_argument)); GET_ARGUMENT_L("rotation_rijkstuv=""Rrijkstuv=""RA=" ,Iiteration_polynome_variable_quelconque_dans_HHC_____faire_la_rotation_de__nombre_argument_A ); GET_ARGUMENT_F("PJulia=""pondPA=""alpha=",ponderation_du_point_argument); GET_ARGUMENT_F("PMandelbrot=""pondPC=""beta=",ponderation_du_point_courant); /* Le 20091211141615, les arguments {"pA=","pC="} ont du etre remplaces, pour des raisons */ /* d'ambiguite, par {"pondPA=","pondPC="}. */ /* */ /* Le 20150305205746 {"PJulia=","PMandelbrot="} ont ete introduits... */ GET_ARGUMENT_F("seuil=",seuil_du_carre_du_module); GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie); GET_ARGUMENT_F("ponderationA1=",ponderation_argumentA1); GET_ARGUMENT_F("ponderationA2=",ponderation_argumentA2); GET_ARGUMENT_F("ponderationA3=",ponderation_argumentA3); GET_ARGUMENT_F("ponderationA4=",ponderation_argumentA4); GET_ARGUMENT_F("ponderationA5=",ponderation_argumentA5); GET_ARGUMENT_F("ponderationA6=",ponderation_argumentA6); GET_ARGUMENT_F("ponderationA7=",ponderation_argumentA7); PROCESS_ARGUMENTS_DE_PARAMETRAGE_DES_ITERATIONS_1; /* Introduit sous cette forme le 20070108143853... */ GET_ARGUMENT_L("F_ITERATION_DANS_HHC_____compatibilite_20110825=""compatibilite_20110825=" ,F_ITERATION_DANS_HHC_____compatibilite_20110825 ); /* Introduit le 20110826123224. */ CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_HYPER_HYPER_COMPLEXES; /* Le 20121103105216, 'CONTROLE_DE_LA_TRANFORMATION_DE_LA_SUITE_DE_POINTS_DANS_HC' a ete */ /* supprime puisqu'il concernait 'HC' et non pas 'HHC'... */ CONTROLE_DU_REPLIEMENT_DE_L_ESPACE_DES_NOMBRES_HYPER_HYPER_COMPLEXES; /* Cette extension a ete introduite le 20121104081143... */ ) ); Test(EST_VRAI(les_fenetres_de_depart_et_d_arrivee_sont_identiques)) 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 begin_nouveau_block Bblock BDEFV(albumF,definition_du_polynome); /* Definition de l'album Argument definissant le polynome... */ CALi(dAloadF_album(definition_du_polynome ,nom_imageA ,nom_postfixe ,premiere_image_de_l_album ,pas_des_images_de_l_album ,nombre_de_chiffres_de_l_album ,FAUX ,UNDEF ,FAUX ) ); /* Chargement de l'album Argument... */ Test(PAS_D_ERREUR(CODE_ERREUR)) Bblock Test(IL_FAUT(FHHCproduit_____utiliser_la_definition_theorique_du_produit)) Bblock Test(I8OU(IFNE_chaine(nom_definition_de_exposa1_Rho,NOM_PIPE) ,IFNE_chaine(nom_definition_de_ponder1_A1_,NOM_PIPE) ,IFNE_chaine(nom_definition_de_ponder1_A2_,NOM_PIPE) ,IFNE_chaine(nom_definition_de_ponder1_A3_,NOM_PIPE) ,IFNE_chaine(nom_definition_de_ponder1_A4_,NOM_PIPE) ,IFNE_chaine(nom_definition_de_ponder1_A5_,NOM_PIPE) ,IFNE_chaine(nom_definition_de_ponder1_A6_,NOM_PIPE) ,IFNE_chaine(nom_definition_de_ponder1_A7_,NOM_PIPE) ) ) Bblock PRINT_ATTENTION("le parametrage du produit non theorique n'a pas de sens ici et est donc ignore"); Eblock ATes Bblock Eblock ETes Eblock ATes Bblock Test(IFET(IL_FAUT(avertir_s_il_y_a_parametrage_du_produit_non_theorique_avec_plus_d_une_image_calculee) ,IFGT(nombre_d_images,UN) ) ) Bblock PRINT_ATTENTION("le parametrage du produit non theorique n'a de sens que pour le calcul d'une image unique"); Eblock ATes Bblock Eblock ETes CALi(IFinitialisation(definition_de_exposa1_Rho,FU)); CALi(IFinitialisation(definition_de_ponder1_A1_,FU)); CALi(IFinitialisation(definition_de_ponder1_A2_,FU)); CALi(IFinitialisation(definition_de_ponder1_A3_,FU)); CALi(IFinitialisation(definition_de_ponder1_A4_,FU)); CALi(IFinitialisation(definition_de_ponder1_A5_,FU)); CALi(IFinitialisation(definition_de_ponder1_A6_,FU)); CALi(IFinitialisation(definition_de_ponder1_A7_,FU)); /* A priori, les valeurs les plus probables sont forcees... */ /* */ /* On notera le 20210318160046 une grosse erreur : les initialisations de {A4,A5,A6,A7} */ /* manquaient depuis toujours... */ Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_exposa1_Rho ,nom_definition_de_exposa1_Rho ) ) ) ) Bblock Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_A1_ ,nom_definition_de_ponder1_A1_ ) ) ) ) Bblock Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_A2_ ,nom_definition_de_ponder1_A2_ ) ) ) ) Bblock Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_A3_ ,nom_definition_de_ponder1_A3_ ) ) ) ) Bblock Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_A4_ ,nom_definition_de_ponder1_A4_ ) ) ) ) Bblock Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_A5_ ,nom_definition_de_ponder1_A5_ ) ) ) ) Bblock Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_A6_ ,nom_definition_de_ponder1_A6_ ) ) ) ) Bblock Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_A7_ ,nom_definition_de_ponder1_A7_ ) ) ) ) Bblock Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; Eblock ETes Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; Eblock ETes Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; Eblock ETes Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; Eblock ETes Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; Eblock ETes Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; Eblock ETes Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; Eblock ETes Eblock ATes Bblock Test__CODE_ERREUR__ERREUR07; Eblock ETes Eblock ETes CALCUL_DU_RAPPORT_DE_REDUCTION; 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... */ 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 octonions lorsqu'elle est variable. */ vTRANSFORMAT_31(composante_K,sCOMPOSANTE_K,numero_de_la_periode_courante,fichier_COMPOSANTE_K); /* Calcul de la composante 'K' des octonions lorsqu'elle est variable. */ vTRANSFORMAT_31(composante_S,sCOMPOSANTE_S,numero_de_la_periode_courante,fichier_COMPOSANTE_S); /* Calcul de la composante 'S' des octonions lorsqu'elle est variable. */ vTRANSFORMAT_31(composante_T,sCOMPOSANTE_T,numero_de_la_periode_courante,fichier_COMPOSANTE_T); /* Calcul de la composante 'T' des octonions lorsqu'elle est variable. */ vTRANSFORMAT_31(composante_U,sCOMPOSANTE_U,numero_de_la_periode_courante,fichier_COMPOSANTE_U); /* Calcul de la composante 'U' des octonions lorsqu'elle est variable. */ vTRANSFORMAT_31(composante_V,sCOMPOSANTE_V,numero_de_la_periode_courante,fichier_COMPOSANTE_V); /* Calcul de la composante 'V' des octonions lorsqu'elle est variable. */ CALS(Iiteration_polynome_variable_quelconque_dans_HHC(ImageR ,definition_du_polynome ,definition_de_exposa1_Rho ,definition_de_ponder1_A1_ ,definition_de_ponder1_A2_ ,definition_de_ponder1_A3_ ,definition_de_ponder1_A4_ ,definition_de_ponder1_A5_ ,definition_de_ponder1_A6_ ,definition_de_ponder1_A7_ ,ponderation_du_point_argument ,ADRESSE(point_argument) ,ponderation_du_point_courant ,ADRESSE(coin_bas_gauche_de_depart_puis_courant) ,ADRESSE(coin_haut_droite_de_depart_puis_courant) ,composante_J ,composante_K ,composante_S ,composante_T ,composante_U ,composante_V ,nombre_d_iterations ,seuil_du_carre_du_module ,visualisation_arguments_de_sortie ,ponderation_argumentA1 ,ponderation_argumentA2 ,ponderation_argumentA3 ,ponderation_argumentA4 ,ponderation_argumentA5 ,ponderation_argumentA6 ,ponderation_argumentA7 ) ); /* Calcul de l'ensemble de Julia dans l'espace hyper-hyper-Complexe. */ GENERATION_D_UN_TYPE_D_IMAGE(VRAI ,nom_imageR ,ImageR ); /* Generation des images demandees. */ Eblock ATes Bblock Eblock ETes REDUCTION_DE_LA_FENETRE_COURANTE; Eblock EDoI Eblock ATes Bblock Eblock ETes EDEFV(albumF,definition_du_polynome); /* Definition de l'album Argument definissant le polynome... */ Eblock end_nouveau_block RETU_Commande; Eblock ECommande