/*************************************************************************************************************************************/ /* */ /* D E F I N I T I O N D E L A S P I R A L E D E B O U T E I L L E S D E J E E N E R : */ /* */ /* */ /* Exemple : */ /* */ /* L'image 'v $xiirs/BKLN.K1.1$LUMIN' */ /* en est un exemple... */ /* */ /* */ /* Author of '$xrs/bKlein.G1$K' : */ /* */ /* Jean-Francois COLONNA (LACTAMME, 20200219110513). */ /* */ /*************************************************************************************************************************************/ /*===================================================================================================================================*/ /*************************************************************************************************************************************/ /* */ /* 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_GAIRE_DE_L_ARITHMETIQUE_ETENDUE_DES_NOMBRES_GLOTTANTS @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/bKlein.G8.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/bKlein.G9.I" /*===================================================================================================================================*/ /* :Debut_listMN_SPIRALE_DE_JEENER_G1: */ /*************************************************************************************************************************************/ /* */ /* D E F I N I T I O N D E L A S P I R A L E D E B O U T E I L L E S D E J E E N E R : */ /* */ /* */ /* Definition de la surface (due a Patrice Jeener, 'v $xrs/bKlein.G1$K') : */ /* */ /* Elle est definie parametriquement */ /* en fonction des deux parametres 'u' */ /* et 'v' : */ /* */ /* W(u,v) = J.sin(K.v) + L */ /* */ /* v */ /* --- */ /* A */ /* F (u,v) = e .[(B.cos(C.v)) + (D.cos(E.v)) + (F.W(u,v).sin(G.v).cos(H.u))] */ /* x */ /* */ /* v */ /* --- */ /* A */ /* F (u,v) = e .[(B.sin(C.v)) - (D.sin(E.v)) + (F.W(u,v).cos(G.v).cos(H.u))] */ /* y */ /* */ /* v */ /* --- */ /* A */ /* F (u,v) = e .W(u,v).sin(I.u) */ /* z */ /* */ /* avec : */ /* */ /* u E [ -pi, +pi ] */ /* */ /* v E [ -2.pi, +2.pi ] */ /* */ /* et : */ /* */ /* A = +10 */ /* B = +8 */ /* C = +1 */ /* D = +1 */ /* E = B */ /* F = -9/16 */ /* G = +7/2 */ /* H = +1 */ /* I = +1 */ /* J = +1 */ /* K = +9 */ /* L = +3/2 */ /* */ /* en general... */ /* */ /* */ /*************************************************************************************************************************************/ /* :Fin_listMN_SPIRALE_DE_JEENER_G1: */ #include xrs/bKlein.G4.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/bKlein.G3.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/bKlein.G1.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 E L A S P I R A L E D E B O U T E I L L E S D E J E E N E R : */ /* */ /*************************************************************************************************************************************/ BCommande(nombre_d_arguments,arguments) /*-----------------------------------------------------------------------------------------------------------------------------------*/ Bblock #include xrs/surfaces.41.I" #include xrs/surfaces.22.I" #include xrs/bKlein.G2.I" /*..............................................................................................................................*/ INITIALISATIONS_GENERALES; /* Initialisations generales faites au tout debut... */ iTRANSFORMAT_31(liste_PARAMETRE_A,PARAMETRE_A_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_B,PARAMETRE_B_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_C,PARAMETRE_C_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_D,PARAMETRE_D_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_E,PARAMETRE_E_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_F,PARAMETRE_F_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_G,PARAMETRE_G_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_H,PARAMETRE_H_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_I,PARAMETRE_I_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_J,PARAMETRE_J_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_K,PARAMETRE_K_IMPLICITE); iTRANSFORMAT_31(liste_PARAMETRE_L,PARAMETRE_L_IMPLICITE); /* Initialisation des parametres {a,...l}. */ #include xrv/champs_5.1A.I" #include xrs/surfaces.13.I" GET_ARGUMENTSv(nombre_d_arguments ,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES; GET_FICHIER_PARAMETRE("PARAMETRE_A=",_A); GET_FICHIER_PARAMETRE("PARAMETRE_B=",_B); GET_FICHIER_PARAMETRE("PARAMETRE_C=",_C); GET_FICHIER_PARAMETRE("PARAMETRE_D=",_D); GET_FICHIER_PARAMETRE("PARAMETRE_E=",_E); GET_FICHIER_PARAMETRE("PARAMETRE_F=",_F); GET_FICHIER_PARAMETRE("PARAMETRE_G=",_G); GET_FICHIER_PARAMETRE("PARAMETRE_H=",_H); GET_FICHIER_PARAMETRE("PARAMETRE_I=",_I); GET_FICHIER_PARAMETRE("PARAMETRE_J=",_J); GET_FICHIER_PARAMETRE("PARAMETRE_K=",_K); GET_FICHIER_PARAMETRE("PARAMETRE_L=",_L); 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_DIFFERENTIATION(BLOC(Bblock EGAL(parametre_a,sPARAMETRE_A(numero_de_la_periode_courante)); EGAL(parametre_b,sPARAMETRE_B(numero_de_la_periode_courante)); EGAL(parametre_c,sPARAMETRE_C(numero_de_la_periode_courante)); EGAL(parametre_d,sPARAMETRE_D(numero_de_la_periode_courante)); EGAL(parametre_e,sPARAMETRE_E(numero_de_la_periode_courante)); EGAL(parametre_f,sPARAMETRE_F(numero_de_la_periode_courante)); EGAL(parametre_g,sPARAMETRE_G(numero_de_la_periode_courante)); EGAL(parametre_h,sPARAMETRE_H(numero_de_la_periode_courante)); EGAL(parametre_i,sPARAMETRE_I(numero_de_la_periode_courante)); EGAL(parametre_j,sPARAMETRE_J(numero_de_la_periode_courante)); EGAL(parametre_k,sPARAMETRE_K(numero_de_la_periode_courante)); EGAL(parametre_l,sPARAMETRE_L(numero_de_la_periode_courante)); Eblock ) /* Calcul de la valeur des parametres {a,...,l}. */ ,Fx_BouteilleDeKlein_G,Fy_BouteilleDeKlein_G,Fz_BouteilleDeKlein_G /* Fonctions definissant parametriquement la surface... */ ,Pxyz_BouteilleDeKlein_G ,BLOC(VIDE;) /* Il n'y a rien de specifique a calculer pour le point courant... */ ); RETU_Commande; Eblock ECommande