/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S V A G U E S 1 D E J E E N E R : */
/* */
/* */
/* Exemple : */
/* */
/* L'image 'v $xiirs/VAGU.11.21' */
/* en est un exemple... */
/* */
/* */
/* Author of '$xrs/vagues.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20200319171642). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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_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/vagues.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/vagues.19.I"
/*===================================================================================================================================*/
/* :Debut_listMN_VAGUES_DE_JEENER_11: */
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S V A G U E S 1 D E J E E N E R : */
/* */
/* */
/* Definition de la surface (due a Patrice Jeener, 'v $xrs/vagues.11$K') : */
/* */
/* Elles sont definies parametriquement */
/* en fonction des deux parametres 'u' et 'v' : */
/* */
/* */
/* Comb(n,m) = Combinaisons de 'n' objets 'm' a 'm', */
/* Prim(f,x) = Primitive de la fonction reelle 'f' par rapport a 'x'. */
/* */
/* */
/* V.U A.Prim(1,v) */
/* F (u,v) = (-1) .[------------- */
/* x 1! */
/* */
/* 2 0 0 2 */
/* B.Comb(2,0).Prim(v ,v).(iu) + C.Comb(2,2).Prim(v ,v).(iu) */
/* - ------------------------------------------------------------- */
/* 3! */
/* */
/* 4 0 2 2 0 4 */
/* D.Comb(4,0).Prim(v ,v).(iu) + E.Comb(4,2).Prim(v ,v).(iu) + F.Comb(4,4).Prim(v ,v).(iu) */
/* + -------------------------------------------------------------------------------------------- */
/* 5! */
/* */
/* 6 0 4 2 2 4 0 6 */
/* G.Comb(6,0).Prim(v ,v).(iu) + H.Comb(6,2).Prim(v ,v).(iu) + I.Comb(6,4).Prim(v ,v).(iu) + J.Comb(6,6).Prim(v ,v).(iu) */
/* - --------------------------------------------------------------------------------------------------------------------------- */
/* 7! */
/* ] */
/* */
/* */
/* 2 2 */
/* 1 log(v + u + L) */
/* F (u,v) = K.---.------------------ */
/* y 2 0! */
/* */
/* 1 0 */
/* M.Comb(1,0).Prim(v ,v).(iu) */
/* - ------------------------------ */
/* 2! */
/* */
/* 3 0 1 2 */
/* N.Comb(3,0).Prim(v ,v).(iu) + P.Comb(3,2).Prim(v ,v).(iu) */
/* + ------------------------------------------------------------- */
/* 4! */
/* */
/* 5 0 3 2 1 4 */
/* Q.Comb(5,0).Prim(v ,v).(iu) + R.Comb(5,2).Prim(v ,v).(iu) + S.Comb(5,4).Prim(v ,v).(iu) */
/* - -------------------------------------------------------------------------------------------- */
/* 6! */
/* */
/* */
/* F (u,v) = T.atan(v,u) + U.pi */
/* z */
/* */
/* avec : */
/* */
/* u E [ -pi/2, +pi/2 ] */
/* */
/* v E [ 0 , +5.pi/4 ] */
/* */
/* et : */
/* */
/* A = +1 */
/* B = +1 */
/* C = +1 */
/* D = +1 */
/* E = +1 */
/* F = +1 */
/* G = +1 */
/* H = +1 */
/* I = +1 */
/* J = +1 */
/* K = +epsilon (afin d'eviter des problemes avec le 'log(...)') */
/* L = +1 */
/* M = +1 */
/* N = +1 */
/* P = +1 */
/* Q = +1 */
/* R = +1 */
/* S = +1 */
/* T = +1 */
/* U E [Um,UM] = [0,0] (segment [0,3] par exemple...) */
/* V = 1 (V=0 inhibe la bascule 'U' de la fonction 'Fx(...)') */
/* */
/* en general... */
/* */
/* */
/*************************************************************************************************************************************/
/* :Fin_listMN_VAGUES_DE_JEENER_11: */
#include xrs/vagues.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/vagues.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/vagues.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 E S V A G U E S 1 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/vagues.12.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);
iTRANSFORMAT_31(liste_PARAMETRE_M,PARAMETRE_M_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_N,PARAMETRE_N_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_P,PARAMETRE_P_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_Q,PARAMETRE_Q_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_R,PARAMETRE_R_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_S,PARAMETRE_S_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_T,PARAMETRE_T_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_Um,PARAMETRE_Um_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_UM,PARAMETRE_UM_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_V,PARAMETRE_V_IMPLICITE);
/* Initialisation des parametres {a,...v}. */
#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);
GET_FICHIER_PARAMETRE("PARAMETRE_M=",_M);
GET_FICHIER_PARAMETRE("PARAMETRE_N=",_N);
GET_FICHIER_PARAMETRE("PARAMETRE_P=",_P);
GET_FICHIER_PARAMETRE("PARAMETRE_Q=",_Q);
GET_FICHIER_PARAMETRE("PARAMETRE_R=",_R);
GET_FICHIER_PARAMETRE("PARAMETRE_S=",_S);
GET_FICHIER_PARAMETRE("PARAMETRE_T=",_T);
GET_FICHIER_PARAMETRE("PARAMETRE_Um=",_Um);
GET_FICHIER_PARAMETRE("PARAMETRE_UM=",_UM);
GET_FICHIER_PARAMETRE("PARAMETRE_V=",_V);
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_GENERALE_DE_LA_SURFACE(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));
EGAL(parametre_m,sPARAMETRE_M(numero_de_la_periode_courante));
EGAL(parametre_n,sPARAMETRE_N(numero_de_la_periode_courante));
EGAL(parametre_p,sPARAMETRE_P(numero_de_la_periode_courante));
EGAL(parametre_q,sPARAMETRE_Q(numero_de_la_periode_courante));
EGAL(parametre_r,sPARAMETRE_R(numero_de_la_periode_courante));
EGAL(parametre_s,sPARAMETRE_S(numero_de_la_periode_courante));
EGAL(parametre_t,sPARAMETRE_T(numero_de_la_periode_courante));
EGAL(parametre_um,sPARAMETRE_Um(numero_de_la_periode_courante));
EGAL(parametre_uM,sPARAMETRE_UM(numero_de_la_periode_courante));
EGAL(parametre_v,sPARAMETRE_V(numero_de_la_periode_courante));
/* Calcul de la valeur des parametres {a,...,v}. */
Test(N_EST_PAS_ENTIER(parametre_um))
Bblock
PRINT_ATTENTION("le parametre 'um' n'est pas entier");
Eblock
ATes
Bblock
Eblock
ETes
Test(N_EST_PAS_ENTIER(parametre_uM))
Bblock
PRINT_ATTENTION("le parametre 'uM' n'est pas entier");
Eblock
ATes
Bblock
Eblock
ETes
Eblock
)
,BLOC(Bblock
DEFV(Float,INIT(facteur_de_la_translation_de_uv,FZERO));
/* Afin de tester le changement des associations {u,v} <--> {1,2}. */
DoIn(parametre_u,parametre_um,parametre_uM,FU)
Bblock
EGAL(facteur_de_la_translation_de_u_effectif_lors_d_un_lister
,facteur_de_la_translation_de_uv
);
EGAL(facteur_de_la_translation_de_v_effectif_lors_d_un_lister
,facteur_de_la_translation_de_uv
);
GENERATION_D_UNE_IMAGE_DE_LA_SURFACE_SANS_DIFFERENTIATION(Fx_VaguesDeJeener_1
,Fy_VaguesDeJeener_1
,Fz_VaguesDeJeener_1
,BLOC(VIDE;)
)
/* Il n'y a rien de specifique a calculer pour le point courant... */
INCR(facteur_de_la_translation_de_uv,FU);
Eblock
EDoI
Eblock
)
,Fx_VaguesDeJeener_1,Fy_VaguesDeJeener_1,Fz_VaguesDeJeener_1
/* Fonctions definissant parametriquement la surface... */
,Pxyz_VaguesDeJeener_1
);
/* On notera que c'est 'GENERATION_D_UNE_IMAGE_DE_LA_SURFACE_SANS_DIFFERENTIATION(...)' qui */
/* est utilies et non 'GENERATION_D_UNE_IMAGE_DE_LA_SURFACE_AVEC_DIFFERENTIATION(...)' a */
/* cause de la "lourdeur" des fonctions {Fx,Fy,Fz}... */
RETU_Commande;
Eblock
ECommande