####################################################################################################################################### # # # G E N E R A T I O N D ' U N E S P H E R E F R A C T A L E # # P A R L A M E T H O D E D E S P R O J E C T O R S : # # # # # # Author of '$xiirs/.PROJ.21.1.$U' : # # # # Jean-Francois Colonna (LACTAMME, 20041111081718). # # # ####################################################################################################################################### $Z setParam _____Rayon 1.0 $Z setParam _____Zoom 0.8 $Z setParam _____RotationX $pis2 $Z setParam _____RotationY 0 $Z # Parametre introduit le 20041125105450... # $Z setParam _____RotationZ 0 $Z # Parametre introduit le 20041125105450... # $Z setParam _____Arguments "$K_VIDE" $Z # Parametre introduit le 20060130140839 pour permettre, par exemple, de generer les normales. # $Z setParam _____Redirection "$K_VIDE" $Z # Parametre introduit le 20191015172410... # $Z if ($?_____imagesRFx == $EXIST) then $Z # Test introduit le 20041113194546... # $Z $xci/acces$X \ $Z A=$xTV/Fx \ $Z standard=FAUX \ $Z zero=FAUX \ $Z R=$_____imagesRFx \ $Z $formatI $Z # La memorisation de '$_____imagesRFx' a ete introduite le 20041112125655... # $Z else $Z endif $Z if ($?_____imagesRFy == $EXIST) then $Z # Test introduit le 20041113194546... # $Z $xci/acces$X \ $Z A=$xTV/Fy \ $Z standard=FAUX \ $Z zero=FAUX \ $Z R=$_____imagesRFy \ $Z $formatI $Z # La memorisation de '$_____imagesRFy' a ete introduite le 20041112125655... # $Z else $Z endif $Z if ($?_____imagesRFz == $EXIST) then $Z # Test introduit le 20041113194546... # $Z $xci/acces$X \ $Z A=$xTV/Fz \ $Z standard=FAUX \ $Z zero=FAUX \ $Z R=$_____imagesRFz \ $Z $formatI $Z # La memorisation de '$_____imagesRFz' a ete introduite le 20041112125655... # $Z else $Z endif $Z if ($?_____imagesRfractalXY == $EXIST) then $Z $xci/acces$X \ $Z A=$xTV/FRACTAL.XY \ $Z standard=FAUX \ $Z zero=FAUX \ $Z R=$_____imagesRfractalXY \ $Z $formatI $Z # La memorisation de '$_____imagesRfractalXY' a ete introduite le 20041114081922... # $Z else $Z endif $Z if ($?_____imagesRfractalZ == $EXIST) then $Z $xci/acces$X \ $Z A=$xTV/FRACTAL.Z \ $Z standard=FAUX \ $Z zero=FAUX \ $Z R=$_____imagesRfractalZ \ $Z $formatI $Z # La memorisation de '$_____imagesRfractalZ' a ete introduite le 20041114081922... # $Z else $Z endif $Z set Parametres="$K_VIDE" $Z set Parametres="$Parametres"" np=1" $Z set Parametres="$Parametres"" brume=FAUX" $Z set Parametres="$Parametres"" pu=0.02 pv=0.02" $Z set Parametres="$Parametres"" eu=20 ev=20" $Z set Parametres="$Parametres"" dpu=40 dpv=40" $Z set Parametres="$Parametres"" ZOOM=$_____Zoom" $Z set Parametres="$Parametres"" ROTATION_OX=$_____RotationX" $Z set Parametres="$Parametres"" ROTATION_OY=$_____RotationY" $Z set Parametres="$Parametres"" ROTATION_OZ=$_____RotationZ" $Z set Parametres="$Parametres"" Lz=100" $Z set Parametres="$Parametres"" rayon_de_visualisation=0.015" $Z set Parametres="$Parametres"" N_AU_CARRE=VRAI" $Z set Parametres="$Parametres"" edc=VRAI" $Z eval $xrs/project2D.11$X \ $Z aFx=$_____Rayon Fx=$xTV/Fx \ $Z aFy=$_____Rayon Fy=$xTV/Fy \ $Z aFz=$_____Rayon Fz=$xTV/Fz \ $Z mu=0 Mu=0.5 \ $Z mv=0 Mv=1.0 \ $Z $Parametres \ $Z chiffres=0 \ $Z R=$_____imagesR \ $Z $_____Arguments \ $Z $formatI \ $Z $_____Redirection $Z # Rappelons que la sphere "unite" est definie par : # $Z # # $Z # F (u,v) = sin(u).cos(v) # $Z # x # $Z # # $Z # F (u,v) = sin(u).sin(v) # $Z # y # $Z # # $Z # F (u,v) = cos(u) # $Z # z # $Z # # $Z # avec : # $Z # # $Z # u E [ 0 , pi ] # $Z # # $Z # v E [ 0 , 2.pi ] # $Z # # $Z # On voit que 'u' n'est definie que dans [0,pi], c'est pourquoi, on trouve ci-dessus : # $Z # # $Z # mu=0 Mu=0.5 # $Z # # $Z # de facon a n'exploiter que la moitie inferieure (partie gauche des images) des deux champs # $Z # '$xTV/SINUS.u' et '$xTV/COSINUS.u' qui ont donc ete calcules dans [0,2.pi]. On ne peut # $Z # malheureusement generer ces deux derniers dans [0,pi] a cause de la derivation numerique # $Z # 'v $xrs/surfaces.12$I DERIVATION_PARTIELLE' qui peut provoquer, via les 'ADD2(...)' et # $Z # les 'SOUS(...)', des sorties hors du champ calcule. Etant donne que '$xrs/project2D.11$X' # $Z # est par defaut en mode "F?_periodiser_?=VRAI", il est donc essentiel que les deux champs # $Z # champs '$xTV/SINUS.u' et '$xTV/COSINUS.u' soient reellement periodiques, ce qui empeche # $Z # donc de les calculer dans [0,pi], alors que dans [0,2.pi] ils sont bien periodiques... # $Z # Evidemment, on trouve ci-dessus : # $Z # # $Z # mv=0 Mv=1.0 # $Z # # $Z # 'v' etant elle, contrairement a 'u', bien definie dans [0,2.pi]... # $Z # # $Z # Le 'eval' a ete introduit le 20191015172410 au cas ou il '$_____Redirection' est non vide. # #20191017093007____:$Z $DELETE $xTV/FRACTAL.1 # #20191017093007____:$Z $DELETE $xTV/FRACTAL.2 # #20191017093007____:$Z $DELETE $xTV/FRACTAL.3 # #20191017093007____:$Z $DELETE $xTV/FRACTAL.4 # #20191017093007____:$Z $DELETE $xTV/FRACTAL.5 # #20191017093007____:$Z $DELETE $xTV/FRACTAL.XY # #20191017093007____:$Z $DELETE $xTV/SINUS.u # #20191017093007____:$Z $DELETE $xTV/SINUS.v # #20191017093007____:$Z $DELETE $xTV/COSINUS.u # #20191017093007____:$Z $DELETE $xTV/COSINUS.v # #20191017093007____:$Z $DELETE $xTV/FRACTAL.Z # #20191017093007____:$Z $DELETE $xTV/Fx # #20191017093007____:$Z $DELETE $xTV/Fy # #20191017093007____:$Z $DELETE $xTV/Fz #