/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N P R O J E C T O R B I D I M E N S I O N N E L : */
/* */
/* */
/* Author of '$xrs/project2D.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20041107113716). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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_____MODULE_NON_DERIVABLE
@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
/* Introduit le 20041229181213 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 20111209172032... */
#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/project2D.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/project2D.19.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N P R O J E C T O R B I D I M E N S I O N N E L : */
/* */
/* */
/* Definition de la surface ('v $xrs/project2D.11$K') : */
/* */
/* Le projector bidimensionnel 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') : */
/* */
/* F (u,v) = A .ImageX(u,v) + T */
/* x x x */
/* */
/* F (u,v) = A .ImageY(u,v) + T */
/* y y y */
/* */
/* F (u,v) = A .ImageZ(u,v) + T */
/* z z z */
/* */
/* (ou 'Image?(u,v)' sont des images flottantes indexees */
/* par (u,v) -coordonnees flottantes {Xf,Yf} respectivement- */
/* a l'aide d'une procedure d'interpolation bicubique) */
/* */
/* avec : */
/* */
/* u E [ 0 , 1 ] (--O--> Xf) */
/* */
/* v E [ 0 , 1 ] (--O--> Yf) */
/* */
/* */
/*************************************************************************************************************************************/
#include xrs/project2D.12.I"
#include xrs/project2D.14.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S T R O I S I M A G E S D E F I N I S S A N T */
/* L E S T R O I S P R O J E C T I O N S : */
/* */
/*************************************************************************************************************************************/
#define NORMALISER_LES_COMPOSANTES \
FAUX
DEFV(Local,DEFV(Logical,INIT(normaliser_les_composantes,NORMALISER_LES_COMPOSANTES)));
/* Introduit le 20091023210137 pour etre utilise lors d'un appel a 'v $xrs/SurfParti.01$Z'. */
#include xrs/project2D.1A.I"
/* Introduit sous cette forme le 20150604143557... */
DEFV(Local,DEFV(Float,INIT(Image_Fx_amplificateur,AMPLIFICATEUR_DES_IMAGES)));
DEFV(Local,DEFV(Float,INIT(Image_Fy_amplificateur,AMPLIFICATEUR_DES_IMAGES)));
DEFV(Local,DEFV(Float,INIT(Image_Fz_amplificateur,AMPLIFICATEUR_DES_IMAGES)));
/* Amplificateurs pour {Fx,Fy,Fz}. */
/* */
/* Reintroduit le 20150605143901 hors de 'ENTREE_DE_LA_DEFINITION_DES_IMAGES_Fx_Fy_Fz' */
/* afin de faciliter l'usage 'v $xrv/particule.10$K xrs.project2D.11.I'... */
DEFV(Local,DEFV(Float,INIT(Image_Fx_translateur,TRANSLATEUR_DES_IMAGES)));
DEFV(Local,DEFV(Float,INIT(Image_Fy_translateur,TRANSLATEUR_DES_IMAGES)));
DEFV(Local,DEFV(Float,INIT(Image_Fz_translateur,TRANSLATEUR_DES_IMAGES)));
/* Translateurs pour {Fx,Fy,Fz}. */
/* */
/* Reintroduit le 20150605143901 hors de 'ENTREE_DE_LA_DEFINITION_DES_IMAGES_Fx_Fy_Fz' */
/* afin de faciliter l'usage 'v $xrv/particule.10$K xrs.project2D.11.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"
/* Le 20060505173308 en calculant 'v $xiirs/SPHE.72' et 'v $xiirs/SPHE.73.2D', j'ai note */
/* une difference de comportement au niveau du parametre "eda=" (ou de son "complement" */
/* "edc=") venant du fait que le programme 'v $xrs/project2D.11$K' ne definit pas */
/* 'RENORMALISER_ARBITRAIREMENT_LES_DIFFERENTIELLES' et c'est donc la valeur par */
/* defaut ('VRAI') qui est utilisee. Alors qu'en ce qui concerne 'v $xrs/project3D.11$K', */
/* celui-ci lui donne la valeur 'FAUX'... */
#define EDITER_LES_EXTREMA_S_ILS_SONT_HORS_PREVISIONS \
FAUX \
/* Afin de ne pas editer par defaut a la fin la valeur des extrema de {x,y,z} et {dx,dy,dz} */ \
/* lorsqu'ils sont en dehors de {?min_ESPACE,?max_ESPACE} et de {d?min_ESPACE,d?max_ESPACE} */ \
/* respectivement. Ceci a ete introduit le 20041112090812. */
#include xrs/project2D.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/project2D.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 R O J E C T O R B I D I M E N S I O N N E L : */
/* */
/*************************************************************************************************************************************/
/* Le 20041120122429, je note qu'il est impossible de mettre ici (comme je l'esperai */
/* afin d'optimiser les 'LIO9(...)' de 'v $xrs/surfaces.12$I LIO9' qui sont un peu */
/* "lourds" au niveau de la compilation...) : */
/* */
/* #include xrs/surfaces.31.I" */
/* GENERE_LES_DIFFERENTIELLES_DES_FONCTIONS_F(Fx_Projector2D_1 */
/* ,Fy_Projector2D_1 */
/* ,Fz_Projector2D_1 */
/* ); */
/* */
/* car, en effet, les fonctions {Fx_Projector2D_1,Fy_Projector2D_1,Fz_Projector2D_1} */
/* referencent {Image_Fx,Image_Fy,Image_Fz} qui ne sont pas encore definies. */
/* */
/* Le fichier '$xrs/surfaces.31$I' est donc mis dans la categorie '$INUTILE' (afin de ne */
/* pas le perdre car il est malgre tout assez spactaculaire et que cela pourra peut-etre */
/* servir un jour... */
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/surfaces.41.I"
/* Introduit le 20150528140108... */
#include xrs/project2D.1B.I"
/* Introduit sous cette forme le 20150604143557... */
#include xrs/surfaces.22.I"
/* Introduit le 20040512103205 a des fins de faisabilite, avant de l'introduire dans tous */
/* les '$K' de '$xrs' apres cette date... */
/*..............................................................................................................................*/
EGAL(AXES_COORDONNEES_FERMES_OUVERTS_____compatibilite_20090929,VRAI);
/* Introduit le 20101005092845 car, en effet, sans cela les images faites anterieurement */
/* au 20090929092932 n'etaient pas reproductibles (des defauts importants apparaissaient */
/* en particulier a la peripherie des "objets"). */
/* */
/* ATTENTION, lors de la generation de 'v $xiirs/PROK.P3.0', il a fallu mettre cette option */
/* a 'FAUX' ('v $xiirs/$Fnota COORDONNEES_____compatibilite_20090929.FAUX'). Cela a conduit */
/* a l'introduction 'v $xcc/LFonctionR$vv$D/ArgumentsUse$vv$gre _____compatibilite_'... */
INITIALISATIONS_GENERALES;
/* Initialisations generales faites au tout debut... */
#include xrv/champs_5.1A.I"
#define DEFINITION_DES_INTERPOLATIONS_Fx_Fy_Fz
/* Afin de generer les 'GET_ARGUMENT_L(...)' de compatibilite (introduit le 20150528104906). */
#include xrs/surfaces.13.I"
#include xrs/project2D.1C.I"
/* Introduit le 20150604143557... */
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES;
ENTREE_DE_LA_DEFINITION_DES_IMAGES_Fx_Fy_Fz;
/* Introduit sous cette forme le 20150604143557... */
GET_ARGUMENT_F("aFx=",Image_Fx_amplificateur);
GET_ARGUMENT_F("tFx=""bFx=",Image_Fx_translateur);
GET_ARGUMENT_F("aFy=",Image_Fy_amplificateur);
GET_ARGUMENT_F("tFy=""bFy=",Image_Fy_translateur);
GET_ARGUMENT_F("aFz=",Image_Fz_amplificateur);
GET_ARGUMENT_F("tFz=""bFz=",Image_Fz_translateur);
/* Reintroduit le 20150605143901 hors de 'ENTREE_DE_LA_DEFINITION_DES_IMAGES_Fx_Fy_Fz' */
/* afin de faciliter l'usage 'v $xrv/particule.10$K xrs.project2D.11.I'... */
GET_ARGUMENT_L("lineaire=""bilineaire=""lineaire2D=""bilineaire2D="
,FFload_point_coordonnees_01_____utiliser_l_interpolation_bilineaire
);
GET_ARGUMENT_N("cubique=""bicubique=""cubique2D=""bicubique2D="
,FFload_point_coordonnees_01_____utiliser_l_interpolation_bilineaire
);
/* ATTENTION : 'FFload_point_coordonnees_01_____utiliser_l_interpolation_bilineaire' est */
/* aussi reference dans 'v $xrs/surfaces.13$I tbilineaire=' (note le 20060602093328...). */
/* */
/* Introduit hors de 'ENTREE_DE_LA_DEFINITION_DES_IMAGES_Fx_Fy_Fz' afin d'eviter des doubles */
/* definitions avec 'v $xrv/particule.10$K ENTREE_DE_LA_DEFINITION_DES_IMAGES_Fx_Fy_Fz'. */
GET_ARGUMENT_L("normaliser_composantes=""normaliser=""norm=",normaliser_les_composantes);
GET_ARGUMENT_F("origine=""o=",IFnormalisation_____niveau_origine___de_normalisation);
GET_ARGUMENT_F("extremite=""e=",IFnormalisation_____niveau_extremite_de_normalisation);
GET_ARGUMENT_L("conserver_zero=""cz=",IFnormalisation_____conserver_le_zero);
/* Arguments introduits le 20091023210137... */
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_TRANSFORMATION_DE_u_v;
/* Introduit le 20060117145814. */
PROCESS_ARGUMENTS_DE_PERMUTATION_DES_COORDONNEES_ET_DES_DIFFERENTIELLES;
/* Introduit le 20050615152759. */
PROCESS_ARGUMENTS_DE_PROJECTION_DE_LA_SURFACE;
/* Introduit le 20050714092926 (depuis que 'v $xrs/surfaces.14$I 20050713143505'...). */
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_TEXTURE;
PROCESS_ARGUMENTS_DE_VISUALISATION_DE_LA_NORMALE;
PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;
/* On notera le 20170111114619 qu'en ajustant specifiquement 'u' ou 'v', on peut ainsi */
/* tracer des courbes. En s'inspirant de 'v $xiirs/PROJ.11', en faisant : */
/* */
/* mu=0 Mu=1.0 */
/* mv=0.5 Mv=0.5 mv=Mv */
/* */
/* on obtient l'equateur d'une sphere, alors qu'avec : */
/* */
/* mu=0.25 Mu=0.25 mu=Mu */
/* mv=0 Mv=1.0 */
/* */
/* on obtient un meridien. ATTENTION : on n'oubliera pas d'ajuster si besoin est les */
/* differentes rotations... */
/* */
/* Au passage, pour obtenir une sphere complete, il faut prendre : */
/* */
/* mu=0 Mu=0.5 */
/* mv=0 Mv=1.0 */
/* */
/* On notera aussi que l'on peut exploiter une droite {u,v}. Ainsi, avec : */
/* */
/* Aduv=+16.0 */
/* Bduv=-1.0 */
/* Cduv=0 */
/* */
/* et les options : */
/* */
/* maillage_uv=FAUX droite_uv=VRAI */
/* */
/* on obtient une "droite" qui s'enroule Aduv/(Mv/Mu) (= 16.0/(1.0/0.5) = 8) fois autour */
/* de la sphere, d'un pole a l'autre... */
PROCESS_ARGUMENTS_DE_DEFINITION_D_UNE_DROITE_DU_PLAN_UV;
/* Introduit le 20060831113702... */
PROCESS_ARGUMENTS_DE_DIFFERENTIATION_PAR_DIFFERENCES;
/* Introduit le 20041118144233... */
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
BDEFV(imageF,Image_Fx_Argument);
BDEFV(imageF,Image_Fy_Argument);
BDEFV(imageF,Image_Fz_Argument);
/* Definition des trois projections Argument... */
Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(Image_Fx_Argument,nom_imageFx))))
Bblock
Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(Image_Fy_Argument,nom_imageFy))))
Bblock
Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(Image_Fz_Argument,nom_imageFz))))
Bblock
BDEFV(imageF,Image_Fx);
BDEFV(imageF,Image_Fy);
BDEFV(imageF,Image_Fz);
/* Definition des trois projections apres une eventuelle normalisation... */
Test(IL_FAUT(normaliser_les_composantes))
/* Possibilite introduite le 20091023210137... */
Bblock
CALS(IFnormalisation_automatique(Image_Fx,Image_Fx_Argument));
CALS(IFnormalisation_automatique(Image_Fy,Image_Fy_Argument));
CALS(IFnormalisation_automatique(Image_Fz,Image_Fz_Argument));
Eblock
ATes
Bblock
CALS(IFmove(Image_Fx,Image_Fx_Argument));
CALS(IFmove(Image_Fy,Image_Fy_Argument));
CALS(IFmove(Image_Fz,Image_Fz_Argument));
Eblock
ETes
VISUALISATION_DE_LA_SURFACE_AVEC_DIFFERENCES____(BLOC(VIDE;)
,Fx_Projector2D_1,Fy_Projector2D_1,Fz_Projector2D_1
/* Fonctions definissant parametriquement la surface... */
,Pxyz_Projector2D_1
,BLOC(VIDE;)
/* Il n'y a rien de specifique a calculer pour le point courant... */
);
/* Nota : les noms '*_Projector_1' ont ete changes en '*_Projector_2' le 20050619170241, */
/* puis en '*_Projector2D_1' le 20050620102241... */
EDEFV(imageF,Image_Fz);
EDEFV(imageF,Image_Fy);
EDEFV(imageF,Image_Fx);
/* Definition des trois projections apres une eventuelle normalisation... */
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
EDEFV(imageF,Image_Fz_Argument);
EDEFV(imageF,Image_Fy_Argument);
EDEFV(imageF,Image_Fx_Argument);
/* Definition des trois projections Argument... */
Eblock
end_nouveau_block
RETU_Commande;
Eblock
ECommande