#!/bin/csh

#######################################################################################################################################
#                                                                                                                                     #
#         G E N E R A T I O N   D ' I M A G E S   S U R   L E S   F L U C T U A T I O N S   D U   V I D E  :                          #
#                                                                                                                                     #
#                                                                                                                                     #
#         Utilisation :                                                                                                               #
#                                                                                                                                     #
#                   $xrq/vide.L3$Z                                                                                                    #
#                                                                                                                                     #
#                                                                                                                                     #
#              *     * * * * * *   * * * * * *   * * * * * *   *         *   * * * * * *   *   * * * * * *   *         *              #
#                         *             *        *             **        *        *        *   *         *   **        *              #
#             * *         *             *        *             * *       *        *        *   *         *   * *       *              #
#                         *             *        *             *  *      *        *        *   *         *   *  *      *              #
#            *   *        *             *        *             *   *     *        *        *   *         *   *   *     *              #
#                         *             *        * * *         *    *    *        *        *   *         *   *    *    *              #
#           * * * *       *             *        *             *     *   *        *        *   *         *   *     *   *              #
#                         *             *        *             *      *  *        *        *   *         *   *      *  *              #
#          *       *      *             *        *             *       * *        *        *   *         *   *       * *              #
#                         *             *        *             *        **        *        *   *         *   *        **              #
#         *         *     *             *        * * * * * *   *         *        *        *   * * * * * *   *         *              #
#                                                                                                                                     #
#                                                                                                                                     #
#         ATTENTION :                                                                                                                 #
#                                                                                                                                     #
#                     Les parametres ici definis ont ete utilises                                                                     #
#                   le 20031106103500 pour generer les images                                                                         #
#                   'v $xiirq/VIDE.11' et 'v $xiirq/VIDE.12'.                                                                         #
#                   Il convient donc de ne pas les modifier...                                                                        #
#                   Il en est de meme le 20031107090054 en ce qui                                                                     #
#                   concerne la sequence 'v _____xivPdf_01_2/015152_015351'.                                                          #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xrq/vide.L3$Z' :                                                                                                #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, 20031103141146).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

source    $xrq/proton.K2.2$Y
                                        # Definition du choix des machines a utiliser pour generer le film.                           #
source    $xrq/nucleon.L4.1$Y
                                        # Preparation eventuelle de l'enregistrement...                                               #
source    $xrq/nucleon.L4.2$Y
                                        # Quelques parametres...                                                                      #

set       imageO=1
                                        # Definition du numero de la premiere image.                                                  #
set       imageV=$imageO
                                        # Definition de la premiere image a visualiser.                                               #
#20031104153210____:set       Nimages=100                                                                                             #
set       Nimages=200
                                        # Definition du nombre d'images a generer.                                                    #
set       simultanees=$INFINI
                                        # Nombre d'images simultanees autorisees.                                                     #

set       positionnement="p0=$imageO pV=$imageV"

set       convolution="noyau_variable=FAUX cNOIR=FAUX convolution=1 facteur=0.100"
set       accumulation="trainees=VRAI renormaliser=VRAI attenuation_des_trainees=0.40 attenuation_des_images=0.60 mode=0"
                                        # ATTENTION : il ne faut pas que les trainees soient trop importantes a cause de l'effet      #
                                        # d'oscillation horizontale, car alors cela induit un effet visuel de persistance penible...  #

set       evenements="graine=1789"
set       evenements="$evenements ""ralentissement=0.1100 particules=600"
set       evenements="$evenements ""auto_reguler=VRAI Eregulation=1 Cregulation=40.0"
set       evenements="$evenements ""nefqv=2400"
set       evenements="$evenements ""majorer=VRAI"
#20031104153210____:set       evenements="$evenements ""T=5.211938987471126e-26"                                                      #
set       evenements="$evenements ""T=2.605969493735563e-25"
                                        # Definition des parametres fondamentaux.                                                     #

set       imageN = $imageO
@         imageN = $imageN + $Nimages
@         imageN = $imageN - 1
                                        # Definition du numero de la derniere image.                                                  #

set       stereo0="0.00"
set       stereoD="0.04"
set       stereoG="-0.04"

set       stereo=$stereo0
set       Stereo="stereo=$stereo"
                                        # Choix du mode stereoscopique.                                                               #
                                        #                                                                                             #
                                        #                   'stereo0' --> absent,                                                     #
                                        #                   'stereoD' --> DROITE,                                                     #
                                        #                   'stereoG' --> GAUCHE.                                                     #
                                        #                                                                                             #

set       InomT="VIDE"
                                        # On choisit un nom court afin de ne pas avoir de problemes avec la                           #
                                        # constante 'LONGUEUR_MAXIMALE_D_UN_NOM_RELATIF_DE_FICHIER'.                                  #
switch    ("$stereo")
                                        # Definition du nom relatif des images.                                                       #
     case "$stereo0":
          breaksw
     case "$stereoD":
          set       InomT="$InomT$DROITE"
          breaksw
     case "$stereoG":
          set       InomT="$InomT$GAUCHE"
          breaksw
     default:
          EROR      "Mode stereoscopique non reconnu."
          breaksw
endsw

source    $xrq/proton.K2.1$Y
                                        # Definition du choix du directory de generation des images.                                  #

source    $xrq/nucleon.L2.3$Y
                                        # Pour fixer ou modifier des parametres de derniere minute...                                 #

set       images="$DnomT$InomT"

set       generation="R=$images simultanees=$simultanees synchroniser=FAUX trace=FAUX messages=FAUX"
#20031102160150____:set       generation="$generation ""parametres=VRAI"                                                              #
set       generation="$generation ""parametres=FAUX"
set       generation="$generation ""chronometrer=FAUX"
set       generation="$generation ""randomiser=FAUX interieur=1.0 exterieur=1.0 seuil=0.9 vqR=VRAI vaqR=VRAI"
set       generation="$generation ""couronne=0"
set       generation="$generation ""rqR=0.0280e-15 rqV=0.0280e-15 rg=0.0093e-15 nnsG=16 npsG=16 dvG=0.1"
#20031104153456____:set       generation="$generation ""N_AU_CARRE=VRAI"                                                              #
set       generation="$generation ""N_AU_CARRE=FAUX"
#20031105124011____:set       generation="$generation ""binariser_couleur=FAUX"                                                       #
set       generation="$generation ""binariser_couleur=VRAI"
                                        # Definition des arguments d'entrees-sorties ; plusieurs possibilites :                       #
                                        #                                                                                             #
                                        #         nnsG=16 npsG=16 dvG=0.1     : bien esthetiquement, mais tres lent,                  #
                                        #         nnsG=8  npsG=8  dvG=0.2                                                             #
                                        #         nnsG=4  npsG=4  dvG=0.4     : quatre fois plus rapide, et degradation acceptable... #
                                        #                                                                                             #

set       generation="$generation ""queue=VRAI"
                                        # Introduit le 20031101102343 afin de generer des noms "standards" a priori...                #

@         npA = $imageO
                                        # Definition de 'npA' (numero de la premiere periode),                                        #

@         npX = $imageN
                                        # Definition de 'npX' (numero de la derniere periode).                                        #

@         npB = $npX * 7
@         npB = $npB / 17

@         npC = $npX * 17
@         npC = $npC / 17

if        ($npB <= $npA) then
          @         npB = $npA + 1
else
endif
if        ($npC <= $npB) then
          @         npC = $npB + 1
else
endif

alias     valeurs   'set       valD=(\!:1) ; set valA=(\!:2) ; set Fval="\!:3"                                                    ;    \\
                     set       D1=1   ; set D2=`expr $D1 + 1` ; set D3=`expr $D2 + 1`                                             ;    \\
                     set       A1=$D1 ; set A2=$D2            ; set A3=$D3                                                        ;    \\
                     set       Depart  = "D=$valD[$D1] vD=$valD[$D2] dD=$valD[$D3]"                                               ;    \\
                     set       Arrivee = "A=$valA[$A1] vA=$valA[$A2] dA=$valA[$A3]"                                               ;    \\
                     $xci/valeurs_inte$X $Depart $Arrivee cubique=VRAI                                                  >>! $Fval ;    \\
                     unset     valD valA Fval D1 D2 D3 A1 A2 A3 Depart Arrivee'
                                        # Generation incrementale des fichiers des transformations par interpolation cubique.         #

set       Rayon = 0.50
                                        # Argument :                                                                                  #
                                        #                                                                                             #
                                        #                   R                                                                         #
                                        #                                                                                             #
set       Expo = -0.09543
                                        # Argument :                                                                                  #
                                        #                                                                                             #
                                        #                   d = -3 / (10.pi)                                                          #
                                        #                   d = -0.09543                                                              #
                                        #                                                                                             #
set       Beta = 1.00
                                        # Argument :                                                                                  #
                                        #                                                                                             #
                                        #                   b                                                                         #
                                        #                                                                                             #

set       luminance = "$K_VIDE"
set       luminance = "$luminance ""$K_VIDE"

set       translate = "$K_VIDE"
set       translate = "$translate ""$K_VIDE"

FileTmpB  Fzoom
set       zoom1     = ($npA 1.000000 00.000000)
set       zoom2     = ($npC 1.000000 00.000000)
valeurs   $zoom1    $zoom2    $Fzoom
set       zoom      = "$K_VIDE"
set       zoom      = "$zoom ""ZOOM=$Fzoom"

set       rot_OX    = "$K_VIDE"
set       rot_OX    = "$rot_OX ""$K_VIDE"

set       rot_OY    = "$K_VIDE"
set       rot_OY    = "$rot_OY ""$K_VIDE"

set       rot_OZ    = "$K_VIDE"
set       rot_OZ    = "$rot_OZ ""$K_VIDE"

set       Lumiere   = "$K_VIDE"
set       Lumiere   = "$Lumiere ""Lx=0"
set       Lumiere   = "$Lumiere ""Ly=0"
set       Lumiere   = "$Lumiere ""Lz=2"
                                        # Definition de la source lumineuse, telle que pratiquement toutes les particules soient      #
                                        # derriere elle par rapport a l'observateur...                                                #

set       transformation="$K_VIDE"
set       transformation="$transformation $luminance"
set       transformation="$transformation $translate"
set       transformation="$transformation $zoom"
set       transformation="$transformation $rot_OX $rot_OY $rot_OZ"
                                        # Definition de la suite des transformations.                                                 #

Pal
                                        # Definition du format des images 'Pal'...                                                    #

$DELETE_VALIDE      $images*                                                                                            >& $nul

nohup     $xrq/vide.L0$X      np=$Nimages $evenements $transformation $convolution $accumulation $positionnement $Lumiere               \
                              $Stereo $generation                                                                                       \
                                                                                                                        $formatI
                                        # Simulation...                                                                               #

FileTmpE  Fzoom
                                        # Et elimination des fichiers des transformations...                                          #



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