/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N */
/* D ' U N E V A R I E T E D E H O R N E R " C I R C U L A I R E " D U P R E M I E R D E G R E : */
/* */
/* */
/* Author of '$xrs/Horner61.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20201122115610). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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_3_03(...)' 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/Horner61.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/Horner61.19.I"
/*===================================================================================================================================*/
/* :Debut_listMN_HORNER61_11: */
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N */
/* D ' U N E V A R I E T E D E H O R N E R " C I R C U L A I R E " D U P R E M I E R D E G R E : */
/* */
/* */
/* Definition de la variete ('v $xrs/Horner61.11$K') : */
/* */
/* Elle est definie parametriquement */
/* en fonction des trois parametres 'u', */
/* 'v' et 'W' : */
/* */
/* F (u,v,w) = RX.HORNER_6_01(cos(u),sin(u),cos(v),sin(v),cos(w),sin(w),{pXijklm}) */
/* x */
/* */
/* F (u,v,w) = RY.HORNER_6_01(cos(u),sin(u),cos(v),sin(v),cos(w),sin(w),{pYijklm}) */
/* y */
/* */
/* F (u,v,w) = RZ.HORNER_6_01(cos(u),sin(u),cos(v),sin(v),cos(w),sin(w),{pZijklm}) */
/* z */
/* */
/* ('v $ximd/operator.1$FON HORNER_6_01') avec : */
/* */
/* {pXijklm} un ensemble de coefficients (i,j,k,l,m E {0,1}) */
/* {pYijklm} un ensemble de coefficients (i,j,k,l,m E {0,1}) */
/* {pZijklm} un ensemble de coefficients (i,j,k,l,m E {0,1}) */
/* */
/* et : */
/* u E [ 0 , 2.p ] */
/* */
/* v E [ 0 , 2.p ] */
/* */
/* w E [ 0 , 2.p ] */
/* */
/* (ou 'p' designe 'pi'). */
/* */
/* */
/*************************************************************************************************************************************/
/* :Fin_listMN_HORNER61_11: */
#define SI_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_FLOTTANTS_EST_ACTIVABLE_ALORS_L_UTILISER_AUSSI_POUR_LES_DERIVEES
#include xrs/Horner61.12.I"
#include xrs/Horner61.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/volumes.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 xci/sequence.01.I"
/* ATTENTION, on definit ainsi le symbole 'DERNIERE_IMAGE' qui ne sert a rien ici, puisque */
/* c'est en effet 'Zmax' qui joue ce role... */
#include xrs/Horner61.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/Horner61.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 V A R I E T E D E H O R N E R " C I R C U L A I R E " D U P R E M I E R D E G R E : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/volumes.41.I"
#include xrs/volumes.22.I"
/*..............................................................................................................................*/
INITIALISATIONS_GENERALES;
/* Initialisations generales faites au tout debut... */
#include xrv/champs_5.1A.I"
#include xrs/volumes.13.I"
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_ARGUMENT_F("RX=",parametre_RX);
GET_ARGUMENT_F("pX111111=",pX_a111111);
GET_ARGUMENT_F("pX111110=",pX_a111110);
GET_ARGUMENT_F("pX111101=",pX_a111101);
GET_ARGUMENT_F("pX111100=",pX_a111100);
GET_ARGUMENT_F("pX111011=",pX_a111011);
GET_ARGUMENT_F("pX111010=",pX_a111010);
GET_ARGUMENT_F("pX111001=",pX_a111001);
GET_ARGUMENT_F("pX111000=",pX_a111000);
GET_ARGUMENT_F("pX110111=",pX_a110111);
GET_ARGUMENT_F("pX110110=",pX_a110110);
GET_ARGUMENT_F("pX110101=",pX_a110101);
GET_ARGUMENT_F("pX110100=",pX_a110100);
GET_ARGUMENT_F("pX110011=",pX_a110011);
GET_ARGUMENT_F("pX110010=",pX_a110010);
GET_ARGUMENT_F("pX110001=",pX_a110001);
GET_ARGUMENT_F("pX110000=",pX_a110000);
GET_ARGUMENT_F("pX101111=",pX_a101111);
GET_ARGUMENT_F("pX101110=",pX_a101110);
GET_ARGUMENT_F("pX101101=",pX_a101101);
GET_ARGUMENT_F("pX101100=",pX_a101100);
GET_ARGUMENT_F("pX101011=",pX_a101011);
GET_ARGUMENT_F("pX101010=",pX_a101010);
GET_ARGUMENT_F("pX101001=",pX_a101001);
GET_ARGUMENT_F("pX101000=",pX_a101000);
GET_ARGUMENT_F("pX100111=",pX_a100111);
GET_ARGUMENT_F("pX100110=",pX_a100110);
GET_ARGUMENT_F("pX100101=",pX_a100101);
GET_ARGUMENT_F("pX100100=",pX_a100100);
GET_ARGUMENT_F("pX100011=",pX_a100011);
GET_ARGUMENT_F("pX100010=",pX_a100010);
GET_ARGUMENT_F("pX100001=",pX_a100001);
GET_ARGUMENT_F("pX100000=",pX_a100000);
GET_ARGUMENT_F("pX011111=",pX_a011111);
GET_ARGUMENT_F("pX011110=",pX_a011110);
GET_ARGUMENT_F("pX011101=",pX_a011101);
GET_ARGUMENT_F("pX011100=",pX_a011100);
GET_ARGUMENT_F("pX011011=",pX_a011011);
GET_ARGUMENT_F("pX011010=",pX_a011010);
GET_ARGUMENT_F("pX011001=",pX_a011001);
GET_ARGUMENT_F("pX011000=",pX_a011000);
GET_ARGUMENT_F("pX010111=",pX_a010111);
GET_ARGUMENT_F("pX010110=",pX_a010110);
GET_ARGUMENT_F("pX010101=",pX_a010101);
GET_ARGUMENT_F("pX010100=",pX_a010100);
GET_ARGUMENT_F("pX010011=",pX_a010011);
GET_ARGUMENT_F("pX010010=",pX_a010010);
GET_ARGUMENT_F("pX010001=",pX_a010001);
GET_ARGUMENT_F("pX010000=",pX_a010000);
GET_ARGUMENT_F("pX001111=",pX_a001111);
GET_ARGUMENT_F("pX001110=",pX_a001110);
GET_ARGUMENT_F("pX001101=",pX_a001101);
GET_ARGUMENT_F("pX001100=",pX_a001100);
GET_ARGUMENT_F("pX001011=",pX_a001011);
GET_ARGUMENT_F("pX001010=",pX_a001010);
GET_ARGUMENT_F("pX001001=",pX_a001001);
GET_ARGUMENT_F("pX001000=",pX_a001000);
GET_ARGUMENT_F("pX000111=",pX_a000111);
GET_ARGUMENT_F("pX000110=",pX_a000110);
GET_ARGUMENT_F("pX000101=",pX_a000101);
GET_ARGUMENT_F("pX000100=",pX_a000100);
GET_ARGUMENT_F("pX000011=",pX_a000011);
GET_ARGUMENT_F("pX000010=",pX_a000010);
GET_ARGUMENT_F("pX000001=",pX_a000001);
GET_ARGUMENT_F("pX000000=",pX_a000000);
GET_ARGUMENT_F("RY=",parametre_RY);
GET_ARGUMENT_F("pY111111=",pY_a111111);
GET_ARGUMENT_F("pY111110=",pY_a111110);
GET_ARGUMENT_F("pY111101=",pY_a111101);
GET_ARGUMENT_F("pY111100=",pY_a111100);
GET_ARGUMENT_F("pY111011=",pY_a111011);
GET_ARGUMENT_F("pY111010=",pY_a111010);
GET_ARGUMENT_F("pY111001=",pY_a111001);
GET_ARGUMENT_F("pY111000=",pY_a111000);
GET_ARGUMENT_F("pY110111=",pY_a110111);
GET_ARGUMENT_F("pY110110=",pY_a110110);
GET_ARGUMENT_F("pY110101=",pY_a110101);
GET_ARGUMENT_F("pY110100=",pY_a110100);
GET_ARGUMENT_F("pY110011=",pY_a110011);
GET_ARGUMENT_F("pY110010=",pY_a110010);
GET_ARGUMENT_F("pY110001=",pY_a110001);
GET_ARGUMENT_F("pY110000=",pY_a110000);
GET_ARGUMENT_F("pY101111=",pY_a101111);
GET_ARGUMENT_F("pY101110=",pY_a101110);
GET_ARGUMENT_F("pY101101=",pY_a101101);
GET_ARGUMENT_F("pY101100=",pY_a101100);
GET_ARGUMENT_F("pY101011=",pY_a101011);
GET_ARGUMENT_F("pY101010=",pY_a101010);
GET_ARGUMENT_F("pY101001=",pY_a101001);
GET_ARGUMENT_F("pY101000=",pY_a101000);
GET_ARGUMENT_F("pY100111=",pY_a100111);
GET_ARGUMENT_F("pY100110=",pY_a100110);
GET_ARGUMENT_F("pY100101=",pY_a100101);
GET_ARGUMENT_F("pY100100=",pY_a100100);
GET_ARGUMENT_F("pY100011=",pY_a100011);
GET_ARGUMENT_F("pY100010=",pY_a100010);
GET_ARGUMENT_F("pY100001=",pY_a100001);
GET_ARGUMENT_F("pY100000=",pY_a100000);
GET_ARGUMENT_F("pY011111=",pY_a011111);
GET_ARGUMENT_F("pY011110=",pY_a011110);
GET_ARGUMENT_F("pY011101=",pY_a011101);
GET_ARGUMENT_F("pY011100=",pY_a011100);
GET_ARGUMENT_F("pY011011=",pY_a011011);
GET_ARGUMENT_F("pY011010=",pY_a011010);
GET_ARGUMENT_F("pY011001=",pY_a011001);
GET_ARGUMENT_F("pY011000=",pY_a011000);
GET_ARGUMENT_F("pY010111=",pY_a010111);
GET_ARGUMENT_F("pY010110=",pY_a010110);
GET_ARGUMENT_F("pY010101=",pY_a010101);
GET_ARGUMENT_F("pY010100=",pY_a010100);
GET_ARGUMENT_F("pY010011=",pY_a010011);
GET_ARGUMENT_F("pY010010=",pY_a010010);
GET_ARGUMENT_F("pY010001=",pY_a010001);
GET_ARGUMENT_F("pY010000=",pY_a010000);
GET_ARGUMENT_F("pY001111=",pY_a001111);
GET_ARGUMENT_F("pY001110=",pY_a001110);
GET_ARGUMENT_F("pY001101=",pY_a001101);
GET_ARGUMENT_F("pY001100=",pY_a001100);
GET_ARGUMENT_F("pY001011=",pY_a001011);
GET_ARGUMENT_F("pY001010=",pY_a001010);
GET_ARGUMENT_F("pY001001=",pY_a001001);
GET_ARGUMENT_F("pY001000=",pY_a001000);
GET_ARGUMENT_F("pY000111=",pY_a000111);
GET_ARGUMENT_F("pY000110=",pY_a000110);
GET_ARGUMENT_F("pY000101=",pY_a000101);
GET_ARGUMENT_F("pY000100=",pY_a000100);
GET_ARGUMENT_F("pY000011=",pY_a000011);
GET_ARGUMENT_F("pY000010=",pY_a000010);
GET_ARGUMENT_F("pY000001=",pY_a000001);
GET_ARGUMENT_F("pY000000=",pY_a000000);
GET_ARGUMENT_F("RZ=",parametre_RZ);
GET_ARGUMENT_F("pZ111111=",pZ_a111111);
GET_ARGUMENT_F("pZ111110=",pZ_a111110);
GET_ARGUMENT_F("pZ111101=",pZ_a111101);
GET_ARGUMENT_F("pZ111100=",pZ_a111100);
GET_ARGUMENT_F("pZ111011=",pZ_a111011);
GET_ARGUMENT_F("pZ111010=",pZ_a111010);
GET_ARGUMENT_F("pZ111001=",pZ_a111001);
GET_ARGUMENT_F("pZ111000=",pZ_a111000);
GET_ARGUMENT_F("pZ110111=",pZ_a110111);
GET_ARGUMENT_F("pZ110110=",pZ_a110110);
GET_ARGUMENT_F("pZ110101=",pZ_a110101);
GET_ARGUMENT_F("pZ110100=",pZ_a110100);
GET_ARGUMENT_F("pZ110011=",pZ_a110011);
GET_ARGUMENT_F("pZ110010=",pZ_a110010);
GET_ARGUMENT_F("pZ110001=",pZ_a110001);
GET_ARGUMENT_F("pZ110000=",pZ_a110000);
GET_ARGUMENT_F("pZ101111=",pZ_a101111);
GET_ARGUMENT_F("pZ101110=",pZ_a101110);
GET_ARGUMENT_F("pZ101101=",pZ_a101101);
GET_ARGUMENT_F("pZ101100=",pZ_a101100);
GET_ARGUMENT_F("pZ101011=",pZ_a101011);
GET_ARGUMENT_F("pZ101010=",pZ_a101010);
GET_ARGUMENT_F("pZ101001=",pZ_a101001);
GET_ARGUMENT_F("pZ101000=",pZ_a101000);
GET_ARGUMENT_F("pZ100111=",pZ_a100111);
GET_ARGUMENT_F("pZ100110=",pZ_a100110);
GET_ARGUMENT_F("pZ100101=",pZ_a100101);
GET_ARGUMENT_F("pZ100100=",pZ_a100100);
GET_ARGUMENT_F("pZ100011=",pZ_a100011);
GET_ARGUMENT_F("pZ100010=",pZ_a100010);
GET_ARGUMENT_F("pZ100001=",pZ_a100001);
GET_ARGUMENT_F("pZ100000=",pZ_a100000);
GET_ARGUMENT_F("pZ011111=",pZ_a011111);
GET_ARGUMENT_F("pZ011110=",pZ_a011110);
GET_ARGUMENT_F("pZ011101=",pZ_a011101);
GET_ARGUMENT_F("pZ011100=",pZ_a011100);
GET_ARGUMENT_F("pZ011011=",pZ_a011011);
GET_ARGUMENT_F("pZ011010=",pZ_a011010);
GET_ARGUMENT_F("pZ011001=",pZ_a011001);
GET_ARGUMENT_F("pZ011000=",pZ_a011000);
GET_ARGUMENT_F("pZ010111=",pZ_a010111);
GET_ARGUMENT_F("pZ010110=",pZ_a010110);
GET_ARGUMENT_F("pZ010101=",pZ_a010101);
GET_ARGUMENT_F("pZ010100=",pZ_a010100);
GET_ARGUMENT_F("pZ010011=",pZ_a010011);
GET_ARGUMENT_F("pZ010010=",pZ_a010010);
GET_ARGUMENT_F("pZ010001=",pZ_a010001);
GET_ARGUMENT_F("pZ010000=",pZ_a010000);
GET_ARGUMENT_F("pZ001111=",pZ_a001111);
GET_ARGUMENT_F("pZ001110=",pZ_a001110);
GET_ARGUMENT_F("pZ001101=",pZ_a001101);
GET_ARGUMENT_F("pZ001100=",pZ_a001100);
GET_ARGUMENT_F("pZ001011=",pZ_a001011);
GET_ARGUMENT_F("pZ001010=",pZ_a001010);
GET_ARGUMENT_F("pZ001001=",pZ_a001001);
GET_ARGUMENT_F("pZ001000=",pZ_a001000);
GET_ARGUMENT_F("pZ000111=",pZ_a000111);
GET_ARGUMENT_F("pZ000110=",pZ_a000110);
GET_ARGUMENT_F("pZ000101=",pZ_a000101);
GET_ARGUMENT_F("pZ000100=",pZ_a000100);
GET_ARGUMENT_F("pZ000011=",pZ_a000011);
GET_ARGUMENT_F("pZ000010=",pZ_a000010);
GET_ARGUMENT_F("pZ000001=",pZ_a000001);
GET_ARGUMENT_F("pZ000000=",pZ_a000000);
PROCESS_ARGUMENTS_DE_DERIVATION_PARTIELLE_NUMERIQUE;
PROCESS_ARGUMENTS_DE_COLORIAGE;
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v_w;
/* Introduit le 20221221185313... */
PROCESS_ARGUMENTS_DE_PROJECTION_DU_VOLUME;
PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;
)
);
#include xrs/volumes.12.I"
/* Visualisation et evolution de la variete au cours du temps... */
VISUALISATION_DU_VOLUME_AVEC_DIFFERENCES____(BLOC(BLOC(VIDE;)
)
,Fx_Horner61_1,Fy_Horner61_1,Fz_Horner61_1
/* Fonctions definissant parametriquement la variete... */
,Pxyz_Horner61_1
,BLOC(VIDE;)
/* Il n'y a rien de specifique a calculer pour le point courant... */
);
RETU_Commande;
Eblock
ECommande
Copyright © Jean-François COLONNA, 2020-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2020-2024.