#######################################################################################################################################
#                                                                                                                                     #
#         G E N E R A T I O N   D ' U N   M I L I E U   E V O L U T I F                                                               #
#         A   C O U R B U R E   N E G A T I V E   P U I S   P O S I T I V E                                                           #
#         E T   A N T I - A L I A S E  :                                                                                              #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xivP/disk.000000009/face.2/.REFL.E.12.3.$U' :                                                                                  #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$Z                  setParam            _____StabiliteDebut           0.25
$Z                                      # Parametre introduit le 20021218142004.                                                      #

$Z                  setParam            _____StabiliteFin             0.75
$Z                                      # Parametre introduit le 20021218142004.                                                      #

$Z                  $DELETE             $xTV/INTER_0_1.1

$Z                  set                 PremiereDebut=$Premiere
$Z                  set                 DerniereDebut=`$xcg/MUL2.01$X a=$_____StabiliteDebut b=$Derniere entier=VRAI signe="$K_VIDE"`
$Z                  set                 PremiereMilieu=`$xcg/ADD2.01$X a=$DerniereDebut b=1 entier=VRAI signe="$K_VIDE"`
$Z                  set                 PremiereFin=`$xcg/MUL2.01$X a=$_____StabiliteFin b=$Derniere entier=VRAI signe="$K_VIDE"`
$Z                  set                 DerniereMilieu=`$xcg/SOUS.01$X a=$PremiereFin b=1 entier=VRAI signe="$K_VIDE"`
$Z                  set                 DerniereFin=$Derniere

$Z                  $xci/valeurs_inte$X premiere=$PremiereDebut  derniere=$DerniereDebut                                                \
$Z                                      vD=0 vA=0                                                                                       \
$Z                                      cubique=FAUX                                                                                    \
$Z                                                                                                            >>! $xTV/INTER_0_1.1
$Z                  $xci/valeurs_inte$X premiere=$PremiereMilieu derniere=$DerniereMilieu                                               \
$Z                                      vD=0 vA=1                                                                                       \
$Z                                      cubique=FAUX                                                                                    \
$Z                                                                                                            >>! $xTV/INTER_0_1.1
$Z                  $xci/valeurs_inte$X premiere=$PremiereFin    derniere=$DerniereFin                                                  \
$Z                                      vD=1 vA=1                                                                                       \
$Z                                      cubique=FAUX                                                                                    \
$Z                                                                                                            >>! $xTV/INTER_0_1.1
$Z                                      # Cette facilite introduite le 20021218142004 permet d'avoir au debut et a la fin de la       #
$Z                                      # simulation une geometrie fixe...                                                            #

$Z                  SETParam            _____EInterpole               5.3

$Z                  $DELETE             $xTV/INTER_0_1.2

$Z                  $xrv/PUIX.01$X                                                                                                      \
$Z                                      ne=0                                                                                            \
$Z                                      fichier=$xTV/INTER_0_1.1                                                                        \
$Z                                      exposant=$_____EInterpole                                                                       \
$Z                                                                                                            >   $xTV/INTER_0_1.2

$Z                  SETParam            _____A_tc                     98.5

$Z                  SETParam            _____B_tc                     1.5

$Z                  $DELETE             $xTV/PARAMETRE_tc

$Z                  $xrv/AXPB.01$X                                                                                                      \
$Z                                      ne=0                                                                                            \
$Z                                      fichier=$xTV/INTER_0_1.2                                                                        \
$Z                                      a=$_____A_tc                                                                                    \
$Z                                      b=$_____B_tc                                                                                    \
$Z                                                                                                            >   $xTV/PARAMETRE_tc

$Z                  SETParam            _____A_Exy                    23.64

$Z                  SETParam            _____B_Exy                    0.7

$Z                  $DELETE             $xTV/PARAMETRE_Exy

$Z                  $xrv/AXPB.01$X                                                                                                      \
$Z                                      ne=0                                                                                            \
$Z                                      fichier=$xTV/INTER_0_1.2                                                                        \
$Z                                      a=$_____A_Exy                                                                                   \
$Z                                      b=$_____B_Exy                                                                                   \
$Z                                                                                                            >   $xTV/PARAMETRE_Exy

$Z                  $DELETE             $xTV/PARAMETRES.1

$Z                  $PAST                                                                                                               \
$Z                                      $xTV/PARAMETRE_tc                                                                               \
$Z                                      $xTV/PARAMETRE_Exy                                                                              \
$Z                                      $xTV/PARAMETRE_Exy                                                                        |     \
$Z                  $R                                                                                                                  \
$Z                                      "$K_TAB" "$K_BLANC"                                                                       |     \
$Z                  $AW                                                                                                                 \
$Z                                      ' { print "tc=" $1 " Ex=" $2 " Ey=" $3 } '                                                |     \
$Z                  $SE                                                                                                                 \
$Z                                      -e 's+^+$xci/trefle$X +'                                                                        \
$Z                                                                                                            >   $xTV/PARAMETRES.1

$Z                  SETParam            _____STrefle                  0.985

$Z                  $DELETE             $xTV/PARAMETRES.2

$Z                  $xci/genere$X       RVB=FAUX en_tete=FAUX                                                                           \
$Z                                      premiere=$Premiere derniere=$Derniere                                                           \
$Z                                      C='                                                                                           \\\
$Z                                                   standard=FAUX                                                                    \\\
$Z                                                                                                                      $formatI  |   \\\
$Z                                         $xci/normalise.01$X                                                                        \\\
$Z                                                                                                                      $formatI  |   \\\
$Z                                         $xci/passe_bande$X                                                                         \\\
$Z                                                   standard=FAUX                                                                    \\\
$Z                                                   Finf=$_____STrefle Fniveau_inferieur=$_____STrefle                               \\\
$Z                                                   Fsup=1             Fniveau_superieur=1                                           \\\
$Z                                                                                                                      $formatI  |   \\\
$Z                                         $xci/acces$X                                                                               \\\
$Z                                                   standard=FAUX                                                                    \\\
$Z                                                   zero=FAUX                                                                        \\\
$Z                                                   R=$xTV/BOITE.2.%04d                                                              \\\
$Z                                                                                                                      $formatI'       \
$Z                                                                                                            >   $xTV/PARAMETRES.2

$Z                  $DELETE             $xTV/PARAMETRES$Y

$Z                  $CA                 $xcs/csh$Y                                                                                      \
$Z                                                                                                            >   $xTV/PARAMETRES$Y
$Z                  $PAST                                                                                                               \
$Z                                      $xTV/PARAMETRES.1                                                                               \
$Z                                      $xTV/PARAMETRES.2                                                                         |     \
$Z                  $R                                                                                                                  \
$Z                                      "$K_TAB" "$K_BLANC"                                                                             \
$Z                                                                                                            >>! $xTV/PARAMETRES$Y

$Z                  source              $xTV/PARAMETRES$Y

$c                  #include  <stdio.h>
$c
$c                  extern    int       atoi();
$c                  extern    double    atof();
$c                  extern    char      *getenv();
$c
$c                  #define   NA0       (atoi(getenv("PremiereA")))
$c                  #define   NA        (atoi(getenv("DerniereA")))
$c
$c                  int       main()
$c                            {
$c                            int       n;
$c
$c                            for       (n=NA0 ; n<=NA ; n++)
$c                                      {
$c                                      printf("$xci/vraies_C$X
$c                                                                  A=$xTV/BOITE.2.%04d
$c                                                                  p=$xiP/gris.11.1s2
$c                                                                  RVB=VRAI
$c                                                                  R=$xTV/BOITE.3.%04d
$c                                                                                                $formatI\n"
$c
$c                                            ,n,n
$c                                             );
$c                                      /* Definition du decor "dynamique" dans '$xTV/BOITE.3.????.?'.                               */
$c
$c                                      printf("$xci/seuil$X
$c                                                                  A=$xTV/BOITE.2.%04d
$c                                                                                                $formatI  |
$c                                              $xci/complement$X
$c                                                                  R=$xTV/BOITE.2.%04d
$c                                                                                                $formatI\n"
$c
$c                                            ,n,n
$c                                             );
$c                                      /* Et enfin generation finale du niveau dynamique...                                         */
$c                                      }
$c                            }

$Z                  $DELETE             $xTV/INTER_0_1.1
$Z                  $DELETE             $xTV/INTER_0_1.2
$Z                  $DELETE             $xTV/PARAMETRE_tc
$Z                  $DELETE             $xTV/PARAMETRE_Exy
$Z                  $DELETE             $xTV/PARAMETRES.1
$Z                  $DELETE             $xTV/PARAMETRES.2
$Z                  $DELETE             $xTV/PARAMETRES$Y



Copyright © Jean-François Colonna, 2019-2021.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 2019-2021.