/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N E I N T E R P O L A T I O N E N T R E */
/* U N E L L I P O I D E E T L A S U R F A C E D E B O Y : */
/* */
/* */
/* Exemple : */
/* */
/* L'image 'v $xias/EBOY.11$R16' */
/* en est un exemple... */
/* */
/* */
/* Nota : */
/* */
/* Ce programme realise l'interpolation */
/* entre un ellipsoide et la surface de Boy. */
/* On notera que l'on utilise l'ellipsoide */
/* et non pas la sphere a cause du parametre */
/* 'PARAMETRE_R' qui existe malheureusement */
/* a la fois dans les fichiers '$xrs/sphere.1*' */
/* et '$xrs/sBoy.1*'... */
/* */
/* */
/* Author of '$xrs/elli_sBoy.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 19990121144036). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 20080126104953, mais le 20080126173152 j'ai du y renoncer, */
/* les fichiers generes etant monstrueux (plus de 10 Go). */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 20041229181027 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 20130306152417... */
#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/sBoy.18.I"
/* Introduit le 20030130120907 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/sBoy.19.I"
/* Introduit le 20030130120907 a cause de 'v $xcc/cpp$Z INCLUDES_1.sed'. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E L ' I N T E R P O L A T I O N E N T R E */
/* L E P L A N E T L A S U R F A C E D E B O Y : */
/* */
/*************************************************************************************************************************************/
#include xrs/sBoy.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/ellipsoid.12.I"
/* ATTENTION, il est imperatif que l'ordre soit 'v $xrs/sBoy.13$I', puis */
/* 'v $xrs/ellipsoid.13$I' a cause de la definition de 'RAYON_DE_VISUALISATION' qui est */
/* faite par le premier de ces fichiers, via 'v $xrs/surfaces.11$I'... */
#include xrs/sBoy.13.I"
/* ATTENTION, il est imperatif que l'ordre soit 'v $xrs/sBoy.13$I', puis */
/* 'v $xrs/ellipsoid.13$I' a cause de la definition de 'RAYON_DE_VISUALISATION' qui est */
/* faite par le premier de ces fichiers, via 'v $xrs/surfaces.11$I'... */
#include xrs/ellipsoid.13.I"
/* ATTENTION, il est imperatif que l'ordre soit 'v $xrs/sBoy.13$I', puis */
/* 'v $xrs/ellipsoid.13$I' a cause de la definition de 'RAYON_DE_VISUALISATION' qui est */
/* faite par le premier de ces fichiers, via 'v $xrs/surfaces.11$I'... */
#include xrs/elli_sBoy.13.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E L ' I N T E R P O L A T I O N E N T R E */
/* L E P L A N E T L A S U R F A C E D E B O Y : */
/* */
/*************************************************************************************************************************************/
#include xrs/sBoy.12.I"
#include xrs/elli_sBoy.12.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/ellipsoid.11.I"
#include xrs/sBoy.11.I"
#include xrs/elli_sBoy.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 I N T E R P O L A T I O N E N T R E */
/* U N E L L I P O I D E E T L A S U R F A C E D E B O Y : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/surfaces.41.I"
/* Introduit le 20150528135749... */
#include xrs/surfaces.22.I"
/*..............................................................................................................................*/
INITIALISATIONS_GENERALES;
/* Initialisations generales faites au tout debut... */
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'. */
iTRANSFORMAT_31(liste_PARAMETRE_R,PARAMETRE_R_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_C0,PARAMETRE_C0_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_C1,PARAMETRE_C1_IMPLICITE);
iTRANSFORMAT_31(liste_PARAMETRE_C2,PARAMETRE_C2_IMPLICITE);
/* Initialisation des parametres 'R' et '{C0,C1,C2}'. */
iTRANSFORMAT_31(liste_PONDERATION_DE_L_ELLIPSOIDE,PONDERATION_DE_L_ELLIPSOIDE_IMPLICITE);
iTRANSFORMAT_31(liste_PONDERATION_DE_LA_SURFACE_DE_BOY,PONDERATION_DE_LA_SURFACE_DE_BOY_IMPLICITE);
/* Initialisation des ponderations de l'ellipsoide et de la surface de Boy. */
#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_R=",_R);
GET_FICHIER_PARAMETRE("PARAMETRE_C0=",_C0);
GET_FICHIER_PARAMETRE("PARAMETRE_C1=",_C1);
GET_FICHIER_PARAMETRE("PARAMETRE_C2=",_C2);
PROCESS_ARGUMENT_FICHIER("PONDERATION_DE_L_ELLIPSOIDE="
,fichier_PONDERATION_DE_L_ELLIPSOIDE
,liste_PONDERATION_DE_L_ELLIPSOIDE
,PONDERATION_DE_L_ELLIPSOIDE_IMPLICITE
,gTRANSFORMAT_31
);
PROCESS_ARGUMENT_FICHIER("PONDERATION_DE_LA_SURFACE_DE_BOY="
,fichier_PONDERATION_DE_LA_SURFACE_DE_BOY
,liste_PONDERATION_DE_LA_SURFACE_DE_BOY
,PONDERATION_DE_LA_SURFACE_DE_BOY_IMPLICITE
,gTRANSFORMAT_31
);
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;
/* Introduit le 20060117145618. */
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
/* Introduit le 20050615091756. */
PROCESS_ARGUMENTS_DE_PROJECTION_DE_LA_SURFACE;
/* Introduit le 20041126095232. */
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 20060831113542... */
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));
/* 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'. */
EGAL(rayon_surface_de_boy,sPARAMETRE_R(numero_de_la_periode_courante));
EGAL(coefficient0_surface_de_boy
,sPARAMETRE_C0(numero_de_la_periode_courante)
);
EGAL(coefficient1_surface_de_boy
,sPARAMETRE_C1(numero_de_la_periode_courante)
);
EGAL(coefficient2_surface_de_boy
,sPARAMETRE_C2(numero_de_la_periode_courante)
);
Eblock
)
/* Calcul de la valeur des parametres 'R' et '{C0,C1,C2}'. */
EGAL(ponderation_du_Ellipsoide
,sPONDERATION_DE_L_ELLIPSOIDE(numero_de_la_periode_courante)
);
EGAL(ponderation_de_la_surface_de_Boy
,sPONDERATION_DE_LA_SURFACE_DE_BOY(numero_de_la_periode_courante)
);
/* Calcul des parametres de ponderation de l'ellipsoide et de la surface de Boy. */
,Fx_Ellipsoide_SurfaceDeBoy_1
,Fy_Ellipsoide_SurfaceDeBoy_1
,Fz_Ellipsoide_SurfaceDeBoy_1
/* Fonctions definissant parametriquement la surface... */
,Pxyz_Ellipsoide_SurfaceDeBoy_1
,BLOC(VIDE;)
/* Il n'y a rien de specifique a calculer pour le point courant... */
);
RETU_Commande;
Eblock
ECommande