/*************************************************************************************************************************************/ /* */ /* D E F I N I T I O N S R E L A T I V E S A L A D E R I V A T I O N F O R M E L L E */ /* D ' U N E L L I P S O I D E : */ /* */ /* */ /* Author of '$xrs/ellipsoid.11$I' : */ /* */ /* Jean-Francois Colonna (LACTAMME, 1993??????????). */ /* */ /*************************************************************************************************************************************/ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* D E F I N I T I O N D E S T R O I S F O N C T I O N S ' F ' : */ /* */ /* */ /* Definition ('v $xrs/ellipsoid.11$K') : */ /* */ /* L'ellipsoide est definie parametriquement */ /* en fonction des deux parametres 'u' */ /* (appele aussi 'distance polaire' ou */ /* 'theta' ou encore 'latitude') et 'v' */ /* (appele aussi 'longitude' ou 'phi') : */ /* */ /* F (u,v) = a.sin(u).cos(v) */ /* x */ /* */ /* F (u,v) = b.sin(u).sin(v) */ /* y */ /* */ /* F (u,v) = c.cos(u) */ /* z */ /* */ /* avec (voir l'Aide-Memoire de Mathematiques, */ /* pages 321 et 322) : */ /* */ /* u E [ 0 , p ] */ /* */ /* v E [ 0 , 2.p ] */ /* */ /* (ou 'p' designe 'pi'). */ /* */ /* */ /*************************************************************************************************************************************/ /* ATTENTION, il est imperatif que les fonctions derivables formellement soient definies */ /* dans un fichier a part afin qu'elles ne soient pas l'objet d'une double definition. En */ /* effet, on trouve dans 'v $xcc/cpp$Z' : */ /* */ /* $CA $module$w */ /* | $xcp/cpp$X ... -c$PASSE_D -e$PASSE_5 */ /* > $fichierR */ /* */ /* Ainsi, si par exemple la definition de 'Fx(u,v)' etait faite localement dans le */ /* fichier '$xrf/recursif.11$K', cette definition apparaitrait deux fois dans le fichier */ /* '$module$w' : une premiere fois en tete, provenant de '$PASSE_D' ou l'on ne trouve que */ /* les definitions du type '-define ...', puis une deuxieme fois a la suite, la ou se trouve */ /* concatene au fichier '$module$W' des '-define ...' le fichier '$fichierA' a compiler... */ -define Fx_Ellipsoide_1(u,v) \ -_-_-_- ADD2(Xcartesienne_3D(parametre_a \ -_-_-_- ,v \ -_-_-_- ,u \ -_-_-_- ) \ -_-_-_- ,FZERO \ -_-_-_- ) /* Definition de la fonction F (u,v). */ /* x */ -define Fy_Ellipsoide_1(u,v) \ -_-_-_- ADD2(Ycartesienne_3D(parametre_b \ -_-_-_- ,v \ -_-_-_- ,u \ -_-_-_- ) \ -_-_-_- ,FZERO \ -_-_-_- ) /* Definition de la fonction F (u,v). */ /* y */ -define Fz_Ellipsoide_1(u,v) \ -_-_-_- ADD2(Zcartesienne_3D(parametre_c \ -_-_-_- ,v \ -_-_-_- ,u \ -_-_-_- ) \ -_-_-_- ,FZERO \ -_-_-_- ) /* Definition de la fonction F (u,v). */ /* z */ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* I N I T I A L I S A T I O N S R E L A T I V E S A L A P S E U D O - P R O J E C T I O N : */ /* */ /*************************************************************************************************************************************/ -define Pxyz_Ellipsoide_1 \ -_-_-_- Bblock \ -_-_-_- EGAL(minimum_de_u,FZERO); \ -_-_-_- EGAL(maximum_de_u,CERCLE_TRIGONOMETRIQUE); \ -_-_-_- EGAL(minimum_de_v,FZERO); \ -_-_-_- EGAL(maximum_de_v,CERCLE_TRIGONOMETRIQUE); \ -_-_-_- /* Afin d'assurer la periodicite du champ... */ \ -_-_-_- iTRANSFORMAT_31(liste_PARAMETRE_A,FU); \ -_-_-_- iTRANSFORMAT_31(liste_PARAMETRE_B,FU); \ -_-_-_- iTRANSFORMAT_31(liste_PARAMETRE_C,FU); \ -_-_-_- /* Afin de definir un ellipsoide "unite"... */ \ -_-_-_- Eblock /* Initialisations specifiques a cette surface destinees a permettre la reinjection des */ /* trois pseudo-projections {Projection_de_Fx,Projection_de_Fy,Projection_de_Fz} dans */ /* 'v $xrs/project2D.11$K' (introduit le 20050203093002)... */