####################################################################################################################################### # # # 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