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

$Z                  SET                 Premiere=1
$Z                  SET                 Derniere=128

$Z                  setParam            _____Gauss                    128

$Z                  SETParam            _____SEUIL_0_1                0.1

$Z                  SETENV              SEUIL_0_1           $_____SEUIL_0_1

$Z                  SETParam            _____SEUIL_0_N                0.9

$Z                  SETENV              SEUIL_0_N           $_____SEUIL_0_N

$Z                  SETParam            _____SEUIL_1_1                253

$Z                  SETENV              SEUIL_1_1           $_____SEUIL_1_1

$Z                  SETParam            _____SEUIL_1_N                253

$Z                  SETENV              SEUIL_1_N           $_____SEUIL_1_N

$Z                  SETParam            _____INF_____1                140

$Z                  SETENV              INF_____1           $_____INF_____1

$Z                  SETParam            _____INF_____N                110

$Z                  SETENV              INF_____N           $_____INF_____N

$Z                  SETParam            _____SUP_____1                150

$Z                  SETENV              SUP_____1           $_____SUP_____1

$Z                  SETParam            _____SUP_____N                120

$Z                  SETENV              SUP_____N           $_____SUP_____N

$Z                  SETParam            _____SEUIL_2_1                200

$Z                  SETENV              SEUIL_2_1           $_____SEUIL_2_1

$Z                  SETParam            _____SEUIL_2_N                200

$Z                  SETENV              SEUIL_2_N           $_____SEUIL_2_N

$X                  echo                "neutre"
$X                                      # commande destinee uniquement a mettre a jour l'environnement 'env'...                       #

$Z                  Sdu

$Z                  $xci/gauss$X                                                                                                        \
$Z                                      c=$_____Gauss                                                                                   \
$Z                                      R=$xTV/RANDOM.12                                                                                \
$Z                                                                                                                      $formatI
$Z                  $xci/gauss$X                                                                                                        \
$Z                                      R=$xTV/GAUSS                                                                                    \
$Z                                                                                                                      $formatI

$c                  #include  <stdio.h>
$c
$c                  extern    int       atoi();
$c                  extern    double    atof();
$c                  extern    char      *getenv();
$c
$c                  #define   N0        Premiere
$c                  #define   N         Derniere
$c
$c                  #define   SEUIL_0_1 (atof(getenv("SEUIL_0_1")))
$c                  #define   SEUIL_0_N (atof(getenv("SEUIL_0_N")))
$c
$c                  #define   SEUIL_1_1 (atoi(getenv("SEUIL_1_1")))
$c                  #define   SEUIL_1_N (atoi(getenv("SEUIL_1_N")))
$c
$c                  #define   INF_____1 (atoi(getenv("INF_____1")))
$c                  #define   INF_____N (atoi(getenv("INF_____N")))
$c
$c                  #define   SUP_____1 (atoi(getenv("SUP_____1")))
$c                  #define   SUP_____N (atoi(getenv("SUP_____N")))
$c
$c                  #define   SEUIL_2_1 (atoi(getenv("SEUIL_2_1")))
$c                  #define   SEUIL_2_N (atoi(getenv("SEUIL_2_N")))
$c
$c                  #define   BARd(a,b) ((((((a)*((double)(N-n))))+((b)*((double)(n-N0)))))/((double)(N-N0)))
$c                  #define   BARi(a,b) (int)(((double)((((a)*(N-n)))+((b)*(n-N0))))/((double)(N-N0)))
$c
$c                  int       main()
$c                            {
$c                            int       n;
$c
$c                            for       (n=N0 ; n<=N ; n++)
$c                                      {
$c                                      double    seuil_0_1=BARd(SEUIL_0_1,SEUIL_0_N);
$c                                      int       seuil_1_1=BARi(SEUIL_1_1,SEUIL_1_N);
$c                                      int       inf______=BARi(INF_____1,INF_____N);
$c                                      int       sup______=BARi(SUP_____1,SUP_____N);
$c                                      int       seuil_2_1=BARi(SEUIL_2_1,SEUIL_2_N);
$c
$c                                      printf("$xci/random.02$X
$c                                                          seuil=%f
$c                                                                                                  $formatI  |
$c                                              $xci/filtre.01$X
$c                                                          N=$xTV/RANDOM.12
$c                                                          seuil=%d
$c                                                                                                  $formatI  |
$c                                              $xci/passe_bande$X
$c                                                          inf=%d
$c                                                          sup=%d
$c                                                          R=$xTV/RANDOM.34
$c                                                                                                  $formatI\n"
$c                                            ,seuil_0_1
$c                                            ,seuil_1_1
$c                                            ,inf______
$c                                            ,sup______
$c                                             );
$c
$c                                      printf("$xci/Sy$X
$c                                                          A=$xTV/RANDOM.34
$c                                                                                                  $formatI  |
$c                                              $xci/minmax$X
$c                                                          A1=$xTV/RANDOM.34
$c                                                          R=$xTV/RANDOM.45
$c                                                                                                  $formatI\n"
$c                                             );
$c                                      printf("$xci/Sx$X
$c                                                          A=$xTV/RANDOM.45
$c                                                                                                  $formatI  |
$c                                              $xci/minmax$X
$c                                                          A1=$xTV/RANDOM.45
$c                                                                                                  $formatI  |
$c                                              $xci/acces$X
$c                                                          r=VRAI
$c                                                                                                  $formatI  |
$c                                              $xci/filtre.01$X
$c                                                          N=$xTV/GAUSS
$c                                                          seuil=%d
$c                                                                                                  $formatI  |
$c                                              $xci/format.01$X
$c                                                                                                  $formatI
$c                                                          mode=7
$c                                                          R=$_____imagesR.%04d
$c                                                                                                  $formatR_Pal\n"
$c                                            ,seuil_2_1
$c                                            ,n
$c                                             );
$c                                      }
$c                            }

$Z        Pal



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