/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N */
/* D ' U N E S U R F A C E D E H O R N E R " L I N E A I R E " D U Q U A T R I E M E D E G R E : */
/* */
/* */
/* Author of '$xrs/Horner24.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20201120101627). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* I N T E R F A C E ' listG ' : */
/* */
/* */
/* :Debut_listG: */
/* :Fin_listG: */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D I R E C T I V E S S P E C I F I Q U E S D E C O M P I L A T I O N : */
/* */
/*************************************************************************************************************************************/
@define PRAGMA_CPP_____MODULE_NON_DERIVABLE
#nodefine PRAGMA_CPP_____AUTORISER_LE_GooF
#nodefine PRAGMA_CPP_____SI_LE_GooF_EST_ACTIVABLE_ALORS_FAIRE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_FLOTTANTS
/* L'arithmetique etendue a peu d'interet avec la procedure 'HORNER_2_04(...)' car, en */
/* effet, elle ne contient que des 'ADD2(...)'s et des 'MUL2(...)'s ; changer leur */
/* definition ne provoquerait donc pas des petites modifications subtiles..., tout en etant */
/* interminable au niveau pre-processing avec 'GooF'... */
@define PRAGMA_CL_____MODULE_NON_OPTIMISABLE
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* F I C H I E R S D ' I N C L U D E S : */
/* */
/*************************************************************************************************************************************/
#include INCLUDES_BASE
#include image_image_IMAGESF_EXT
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N S D E B A S E E T U N I V E R S E L L E S : */
/* */
/*************************************************************************************************************************************/
#include xrv/dimensionnement.01.vv.I"
#include xrk/attractor.11.I"
#include xrv/particule.41.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 3 */
/* D E F I N I T I O N D E L ' E S P A C E P H Y S I Q U E D A N S R ( D E B U T ) : */
/* */
/* */
/* Nota : */
/* */
/* Les extrema des coordonnees {x,y,z} */
/* ainsi que ceux de leurs differentielles */
/* {dx,dy,dz} sont fixees un peu arbitrairement */
/* et sans etre parametrees. */
/* */
/* */
/*************************************************************************************************************************************/
#include xrs/Horner24.18.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 3 */
/* D E F I N I T I O N D E L ' E S P A C E P H Y S I Q U E D A N S R ( D E B U T ) : */
/* */
/*************************************************************************************************************************************/
#include xrk/attractor.12.I"
#include xrs/Horner24.19.I"
/*===================================================================================================================================*/
/* :Debut_listMN_HORNER24_11: */
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N */
/* D ' U N E S U R F A C E D E H O R N E R " L I N E A I R E " D U Q U A T R I E M E D E G R E : */
/* */
/* */
/* Definition de la surface ('v $xrs/Horner24.11$K') : */
/* */
/* Elle 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) = RX.HORNER_2_04(u,v,{pXij}) */
/* x */
/* */
/* F (u,v) = RY.HORNER_2_04(u,v,{pYij}) */
/* y */
/* */
/* F (u,v) = RZ.HORNER_2_04(u,v,{pZij}) */
/* z */
/* */
/* ('v $ximd/operator.1$FON HORNER_2_04') avec : */
/* */
/* {pXij} un ensemble de coefficients (i,j E {0,1,2,3,4}) */
/* {pYij} un ensemble de coefficients (i,j E {0,1,2,3,4}) */
/* {pZij} un ensemble de coefficients (i,j E {0,1,2,3,4}) */
/* */
/* et : */
/* */
/* avec : */
/* */
/* u E [ -1 , +1 ] */
/* */
/* v E [ -1 , +1 ] */
/* */
/* */
/*************************************************************************************************************************************/
/* :Fin_listMN_HORNER24_11: */
#define SI_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_FLOTTANTS_EST_ACTIVABLE_ALORS_L_UTILISER_AUSSI_POUR_LES_DERIVEES
#include xrs/Horner24.12.I"
#include xrs/Horner24.14.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* F O N C T I O N S D E V I S U A L I S A T I O N E T D ' I N T E R P O L A T I O N : */
/* */
/*************************************************************************************************************************************/
#include xrs/surfaces.21.I"
#define EDITER_LES_EXTREMA_S_ILS_SONT_HORS_PREVISIONS \
FAUX \
/* En consequence de quoi, le message d'erreur relatif aux extrema de {x,y,z} et/ou de */ \
/* {dx,dy,dz} n'a plus de sens. */
#include xrs/Horner24.13.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 ' : */
/* */
/*************************************************************************************************************************************/
#include xrs/Horner24.11.I"
/* 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 'RAYON(h,theta)' 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... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* P O U R D E S R A I S O N S D E C O M P A T I B I L I T E : */
/* */
/*************************************************************************************************************************************/
#include xrk/integr.1B.vv.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S I N I T I A L I S A T I O N S : */
/* */
/*************************************************************************************************************************************/
#include xrk/attractor.18.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N */
/* D ' U N E S U R F A C E D E H O R N E R " L I N E A I R E " D U Q U A T R I E M E D E G R E : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/surfaces.41.I"
#include xrs/surfaces.22.I"
/*..............................................................................................................................*/
INITIALISATIONS_GENERALES;
/* Initialisations generales faites au tout debut... */
#include xrv/champs_5.1A.I"
#include xrs/surfaces.13.I"
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_ARGUMENT_F("RX=",parametre_RX);
GET_ARGUMENT_F("pX44=",pX_a44);
GET_ARGUMENT_F("pX43=",pX_a43);
GET_ARGUMENT_F("pX42=",pX_a42);
GET_ARGUMENT_F("pX41=",pX_a41);
GET_ARGUMENT_F("pX40=",pX_a40);
GET_ARGUMENT_F("pX34=",pX_a34);
GET_ARGUMENT_F("pX33=",pX_a33);
GET_ARGUMENT_F("pX32=",pX_a32);
GET_ARGUMENT_F("pX31=",pX_a31);
GET_ARGUMENT_F("pX30=",pX_a30);
GET_ARGUMENT_F("pX24=",pX_a24);
GET_ARGUMENT_F("pX23=",pX_a23);
GET_ARGUMENT_F("pX22=",pX_a22);
GET_ARGUMENT_F("pX21=",pX_a21);
GET_ARGUMENT_F("pX20=",pX_a20);
GET_ARGUMENT_F("pX14=",pX_a14);
GET_ARGUMENT_F("pX13=",pX_a13);
GET_ARGUMENT_F("pX12=",pX_a12);
GET_ARGUMENT_F("pX11=",pX_a11);
GET_ARGUMENT_F("pX10=",pX_a10);
GET_ARGUMENT_F("pX04=",pX_a04);
GET_ARGUMENT_F("pX03=",pX_a03);
GET_ARGUMENT_F("pX02=",pX_a02);
GET_ARGUMENT_F("pX01=",pX_a01);
GET_ARGUMENT_F("pX00=",pX_a00);
GET_ARGUMENT_F("RY=",parametre_RY);
GET_ARGUMENT_F("pY44=",pY_a44);
GET_ARGUMENT_F("pY43=",pY_a43);
GET_ARGUMENT_F("pY42=",pY_a42);
GET_ARGUMENT_F("pY41=",pY_a41);
GET_ARGUMENT_F("pY40=",pY_a40);
GET_ARGUMENT_F("pY34=",pY_a34);
GET_ARGUMENT_F("pY33=",pY_a33);
GET_ARGUMENT_F("pY32=",pY_a32);
GET_ARGUMENT_F("pY31=",pY_a31);
GET_ARGUMENT_F("pY30=",pY_a30);
GET_ARGUMENT_F("pY24=",pY_a24);
GET_ARGUMENT_F("pY23=",pY_a23);
GET_ARGUMENT_F("pY22=",pY_a22);
GET_ARGUMENT_F("pY21=",pY_a21);
GET_ARGUMENT_F("pY20=",pY_a20);
GET_ARGUMENT_F("pY14=",pY_a14);
GET_ARGUMENT_F("pY13=",pY_a13);
GET_ARGUMENT_F("pY12=",pY_a12);
GET_ARGUMENT_F("pY11=",pY_a11);
GET_ARGUMENT_F("pY10=",pY_a10);
GET_ARGUMENT_F("pY04=",pY_a04);
GET_ARGUMENT_F("pY03=",pY_a03);
GET_ARGUMENT_F("pY02=",pY_a02);
GET_ARGUMENT_F("pY01=",pY_a01);
GET_ARGUMENT_F("pY00=",pY_a00);
GET_ARGUMENT_F("RZ=",parametre_RZ);
GET_ARGUMENT_F("pZ44=",pZ_a44);
GET_ARGUMENT_F("pZ43=",pZ_a43);
GET_ARGUMENT_F("pZ42=",pZ_a42);
GET_ARGUMENT_F("pZ41=",pZ_a41);
GET_ARGUMENT_F("pZ40=",pZ_a40);
GET_ARGUMENT_F("pZ34=",pZ_a34);
GET_ARGUMENT_F("pZ33=",pZ_a33);
GET_ARGUMENT_F("pZ32=",pZ_a32);
GET_ARGUMENT_F("pZ31=",pZ_a31);
GET_ARGUMENT_F("pZ30=",pZ_a30);
GET_ARGUMENT_F("pZ24=",pZ_a24);
GET_ARGUMENT_F("pZ23=",pZ_a23);
GET_ARGUMENT_F("pZ22=",pZ_a22);
GET_ARGUMENT_F("pZ21=",pZ_a21);
GET_ARGUMENT_F("pZ20=",pZ_a20);
GET_ARGUMENT_F("pZ14=",pZ_a14);
GET_ARGUMENT_F("pZ13=",pZ_a13);
GET_ARGUMENT_F("pZ12=",pZ_a12);
GET_ARGUMENT_F("pZ11=",pZ_a11);
GET_ARGUMENT_F("pZ10=",pZ_a10);
GET_ARGUMENT_F("pZ04=",pZ_a04);
GET_ARGUMENT_F("pZ03=",pZ_a03);
GET_ARGUMENT_F("pZ02=",pZ_a02);
GET_ARGUMENT_F("pZ01=",pZ_a01);
GET_ARGUMENT_F("pZ00=",pZ_a00);
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
PROCESS_ARGUMENTS_DE_PROJECTION_DE_LA_SURFACE;
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_TEXTURE;
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_NORMALE;
PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;
PROCESS_ARGUMENTS_DE_DEFINITION_D_UNE_DROITE_DU_PLAN_UV;
PROCESS_ARGUMENTS_D_ACTIVATION_DE_L_ARITHMETIQUE_ETENDUE_ET_NON_BASIQUE_POUR_Fx_Fy_Fz;
PROCESS_ARGUMENTS_DE_GESTION_DE_LA_RELATIVITE_RESTREINTE;
)
);
#include xrs/surfaces.12.I"
/* Visualisation et evolution de la surface au cours du temps... */
VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENCES____(BLOC(BLOC(VIDE;)
)
,Fx_Horner24_1,Fy_Horner24_1,Fz_Horner24_1
/* Fonctions definissant parametriquement la surface... */
,Pxyz_Horner24_1
,BLOC(VIDE;)
/* Il n'y a rien de specifique a calculer pour le point courant... */
);
RETU_Commande;
Eblock
ECommande