/*************************************************************************************************************************************/
/* */
/* M O D E L E B I D I M E N S I O N N E L D E C O T E S A V E C " E T C H I N G " : */
/* */
/* */
/* Author of '$xci/CoastL_2D.21$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20020519175926). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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_CPP_____POUR_LES_Istore_Iupdate_NE_PAS_VERIFIER_LES_MULTIPLES_UTILISATIONS_DE_L_IMAGE_____IFmageR
/* Introduit le 20060124113238. */
@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
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* P A R A M E T R E S : */
/* */
/*************************************************************************************************************************************/
#include xci/sequence.01.I"
#define EDITER_LES_PARAMETRES_DU_NOYAU \
FAUX \
/* Faut-il editer les parametres du noyau et sortir immediatement ('VRAI') ou exploiter */ \
/* le modele d'etching de cotes bidimensionnel ('FAUX') ? */
#define LISTER_LE_CONTENU_DU_NOYAU \
FAUX \
/* Faut-il lister le contenu du noyau ('VRAI') ou pas ('FAUX') ? */
#define GENERER_LES_CONDITIONS_INITIALES \
VRAI \
/* Doit-on generer les conditions initiales en tant que premiere image ('VRAI') ou bien */ \
/* le premier resultat du modele ('FAUX') ? */
#define IMPOSER_LES_EXTREMA_DES_CONDITIONS_INITIALES \
FAUX
#define MINIMUM_IMPOSE_POUR_LES_CONDITIONS_INITIALES \
NIVEAU_DE_LA_TERRE_DANS_Ietche__CoastLines_2D
#define MAXIMUM_IMPOSE_POUR_LES_CONDITIONS_INITIALES \
NIVEAU_DE_LA___MER_DANS_Ietche__CoastLines_2D
/* Faut-il imposer ('VRAI') ou pas ('FAUX') les extrema des conditions initiales, et si */
/* 'IL_FAUT(imposer_les_extrema_des_conditions_initiales)', valeurs implicites ? */
#define EDITER_LES_EXTREMA_DES_CONDITIONS_INITIALES \
FAUX \
/* Faut-il editer les extrema des des conditions initiales, ('VRAI') ou pas ('FAUX') ? */
#define NOMBRE_MINIMAL_D_ITERATIONS \
ZERO \
/* Nombre minimal d'iterations utilise dans le cas ou c'est le temps ou le pas de temps */ \
/* qui est utilise pour stopper les iterations. */
#define NOMBRE_MAXIMAL_D_ITERATIONS \
UN \
/* Nombre maximal d'iterations. */
#define SEUIL_DU_PAS_DE_TEMPS \
F_INFINI \
/* Seuil du pas de temps au-dela duquel on arrete d'iterer pour generer une image. Ce */ \
/* dispositif a ete introduit le 20020711162157 afin de tenter de visualiser le phenomene */ \
/* d'avalanche. Ainsi, on devrait regrouper dans une meme image tous les petits de temps */ \
/* et laisser plutot isoles les grands pas de temps. */
#define SEUIL_DU_CUMUL_DES_PAS_DE_TEMPS \
F_INFINI \
/* De meme pour le cumul des pas de temps reinitialise au debut de chaque nouvelle image. */
#define ARRETER_D_ITERER_S_IL_N_Y_A_PAS_EU_DE_POINTS_ETCHES \
FAUX \
/* Faut-il s'arreter ('VRAI') ou continuer d'iterer ('FAUX') lorsqu'a l'iteration courante */ \
/* aucun point n'a ete "etche" ? */
#define DECREMENTER_LE_MARQUEUR_DES_POINTS_ETCHES_APRES_L_ETCHING \
VRAI \
/* Pour decrementer 'Ietche__CoastLines_2D_____marqueur_des_points_etches' apres chaque */ \
/* iteration ('VRAI') ou avant ('FAUX'). En decrementant avant, on garantit ainsi que le */ \
/* niveau de marquage de la mer est strictement superieur a celui de la cote. Ceci a ete */ \
/* introduit le 20040909133217... */
#define DECREMENT_DU_MARQUEUR_DES_POINTS_ETCHES \
FZERO \
/* Decrement de 'Ietche__CoastLines_2D_____marqueur_des_points_etches' avant ou apres */ \
/* chaque iteration... */
#define NOMBRE_MAXIMAL_D_ETCHERS \
FLOT(CENT_MILLE) \
/* Nombre maximal (et initial) d'"etcher"s. */
#define CONSTANTE_DE_BOLTZMANN_ETCHE______COTES_2D \
FU \
/* Afin de simplifier l'utilisation et eviter les echelles physiques. */
#define TEMPERATURE_GLOBALE \
FU \
/* Temperature globale. */
#define GRAINE_INITIALE_DE_SELECTION_DES_POINTS_A_ETCHER \
TROIS \
/* Pour choisir les points a "etcher". */
#define EDITER_QUELQUES_NOMBRES_UTILES \
FAUX \
/* Faut-il editer le nombre d'"etchers" avant "etching" et le nombre de points "etches" */ \
/* apres "etching" ('VRAI') ou pas ('FAUX') ? */
#define NOMBRE_DE_CHIFFRES_A_EDITER \
NOMBRE_DE_CHIFFRES \
/* Nombre de chiffres a editer pour les numeros d'image. */
#define CONVERTIR_EN_IMAGES_STANDARDS \
VRAI \
/* Faut-il convertir les images Resultats en images standards ('VRAI') ou les laisser */ \
/* sous forme d'images flottantes ('FAUX') ? */
#define EDITER_L_EVENTUELLE_GRAINE_INITIALE_SUIVANTE \
FAUX \
/* Doit-on editer ('VRAI') ou pas ('FAUX') la graine a re-injecter dans le modele pour */ \
/* poursuivre de facon coherente la generation ? */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* M A C R O S U T I L E S : */
/* */
/*************************************************************************************************************************************/
#include xci/Verhulst.01.I"
/* Introduit le 20040916112632. */
#define INCREMENT_DES_GRAINES \
CENT_MILLE \
/* Pour incrementer les graines du generateur aleatoire apres chaque iteration, en */ \
/* particulier pour 'v $xiii/di_image$FON valeur_aleatoire_de_selection_du_point_a_etcher'. */ \
/* ATTENTION : avant le 20020603140431, 'INCREMENT_DES_GRAINES' valait 'I'. Mais a cause */ \
/* de la modification 'v $xiii/di_image$FON ADD2.graine_de_selection_des_points_a_etcher', */ \
/* il est essentiel que 'graine_de_selection_des_points_a_etcher' progresse (de +1 via */ \
/* 'index_valeur_aleatoire' dans '$xiii/di_image$FON') differemment que progresse */ \
/* 'graine_initiale_de_selection_des_points_a_etcher' (de 'INCREMENT_DES_GRAINES' ici) car */ \
/* sinon on retrouverait a chaque iteration quasiment toutes le valeurs aleatoires trouvees */ \
/* a l'iteration precedente decalees d'un "cran"... */
#define TITRE_graine_initiale_de_selection_des_points_a_etcher \
"graine=" \
/* Titre utile. */
#define NOMBRE_DE_CHIFFRES_NECESSAIRES(dimension) \
ADD2(NOMBRE_DE_CHIFFRES_DECIMAUX(dimension),UN) \
/* Choix du nombre de chiffres necessaires pour editer une dimension. On notera qu'il y */ \
/* a deux 'TRPU(...)', le deuxieme etant destine au signe qui compte pour un caractere. */ \
/* */ \
/* La procedure 'NOMBRE_DE_CHIFFRES_DECIMAUX(...)' fut introduite le 20051210183717. Les */ \
/* deux 'TRPU(...)' furent donc remplaces par 'un 'ADD2(...,UN)'... */
#define DECREMENTATION_DU_MARQUEUR_DES_POINTS_ETCHES(decrementer) \
Bblock \
Test(IL_FAUT(decrementer)) \
Bblock \
Test(IL_FAUT(Ietche__CoastLines_2D_____marquer_les_points_etches_arbitrairement)) \
Bblock \
DECR(Ietche__CoastLines_2D_____marqueur_des_points_etches \
,decrement_du_marqueur_des_points_etches \
); \
/* Gestion du marqueur des points etches lorsqu'il est utile... */ \
Eblock \
ATes \
Bblock \
Eblock \
ETes \
Eblock \
ATes \
Bblock \
Eblock \
ETes \
Eblock \
/* Decrementation du marqueur des points etches. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* M O D E L E B I D I M E N S I O N N E L D E C O T E S A V E C " E T C H I N G " : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_imageA),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_imageN),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_imageT),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_imageF),NOM_PIPE));
DEFV(Int,INIT(premiere_image,PREMIERE_IMAGE));
/* Numero de la premiere image, */
DEFV(Int,INIT(derniere_image,DERNIERE_IMAGE));
/* Numero de la derniere image. */
DEFV(Int,INIT(pas_des_images,PAS_DES_IMAGES));
/* Pas de passage d'un numero d'image a une autre. */
DEFV(Int,INIT(nombre_de_chiffres,NOMBRE_DE_CHIFFRES));
/* Nombre de chiffres codant le numero des images de la serie... */
DEFV(Int,INIT(numero_d_image,UNDEF));
/* Numero de l'image courante. */
DEFV(Logical,INIT(editer_les_parametres_du_noyau,EDITER_LES_PARAMETRES_DU_NOYAU));
/* Faut-il editer les parametres du noyau et sortir immediatement ('VRAI') ou exploiter */
/* le modele d'etching de cotes bidimensionnel ('FAUX') ? */
DEFV(Logical,INIT(lister_le_contenu_du_noyau,LISTER_LE_CONTENU_DU_NOYAU));
/* Faut-il lister le contenu du noyau ('VRAI') ou pas ('FAUX') ? */
DEFV(Logical,INIT(generer_les_conditions_initiales,GENERER_LES_CONDITIONS_INITIALES));
/* Doit-on generer les conditions initiales en tant que premiere image ('VRAI') ou bien */
/* le premier resultat du modele ('FAUX') ? */
DEFV(Logical,INIT(imposer_les_extrema_des_conditions_initiales,IMPOSER_LES_EXTREMA_DES_CONDITIONS_INITIALES));
DEFV(genere_Float,INIT(minimum_impose_pour_les_conditions_initiales,MINIMUM_IMPOSE_POUR_LES_CONDITIONS_INITIALES));
DEFV(genere_Float,INIT(maximum_impose_pour_les_conditions_initiales,MAXIMUM_IMPOSE_POUR_LES_CONDITIONS_INITIALES));
/* Faut-il imposer ('VRAI') ou pas ('FAUX') les extrema des conditions initiales, et si */
/* 'IL_FAUT(imposer_les_extrema_des_conditions_initiales)', valeurs implicites ? */
DEFV(Logical,INIT(editer_les_extrema_des_conditions_initiales,EDITER_LES_EXTREMA_DES_CONDITIONS_INITIALES));
/* Faut-il editer les extrema des des conditions initiales, ('VRAI') ou pas ('FAUX') ? */
DEFV(Positive,INIT(nombre_minimal_d_iterations,NOMBRE_MINIMAL_D_ITERATIONS));
/* Nombre minimal d'iterations utilise dans le cas ou c'est le temps ou le pas de temps */
/* qui est utilise pour stopper les iterations. */
DEFV(Positive,INIT(nombre_maximal_d_iterations,NOMBRE_MAXIMAL_D_ITERATIONS));
/* Nombre maximal d'iterations. */
DEFV(Float,INIT(seuil_du_pas_de_temps,SEUIL_DU_PAS_DE_TEMPS));
/* Seuil du pas de temps au-dela duquel on arrete d'iterer pour generer une image. Ce */
/* dispositif a ete introduit le 20020711162157 afin de tenter de visualiser le phenomene */
/* d'avalanche. Ainsi, on devrait regrouper dans une meme image tous les petits de temps */
/* et laisser plutot isoles les grands pas de temps. */
DEFV(Float,INIT(seuil_du_cumul_des_pas_de_temps,SEUIL_DU_CUMUL_DES_PAS_DE_TEMPS));
/* De meme pour le cumul des pas de temps reinitialise au debut de chaque nouvelle image. */
DEFV(Logical,INIT(arreter_d_iterer_s_il_n_y_a_pas_eu_de_points_etches,ARRETER_D_ITERER_S_IL_N_Y_A_PAS_EU_DE_POINTS_ETCHES));
/* Faut-il s'arreter ('VRAI') ou continuer d'iterer ('FAUX') lorsqu'a l'iteration courante */
/* aucun point n'a ete "etche" ? */
DEFV(Logical,INIT(decrementer_le_marqueur_des_points_etches_apres_l_etching
,DECREMENTER_LE_MARQUEUR_DES_POINTS_ETCHES_APRES_L_ETCHING
)
);
/* Pour decrementer 'Ietche__CoastLines_2D_____marqueur_des_points_etches' apres chaque */
/* iteration ('VRAI') ou avant ('FAUX'). En decrementant avant, on garantit ainsi que le */
/* niveau de marquage de la mer est strictement superieur a celui de la cote. Ceci a ete */
/* introduit le 20040909133217... */
DEFV(Float,INIT(decrement_du_marqueur_des_points_etches,DECREMENT_DU_MARQUEUR_DES_POINTS_ETCHES));
/* Decrement de 'Ietche__CoastLines_2D_____marqueur_des_points_etches' avant ou apres */
/* chaque iteration... */
DEFV(Float,INIT(nombre_maximal_d_etchers,NOMBRE_MAXIMAL_D_ETCHERS));
/* Ombre maximal (et initial) d'"etcher"s. */
DEFV(Int,INIT(graine_initiale_de_selection_des_points_a_etcher,GRAINE_INITIALE_DE_SELECTION_DES_POINTS_A_ETCHER));
/* Pour choisir les points a "etcher". */
DEFV(Float,INIT(temperature_globale,TEMPERATURE_GLOBALE));
/* Temperature globale. */
DEFV(Logical,INIT(editer_quelques_nombres_utiles,EDITER_QUELQUES_NOMBRES_UTILES));
/* Faut-il editer le nombre d'"etchers" avant "etching" et le nombre de points "etches" */
/* apres "etching" ('VRAI') ou pas ('FAUX') ? */
DEFV(Int,INIT(nombre_de_chiffres_a_editer,NOMBRE_DE_CHIFFRES_A_EDITER));
/* Nombre de chiffres a editer pour les numeros. */
DEFV(Logical,INIT(convertir_en_images_standards,CONVERTIR_EN_IMAGES_STANDARDS));
/* Faut-il convertir les images Resultats en images standards ('VRAI') ou les laisser */
/* sous forme d'images flottantes ('FAUX') ? */
DEFV(Logical,INIT(editer_l_eventuelle_graine_initiale_suivante,EDITER_L_EVENTUELLE_GRAINE_INITIALE_SUIVANTE));
/* Doit-on editer ('VRAI') ou pas ('FAUX') la graine a re-injecter dans le modele pour */
/* poursuivre de facon coherente la generation ? */
/*..............................................................................................................................*/
EGAL(Ietche__CoastLines_2D_____constante_de_Boltzmann,CONSTANTE_DE_BOLTZMANN_ETCHE______COTES_2D);
GET_ARGUMENTSi(nombre_d_arguments
,BLOC(GET_ARGUMENT_C("imageA=""A=",nom_imageA);
GET_ARGUMENT_C("imageN=""N=",nom_imageN);
GET_ARGUMENT_C("imageT=""T=",nom_imageT);
GET_ARGUMENT_C("imageR=""R=",nom_imageR);
GET_ARGUMENT_C("imageF=""F=",nom_imageF);
GET_ARGUMENT_I("premiere=",premiere_image);
GET_ARGUMENT_I("derniere=",derniere_image);
GET_ARGUMENT_I("pas=",pas_des_images);
GET_ARGUMENT_I("chiffres=",nombre_de_chiffres);
GET_ARGUMENT_L("initiales=",generer_les_conditions_initiales);
GET_ARGUMENT_L("periodiser_X=""periodX=",Ietche__CoastLines_2D_____periodiser_X);
GET_ARGUMENT_L("periodiser_Y=""periodY=",Ietche__CoastLines_2D_____periodiser_Y);
GET_ARGUMENT_L("symetriser_X=""symetrX=",Ietche__CoastLines_2D_____symetriser_X);
GET_ARGUMENT_L("symetriser_Y=""symetrY=",Ietche__CoastLines_2D_____symetriser_Y);
GET_ARGUMENT_L("prolonger_X=""prolX=",Ietche__CoastLines_2D_____prolonger_X);
GET_ARGUMENT_L("prolonger_Y=""prolY=",Ietche__CoastLines_2D_____prolonger_Y);
GET_ARGUMENT_F("niveau_hors_image=""niveau=",Ietche__CoastLines_2D_____niveau_hors_image);
GET_ARGUMENT_L("imposer_extrema=""extrema=",imposer_les_extrema_des_conditions_initiales);
GET_ARGUMENT_F("minimum=",minimum_impose_pour_les_conditions_initiales);
GET_ARGUMENT_F("maximum=",maximum_impose_pour_les_conditions_initiales);
GET_ARGUMENT_L("editer_extrema=",editer_les_extrema_des_conditions_initiales);
GET_ARGUMENT_L("parallele=",Ietche__CoastLines_2D_____explorer_parallelement_la_cote);
GET_ARGUMENT_I("nombre_maximal=""nmpe=",Ietche__CoastLines_2D_____nombre_maximal_de_points_a_etcher);
GET_ARGUMENT_I("minimum_iterations=",nombre_minimal_d_iterations);
GET_ARGUMENT_I("ni=""maximum_iterations=""iterations=",nombre_maximal_d_iterations);
/* Le 20050623141115, "n=" a ete remplace par "ni=" (double definition...). */
GET_ARGUMENT_F("seuil_dt=",seuil_du_pas_de_temps);
GET_ARGUMENT_F("seuil_cdt=",seuil_du_cumul_des_pas_de_temps);
GET_ARGUMENT_L("arreter=""arret=",arreter_d_iterer_s_il_n_y_a_pas_eu_de_points_etches);
GET_ARGUMENT_L("noyau=",Ietche__CoastLines_2D_____utiliser_un_noyau);
GET_ARGUMENT_I("demi_dimension=""dd=",Ietche__CoastLines_2D_____demi_dimension_effective_du_noyau);
GET_ARGUMENT_L("editer_noyau=",editer_les_parametres_du_noyau);
GET_ARGUMENT_L("lister_noyau=",lister_le_contenu_du_noyau);
GET_ARGUMENT_F("sX__Y=""sXY=",Ietche__CoastLines_2D_____ponderation_du_point_sX__Y);
GET_ARGUMENT_F("_X_sY=""XsY=",Ietche__CoastLines_2D_____ponderation_du_point__X_sY);
GET_ARGUMENT_F("pX__Y=""pXY=",Ietche__CoastLines_2D_____ponderation_du_point_pX__Y);
GET_ARGUMENT_F("_X_pY=""XpY=",Ietche__CoastLines_2D_____ponderation_du_point__X_pY);
GET_ARGUMENT_I("dx=""dX=",Ietche__CoastLines_2D_____delta_X);
GET_ARGUMENT_I("dy=""dY=",Ietche__CoastLines_2D_____delta_Y);
GET_ARGUMENT_L("marquer_taux=",Ietche__CoastLines_2D_____marquer_les_points_etches_avec_le_taux);
GET_ARGUMENT_L("marquer_temps=",Ietche__CoastLines_2D_____marquer_les_points_etches_avec_le_temps_simule);
GET_ARGUMENT_L("marquer_pas_de_temps="
,Ietche__CoastLines_2D_____marquer_les_points_etches_avec_le_pas_de_temps_simule
);
GET_ARGUMENT_L("marquer_arbitraire=",Ietche__CoastLines_2D_____marquer_les_points_etches_arbitrairement);
GET_ARGUMENT_F("points_etches=""mpe=",Ietche__CoastLines_2D_____marqueur_des_points_etches);
GET_ARGUMENT_L("decrementer_apres=""dap=",decrementer_le_marqueur_des_points_etches_apres_l_etching);
GET_ARGUMENT_N("decrementer_avant=""dav=",decrementer_le_marqueur_des_points_etches_apres_l_etching);
GET_ARGUMENT_F("dmpe=""decrement_points_etches=",decrement_du_marqueur_des_points_etches);
GET_ARGUMENT_I(TITRE_graine_initiale_de_selection_des_points_a_etcher
,graine_initiale_de_selection_des_points_a_etcher
);
GET_ARGUMENT_L("editer_graine=",editer_l_eventuelle_graine_initiale_suivante);
GET_ARGUMENT_F("etchers=""Metchers=",nombre_maximal_d_etchers);
GET_ARGUMENT_F("Cetchers=",Ietche__CoastLines_2D_____nombre_courant_d_etchers);
GET_ARGUMENT_F("concentration=""c=",Ietche__CoastLines_2D_____concentration_initiale_d_etchers);
GET_ARGUMENT_F("frequence=""f=",Ietche__CoastLines_2D_____frequence_de_la_reaction);
/* ATTENTION : ces deux parametres sont inutiles ('v $xiii/di_image$FON 20020604114847'). */
GET_ARGUMENT_F("energie_minimale=""em=",Ietche__CoastLines_2D_____energie_minimale);
GET_ARGUMENT_F("energie_maximale=""eM=",Ietche__CoastLines_2D_____energie_maximale);
GET_ARGUMENT_F("temperature=""t=",temperature_globale);
GET_ARGUMENT_L("editer_nombres=",editer_quelques_nombres_utiles);
GET_ARGUMENT_I("nombre=""nchiffres=""nc=",nombre_de_chiffres_a_editer);
/* Le 20050623141115, "chiffres=" a ete remplace par "nchiffres=" (double definition...) et */
/* "n=" a ete remplace par "nc=" (double definition...). */
GET_ARGUMENT_L("editer_coordonnees="
,Ietche__CoastLines_2D_____editer_les_coordonnees_des_points_parmi_les_nombres_utiles
);
GET_ARGUMENT_L("convertir=",convertir_en_images_standards);
)
);
Test(IL_FAUT(editer_les_parametres_du_noyau))
Bblock
CAL3(Prme4("Xmin=%d Xmax=%d Ymin=%d Ymax=%d\n"
,PINTE(Ietche__CoastLines_2D_____XYmin),PINTE(Ietche__CoastLines_2D_____XYmax)
,PINTE(Ietche__CoastLines_2D_____XYmin),PINTE(Ietche__CoastLines_2D_____XYmax)
)
);
CAL3(Prme4("XminR=%d XmaxR=%d YminR=%d YmaxR=%d\n"
,PINTE(Ietche__CoastLines_2D_____XYmin),PINTE(Ietche__CoastLines_2D_____XYmax)
,PINTE(Ietche__CoastLines_2D_____XYmin),PINTE(Ietche__CoastLines_2D_____XYmax)
)
);
CALS(Fsauts_de_lignes(UN));
CAL3(Prme4("XYmaxNe %d %d %d %d\n"
,PINTE(Ietche__CoastLines_2D_____XYmin),PINTE(Ietche__CoastLines_2D_____XYmax)
,PINTE(Ietche__CoastLines_2D_____XYmin),PINTE(Ietche__CoastLines_2D_____XYmax)
)
);
/* ATTENTION : on ne devra pas etre surpris par les valeurs ainsi obtenues. Par exemple, on */
/* obtiendra : */
/* */
/* XYmax 128 128 */
/* */
/* avec la valeur 64 pour 'DEMI_DIMENSION_MAXIMALE_DU_NOYAU_DANS_Ietche__CoastLines_2D', en */
/* se rappelant que l'on passe a 'DimNo_Ietche__CoastLines_2D' par un 'DOUP(...)'. Les */
/* extrema ainsi obtenus (128) different d'une unite de ceux fournis par 'Suq' (127) et */
/* les parites des deux dimensions respectives (129 et 128) sont elles-aussi differentes... */
/* */
/* ATTENTION : passage de 'XYmax' a 'XYmaxNe' le 20020909112942. */
/* */
/* Pour l'usage de 'PINTE(...)' voir 'v $xil/defi_K2$vv$DEF 20230311135805' qui a ete */
/* introduit pour ces impressions le 20230311144025... */
Eblock
ATes
Bblock
DEFV(Int,INIT(graine_courante_de_selection_des_points_a_etcher
,graine_initiale_de_selection_des_points_a_etcher
)
);
/* Pour choisir les points a "etcher". */
CALi(IFinitialisation(IFmageR,FZERO));
/* Initialisation de l'image Resultat flottante (dans tous les cas car elle est calculee */
/* meme si 'IL_FAUT(convertir_en_images_standards)'). */
Test(IL_FAUT(convertir_en_images_standards))
Bblock
CALi(Inoir(ImageR));
/* Initialisation de l'image Resultat standard si elle est utile... */
Eblock
ATes
Bblock
Eblock
ETes
Test(IFNE_chaine(nom_imageN,NOM_PIPE))
Bblock
PUSH_DIMENSIONS_2D;
SET_DIMENSIONS_2D_SANS_VALIDATION(Ietche__CoastLines_2D_____XYmin,Ietche__CoastLines_2D_____XYmax
,Ietche__CoastLines_2D_____XYmin,Ietche__CoastLines_2D_____XYmax
);
/* Le noyau va etre recupere comme une "petite" image dont les dimensions sont definies */
/* par {Ietche__CoastLines_2D_____XYmin,Ietche__CoastLines_2D_____XYmax} et ce */
/* independemment de 'Ietche__CoastLines_2D_____demi_dimension_effective_du_noyau'... */
Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(PIetche__CoastLines_2D_____noyau,nom_imageN))))
/* Chargement de l'image du noyau... */
/* */
/* ATTENTION : on notera que l'image 'nom_imageN' doit avoir comme taille la taille */
/* maximale du noyau correspondant a 'v $xiii/di_image$DEF DimNo_Ietche__CoastLines_2D', */
/* et non pas celle qui correspond a la taille effective de ce meme noyau */
/* 'v $xiii/di_image$FON Ietche__CoastLines_2D_____demi_dimension_effective_du_noyau'. */
/* */
/* ATTENTION : on notera de plus que l'on passe par 'PIetche__CoastLines_2D_____noyau' car, */
/* en effet, 'Ietche__CoastLines_2D_____noyau' n'est pas un vrai 'albumF' puisqu'il n'est */
/* cree par une allocation dynamique et variable... */
Bblock
EGAL(Ietche__CoastLines_2D_____initialiser_le_noyau,FAUX);
/* Evidemment, l'initialisation est faite... */
Test(IL_FAUT(lister_le_contenu_du_noyau))
Bblock
begin_imageQ(DoIn,Ietche__CoastLines_2D_____XYmin_effectif,Ietche__CoastLines_2D_____XYmax_effectif,PasY
,DoIn,Ietche__CoastLines_2D_____XYmin_effectif,Ietche__CoastLines_2D_____XYmax_effectif,PasX
)
Bblock
CAL3(Prme5("noyau(%+0*d,%+0*d)=%+.^^^\n"
,NOMBRE_DE_CHIFFRES_NECESSAIRES(DEMI_DIMENSION_MAXIMALE_DU_NOYAU_DANS_Ietche__CoastLines_2D),X
,NOMBRE_DE_CHIFFRES_NECESSAIRES(DEMI_DIMENSION_MAXIMALE_DU_NOYAU_DANS_Ietche__CoastLines_2D),Y
,ACCES_NOYAU_DANS_Ietche__CoastLines_2D(X,Y)
)
);
/* Le 20060105153658, le format "16g" est passe a "^^g" pour plus de souplesse... */
/* */
/* Le 20091123123102, le format "^^g" est passe a "^^^" pour plus de souplesse... */
Eblock
end_imageQ(EDoI,EDoI)
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
INITIALISATION_EVENTUELLE_DU_NOYAU_DANS_Ietche__CoastLines_2D;
/* Initialisation du noyau en cas d'erreur... */
Eblock
ETes
PULL_DIMENSIONS_2D;
/* Restauration des dimensions de travail... */
Eblock
ATes
Bblock
Eblock
ETes
LOAD_D_UNE_IMAGE_FLOTTANTE_DE_PARAMETRAGE(IFmageA1,nom_imageT,temperature_globale,"temperature globale");
Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(IFmageA,nom_imageA))))
Bblock
DEFV(genere_Float,INIT(niveau_minimum_de_IFmageA,minimum_impose_pour_les_conditions_initiales));
DEFV(genere_Float,INIT(niveau_maximum_de_IFmageA,maximum_impose_pour_les_conditions_initiales));
/* Extrema du champ initial ('IFmageA'). */
Test(IL_NE_FAUT_PAS(imposer_les_extrema_des_conditions_initiales))
Bblock
CALS(IFnivo_extrema(IFmageA
,ADRESSE(niveau_minimum_de_IFmageA)
,ADRESSE(niveau_maximum_de_IFmageA)
)
);
/* Recherche des extrema de 'IFmageA'. */
Eblock
ATes
Bblock
Eblock
ETes
Test(IL_FAUT(editer_les_extrema_des_conditions_initiales))
Bblock
CAL3(Prme1("minimum des conditions initiales=%+.^^^\n",niveau_minimum_de_IFmageA));
CAL3(Prme1("maximum des conditions initiales=%+.^^^\n",niveau_maximum_de_IFmageA));
/* Edition des extrema de 'IFmageA'. */
/* */
/* Le 20060105153658, le format "16g" est passe a "^^g" pour plus de souplesse... */
/* */
/* Le 20091123123102, le format "^^g" est passe a "^^^" pour plus de souplesse... */
Eblock
ATes
Bblock
Eblock
ETes
EGAL(Ietche__CoastLines_2D_____niveau_de_la___mer,niveau_maximum_de_IFmageA);
EGAL(Ietche__CoastLines_2D_____niveau_de_la_terre,niveau_minimum_de_IFmageA);
/* Mise a jour du modele... */
Test(IL_FAUT(Ietche__CoastLines_2D_____marquer_les_points_etches_arbitrairement))
Bblock
Test(IFOU(SONT_DE_MEME_SIGNE(NEUT(Ietche__CoastLines_2D_____marqueur_des_points_etches)
,Ietche__CoastLines_2D_____niveau_de_la_terre
)
,SONT_DE_MEME_SIGNE(SOUS(Ietche__CoastLines_2D_____marqueur_des_points_etches
,MUL2(PRED(MUL2(NBRE(premiere_image,derniere_image)
,nombre_maximal_d_iterations
)
)
,decrement_du_marqueur_des_points_etches
)
)
,Ietche__CoastLines_2D_____niveau_de_la_terre
)
)
)
Bblock
PRINT_ATTENTION("il va y avoir confusion entre le marqueur des points etches et le niveau de la terre");
CAL3(Prme1("niveau de la terre................=%+.^^^\n",Ietche__CoastLines_2D_____niveau_de_la_terre));
/* Le 20060105153658, le format "16g" est passe a "^^g" pour plus de souplesse... */
/* */
/* Le 20091123123102, le format "^^g" est passe a "^^^" pour plus de souplesse... */
CAL3(Prme1("nombre total d'iterations.........=%d\n"
,MUL2(NBRE(premiere_image,derniere_image),nombre_maximal_d_iterations)
)
);
CAL3(Prme1("marqueur initial des points etches=%+.^^^\n"
,NEUT(Ietche__CoastLines_2D_____marqueur_des_points_etches)
)
);
/* Le 20060105153658, le format "16g" est passe a "^^g" pour plus de souplesse... */
/* */
/* Le 20091123123102, le format "^^g" est passe a "^^^" pour plus de souplesse... */
CAL3(Prme1("marqueur final des points etches..=%+.^^^\n"
,SOUS(Ietche__CoastLines_2D_____marqueur_des_points_etches
,MUL2(PRED(MUL2(NBRE(premiere_image,derniere_image),nombre_maximal_d_iterations))
,decrement_du_marqueur_des_points_etches
)
)
)
);
/* Le 20060105153658, le format "16g" est passe a "^^g" pour plus de souplesse... */
/* */
/* Le 20091123123102, le format "^^g" est passe a "^^^" pour plus de souplesse... */
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
DoIn(numero_d_image,premiere_image,derniere_image,pas_des_images)
Bblock
DEFV(Float,INIT(cumul_des_pas_de_temps,FZERO));
/* Afin de pouvoir eventuellement s'arreter lorsque le cumul des pas de temps depasse */
/* son seuil... */
DEFV(CHAR,INIC(POINTERc(nom_imagesR),NOM_PIPE));
Test(IFET(IL_FAUT(generer_les_conditions_initiales),IFEQ(numero_d_image,premiere_image)))
Bblock
CALS(IFmove(IFmageR,IFmageA));
Eblock
ATes
Bblock
DEFV(Logical,INIT(iterer_l_etching,VRAI));
DEFV(Int,INIT(compteur_des_iterations_d_etching,PREMIERE_ITERATION_D_UN_Repe));
/* Afin de simuler un 'Repe(...)' a nombre variable d'iterations... */
Tant(IL_FAUT(iterer_l_etching))
Bblock
Test(IZGT(Ietche__CoastLines_2D_____nombre_courant_d_etchers))
/* 'Ietche__CoastLines_2D_____nombre_de_points_etches_lors_de_l_iteration_courante' n'est */
/* pas teste par rapport a sa nullite afin d'arreter dans ce cas les iterations car, en */
/* effet, il s'agit d'un processus aleatoire et s'il reste des "etchers", rien ne prouve */
/* qu'a l'iteration suivante il n'y aura pas de points "etches" meme s'il n'y en a pas eu a */
/* l'iteration courante. Malgre la remarque precedente, le 20020524090327, l'indicateur */
/* 'arreter_d_iterer_s_il_n_y_a_pas_eu_de_points_etches' a ete introduit... */
Bblock
Test(IL_FAUT(editer_quelques_nombres_utiles))
Bblock
CAL3(Prme3("iteration %d de l'image %0*d :"
,compteur_des_iterations_d_etching
,nombre_de_chiffres_a_editer
,numero_d_image
)
);
Eblock
ATes
Bblock
Eblock
ETes
#define decrementer_apres_l_etching \
decrementer_le_marqueur_des_points_etches_apres_l_etching \
/* Pour raccourcir certaines lignes qui suivent... */
DECREMENTATION_DU_MARQUEUR_DES_POINTS_ETCHES(NOTL(decrementer_apres_l_etching));
CALS(Ietche__CoastLines_2D_a_temperatures_locales(IFmageR
,IFmageA
,nombre_maximal_d_etchers
,graine_courante_de_selection_des_points_a_etcher
,IFmageA1
,editer_quelques_nombres_utiles
)
);
/* Et calcul de la generation suivante, avec : */
/* */
/* temperature_locale = IFmageA1 */
/* */
INCR(cumul_des_pas_de_temps,Ietche__CoastLines_2D_____pas_de_temps_simule);
/* Cumul des pas de temps... */
INCR(graine_courante_de_selection_des_points_a_etcher,INCREMENT_DES_GRAINES);
/* Progression des graines du generateur aleatoire puisqu'etant fonction de {X,Y}, il est */
/* fondamental qu'a l'iteration suivante, il ne redonne pas la meme serie, principalement */
/* pour 'v $xiii/di_image$FON valeur_aleatoire_de_selection_du_point_a_etcher'. */
CALS(IFmove(IFmageA,IFmageR));
/* Preparation de l'iteration suivante... */
INCR(compteur_des_iterations_d_etching,I);
/* Simulation du 'Repe(...)' d'"etching"... */
DECREMENTATION_DU_MARQUEUR_DES_POINTS_ETCHES(NEUL(decrementer_apres_l_etching));
#undef decrementer_apres_l_etching
Test(IFET(IL_FAUT(arreter_d_iterer_s_il_n_y_a_pas_eu_de_points_etches)
,IZEQ(Ietche__CoastLines_2D_____nombre_de_points_etches_lors_de_l_iteration_courante)
)
)
Bblock
EGAL(iterer_l_etching,FAUX);
EGAL(derniere_image,numero_d_image);
/* Il n'y a pas eu de nouveaux points "etches" ce qui provoque dans ce contexte un arret */
/* arret des iterations et de la generation des images... */
Eblock
ATes
Bblock
Test(I3OU(IFGT(compteur_des_iterations_d_etching,nombre_maximal_d_iterations)
,IFET(IFGT(Ietche__CoastLines_2D_____pas_de_temps_simule,seuil_du_pas_de_temps)
,IFGT(compteur_des_iterations_d_etching,nombre_minimal_d_iterations)
)
,IFET(IFGT(cumul_des_pas_de_temps,seuil_du_cumul_des_pas_de_temps)
,IFGT(compteur_des_iterations_d_etching,nombre_minimal_d_iterations)
)
)
)
/* Dans le cas ou c'est le temps ou le pas de temps qui sont utilises pour arreter le */
/* processus, j'ai ajoute le 20020712092746 la possibilite d'exiger d'avoir effectue au */
/* moins 'nombre_minimal_d_iterations' iterations... */
Bblock
EGAL(iterer_l_etching,FAUX);
/* On a fait le nombre d'iterations demandees... */
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ETes
Eblock
ATes
Bblock
EGAL(iterer_l_etching,FAUX);
/* En fait, il est devenu impossible d'iterer... */
Eblock
ETes
Eblock
ETan
Eblock
ETes
EGAL(nom_imagesR
,chain_Aconcaten2_sauf_nom_pipe(nom_imageR
,chain_numero(numero_d_image,nombre_de_chiffres)
)
);
Test(IL_FAUT(convertir_en_images_standards))
Bblock
CALS(Ifloat_std_avec_le_vrai_zero(ImageR,IFmageR,niveau_minimum_de_IFmageA,niveau_maximum_de_IFmageA));
/* Renormalisation avec forcage des extrema en utilisant ceux de l'image initiale 'IFmageA'. */
/* ATTENTION, jusqu'au 20020108124607, c'etait 'Ifloat_std(...)' qui etait utilise, ce qui */
/* n'etait pas tres logique puisqu'il est important de conserver le zero... */
CALi(Iupdate_image(nom_imagesR,ImageR));
Eblock
ATes
Bblock
CALi(IupdateF_image(nom_imagesR,IFmageR));
Eblock
ETes
CALZ_FreCC(nom_imagesR);
Eblock
EDoI
Test(IFET(IFNE_chaine(nom_imageF,NOM_PIPE),IL_FAUT(convertir_en_images_standards)))
Bblock
CALi(IupdateF_image(nom_imageF,IFmageR));
/* Dans le cas ou la sortie est du type "image standard", la possibilite est offerte de */
/* generer les conditions Finales en mode "non standard" afin qu'elles puissent servir */
/* eventuellement de conditions initiales dans un autre calcul... */
Eblock
ATes
Bblock
Eblock
ETes
Test(IL_FAUT(editer_l_eventuelle_graine_initiale_suivante))
Bblock
CAL3(Prme2("%s%d\n"
,TITRE_graine_initiale_de_selection_des_points_a_etcher
,graine_courante_de_selection_des_points_a_etcher
)
);
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
Eblock
ETes
RETU_Commande;
Eblock
ECommande