/*************************************************************************************************************************************/
/* */
/* 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 - C O M P L E X E : */
/* */
/* */
/* Author of '$xrc/polynome.41$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20040510144612). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 le 20100420142438 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 (introduit le 20091121120100) ? */ \
/* 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(-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
/* Definition des composantes 'J' et 'K' des quaternions. */
#define GENERATEUR_REEL \
PARE(0.0)
#define GENERATEUR_IMAGINAIRE \
PARE(1.0)
#define GENERATEUR_JMAGINAIRE \
PARE(0.0)
#define GENERATEUR_KMAGINAIRE \
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_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 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_Theta_ \
IFmageA2
#define definition_de_ponder1_Phi__ \
IFmageA3
#define definition_de_ponder1_Alpha \
IFmageA4
/* Definitions introduites le 20100210134840... */
#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)
/* 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 ' 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 - 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 le 20100420142438 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_Theta_),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_Phi__),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_definition_de_ponder1_Alpha),NOM_PIPE));
/* Parametrage du "faux produit" (introduit le 20100210134840...). */
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 20091121120100) ? */
/* 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(hyper_complexe,point_argument);
/* Nombre 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_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. */
/*..............................................................................................................................*/
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. */
HCinitialisation(point_argument,GENERATEUR_REEL,GENERATEUR_IMAGINAIRE,GENERATEUR_JMAGINAIRE,GENERATEUR_KMAGINAIRE);
/* Initialisation du nombre 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 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 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
);
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_HC_____premier_exposant_du_polynome
);
/* Argument introduit le 20100316173523... */
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
);
/* Argument introduit le 20100420142438... */
GET_ARGUMENT_C("imageAeRho=""AeRho=""AeR=",nom_definition_de_exposa1_Rho__);
GET_ARGUMENT_L("partie_entiere_ApTheta=""peApTheta=""iApT="
,Iiteration_polynome_variable_quelconque_dans_HC_____prendre_la_partie_entiere_de_ponder1_Theta_
);
GET_ARGUMENT_C("imageApTheta=""ApTheta=""ApT=",nom_definition_de_ponder1_Theta_);
GET_ARGUMENT_L("partie_entiere_ApPhi=""peApPhi=""iApP="
,Iiteration_polynome_variable_quelconque_dans_HC_____prendre_la_partie_entiere_de_ponder1_Phi__
);
GET_ARGUMENT_C("imageApPhi=""ApPhi=""ApP=",nom_definition_de_ponder1_Phi__);
GET_ARGUMENT_L("partie_entiere_ApAlpha=""peApAlpha=""iApA="
,Iiteration_polynome_variable_quelconque_dans_HC_____prendre_la_partie_entiere_de_ponder1_Alpha
);
GET_ARGUMENT_C("imageApAlpha=""ApAlpha=""ApA=",nom_definition_de_ponder1_Alpha);
/* Arguments introduits le 20100210134840 et les abbreviations {"AeR=","ApT=","ApP=","ApA="} */
/* le furent le 20100220160327 pour 'v $xiirc/$Fnota Debut_listG_POLY_61_62'. Les arguments */
/* relatifs a la prise de la partie entiere ont ete introduits le 20100221091228... */
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
);
/* Arguments introduits le 20091121120100... */
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("r=""gr=",HReelle(point_argument));
GET_ARGUMENT_F("i=""gi=",HImaginaire(point_argument));
GET_ARGUMENT_F("j=""gj=",HJmaginaire(point_argument));
GET_ARGUMENT_F("k=""gk=",HKmaginaire(point_argument));
GET_ARGUMENT_L("rotation_rijk=""Rrijk=""RA="
,Iiteration_polynome_variable_quelconque_dans_HC_____faire_la_rotation_de__nombre_argument_A
);
/* Arguments introduits le 20100303082630... */
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 20150305205741 {"PJulia=","PMandelrrot="} ont ete introduits... */
GET_ARGUMENT_F("seuil=",seuil_du_carre_du_module);
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 sous cette forme le 20070108143853... */
GET_ARGUMENT_L("F_ITERATION_DANS_HC_____compatibilite_20110825=""compatibilite_20110825="
,F_ITERATION_DANS_HC_____compatibilite_20110825
);
/* Introduit le 20110826123224. */
CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_HYPER_COMPLEXES;
/* Cette extension a ete introduite le 20091119095338... */
CONTROLE_DE_LA_TRANFORMATION_DE_LA_SUITE_DE_POINTS_DANS_HC;
/* Cette extension a ete introduite le 20100513181816... */
CONTROLE_DU_REPLIEMENT_DE_L_ESPACE_DES_NOMBRES_HYPER_COMPLEXES;
/* Cette extension a ete introduite le 20121104080618... */
)
);
Test(EST_VRAI(les_fenetres_de_depart_et_d_arrivee_sont_identiques))
/* Test introduit le 20091121120100... */
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... */
/* */
/* L'argument 'Inv' ('FAUX') a ete introduit le 20091125134859... */
Test(PAS_D_ERREUR(CODE_ERREUR))
Bblock
Test(IL_FAUT(FHCproduit_____utiliser_la_definition_theorique_du_produit))
/* Test introduit le 20100210134840. C'est a cause de lui que les 'Test(...)'s relatifs */
/* aux 'IloadF_image(...)' ci-apres ne peuvent pas englober (et eventuellement exclure */
/* en cas d'erreur -inexistence d'un fichier-) le calcul du polynome... */
Bblock
Test(I4OU(IFNE_chaine(nom_definition_de_exposa1_Rho__,NOM_PIPE)
,IFNE_chaine(nom_definition_de_ponder1_Theta_,NOM_PIPE)
,IFNE_chaine(nom_definition_de_ponder1_Phi__,NOM_PIPE)
,IFNE_chaine(nom_definition_de_ponder1_Alpha,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)
)
)
/* Test complete le 20100420140056... */
/* */
/* On notera le 20100420142438 qu'il est difficile, pour ne pas dire impossible, de faire */
/* ici un test qui regarderait si ce sont des valeurs numeriques ou des noms de fichier qui */
/* ont ete introduits, d'ou cette "facilite"... */
Bblock
PRINT_ATTENTION("le parametrage du produit non theorique n'a de sens que pour le calcul d'une image unique");
/* Message introduit le 20100210141717. En effet, pour plus d'une image, il conviendrait */
/* que ce parametrage (de meme que la definition du polynome en toute generalite) soit */
/* "multidimensionnel"... et varie donc d'une image a l'autre ce qui serait d'une lourdeur */
/* extreme. La solution est que pour calculer 'N' images dans le cas du produit non */
/* theorique, il suffira d'appeler 'N' fois ce programmes avec des definitions differentes */
/* si besoin est... */
Eblock
ATes
Bblock
Eblock
ETes
CALi(IFinitialisation(definition_de_exposa1_Rho__,FU));
CALi(IFinitialisation(definition_de_ponder1_Theta_,FU));
CALi(IFinitialisation(definition_de_ponder1_Phi__,FU));
CALi(IFinitialisation(definition_de_ponder1_Alpha,FU));
/* A priori, les valeurs les plus probables sont forcees... */
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_Theta_
,nom_definition_de_ponder1_Theta_
)
)
)
)
Bblock
Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_Phi__
,nom_definition_de_ponder1_Phi__
)
)
)
)
Bblock
Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(definition_de_ponder1_Alpha
,nom_definition_de_ponder1_Alpha
)
)
)
)
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
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 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(Iiteration_polynome_variable_quelconque_dans_HC(ImageR
,definition_du_polynome
,definition_de_exposa1_Rho__
,definition_de_ponder1_Theta_
,definition_de_ponder1_Phi__
,definition_de_ponder1_Alpha
,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
,nombre_d_iterations
,seuil_du_carre_du_module
,visualisation_arguments_de_sortie
,ponderation_argumentT
,ponderation_argumentP
,ponderation_argumentA
)
);
/* Calcul de l'ensemble de Julia dans l'espace 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