####################################################################################################################################### # # # D E F I N I T I O N D ' U N E T R A N S F O R M A T I O N T R I D I M E N S I O N N E L L E N E U T R E : # # # # # # Author of '$xiirf/.MENH.21.5.1.$U' : # # # # Jean-Francois Colonna (LACTAMME, 20180928141941). # # # ####################################################################################################################################### $Z FilSTmpB FFrAcTaL $c #include <stdio.h> $c $c extern double cos(); $c extern double sin(); $c $c #define N0 cPremiere $c #define N cDerniere $c $c #define PIs2 1.570796326794897 $c $c #define OTX 0.0 $c #define ETX PIs2 $c $c #define OTY PIs2 $c #define ETY 0.0 $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 TX=((ETX*(n-N0))+(OTX*(N-n)))/(N-N0); $c double TY=((ETY*(n-N0))+(OTY*(N-n)))/(N-N0); $c double Z=((EZ*(n-N0))+(OZ*(N-n)))/(N-N0); $c $c printf("echo %c$xci/lineaire$X $c A=%f B=%f C=0 $c standard=FAUX $c Frenormaliser=VRAI $c R=$_____ImageD$COORD_X.%04d $c $formatI%c $c >> $FFrAcTaL\n" $c ,'\'' $c ,+cos(TX),+sin(TX) $c ,n $c ,'\'' $c ); $c $c printf("echo %c$xci/lineaire$X $c A=%f B=%f C=0 $c standard=FAUX $c Frenormaliser=VRAI $c R=$_____ImageD$COORD_Y.%04d $c $formatI%c $c >> $FFrAcTaL\n" $c ,'\'' $c ,-cos(TY),+sin(TY) $c ,n $c ,'\'' $c ); $c $c printf("echo %c$xci/init$X $c valeur=%f $c standard=FAUX $c R=$_____ImageD$COORD_Z.%04d $c $formatI%c $c >> $FFrAcTaL\n" $c ,'\'' $c ,Z $c ,n $c ,'\'' $c ); $c /* Definition a priori d'une deformation neutre {DX=X,DY=Y,DZ=Z}. */ $c } $c } $Z execPar_V $FFrAcTaL $Z FilSTmpE FFrAcTaL