/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N C U B E A R R O N D I : */
/* */
/* */
/* Author of '$xrs/cube.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 19990519113903). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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_CL_____MODULE_NON_OPTIMISABLE
/* Le 'GooF' fut introduit le 20080126104816, mais le 20080126173152 j'ai du y renoncer, */
/* les fichiers generes etant monstrueux (plus de 10 Go). Il fut retabli le 20080127093210 */
/* grace a 'v $xrs/surfaces.12$I 20080127084712'. Cela fut de nouveau supprime le */
/* 20080127094147 car cela ne suffisait pas a faire disparaitre les monstruosite... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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
/* Introduit le 20041229180953 pour 'v $xrs/surfaces.12$I IFmageA4'... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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"
/* Introduit le 20130306152247... */
#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/sphere.18.I"
/* Introduit le 20030130120355 a cause de 'v $xcc/cpp$Z INCLUDES_1.sed'. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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/sphere.19.I"
/* Introduit le 20030130120355 a cause de 'v $xcc/cpp$Z INCLUDES_1.sed'. */
/*===================================================================================================================================*/
/* :Debut_listMN_CUBE_11: */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N C U B E A R R O N D I : */
/* */
/* */
/* Definition de la surface ('v $xrs/cube.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) = R.Dsin(u).Dcos(v) */
/* x */
/* */
/* F (u,v) = R.Dsin(u).Dsin(v) */
/* y */
/* */
/* F (u,v) = R.Dcos(u) */
/* z */
/* */
/* avec : */
/* */
/* sin(1.t) sin(3.t) sin(5.t) sin(7.t) sin(9.t) sin(11.t) */
/* Dsin(t) = K.[---------- - ---------- + ---------- - ---------- + ---------- - -----------] */
/* 2 2 2 2 2 2 */
/* 1 3 5 7 9 11 */
/* */
/* Dcos(t) = Dsin(t + (p/2)) */
/* */
/* avec : */
/* */
/* u E [ 0 , p ] */
/* */
/* v E [ 0 , 2.p ] */
/* */
/* (ou 'p' designe 'pi'). */
/* */
/* */
/*************************************************************************************************************************************/
/* :Fin_listMN_CUBE_11: */
#define PAS_D_APPROXIMATION_DE_LA_FONCTION_EN_DENT_DE_SCIE \
FRA10(FRA10(FU)) \
/* Afin d'approximer la fonction en dent de scie (introduit le 20190201100404...). */
#include xrs/sphere.12.I"
#include xrs/sphere.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"
#include xrs/sphere.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/cube.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 C U B E A R R O N D I : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/surfaces.41.I"
/* Introduit le 20150528135618... */
#include xrs/surfaces.22.I"
DEFV(Float,INIT(pas_d_approximation_de_la_fonction_en_dent_de_scie,PAS_D_APPROXIMATION_DE_LA_FONCTION_EN_DENT_DE_SCIE));
/* Afin d'approximer la fonction en dent de scie (introduit le 20190201100404...). */
DEFV(Float,INIT(d_coefficient_de_renormalisation,d_FCONSTANTES));
DEFV(Float,INIT(coefficient_de_renormalisation,FLOT__UNDEF));
/* Coefficient intervenant dans la definition des trois fonctions {Fx,Fy,Fz} (introduit le */
/* 20060526175304). */
/*..............................................................................................................................*/
INITIALISATIONS_GENERALES;
/* Initialisations generales faites au tout debut... */
iTRANSFORMAT_31(liste_PARAMETRE_R,PARAMETRE_R_IMPLICITE);
/* Initialisation du parametre 'r'. */
#include xrv/champs_5.1A.I"
#include xrs/surfaces.13.I"
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_FICHIER_PARAMETRE("PARAMETRE_R=",_R);
GET_ARGUMENT_F("pas_dent_de_scie=""pas=",pas_d_approximation_de_la_fonction_en_dent_de_scie);
/* Introduit le 20190201100404. */
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;
/* Introduit le 20060117145521. */
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
/* Introduit le 20050615091657. */
PROCESS_ARGUMENTS_DE_PROJECTION_DE_LA_SURFACE;
/* Introduit le 20041126095151. */
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;
/* Introduit le 20060831113514... */
PROCESS_ARGUMENTS_DE_GESTION_DE_LA_RELATIVITE_RESTREINTE;
)
);
#include xrs/surfaces.12.I"
/* Visualisation et evolution de la surface au cours du temps... */
begin_nouveau_block
Bblock
DEFV(Float,INIT(angle_courant,FLOT__UNDEF));
DEFV(Float,INIT(minimum_de_la_fonction_en_dent_scie,F_INFINI));
DEFV(Float,INIT(maximum_de_la_fonction_en_dent_scie,F_MOINS_L_INFINI));
/* Pour le calcul des extrema d'une approximation de la fonction "en dent de scie"... */
DoIn(angle_courant,FZERO,CERCLE_TRIGONOMETRIQUE,pas_d_approximation_de_la_fonction_en_dent_de_scie)
Bblock
DEFV(Float,INIT(fonction_en_dent_scie,FZERO));
INCR(fonction_en_dent_scie,DSINX_1_1(angle_courant));
INCR(fonction_en_dent_scie,DSINX_1_3(angle_courant));
INCR(fonction_en_dent_scie,DSINX_1_5(angle_courant));
INCR(fonction_en_dent_scie,DSINX_1_7(angle_courant));
INCR(fonction_en_dent_scie,DSINX_1_9(angle_courant));
INCR(fonction_en_dent_scie,DSINX_1_B(angle_courant));
/* Calcul d'une approximation de la fonction "en dent de scie"... */
EGAL(minimum_de_la_fonction_en_dent_scie,MIN2(minimum_de_la_fonction_en_dent_scie,fonction_en_dent_scie));
EGAL(maximum_de_la_fonction_en_dent_scie,MAX2(maximum_de_la_fonction_en_dent_scie,fonction_en_dent_scie));
EGAL(coefficient_de_renormalisation
,MAX2(ABSO(minimum_de_la_fonction_en_dent_scie),ABSO(maximum_de_la_fonction_en_dent_scie))
);
/* Calcul du coefficient de renormalisation. A la date du 20060526183806, il vaut */
/* 1.1921275137943794... */
Eblock
EDoI
Eblock
end_nouveau_block
VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENTIATION(BLOC(Bblock
EGAL(parametre_r,sPARAMETRE_R(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'r'. */
Eblock
)
,Fx_Cube_1,Fy_Cube_1,Fz_Cube_1
/* Fonctions definissant parametriquement la surface... */
,Pxyz_Cube_1
,BLOC(VIDE;)
/* Il n'y a rien de specifique a calculer pour le point courant... */
);
RETU_Commande;
Eblock
ECommande