/*************************************************************************************************************************************/ /* */ /* C A L C U L G E N E R A L D E L ' E N S E M B L E D E M A N D E L B R O T */ /* A V E C M I S E E N M O N T A G N E S D E C H A Q U E I M A G E */ /* D A N S L E P L A N C O M P L E X E : */ /* */ /* */ /* Author of '$xrc/mandel.02$K' : */ /* */ /* Jean-Francois COLONNA (LACTAMME, 1988??????????). */ /* */ /*************************************************************************************************************************************/ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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 20070129111928... */ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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 #include image_image_MONTAGNES_EXT /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* P A R A M E T R E S : */ /* */ /*************************************************************************************************************************************/ #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 (vue d'ensemble). */ #define X_ARRIVEE_GAUCHE \ PARE(-0.234254) #define X_ARRIVEE_DROITE \ PARE(-0.234254) #define Y_ARRIVEE_BAS \ PARE(-0.828110) #define Y_ARRIVEE_HAUT \ PARE(-0.828110) /* Definition de la fenetre d'arrivee (dans les filaments). */ #define OMBRES_PORTEES \ VRAI #define LARGEUR_ZONE_PENOMBRE \ PARE(12.0) /* Pour */ #define DEPTH_CUEING \ VRAI #define MIN_DEPTH_CUEING \ FZERO /* generer */ #define AVION FAUX \ /* la surface */ #define ALIAS VRAI \ /* montagneuse. */ #define TRANS_OX 0.0 #define TRANS_OY 0.80000 /* Translation de la montagne. */ #define X_SOURCE 6.4453 #define Y_SOURCE 2.7343 \ /* Position de la source lumineuse. */ #define MODULATION_DU_CHAMP \ FU \ /* Modulation du champ de Mandelbrot. */ #define NOMBRE_D_ITERATIONS \ SOUS(BLANC,NOIR) \ /* Nombre d'iterations demandees. */ #define PRECISION \ PARE(0.00001) \ /* Pour calculer le rapport de reduction. */ #define VISUALISATION_ARGUMENTS_DE_SORTIE \ FAUX \ /* Faut-il visualiser l'argument de sortie ('VRAI') ou les equipotentielles ('FAUX'). */ #include xci/sequence.01.I" /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* M A C R O S U T I L E S : */ /* */ /*************************************************************************************************************************************/ #include xrc/ITERATION.11.I" /* Introduit le 20091122182507... */ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* C A L C U L G E N E R A L D E L ' E N S E M B L E D E M A N D E L B R O T */ /* A V E C M I S E E N M O N T A G N E S D E C H A Q U E I M A G E */ /* D A N S L E P L A N C O M P L E X E : */ /* */ /*************************************************************************************************************************************/ BCommande(nombre_d_arguments,arguments) /*-----------------------------------------------------------------------------------------------------------------------------------*/ Bblock DEFV(complexe,coin_bas_gauche_de_depart_puis_courant); /* Definition du point situe en bas a gauche de la fenetre de depart, puis de la fenetre */ /* courante au cours du zoom. */ DEFV(complexe,coin_haut_droite_de_depart_puis_courant); /* Definition du point situe en haut a droite de la fenetre de depart, puis de la fenetre */ /* courante au cours du zoom. */ DEFV(complexe,coin_bas_gauche_d_arrivee); /* Definition du point situe en bas a gauche de la fenetre d'arrivee, */ DEFV(complexe,coin_haut_droite_d_arrivee); /* Definition du point situe en haut a droite de la fenetre d'arrivee. */ DEFV(Positive,INIT(nombre_d_iterations,NOMBRE_D_ITERATIONS)); /* Nombre maximal d'iterations a effectuer. */ 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(Logical,INIT(visualisation_arguments_de_sortie,VISUALISATION_ARGUMENTS_DE_SORTIE)); /* Indicateur de choix entre la visualisation du nombre d'iterations ('FAUX'), */ /* et l'argument de sortie ('VRAI'). */ DEFV(Logical,INIT(vue_d_avion,AVION)); /* Pour savoir si l'on genere une vue de cote, ou bien une vue d'avion... */ DEFV(Logical,INIT(ombres_portees,OMBRES_PORTEES)); DEFV(Float,INIT(largeur_zone_penombre,LARGEUR_ZONE_PENOMBRE)); /* Parametres de l'ombrage. */ DEFV(Logical,INIT(depth_cueing,DEPTH_CUEING)); DEFV(Float,INIT(min_depth_cueing,MIN_DEPTH_CUEING)); /* Parametres du depth-cueing. */ DEFV(deltaF_2D,translation_montagne); /* Translation de la surface montagneuse. */ DEFV(Float,INIT(modulation_du_champ,MODULATION_DU_CHAMP)); /* Modulation du champ de Mandelbrot. */ DEFV(pointF_2D,source_lumineuse); /* Position de la source lumineuse. */ DEFV(Int,INIT(nombre_d_images,NOMBRE_D_IMAGES)); /* Nombre d'images a generer. */ DEFV(Int,INIT(numero_d_image,PREMIERE_IMAGE_ABSOLUE)); /* Numero de l'image courante (celle-ci n'est pas necessairement generee : voir */ /* le parametre 'PREMIERE_IMAGE'). */ DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE)); /* Nom de la sequence d'images a generer. */ DEFV(Int,INIT(nombre_de_chiffres,NOMBRE_DE_CHIFFRES)); /* Nombre de chiffres codant le numero des images de la sequence a generer. */ DEFV(CHAR,INIT(POINTERc(nom_image),NOM_UNDEF)); /* Nom courant des images. */ /*..............................................................................................................................*/ EGAL(Z_Buffer_____valeur_initiale,VALEUR_INITIALE_DU_Z_BUFFER_POUR_LES_MONTAGNES); /* Introduit le 20060414171417... */ EGAL(Imontagnes_precises_____visualiser_la_falaise_avant_de_la_montagne,VRAI); /* Introduit le 20060414173128... */ 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. */ INITIALISATION_ACCROISSEMENT_2D(translation_montagne,TRANS_OX,TRANS_OY); INITIALISATION_POINT_2D(source_lumineuse,X_SOURCE,Y_SOURCE); /* Definition de la montagne... */ GET_ARGUMENTSv(nombre_d_arguments ,BLOC(GET_ARGUMENT_C("imageR=""R=",nom_imageR); GET_ARGUMENT_I("chiffres=",nombre_de_chiffres); GET_ARGUMENT_I("n=""images=""N=",nombre_d_images); GET_ARGUMENT_I("iterations=",nombre_d_iterations); GET_ARGUMENT_F("exposantM=""exposant=",Imandelbrot_dans_C_____exposant); /* Arguments introduits le 20091118122305... */ /* */ /* Le 20100608132420, passage de 'GET_ARGUMENT_I(...)' a 'GET_ARGUMENT_F(...)'. */ GET_ARGUMENT_F("seuilM=""seuil=",Imandelbrot_dans_C_____seuil_du_carre_du_module); /* Arguments introduits le 20091118184202... */ GET_ARGUMENT_F("xbgd=""xbgD=",Reelle(coin_bas_gauche_de_depart_puis_courant)); GET_ARGUMENT_F("ybgd=""ybgD=",Imaginaire(coin_bas_gauche_de_depart_puis_courant)); GET_ARGUMENT_F("xhdd=""xhdD=",Reelle(coin_haut_droite_de_depart_puis_courant)); GET_ARGUMENT_F("yhdd=""yhdD=",Imaginaire(coin_haut_droite_de_depart_puis_courant)); GET_ARGUMENT_F("xbga=""xbgA=",Reelle(coin_bas_gauche_d_arrivee)); GET_ARGUMENT_F("ybga=""ybgA=",Imaginaire(coin_bas_gauche_d_arrivee)); GET_ARGUMENT_F("xhda=""xhdA=",Reelle(coin_haut_droite_d_arrivee)); GET_ARGUMENT_F("yhda=""yhdA=",Imaginaire(coin_haut_droite_d_arrivee)); GET_ARGUMENT_F("precision=",precision); GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie); GET_ARGUMENT_F("modulation=",modulation_du_champ); GET_ARGUMENT_L("depth-cueing=",depth_cueing); GET_ARGUMENT_L("avion=",vue_d_avion); GET_ARGUMENT_L("Imontagnes_precises_____compatibilite_20061220=""compatibilite_20061220=" ,Imontagnes_precises_____compatibilite_20061220 ); /* Introduit le 20061220110516... */ GET_ARGUMENT_F("fc1=""facteur_cretes_1=" ,Imontagnes_precises_____facteur_1_interpolation_toutes_lignes_de_crete ); GET_ARGUMENT_F("fc2=""facteur_cretes_2=" ,Imontagnes_precises_____facteur_2_interpolation_toutes_lignes_de_crete ); /* Introduit le 20061220120334 et etendu le 20061220135323 a deux facteurs... */ CONTROLE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_COMPLEXES; /* Cette extension a ete introduite le 20070129111928... */ ) ); CALi(Inoir(ImageR)); /* Initialisation de l'image Resultat. */ CALS(Iblanc(ImageA2)); /* A2 : texture a appliquer sur la montagne. */ CALCUL_DU_RAPPORT_DE_REDUCTION; DoIn(numero_d_image,PREMIERE_IMAGE_ABSOLUE,LSTX(PREMIERE_IMAGE_ABSOLUE,nombre_d_images),PAS_DES_IMAGES) Bblock Test(IFGE(numero_d_image,PREMIERE_IMAGE)) Bblock /* Afin de se positionner dans la sequence... */ CALS(Imandelbrot_dans_C(ImageA1 ,ADRESSE(coin_bas_gauche_de_depart_puis_courant) ,ADRESSE(coin_haut_droite_de_depart_puis_courant) ,nombre_d_iterations ,visualisation_arguments_de_sortie ) ); /* A1 : partie courante de l'ensemble de Mandelbrot. */ CALS(Iinit_Z_Buffer()); CALS(Imontagnes(ImageR ,modulation_du_champ,ImageA1 ,ImageA2 ,VRAI ,ADRESSE(translation_montagne) ,ombres_portees,largeur_zone_penombre ,ADRESSE(source_lumineuse) ,depth_cueing,min_depth_cueing ,vue_d_avion ,ALIAS ) ); /* R : montagnes correspondant a la partie courante de l'ensemble de Mandelbrot. */ EGAL(nom_image,chain_Aconcaten2_sauf_nom_pipe(nom_imageR,chain_numero(numero_d_image,nombre_de_chiffres))); CALi(Iupdate_image(nom_image,ImageR)); CALZ_FreCC(nom_image); Eblock ATes Bblock Eblock ETes REDUCTION_DE_LA_FENETRE_COURANTE; Eblock EDoI RETU_Commande; Eblock ECommande