#######################################################################################################################################
# #
# D I F F U S I O N T R I D I M E N S I O N N E L L E " D O U C E " #
# ( C A L C U L D E L A D E N S I T E D E P A R T I C U L E S ) : #
# #
# #
# Author of '$xiirk/.DIFZ.11.2.$U' : #
# #
# Jean-Francois COLONNA (LACTAMME, 20070602090042). #
# #
#######################################################################################################################################
$Z set NPoInTsX=`$CA $FLiStE_XY$COORD_X | $WCl`
$Z set NPoInTsY=`$CA $FLiStE_XY$COORD_Y | $WCl`
$Z set NPoInTsZ=`$CA $FLiStE_XY$COORD_Z | $WCl`
#20220615102636____:$Z set NPoInTs_=`$xcg/MAX3.01$X x=$NPoInTsX y=$NPoInTsY z=$NPoInTsZ` #
$Z set NPoInTs_=`$xcg/MIN3.01$X x=$NPoInTsX y=$NPoInTsY z=$NPoInTsZ`
$c #include <stdio.h>
$c
$c #define N0 _____Premiere
$c #define N _____Derniere
$c
$c #define Z0 _____Z0
$c #define ZN _____ZN
$c
$c int main()
$c {
$c int n;
$c
$c for (n=N0 ; n<=N ; n++)
$c {
$c double coordonnee_Z=((ZN*(n-N0))+(Z0*(N-n)))/(N-N0);
$c
$c printf("if (`$xcg/parallele.01$X verrou=$_____imagesWW.%04d` == $EXIST) then\n"
$c ,n
$c );
$c /* Gestion du parallelisme... */
$c printf(" $xrv/densite.01$X
$c ne=$NPoInTs_
$c LISTE_X=$FDLiStE_XY$COORD_X
$c LISTE_Y=$FDLiStE_XY$COORD_Y
$c LISTE_Z=$FDLiStE_XY$COORD_Z
$c Z=%f
$c torique=$_____EspaceTore
$c toriqueX=$_____EspaceToreX
$c toriqueY=$_____EspaceToreY
$c toriqueZ=$_____EspaceToreZ
$c rectangulaire=$_____Rcarte grand_carre=$_____GCcarte
$c Eexposant=2
$c ignorer_grandes_distances=$_____IGDcarte
$c seuil_grandes_distances=$_____SGDcarte
$c maximum=$_____Scarte
$c facteur=$_____FDcarte Efacteur=$_____FEcarte
$c R=$_____imagesWW.%04d
$c $formatI\n"
$c ,coordonnee_Z
$c ,n
$c );
$c /* Calcul de la densite des particules aleatoires suivant le gradient... */
$c /* */
$c /* Le 20070702131943 '$FLiStE_XY' a ete remplace par '$FDLiStE_XY' afin de permettre de */
$c /* faire des rotations du champ de densite... */
$c
$c printf(" if ($_____Renormaliser == $EXIST) then\n");
$c printf(" $xci/normalise.01$X
$c A=$_____imagesWW.%04d
$c R=$_____imagesWW.%04d
$c $formatI\n"
$c ,n
$c ,n
$c );
$c /* Normalisation de la densite des particules aleatoires si besoin est... */
$c printf(" else\n");
$c printf(" endif\n");
$c printf("else\n");
$c printf("endif\n");
$c }
$c }
#20070605085104____:$Z FileTmpE FLiStE_XY #