/*************************************************************************************************************************************/
/* */
/* 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 V A R I E T E D E H O R N E R " L I N E A I R E " D U T R O I S I E M E D E G R E : */
/* */
/* */
/* Author of '$xrs/Horner33.11$I' : */
/* */
/* Jean-Francois Colonna (LACTAMME, 20201121120853). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 de la variete ('v $xrs/Horner33.11$K') : */
/* */
/* La variete 'lineaire' de Horner */
/* 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,w) = RX.HORNER_3_03(u,v,w,{pXijk}) */
/* x */
/* */
/* F (u,v,w) = RY.HORNER_3_03(u,v,w,{pYijk}) */
/* y */
/* */
/* F (u,v,w) = RZ.HORNER_3_03(u,v,w,{pZijk}) */
/* z */
/* */
/* ('v $ximd/operator.1$FON HORNER_3_03') avec : */
/* */
/* {pXijk} un ensemble de coefficients (i,j,k E {0,1,2,3}) */
/* {pYijk} un ensemble de coefficients (i,j,k E {0,1,2,3}) */
/* {pZijk} un ensemble de coefficients (i,j,k E {0,1,2,3}) */
/* */
/* et : */
/* */
/* avec : */
/* */
/* u E [ -1 , +1 ] */
/* */
/* v E [ -1 , +1 ] */
/* */
/* w E [ -1 , +1 ] */
/* */
/* */
/*************************************************************************************************************************************/
/* 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,w)' 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_Horner33_1(u,v,w) \
MUL2(parametre_RX \
,HORNER_3_03(u,v,w \
,pX_a333,pX_a332,pX_a331,pX_a330,pX_a323,pX_a322,pX_a321,pX_a320 \
,pX_a313,pX_a312,pX_a311,pX_a310,pX_a303,pX_a302,pX_a301,pX_a300 \
,pX_a233,pX_a232,pX_a231,pX_a230,pX_a223,pX_a222,pX_a221,pX_a220 \
,pX_a213,pX_a212,pX_a211,pX_a210,pX_a203,pX_a202,pX_a201,pX_a200 \
,pX_a133,pX_a132,pX_a131,pX_a130,pX_a123,pX_a122,pX_a121,pX_a120 \
,pX_a113,pX_a112,pX_a111,pX_a110,pX_a103,pX_a102,pX_a101,pX_a100 \
,pX_a033,pX_a032,pX_a031,pX_a030,pX_a023,pX_a022,pX_a021,pX_a020 \
,pX_a013,pX_a012,pX_a011,pX_a010,pX_a003,pX_a002,pX_a001,pX_a000 \
) \
)
/* Definition de la fonction F (u,v,w). */
/* x */
#define Fy_Horner33_1(u,v,w) \
MUL2(parametre_RY \
,HORNER_3_03(u,v,w \
,pY_a333,pY_a332,pY_a331,pY_a330,pY_a323,pY_a322,pY_a321,pY_a320 \
,pY_a313,pY_a312,pY_a311,pY_a310,pY_a303,pY_a302,pY_a301,pY_a300 \
,pY_a233,pY_a232,pY_a231,pY_a230,pY_a223,pY_a222,pY_a221,pY_a220 \
,pY_a213,pY_a212,pY_a211,pY_a210,pY_a203,pY_a202,pY_a201,pY_a200 \
,pY_a133,pY_a132,pY_a131,pY_a130,pY_a123,pY_a122,pY_a121,pY_a120 \
,pY_a113,pY_a112,pY_a111,pY_a110,pY_a103,pY_a102,pY_a101,pY_a100 \
,pY_a033,pY_a032,pY_a031,pY_a030,pY_a023,pY_a022,pY_a021,pY_a020 \
,pY_a013,pY_a012,pY_a011,pY_a010,pY_a003,pY_a002,pY_a001,pY_a000 \
) \
)
/* Definition de la fonction F (u,v,w). */
/* y */
#define Fz_Horner33_1(u,v,w) \
MUL2(parametre_RZ \
,HORNER_3_03(u,v,w \
,pZ_a333,pZ_a332,pZ_a331,pZ_a330,pZ_a323,pZ_a322,pZ_a321,pZ_a320 \
,pZ_a313,pZ_a312,pZ_a311,pZ_a310,pZ_a303,pZ_a302,pZ_a301,pZ_a300 \
,pZ_a233,pZ_a232,pZ_a231,pZ_a230,pZ_a223,pZ_a222,pZ_a221,pZ_a220 \
,pZ_a213,pZ_a212,pZ_a211,pZ_a210,pZ_a203,pZ_a202,pZ_a201,pZ_a200 \
,pZ_a133,pZ_a132,pZ_a131,pZ_a130,pZ_a123,pZ_a122,pZ_a121,pZ_a120 \
,pZ_a113,pZ_a112,pZ_a111,pZ_a110,pZ_a103,pZ_a102,pZ_a101,pZ_a100 \
,pZ_a033,pZ_a032,pZ_a031,pZ_a030,pZ_a023,pZ_a022,pZ_a021,pZ_a020 \
,pZ_a013,pZ_a012,pZ_a011,pZ_a010,pZ_a003,pZ_a002,pZ_a001,pZ_a000 \
) \
)
/* Definition de la fonction F (u,v,w). */
/* 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_Horner33_1 \
Bblock \
BLOC(VIDE;); \
Eblock
/* Initialisations specifiques a cette variete destinees a permettre la reinjection des */
/* trois pseudo-projections {Projection_de_Fx,Projection_de_Fy,Projection_de_Fz} dans */
/* 'v $xrs/project2D.11$K'... */