#######################################################################################################################################
#                                                                                                                                     #
#         D E F I N I T I O N   D ' U N   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  :                      #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xiirc/.LYAP.C3.1.$U' :                                                                                          #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, 20151003100656).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$Z                  setParam            _____LyapunovParam            "$K_VIDE"
$Z                                      # Parametre introduit le 20210817115725...                                                    #

$Z                  SETParam            _____DoublerSections          0
$Z                                      # Parametre introduit le 20210428105809...                                                    #

$Z                  SETParam            _____Compatibilite_20151007   $NEXIST
$Z                                      # Parametre introduit le 20151008084241...                                                    #

$Z                  SETParam            _____EditerPeriode            $EXIST
$Z                                      # Parametre introduit le 20151004092554...                                                    #

$Z                  SETParam            _____FTaux                    $NEXIST
$Z                                      # Parametre introduit le 20151014093215...                                                    #

$Z                  setParam            _____Alpha                    `GetParam $xci/valeurs_trig$X alpha`
$Z                                      # Parametre introduit le 20151006094720...                                                    #

$Z                  setParam            _____Gamma                    1
$Z                                      # Parametre introduit le 20151006094235...                                                    #

$Z                  SETParam            cPremiere                     0

$Z                  SETParam            cDerniere                     255

$Z                  setParam            _____R1bg                     3.000
$Z                                      # Parametre introduit le 20151007090422...                                                    #

$Z                  setParam            _____R2bg                     3.500
$Z                                      # Parametre introduit le 20151007090422...                                                    #

$Z                  setParam            _____R1hd                     3.350
$Z                                      # Parametre introduit le 20151007090422...                                                    #

$Z                  setParam            _____R2hd                     3.850
$Z                                      # Parametre introduit le 20151007090422...                                                    #

$Z                  setParam            _____R3d                      4
$Z                                      # Parametre introduit le 20151007090422...                                                    #

$Z                  setParam            _____R3a                      3
$Z                                      # Parametre introduit le 20151007090422...                                                    #

$Z                  SETParam            _____R4                       `GetParam $xrk/lyapunov.22$X R4`
$Z                                      # Parametre introduit le 20151006095149...                                                    #

$Z                  SETParam            _____R4d                      $_____R4
$Z                                      # Parametre introduit le 20151020092338 et n'ayant de sens qu'en mode parallele,              #
$Z                                      # c'est-a-dire si '$_____Compatibilite_20151007=$NEXIST'...                                   #

$Z                  SETParam            _____R4a                      $_____R4
$Z                                      # Parametre introduit le 20151020092338 et n'ayant de sens qu'en mode parallele,              #
$Z                                      # c'est-a-dire si '$_____Compatibilite_20151007=$NEXIST'...                                   #

$Z                  SETParam            _____R5                       `GetParam $xrk/lyapunov.22$X R5`
$Z                                      # Parametre introduit le 20151006095149...                                                    #

$Z                  SETParam            _____R5d                      $_____R5
$Z                                      # Parametre introduit le 20151020092338 et n'ayant de sens qu'en mode parallele,              #
$Z                                      # c'est-a-dire si '$_____Compatibilite_20151007=$NEXIST'...                                   #

$Z                  SETParam            _____R5a                      $_____R5
$Z                                      # Parametre introduit le 20151020092338 et n'ayant de sens qu'en mode parallele,              #
$Z                                      # c'est-a-dire si '$_____Compatibilite_20151007=$NEXIST'...                                   #

$Z                  SETParam            _____R6                       `GetParam $xrk/lyapunov.22$X R6`
$Z                                      # Parametre introduit le 20151007093155...                                                    #

$Z                  SETParam            _____R6d                      $_____R6
$Z                                      # Parametre introduit le 20151020092338 et n'ayant de sens qu'en mode parallele,              #
$Z                                      # c'est-a-dire si '$_____Compatibilite_20151007=$NEXIST'...                                   #

$Z                  SETParam            _____R6a                      $_____R6
$Z                                      # Parametre introduit le 20151020092338 et n'ayant de sens qu'en mode parallele,              #
$Z                                      # c'est-a-dire si '$_____Compatibilite_20151007=$NEXIST'...                                   #

$Z                  SETParam            _____R7                       `GetParam $xrk/lyapunov.22$X R7`
$Z                                      # Parametre introduit le 20151007093155...                                                    #

$Z                  SETParam            _____R7d                      $_____R7
$Z                                      # Parametre introduit le 20151020092338 et n'ayant de sens qu'en mode parallele,              #
$Z                                      # c'est-a-dire si '$_____Compatibilite_20151007=$NEXIST'...                                   #

$Z                  SETParam            _____R7a                      $_____R7
$Z                                      # Parametre introduit le 20151020092338 et n'ayant de sens qu'en mode parallele,              #
$Z                                      # c'est-a-dire si '$_____Compatibilite_20151007=$NEXIST'...                                   #

$Z                  setParam            _____RotationX                0
$Z                                      # ATTENTION : il est imperatif que les angles de rotation soient nuls sous peine de voir      #
$Z                                      # les surfaces se "dechirer", d'ou l'apparition de zones noires...                            #

$Z                  setParam            _____RotationY                0
$Z                                      # Parametre introduit le 20151014184859...                                                    #
$Z                                      #                                                                                             #
$Z                                      # ATTENTION : il est imperatif que les angles de rotation soient nuls sous peine de voir      #
$Z                                      # les surfaces se "dechirer", d'ou l'apparition de zones noires...                            #

$Z                  setParam            _____RotationZ                0
$Z                                      # Parametre introduit le 20151014184859...                                                    #
$Z                                      #                                                                                             #
$Z                                      # ATTENTION : il est imperatif que les angles de rotation soient nuls sous peine de voir      #
$Z                                      # les surfaces se "dechirer", d'ou l'apparition de zones noires...                            #

$Z                  setParam            _____PaletteBrume             $xiP/gris.01

$Z                  @                   Nimages = $cDerniere - $cPremiere
$Z                  @                   Nimages = $Nimages + 1

$Z                  set                 Nimages_Effectif=`calculINS $Nimages/($_____DoublerSections+1)`
$Z                  set                 cPremiere_Effectif=$cPremiere
$Z                  set                 cDerniere_Effectif=`calculINS $cPremiere+(($cDerniere-$cPremiere+1)/($_____DoublerSections+1))-1`
$Z                                      # Definitions introduites le 20210428105809...                                                #

$Z                  $DELETE             $xTV/TAUX

$Z                  if                  ($_____FTaux == $NEXIST) then
$Z                                      # Test introduit le 20151014093215...                                                         #
$Z                                      set                 NTaux=$_____NTaux

$Z                                      $xci/valeurs_trig$X                                                                             \
$Z                                                          p=1 d=$NTaux                                                                \
$Z                                                          aD=0 aA=`calcul 5*$pi`                                                      \
$Z                                                          alpha=$_____Alpha                                                           \
$Z                                                          gamma=$_____Gamma                                                           \
$Z                                                          entier=VRAI                                                                 \
$Z                                                          epsilon=0.2                                                                 \
$Z                                                                                                                      > $xTV/TAUX
$Z                                      # Generation de la periode utilisee de {0,1,2,...}.                                           #
$Z                  else
$Z                                      set                 NTaux=`$CA $_____FTaux | $WCl`

$Z                                      $CA                 $_____FTaux                                                                 \
$Z                                                                                                                      > $xTV/TAUX
$Z                                      # Utilisation d'un fichier pre-calcule...                                                     #
$Z                  endif

$Z                  if                  ($_____EditerPeriode == $EXIST) then
$Z                                      # Possibilite introduite le 20151004092554...                                                 #
$Z                                      Techo               "Definition de la periode utilisee :"

$Z                                      $xrv/AXPB.01$X                                                                                  \
$Z                                                          ne=0                                                                        \
$Z                                                          fichier=$xTV/TAUX                                                           \
$Z                                                          a=1 b=1                                                                     \
$Z                                                          entier=VRAI                                                                 \
$Z                                                          formater=VRAI signe="$K_VIDE"                                         |     \
$Z                                      $SE                 -e "s/^/R/"                                                           |     \
$Z                                      $R                  "$K_NL" ","                                                           |     \
$Z                                      $xcp/ferme_K_NL$X                                                                         |     \
$Z                                      $SE                 -e 's/,$//'                                                           |     \
$Z                                      $xcp/substitue.01$X                                                                             \
$Z                                                          c01="\(................................................\)"                  \
$Z                                                          s01="\1\n"                                                                  \
$Z                                                          c02="\n\n"                                                                  \
$Z                                                          s02="\n"
$Z                                      # Edition de la periode utilisee de {0,1,2} --> {R1,R2,R3}...                                 #

$Z                                      saut
$Z                  else
$Z                  endif

$Z                  set                 _____ImageT1=$xTV/`GetLeaf $_____imagesR`

$Z                  if                  ($?LArguments == $NEXIST) then
$Z                                      # Test introduit le 20151031100954 pour 'v $xiirc/$Fnota _listG_LYAP_O3_HauteDef'...          #
$Z                                      set                 LArguments="$K_VIDE"
$Z                  else
$Z                  endif

$Z                  set                 LArguments="$LArguments"" $_____LyapunovParam"
$Z                  set                 LArguments="$LArguments"" iterations=$NTaux"
$Z                  set                 LArguments="$LArguments"" LISTE_TAUX=$xTV/TAUX"
$Z                  set                 LArguments="$LArguments"" R1bg=$_____R1bg     R2bg=$_____R2bg"
$Z                  set                 LArguments="$LArguments"" R1hd=$_____R1hd     R2hd=$_____R2hd"

$Z                  if                  ($_____Compatibilite_20151007 == $EXIST) then
$Z                                      # Test introduit le 20151008084418...                                                         #
$Z                                      if                  ($_____R4d != $_____R4a) then
$Z                                                          EROR      "Le calcul sequentiel est impossible avec un 'R4' variable."
$Z                                      else
$Z                                      endif

$Z                                      if                  ($_____R5d != $_____R5a) then
$Z                                                          EROR      "Le calcul sequentiel est impossible avec un 'R5' variable."
$Z                                      else
$Z                                      endif

$Z                                      if                  ($_____R6d != $_____R6a) then
$Z                                                          EROR      "Le calcul sequentiel est impossible avec un 'R6' variable."
$Z                                      else
$Z                                      endif

$Z                                      if                  ($_____R7d != $_____R7a) then
$Z                                                          EROR      "Le calcul sequentiel est impossible avec un 'R7' variable."
$Z                                      else
$Z                                      endif

#20210428105809____:$Z                                      $xrk/lyapunov.22$X                                                     \  #
#20210428105809____:$Z                                                          images=$Nimages                                    \  #
#20210428105809____:$Z                                                          $LArguments                                        \  #
#20210428105809____:$Z                                                          R3d=$_____R3d       R3a=$_____R3a                  \  #
#20210428105809____:$Z                                                          R4=$_____R4                                        \  #
#20210428105809____:$Z                                                          R5=$_____R5                                        \  #
#20210428105809____:$Z                                                          R6=$_____R6                                        \  #
#20210428105809____:$Z                                                          R7=$_____R7                                        \  #
#20210428105809____:$Z                                                          negatifs=VRAI                                      \  #
#20210428105809____:$Z                                                          positifs=FAUX                                      \  #
#20210428105809____:$Z                                                          editer_fenetre_courante=FAUX ChiffresFlot=8        \  #
#20210428105809____:$Z                                                          OX_EX_OY_EY=VRAI                                   \  #
#20210428105809____:$Z                                                          R=$_____ImageT1.                                   \  #
#20210428105809____:$Z                                                                                                       $formatI #

$Z                                      $xrk/lyapunov.22$X                                                                              \
$Z                                                          images=$Nimages_Effectif                                                    \
$Z                                                          $LArguments                                                                 \
$Z                                                          R3d=$_____R3d       R3a=$_____R3a                                           \
$Z                                                          R4=$_____R4                                                                 \
$Z                                                          R5=$_____R5                                                                 \
$Z                                                          R6=$_____R6                                                                 \
$Z                                                          R7=$_____R7                                                                 \
$Z                                                          negatifs=VRAI                                                               \
$Z                                                          positifs=FAUX                                                               \
$Z                                                          editer_fenetre_courante=FAUX ChiffresFlot=8                                 \
$Z                                                          OX_EX_OY_EY=VRAI                                                            \
$Z                                                          R=$_____ImageT1.                                                            \
$Z                                                                                                                      $formatI
$Z                                      # Calcul sequentiel de l'espace de Lyapunov...                                                #
$Z                                      #                                                                                             #
$Z                                      # On rappelle le 20151014142833 que "editer_fenetre_courante=FAUX" permet d'editer la         #
$Z                                      # fenetre de calcul reellement utilisee apres la prise en compte du format des images...      #
$Z                                      #                                                                                             #
$Z                                      # Le 20151014144856 a ete ajoute "OX_EX_OY_EY=VRAI" afin de pouvoir recuperer facilement      #
$Z                                      # les bornes du domaine de definition qui sera [OX,EX]x[OYxEY]...                             #
$Z                                      #                                                                                             #
$Z                                      # Le 20190526081829 a ete introduit "ChiffresFlot=8" au cas ou l'on ferait "a la main"        #
$Z                                      #                                                                                             #
$Z                                      #                   SET       _____Compatibilite_20151007=$EXIST                              #
$Z                                      #                                                                                             #
$Z                                      # et :                                                                                        #
$Z                                      #                                                                                             #
$Z                                      #                   editer_fenetre_courante=VRAI                                              #
$Z                                      #                                                                                             #
$Z                                      # afin de connaitre la fenetre effectivement utilisee lorsque celle-ci n'est pas carree       #
$Z                                      # ('v $xiirc/$Fnota Debut_listG_LYAP_J5_HauteDef' par exemple...).                            #
$Z                  else
#20210428105809____:$Z                                      $xrk/lyapunov.22$Z  $_____ImageT1                                      \  #
#20210428105809____:$Z                                                          $cPremiere $cDerniere                              \  #
#20210428105809____:$Z                                                          $_____R3d $_____R3a                                \  #
#20210428105809____:$Z                                                          $_____R4d $_____R4a                                \  #
#20210428105809____:$Z                                                          $_____R5d $_____R5a                                \  #
#20210428105809____:$Z                                                          $_____R6d $_____R6a                                \  #
#20210428105809____:$Z                                                          $_____R7d $_____R7a                                \  #
#20210428105809____:$Z                                                          "$LArguments"                                         #

$Z                                      $xrk/lyapunov.22$Z  $_____ImageT1                                                               \
$Z                                                          $cPremiere_Effectif $cDerniere_Effectif                                     \
$Z                                                          $_____R3d $_____R3a                                                         \
$Z                                                          $_____R4d $_____R4a                                                         \
$Z                                                          $_____R5d $_____R5a                                                         \
$Z                                                          $_____R6d $_____R6a                                                         \
$Z                                                          $_____R7d $_____R7a                                                         \
$Z                                                          "$LArguments"
$Z                                      # Calcul parallele de l'espace de Lyapunov...                                                 #
$Z                  endif

$Z                  $DELETE             $xTV/TAUX

$Z                  set                 _____ImageT2=$_____ImageT1
$Z                  _____include(_____xiirc/.MANF.DC.5._____U)
$Z                                      # Introduit le 20210428105809...                                                              #

$Z        execALB_F $cPremiere $cDerniere                                                                                               \
$Z                  $xci/seuil$X                                                                                                        \
$Z                                      A=$_____ImageT1.%04d                                                                            \
$Z                                      seuil=30                                                                                        \
$Z                                      R=$_____ImageT1.%04d                                                                            \
$Z                                                                                                  $formatI



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