/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   D E   L ' E N S E M B L E   D E   V E R H U L S T   T R I D I M E N S I O N N E L                            */
/*        ( D E U X   D I M E N S I O N S   P O U R   L E   P L A N   C O M P L E X E   E T   U N E                                  */
/*        P O U R   L E   T A U X   D E   C R O I S S A N C E )                                                                      */
/*        D A N S   L E   P L A N   C O M P L E X E  :                                                                               */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrc/verhulst.02$K' :                                                                                           */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 1991??????????).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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  :                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        P A R A M E T R E S  :                                                                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   X_DEPART_GAUCHE                                                                                                               \
                    PARE(0.0)
#define   X_DEPART_DROITE                                                                                                               \
                    PARE(1.0)
#define   Y_DEPART_BAS                                                                                                                  \
                    PARE(-0.5)
#define   Y_DEPART_HAUT                                                                                                                 \
                    PARE(0.5)
                                        /* Definition de la fenetre de visualisation.                                                */

#define   PARTIE_REELLE_DU_PREMIER_TAUX_DE_CROISSANCE                                                                                   \
                    PARE(2.8)
#define   PARTIE_IMAGINAIRE_DU_PREMIER_TAUX_DE_CROISSANCE                                                                               \
                    PARE(0.0)
#define   PARTIE_REELLE_DU_DERNIER_TAUX_DE_CROISSANCE                                                                                   \
                    PARE(3.4)
#define   PARTIE_IMAGINAIRE_DU_DERNIER_TAUX_DE_CROISSANCE                                                                               \
                    PARE(0.0)
                                        /* Definition de la variation du taux de croissance 'R'.                                     */

#define   NOMBRE_D_ITERATIONS                                                                                                           \
                    SOUS(BLANC,NOIR)                                                                                                    \
                                        /* Nombre d'iterations demandees.                                                            */
#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  :                                                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   D E   L ' E N S E M B L E   D E   V E R H U L S T   T R I D I M E N S I O N N E L                            */
/*        ( D E U X   D I M E N S I O N S   P O U R   L E   P L A N   C O M P L E X E   E T   U N E                                  */
/*        P O U R   L E   T A U X   D E   C R O I S S A N C E )  :                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
     DEFV(complexe,coin_bas_gauche);
                                        /* Definition du point situe en bas a gauche de la fenetre.                                  */
     DEFV(complexe,coin_haut_droite);
                                        /* Definition du point situe en haut a droite de la fenetre.                                 */
     DEFV(complexe,premier_taux_de_croissance);
     DEFV(complexe,dernier_taux_de_croissance);
                                        /* Definition de la variation du taux de croissance 'R'.                                     */
     DEFV(complexe,taux_de_croissance);
                                        /* Taux de croissance courant 'R' de l'ensemble.                                             */
     DEFV(Positive,INIT(nombre_d_iterations,NOMBRE_D_ITERATIONS));
                                        /* Nombre maximal d'iterations a effectuer.                                                  */
     DEFV(Logical,INIT(visualisation_arguments_de_sortie,VISUALISATION_ARGUMENTS_DE_SORTIE));
                                        /* Indicateur de choix entre la visualisation du nombre d'iterations ('FAUX'),               */
                                        /* et l'argument de sortie ('VRAI').                                                         */
     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.                                                                   */
     /*..............................................................................................................................*/
     Cinitialisation(coin_bas_gauche,X_DEPART_GAUCHE,Y_DEPART_BAS);
     Cinitialisation(coin_haut_droite,X_DEPART_DROITE,Y_DEPART_HAUT);
                                        /* Definition de la fenetre de depart.                                                       */

     Cinitialisation(premier_taux_de_croissance
                    ,PARTIE_REELLE_DU_PREMIER_TAUX_DE_CROISSANCE,PARTIE_IMAGINAIRE_DU_PREMIER_TAUX_DE_CROISSANCE
                     );
     Cinitialisation(dernier_taux_de_croissance
                    ,PARTIE_REELLE_DU_DERNIER_TAUX_DE_CROISSANCE,PARTIE_IMAGINAIRE_DU_DERNIER_TAUX_DE_CROISSANCE
                     );
                                        /* Definition de la variation implicite du taux de croissance.                               */

     GET_ARGUMENTSi(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("xbgd=""xbgD=",Reelle(coin_bas_gauche));
                         GET_ARGUMENT_F("ybgd=""ybgD=",Imaginaire(coin_bas_gauche));
                         GET_ARGUMENT_F("xhdd=""xhdD=",Reelle(coin_haut_droite));
                         GET_ARGUMENT_F("yhdd=""yhdD=",Imaginaire(coin_haut_droite));
                         GET_ARGUMENT_L("arguments=",visualisation_arguments_de_sortie);
                         GET_ARGUMENT_F("Rtauxd=""RtauxD=""Rrd=""RrD=""RRd=""RRD=",Reelle(premier_taux_de_croissance));
                         GET_ARGUMENT_F("Itauxd=""ItauxD=""Ird=""IrD=""IRd=""IRD=",Imaginaire(premier_taux_de_croissance));
                         GET_ARGUMENT_F("Rtauxa=""RtauxA=""RrA=""Rra=""RRa=""RRA=",Reelle(dernier_taux_de_croissance));
                         GET_ARGUMENT_F("Itauxa=""ItauxA=""Ira=""IrA=""IRa=""IRA=",Imaginaire(dernier_taux_de_croissance));
                         GET_ARGUMENT_L("simple=",Iverhulst_dans_C_____utiliser_la_forme_la_plus_simple);
                                        /* Argument introduit le 20040921112618...                                                   */
                         )
                    );

     CALi(Inoir(ImageR));
                                        /* Initialisation de l'image Resultat.                                                       */

     Cegal(taux_de_croissance,premier_taux_de_croissance);
                                        /* Initialisation du taux de croissance 'R'.                                                 */

     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(Iverhulst_dans_C(ImageR
                                    ,ADRESSE(coin_bas_gauche)
                                    ,ADRESSE(coin_haut_droite)
                                    ,ADRESSE(taux_de_croissance)
                                    ,nombre_d_iterations
                                    ,visualisation_arguments_de_sortie
                                     )
                    );
               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

          INCR(Reelle(taux_de_croissance)
              ,DIVI(SOUS(Reelle(dernier_taux_de_croissance),Reelle(premier_taux_de_croissance))
                   ,TRMU(nombre_d_images)
                    )
               );
          INCR(Imaginaire(taux_de_croissance)
              ,DIVI(SOUS(Imaginaire(dernier_taux_de_croissance),Imaginaire(premier_taux_de_croissance))
                   ,TRMU(nombre_d_images)
                    )
               );
                                        /* Definition du taux de croissance 'R' courant.                                             */
          Eblock
     EDoI

     RETU_Commande;
     Eblock
ECommande



Copyright © Jean-François COLONNA, 2019-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2019-2024.