/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N P S E U D O - T O R E A D E U X T R O U S : */
/* */
/* */
/* Nota : */
/* */
/* Ce programme s'inspire du programme */
/* 'v $xtc/huit.03$c'... */
/* */
/* */
/* Author of '$xrs/huit.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20161202105922). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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_____AUTORISER_LE_GooF
#define PRAGMA_CPP_____SI_LE_GooF_EST_ACTIVABLE_ALORS_FAIRE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_FLOTTANTS
@define PRAGMA_CPP_____EXECUTER_UNE_COMMANDE_AU_PREALABLE "$xi/DerivFormel$Z"
/* Introduit le 20161209153632 ('PRAGMA_CPP_____EXECUTER__xi_DerivFormel_Z') et mis sous */
/* une forme plus generale le 20161211093510... */
@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/huit.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/huit.19.I"
/*===================================================================================================================================*/
/* :Debut_listMN_HUIT_11: */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 surface ('v $xrs/huit.11$K') : */
/* */
/* Le pseudo-tore a deux trous est defini parametriquement */
/* en fonction des deux parametres 'u' */
/* (appele aussi 'distance polaire' ou */
/* 'theta' ou encore 'latitude') et 'v' */
/* (appele aussi 'longitude' ou 'phi') : */
/* */
/* Choix(u,v) = +1 si v < M.(4.pi/D) */
/* = -1 sinon */
/* */
/* __________ */
/* / A+cos(v) */
/* XPlan(u,v) = R.Choix(u,v).\/ ---------- */
/* B */
/* */
/* sin(v) */
/* YPlan(u,v) = R.-------- */
/* C */
/* Definition d'une courbe plane en forme de "8". */
/* */
/* dXPlan(u,v) = d#XPlan(u,v) */
/* dYPlan(u,v) = d#YPlan(u,v) */
/* La tangente T est donc {+dXPlan(u,v),+dYPlan(u,v). */
/* */
/* Psi(u,v) = arctg(-dXPlan(u,v),+dYPlan(u,v)) */
/* La normale N est orthogonale a la tangente T */
/* et est donc {+dYPlan(u,v),-dXPlan(u,v)}. */
/* */
/* Xellipse(u,v) = Ra.cos(u) */
/* Yellipse(u,v) = 0 */
/* Zellipse(u,v) = Rb.sin(u) */
/* Definition d'une ellipse dans le plan {OX,OZ}. */
/* */
/* F (u,v) = [Xellipse(u,v)*cos(psi)] - [Yellipse(u,v)*sin(psi)] + XPlan(u,v) */
/* x */
/* */
/* F (u,v) = [Xellipse(u,v)*sin(psi)] + [Yellipse(u,v)*cos(psi)] + YPlan(u,v) */
/* y */
/* */
/* F (u,v) = Zellipse(u,v) */
/* z */
/* Definition d'une ellipse dont le centre est le point */
/* courant {XPlan(u,v),YPlan(u,v),0} et situee dans le */
/* plan {N,OZ}. */
/* */
/* et (parametres par defaut) : */
/* */
/* M=1 */
/* D=2 */
/* R=1 */
/* A=1 */
/* B=2 */
/* C=2 */
/* */
/* avec : */
/* */
/* u E [ 0 , 2.p ] */
/* */
/* v E [ 0 , 4.p ] */
/* */
/* (ou 'p' designe 'pi'). */
/* */
/* A cause de la derivee 'd#XPlan(u,v)' qui */
/* fait descendre 'A+cos(v)' en denominateur, */
/* il est imperatif que : */
/* */
/* A+cos(v) # 0 */
/* */
/* Si par malheur : */
/* */
/* A+cos(v) = 0 */
/* */
/* il convient d'ajouter un 'epsilon' a 'v' */
/* afin de supprimer cette nullite, d'ou la */
/* fonction 'v_corrige(v)' qui remplacera 'v' */
/* dans les definitions de 'Fx(u,v)' et de */
/* 'Fy(u,v)'. */
/* */
/* Ainsi, la coordonnee 'u' decrit la petite */
/* ellipse (de "rayons" 'Ra' et 'Rb') et 'v' */
/* decrit la courbe en forme de '8'. */
/* */
/* Tout ceci est tres inspire de 'v $xtc/huit.03$c' */
/* en faisant : */
/* */
/* u = phi E [0,2.p] */
/* v = theta E [0,4.p] */
/* */
/* */
/*************************************************************************************************************************************/
/* :Fin_listMN_HUIT_11: */
#include xrs/huit.12.I"
#include xrs/huit.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 RENORMALISER_ARBITRAIREMENT_LES_DIFFERENTIELLES \
FAUX \
/* Pour ne pas renormaliser arbitrairement les differentielles, et ce afin de prendre en */ \
/* compte leur extrema veritables dans 'RENORMALISATION_D_UNE_LISTE_DE_POINTS'... */
#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 (introduit le 20161203101444...). */
#include xrs/huit.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/huit.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 P S E U D O - T O R E A D E U X T R O U S : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/surfaces.41.I"
#include xrs/surfaces.22.I"
/*..............................................................................................................................*/
INITIALISATIONS_GENERALES;
/* Initialisations generales faites au tout debut... */
iTRANSFORMAT_31(liste_PARAMETRE_Ra,PARAMETRE_Ra_IMPLICITE);
/* Initialisation du parametre 'Ra'. */
iTRANSFORMAT_31(liste_PARAMETRE_Rb,PARAMETRE_Rb_IMPLICITE);
/* Initialisation du parametre 'Rb'. */
iTRANSFORMAT_31(liste_PARAMETRE_M,PARAMETRE_M_IMPLICITE);
/* Initialisation du parametre 'M'. */
iTRANSFORMAT_31(liste_PARAMETRE_D,PARAMETRE_D_IMPLICITE);
/* Initialisation du parametre 'D'. */
iTRANSFORMAT_31(liste_PARAMETRE_R,PARAMETRE_R_IMPLICITE);
/* Initialisation du parametre 'R'. */
iTRANSFORMAT_31(liste_PARAMETRE_A,PARAMETRE_A_IMPLICITE);
/* Initialisation du parametre 'A'. */
iTRANSFORMAT_31(liste_PARAMETRE_B,PARAMETRE_B_IMPLICITE);
/* Initialisation du parametre 'B'. */
iTRANSFORMAT_31(liste_PARAMETRE_C,PARAMETRE_C_IMPLICITE);
/* Initialisation du parametre 'C'. */
#include xrv/champs_5.1A.I"
#include xrs/surfaces.13.I"
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES;
GET_FICHIER_PARAMETRE("PARAMETRE_Ra=",_Ra);
GET_FICHIER_PARAMETRE("PARAMETRE_Rb=",_Rb);
GET_FICHIER_PARAMETRE("PARAMETRE_M=",_M);
GET_FICHIER_PARAMETRE("PARAMETRE_D=",_D);
GET_FICHIER_PARAMETRE("PARAMETRE_R=",_R);
GET_FICHIER_PARAMETRE("PARAMETRE_A=",_A);
GET_FICHIER_PARAMETRE("PARAMETRE_B=",_B);
GET_FICHIER_PARAMETRE("PARAMETRE_C=",_C);
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_DIFFERENTIATION_PAR_DIFFERENCES;
/* Introduit le 20161205111315... */
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_DIFFERENTIATION(BLOC(Bblock
EGAL(parametre_Ra,sPARAMETRE_Ra(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'Ra'. */
EGAL(parametre_Rb,sPARAMETRE_Rb(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'Rb'. */
EGAL(parametre_M,sPARAMETRE_M(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'M'. */
EGAL(parametre_D,sPARAMETRE_D(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'D'. */
EGAL(parametre_R,sPARAMETRE_R(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'R'. */
EGAL(parametre_A,sPARAMETRE_A(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'A'. */
EGAL(parametre_B,sPARAMETRE_B(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'B'. */
EGAL(parametre_C,sPARAMETRE_C(numero_de_la_periode_courante));
/* Calcul de la valeur du parametre 'C'. */
Eblock
)
,Fx_Huit_1,Fy_Huit_1,Fz_Huit_1
/* Fonctions definissant parametriquement la surface... */
,Pxyz_Huit_1
,BLOC(VIDE;)
/* Il n'y a rien de specifique a calculer pour le point courant... */
);
/* On notera l'usage de 'VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENCES____(...)' et non */
/* pas de 'a l'usage de 'VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENTIATION(...)' a cause */
/* du fait que la derivee de 'ATAN(...)' n'est pas definie a la date du 20161202111945... */
/* */
/* La derivee de 'ATAN(...)' a ete introduite le 'v $xil/defi_K2$vv$DEF 20161202134722', */
/* mais je ne l'utilise par car, en effet, etant donnee la complexite de la definition */
/* de cette surface (via en particulier des derivees 'd#'), j'ai peur que la compilation */
/* soit trop "lourde"... */
/* */
/* A la date du 20161205184640, j'ai compris pourquoi la procedure */
/* 'VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENTIATION(...)' ne pouvait pas etre */
/* utilisee a cette date ('v $xrs/huit.11$I 20161205165813'). Ce probleme a ete */
/* resolu dans les jours suivants (le 20161208090234...). */
RETU_Commande;
Eblock
ECommande