#!/bin/csh

#######################################################################################################################################
#                                                                                                                                     #
#         C A L C U L   G E N E R A L   D E   L ' E S P A C E   D E   L Y A P U N O V   T R I D I M E N S I O N N E L                 #
#         A V E C   R E N O R M A L I S A T I O N   G L O B A L E   D E S   I M A G E S  :                                            #
#                                                                                                                                     #
#                                                                                                                                     #
#         Utilisation :                                                                                                               #
#                                                                                                                                     #
#              $xrk/lyapunov.22$Z  <ImageR> <Premiere> <Derniere> <R3d> <R3a> <R4d> <R4a> <R5d> <R5a> <R6d> <R6a> <R7d> <R7a> <LArgs> #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xrk/lyapunov.22$Z' :                                                                                            #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, 20151007130140).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

set       ImageR=$1
set       Premiere=$2
set       Derniere=$3
set       R3d=$4
set       R3a=$5
set       R4d=$6
set       R4a=$7
set       R5d=$8
set       R5a=$9
set       R6d=$10
set       R6a=$11
set       R7d=$12
set       R7a=$13
set       LArguments=($14)
                                        # Arguments de '$xrk/lyapunov.22$X'.                                                          #

@         Nimages = $Derniere - $Premiere
@         Nimages = $Nimages + 1

FileTmpB  LiStE_PlAnS

set       R3_Courant=$R3d
set       R4_Courant=$R4d
set       R5_Courant=$R5d
set       R6_Courant=$R6d
set       R7_Courant=$R7d

set       NuMeRo=$Premiere

set       ArGuMeNtS="$K_VIDE"
set       ArGuMeNtS="$ArGuMeNtS"" standard=FAUX"
set       ArGuMeNtS="$ArGuMeNtS"" negatifs=VRAI"
set       ArGuMeNtS="$ArGuMeNtS"" positifs=FAUX"
set       ArGuMeNtS="$ArGuMeNtS"" editer_fenetre_courante=FAUX"
set       ArGuMeNtS="$ArGuMeNtS"" NombreMaximalChiffresNumeros=0"

while     ($NuMeRo <= $Derniere)
          set       ImAgE="$ImageR""."`$xci/nombres$X premiere=$NuMeRo derniere=$NuMeRo`
          set       Rs_Courant="$K_VIDE"
          set       Rs_Courant="$Rs_Courant"" R3d=$R3_Courant R3a=$R3_Courant"
          set       Rs_Courant="$Rs_Courant"" R4=$R4_Courant"
          set       Rs_Courant="$Rs_Courant"" R5=$R5_Courant"
          set       Rs_Courant="$Rs_Courant"" R6=$R6_Courant"
          set       Rs_Courant="$Rs_Courant"" R7=$R7_Courant"

#20151020093713____:echo "$xrk/lyapunov.22$X  images=1 R3d=$R3_Courant R3a=$R3_Courant $ArGuMeNtS $LArguments R=$ImAgE $formatI"   \  #
#20151020093713____:                                                                                               >>! $LiStE_PlAnS   #

          echo      "$xrk/lyapunov.22$X  images=1 $Rs_Courant $ArGuMeNtS $LArguments R=$ImAgE $formatI"                                 \
                                                                                                              >>! $LiStE_PlAnS
                                        # Generation des commandes de construction de l'espace de Lyapunov...                         #

          if        ($R3a != $R3d) then
                                        # Optimisation introduite le 20151102111740 car, en effet, 'calcul' prend beaucoup de temps.  #
                    set       R3_Courant=`calcul $R3_Courant+(($R3a-$R3d)/($Nimages-1))`
          else
          endif

          if        ($R4a != $R4d) then
                                        # Optimisation introduite le 20151102111740 car, en effet, 'calcul' prend beaucoup de temps.  #
                    set       R4_Courant=`calcul $R4_Courant+(($R4a-$R4d)/($Nimages-1))`
          else
          endif

          if        ($R5a != $R5d) then
                                        # Optimisation introduite le 20151102111740 car, en effet, 'calcul' prend beaucoup de temps.  #
                    set       R5_Courant=`calcul $R5_Courant+(($R5a-$R5d)/($Nimages-1))`
          else
          endif

          if        ($R6a != $R6d) then
                                        # Optimisation introduite le 20151102111740 car, en effet, 'calcul' prend beaucoup de temps.  #
                    set       R6_Courant=`calcul $R6_Courant+(($R6a-$R6d)/($Nimages-1))`
          else
          endif

          if        ($R7a != $R7d) then
                                        # Optimisation introduite le 20151102111740 car, en effet, 'calcul' prend beaucoup de temps.  #
                    set       R7_Courant=`calcul $R7_Courant+(($R7a-$R7d)/($Nimages-1))`
          else
          endif

          @         NuMeRo = $NuMeRo + 1
                                        # Progression dans l'espace de Lyapunov...                                                    #
end

setenv    FindExec_MessagesBP $NEXIST
$xcg/ExecParallele.01$Z       $LiStE_PlAnS
unsetenv  FindExec_MessagesBP
                                        # Calcul parallele de l'espace de Lyapunov...                                                 #

FileTmpE  LiStE_PlAnS

set       MiNiMuM_Global=+$INFINI
set       MaXiMuM_Global=-$INFINI

set       NuMeRo=$Premiere

while     ($NuMeRo <= $Derniere)
          set       ImAgE="$ImageR""."`$xci/nombres$X premiere=$NuMeRo derniere=$NuMeRo`
          set       ExTrEmA=`$xci/extrema$X standard=FAUX A=$ImAgE Prme=VRAI`

          set       MiNiMuM_Global=`$xcg/MIN2.01$X a=$ExTrEmA[1] b=$MiNiMuM_Global`
          set       MaXiMuM_Global=`$xcg/MAX2.01$X a=$ExTrEmA[2] b=$MaXiMuM_Global`
                                        # Recherche des extrema de l'espace de Lyapunov...                                            #

          @         NuMeRo = $NuMeRo + 1
                                        # Progression dans l'espace de Lyapunov...                                                    #
end

set       execALBUM_Para=FAUX
execALBUM $Premiere $Derniere                                                                                                           \
          $xci/acces$X        A=$ImageR.%04d                                                                                            \
                              standard=FAUX zero=FAUX                                                                                   \
                              extrema=FAUX                                                                                              \
                              minimum=$MiNiMuM_Global                                                                                   \
                              maximum=$MaXiMuM_Global                                                                                   \
                              R=$ImageR.%04d                                                                                            \
                                                                                                              $formatI
                                        # Renormalisation globale...                                                                  #



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