#######################################################################################################################################
#                                                                                                                                     #
#         G E N E R A T I O N   D ' U N E   S T R U C T U R E   F R A C T A L E   T R I D I M E N S I O N N E L L E  :                #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xiirf/.FRB3.32.2.1.$U' :                                                                                        #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, 20120802091040).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$Z                  set                 D_____ImageT2=`GetPath $_____ImageT2`
$Z                  Grelatif            $D_____ImageT2 $xp

$Z                  IfBasic             ($nomR != $D_____ImageT2) ThenBasic
$Z                                      # Test introduit le 20190922092217 : lorsque '$_____ImageT2' est dans '$xp', on fait          #
$Z                                      # l'hypothese que la sequence fractale '$_____ImageT2' existe deja...                         #
$Z                  ElseBasic
$Z                                      # Cas ou '$_____ImageT2' n'est pas dans '$xp' : il faut calculer la sequence fractale         #
$Z                                      # '$_____ImageT2'...                                                                          #
$Z                            ATTENTION "Peut-etre utiliser $K_QS'_____ImageT2=$xiff/FRB3.32.0125x0125x0125$D/SECTIONS.1'$K_QS \?"
$Z                                      # Introduit le 20190924125247...                                                              #

$Z                  SETParam            _____MailleX                  0.125
$Z                                      # Parametre introduit le 20150308095530...                                                    #

$Z                  SETParam            _____MailleY                  0.125
$Z                                      # Parametre introduit le 20150308095530...                                                    #

$Z                  SETParam            _____MailleZ                  0.125
$Z                                      # Parametre introduit le 20150308095530...                                                    #

$Z                  SETParam            _____Iterations               16
$Z                                      # Parametre introduit le 20150308095530...                                                    #

$Z                                      FilSTmpB            FFrAcTaL

$c                                      #include  <stdio.h>
$c
$c                                      #define   N0        cPremiere
$c                                      #define   N         cDerniere
$c
$c                                      #define   OZ        0.0
$c                                      #define   EZ        1.0
$c
$c                                      int       main()
$c                                                {
$c                                                int       n;
$c
$c                                                for       (n=N0 ; n<=N ; n++)
$c                                                          {
$c                                                          double    Z=((EZ*(n-N0))+(OZ*(N-n)))/(N-N0);
$c
$c                                                          printf("echo      %c$xci/fract_3D.01$X
$c                                                                                        mx=%f
$c                                                                                        my=%f
$c                                                                                        mz=%f
$c                                                                                        Z=%.16f
$c                                                                                        iterations=%d
$c                                                                                        standard=VRAI
$c                                                                                        R=$_____ImageT2.%04d
$c                                                                                                  $formatI%c
$c                                                                                             >> $FFrAcTaL\n"
$c                                                                ,'\''
$c                                                                ,_____MailleX
$c                                                                ,_____MailleY
$c                                                                ,_____MailleZ
$c                                                                ,Z
$c                                                                ,_____Iterations
$c                                                                ,n
$c                                                                ,'\''
$c                                                                 );
$c                                      /* Aux environs du 20150309091210 je note qu'au format 'Pdu' et avec l'ensemble des          */
$c                                      /* parametres par defaut, sur '$CMAP28 -ex "porte-brancion"-), le calcul de la premiere      */
$c                                      /* image '$_____ImageT2.0001' demande 747.153 secondes (soit entre 12 et 13 minutes)...      */
$c                                                          }
$c                                                }

$Z                                      execPar_V           $FFrAcTaL

$Z                                      FilSTmpE            FFrAcTaL
$Z                  EndifBasic

$Z                  unset               nomR
$Z                  unset               D_____ImageT2



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