#######################################################################################################################################
#                                                                                                                                     #
#         C O U R B E   D E   H I L B E R T   B I D I M E N S I O N N E L L E   E T   E N S E M B L E S   D E   J U L I A  :          #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xiirv/.JULK.I2.1.$U' :                                                                                          #
#                                                                                                                                     #
#                     Jean-Francois COLONNA (LACTAMME, 20221205101654).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$Z                  setParam            _____CadreJuliaX              4

$Z                  setParam            _____CadreJuliaY              4

#20221205171517____:$Z                  setParam            _____FacteurX                 1                                           #
$Z                  setParam            _____FacteurX                 2

#20221205171517____:$Z                  setParam            _____FacteurY                 1                                           #
$Z                  setParam            _____FacteurY                 2

#20221205171517____:$Z                  setParam            _____Iterations               20                                          #
$Z                  setParam            _____Iterations               40

#20221205143326____:$Z                  setParam            _____SeuilInf                 4                                           #
#20221205171517____:$Z                  setParam            _____SeuilInf                 15                                          #
#20221205174240____:$Z                  setParam            _____SeuilInf                 24                                          #
$Z                  setParam            _____SeuilInf                 16

$Z                  setParam            _____SeuilSup                 $GRIS_8

#20221207103426____:$Z                  setParam            _____Attenuation              0.15                                        #
$Z                  setParam            _____Attenuation              0.25

$Z                  set                 LiStE_X=`$CA $FiLeTmP_HILBERT_all$COORD_X`
$Z                  set                 LiStE_Y=`$CA $FiLeTmP_HILBERT_all$COORD_Y`

$Z                  set                 NPoInTs_X=$#LiStE_X
$Z                  set                 NPoInTs_Y=$#LiStE_Y

$Z                  set                 NPoInTs__=`$xcg/MIN2.01$X x=$NPoInTs_X y=$NPoInTs_Y signe="$K_VIDE"`

$Z                  set                 CaRrE_1_X=0
$Z                  set                 CaRrE_1_Y=0

$Z                  set                 CaRrE_2_X=1
$Z                  set                 CaRrE_2_Y=1
$Z                                      # Definition du carre parcouru par la courbe de Hilbert...                                    #

$Z                  set                 CeNtRe_X=`$xcg/MOYE.01$X a=$CaRrE_1_X b=$CaRrE_2_X signe="$K_VIDE"`
$Z                  set                 CeNtRe_Y=`$xcg/MOYE.01$X a=$CaRrE_1_Y b=$CaRrE_2_Y signe="$K_VIDE"`

$Z                  set                 index=1

$Z                  $xci/init$X         niveau=$NOIR                                                                                    \
$Z                                      R=$xTV/JULIA.cumul                                                                              \
$Z                                                                                                  $formatI

$Z                  while               ($index <= $NPoInTs__)
$Z                                      SformatI

#20221205143326____:$Z                                      Suq                 `calculINS pow(2,-($_____Profondeur-4))`              #

#20221210094037____:$Z                                      set            FaCtEuR_FoRmAt=`calculINS pow(2,-($_____Profondeur-4))`    #
$Z                                      set                 FaCtEuR_FoRmAt=`calculINS pow(2,(6-$_____Profondeur))`
$Z                                      # Plus la courbe '$_____CourbeFond' est "tortueuse", plus les enesembles de Julia calcules    #
$Z                                      # ci-apres seront petits...                                                                   #
$Z                                      #                                                                                             #
$Z                                      # On rappelle le 20221205143941 que le format de l'image '$_____ImageR' est 'Sdu' et que      #
$Z                                      # le format 'Sud' est choisi parce qu'il est le quart de 'Sdu'...                             #
#20221205165144____:$Z                                      Sud                 $FaCtEuR_FoRmAt                                       #

#20221210094037____:$Z                                      if        ($FaCtEuR_FoRmAt > 0) then                                      #
$Z                                      # Test introduit le 20221209100956 pour 'v $xiirv/JULK.I5', image pour laquelle on a          #
$Z                                      # '$_____Profondeur=5', d'ou '$FaCtEuR_FoRmAt=0'...                                           #

#20221210094037____:$Z                                                switch    ("$formatI")                                          #
#20221210094037____:$Z                                                          case "$format_Std":                                   #
#20221210094037____:$Z                                                                    Suq                 $FaCtEuR_FoRmAt         #
#20221210094037____:$Z                                                                    breaksw                                     #
#20221210094037____:$Z                                                          case "$format_Sdu":                                   #
#20221210094037____:$Z                                                                    Sud                 $FaCtEuR_FoRmAt         #
#20221210094037____:$Z                                                                    breaksw                                     #
#20221210094037____:$Z                                                          case "$format_Squ":                                   #
#20221210094037____:$Z                                                                    Std                 $FaCtEuR_FoRmAt         #
#20221210094037____:$Z                                                                    breaksw                                     #
#20221210094037____:$Z                                                          default:                                              #
#20221210094037____:$Z                                           EROR "['"'$Falias_wm'"'] Le format '$formatI' n'est pas reconnu."    #
#20221210094037____:$Z                                                                    breaksw                                     #
#20221210094037____:$Z                                                endsw                                                           #

#20221210094037____:$Z                                      else                                                                      #

#20221210094037____:$Z                                                switch    ("$formatI")                                          #
#20221210094037____:$Z                                                          case "$format_Std":                                   #
#20221210094037____:$Z                                                                    Suh                 1                       #
#20221210094037____:$Z                                                                    breaksw                                     #
#20221210094037____:$Z                                                          case "$format_Sdu":                                   #
#20221210094037____:$Z                                                                    Suq                 1                       #
#20221210094037____:$Z                                                                    breaksw                                     #
#20221210094037____:$Z                                                          case "$format_Squ":                                   #
#20221210094037____:$Z                                                                    Sud                 1                       #
#20221210094037____:$Z                                                                    breaksw                                     #
#20221210094037____:$Z                                                          default:                                              #
#20221210094037____:$Z                                           EROR "['"'$Falias_wm'"'] Le format '$formatI' n'est pas reconnu."    #
#20221210094037____:$Z                                                                    breaksw                                     #
#20221210094037____:$Z                                                endsw                                                           #
#20221210094037____:$Z                                      endif                                                                     #

$Z                                      if        ($FaCtEuR_FoRmAt > 0) then

$Z                                                switch    ("$formatI")
$Z                                                          case "$format_Std":
$Z                                                                    Sus                 $FaCtEuR_FoRmAt
$Z                                                                    breaksw
$Z                                                          case "$format_Sdu":
$Z                                                                    Suh                 $FaCtEuR_FoRmAt
$Z                                                                    breaksw
$Z                                                          case "$format_Squ":
$Z                                                                    Suq                 $FaCtEuR_FoRmAt
$Z                                                                    breaksw
$Z                                                          default:
$Z                                                                    EROR "['"'$Falias_wm'"'] Le format '$formatI' n'est pas reconnu."
$Z                                                                    breaksw
$Z                                                endsw

$Z                                      else
$Z                                                EROR "['"'$Falias_wm'"'] Le facteur de format '$FaCtEuR_FoRmAt' est invalide."
$Z                                      endif

$Z                                      set                 CoOrDoNnEe_X=$LiStE_X[$index]
$Z                                      set                 CoOrDoNnEe_Y=$LiStE_Y[$index]

$Z                                      set                 TrAnSlAtIoN_X=`calcul $CoOrDoNnEe_X-$CeNtRe_X`
$Z                                      set                 TrAnSlAtIoN_Y=`calcul $CoOrDoNnEe_Y-$CeNtRe_Y`

$Z                                      set                 GeNeRaTeUr_X=`calcul $_____FacteurX*$TrAnSlAtIoN_X`
$Z                                      set                 GeNeRaTeUr_Y=`calcul $_____FacteurY*$TrAnSlAtIoN_Y`

#20221206112807____:$Z                                      $xrc/julia.01$X     N=1                                                \  #
#20221206112807____:$Z                                                          iterations=$_____Iterations                        \  #
#20221206112807____:$Z                                                          gr=$GeNeRaTeUr_X gi=$GeNeRaTeUr_Y                  \  #
#20221206112807____:$Z                                                          xbg=-$_____CadreJuliaX ybg=-$_____CadreJuliaY      \  #
#20221206112807____:$Z                                                          xhd=+$_____CadreJuliaX yhd=+$_____CadreJuliaY      \  #
#20221206112807____:$Z                                                                                             $formatI  |     \  #
#20221206112807____:$Z                                      $xci/passe_bande$X  inf=$_____SeuilInf sup=$_____SeuilSup              \  #
#20221206112807____:$Z                                                                                             $formatI  |     \  #
#20221206112807____:$Z                                      $xci/acces$X        renormaliser=VRAI                                  \  #
#20221206112807____:$Z                                                          R=$xTV/JULIA.1                                     \  #
#20221206112807____:$Z                                                                                             $formatI           #

$Z                                      $xrc/julia.01$X     N=1                                                                         \
$Z                                                          iterations=$_____Iterations                                                 \
$Z                                                          gr=$GeNeRaTeUr_X gi=$GeNeRaTeUr_Y                                           \
$Z                                                          xbg=-$_____CadreJuliaX ybg=-$_____CadreJuliaY                               \
$Z                                                          xhd=+$_____CadreJuliaX yhd=+$_____CadreJuliaY                               \
$Z                                                                                                                      $formatI  |     \
$Z                                      $xci/passe_bande$X  inf=$_____SeuilInf sup=$_____SeuilSup                                       \
$Z                                                          R=$xTV/JULIA.1                                                              \
$Z                                                                                                                      $formatI

$Z                                      $xci/format.01$X    A=$xTV/JULIA.1                                                              \
$Z                                                          mode=7                                                                      \
$Z                                                                                                                      $formatI        \
$Z                                                          R=$xTV/JULIA.1                                                              \
$Z                                                                                                                      $sformatIR

$Z                                      RformatI

$Z                                      $xci/scroll$X       A=$xTV/JULIA.1                                                              \
$Z                                                          trX=$TrAnSlAtIoN_X toreX=FAUX                                               \
$Z                                                          trY=$TrAnSlAtIoN_Y toreY=FAUX                                               \
$Z                                                                                                                      $formatI  |     \
$Z                                      $xci/or_02$X        A1=$xTV/JULIA.cumul                                                         \
$Z                                                          R=$xTV/JULIA.cumul                                                          \
$Z                                                                                                                      $formatI

$Z                                      @                   index = $index + 1
$Z                  end

$Z                  $xci/vraies_C$X     A=$xTV/JULIA.cumul                                                                              \
$Z                                      p=$_____Palette                                                                                 \
$Z                                      R=$xTV/JULIA.cumul                                                                              \
$Z                                                                                                  $formatI

#20221207103426____:$Z        execZ     $xci/cache.31$Z     $_____CourbeFond $xTV/JULIA.cumul $_____ImageR                            #
#20221207113656____:$Z        execZ     $xci/cache.31$Z     $xTV/JULIA.cumul $_____CourbeFond $_____ImageR                            #
$Z                                      # On notera que l'on ne peut pas utiliser 'v $xci/cache.41$Z' car, en effet, celui-ci         #
$Z                                      # introduit des points NOIRs a l'interieur des ensembles de Julia...                          #

#20221207113656____:$Z        execRVB   $xci/convol.01$X    A=$_____ImageR%s                                                       \  #
#20221207113656____:$Z                                      R=$_____ImageR%s                                                       \  #
#20221207113656____:$Z                                                                                                  $formatI      #

$Z                  set                 MaNdElBrOt_1_X=`calcul $_____FacteurX*($CaRrE_1_X-$CeNtRe_X)`
$Z                  set                 MaNdElBrOt_1_Y=`calcul $_____FacteurY*($CaRrE_1_Y-$CeNtRe_Y)`

$Z                  set                 MaNdElBrOt_2_X=`calcul $_____FacteurX*($CaRrE_2_X-$CeNtRe_X)`
$Z                  set                 MaNdElBrOt_2_Y=`calcul $_____FacteurY*($CaRrE_2_Y-$CeNtRe_Y)`

$Z                  $xrc/mandel.01$X    N=1                                                                                             \
$Z                                      xbg=$MaNdElBrOt_1_X ybg=$MaNdElBrOt_1_Y                                                         \
$Z                                      xhd=$MaNdElBrOt_2_X yhd=$MaNdElBrOt_2_Y                                                         \
$Z                                                                                                  $formatI                      |     \
$Z                  $xci/scale$X        a=$_____Attenuation                                                                             \
$Z                                      R=$xTV/MANDELBROT                                                                               \
$Z                                                                                                  $formatI

$Z        execRVB   $xci/neutre$X       A=$xTV/MANDELBROT                                                                               \
$Z                                      R=$xTV/MANDELBROT%s                                                                             \
$Z                                                                                                  $formatI

#20221207115401____:$Z        execZ     $xci/cache.31$Z     $xTV/MANDELBROT $_____ImageR $_____ImageR                                 #
$Z        execZ     $xci/cache.31$Z     $xTV/MANDELBROT $xTV/JULIA.cumul $_____ImageR
$Z                                      # On notera que l'on ne peut pas utiliser 'v $xci/cache.41$Z' car, en effet, celui-ci         #
$Z                                      # introduit des points NOIRs a l'interieur des ensembles de Julia...                          #

$Z                  FileTmpE            FiLeTmP_HILBERT_all



Copyright © Jean-François COLONNA, 2022-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2022-2024.