/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S F O N C T I O N S P R O P R E S D E L ' A T O M E D ' H Y D R O G E N E */
/* E T V I S U A L I S A T I O N S O U S F O R M E D E S P H E R E S T R I D I M E N S I O N N E L L E S : */
/* */
/* */
/* .:. */
/* -+..-o*-.+oo.:+*o..-- */
/* .-.:+oo:+#-o.:--:*:**--o*+.-. */
/* .+oo. -+*. .:-++o++: .o-- .+**- */
/* .-oo-.. :o+. .:-++o++- .:+: ..:o*o- */
/* .:--.. :o+ ..::---:. :+. ..---. */
/* .-+.-oo- .+o. ..... .-+: :+oo-+*. */
/* .--... ..-+.-:.-.:-:-+.. .::.-+- */
/* :: .:-+- ....:.... :+o+. -+ */
/* .. ....:-+- .. . .. :+o+..-: .: */
/* ..:-:.-- .:-::-+- */
/* .:.:-. .:-+- :--. ... .:-- :-++: .+o+++ */
/* . .. ... .:--- .:-+- .---: ... .. .. */
/* .-- .. ... ..::. ... .. +o. */
/* ..-. .... ..:. .o+: */
/* .-- ... ... ... -o. */
/* .. .. ....... . .....:. .. .. */
/* ..:..- ....:: :+- :-:.:. ++:oo. */
/* .... . .... */
/* . . :. .... .... . ... .... .+.: . */
/* .::: . .... ... ... ... .::. . :oo- */
/* .. .. .-... ...:::. .:.+: .. :+ . */
/* .: .:: ... .. ...:... .: .-: -o+ +: */
/* ..: ... ... ... ..: .-- .:-- */
/* ...:-- ..:.:.::. -oo::. . */
/* . .- .... :: .:+-. :- ..:..+: . */
/* .-..:::.:--..:-:.+. */
/* .. .:: .. */
/* */
/* */
/* (voir l'image 'v $xiirq/HYDR.21.86') */
/* */
/* */
/* Author of '$xrq/hydrogene.51$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 1993??????????). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 20061210113449 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 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 xrq/hydrogene.58.I"
/* Introduit le 20030130112017 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 xrq/hydrogene.59.I"
/* Introduit le 20030130112017 a cause de 'v $xcc/cpp$Z INCLUDES_1.sed'. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D ' U N E S P H E R E : */
/* */
/* */
/* Definition : */
/* */
/* Elle est definie parametriquement */
/* en fonction des deux parametres 'u' */
/* et 'v' : */
/* */
/* F (u,v) = r.cos(u).cos(v) */
/* x */
/* */
/* F (u,v) = r.cos(u).sin(v) */
/* y */
/* */
/* F (u,v) = r.sin(u) */
/* z */
/* */
/* avec : */
/* */
/* p p */
/* u E [ - --- , + --- ] */
/* 2 2 */
/* */
/* v E [ 0 , 2.p ] */
/* */
/* (ou 'p' designe 'pi'). */
/* */
/* */
/*************************************************************************************************************************************/
#include xrs/sphere.12.I"
#include xrq/Hspherik.14.I"
/* Introduit le 20060611085033 suite a 'v $xrs/surfaces.12$I 20060602134331'. */
#include xrs/sphere.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"
/* A cause de 'v $xrs/surfaces.51$I definir_arbitrairement_les_coordonnees_u_et_v' ceci */
/* fut introduit le 20221224163321... */
#define PROFONDEUR_MINIMALE \
COORDONNEE_BARYCENTRIQUE_MINIMALE
#define PROFONDEUR_MAXIMALE \
COORDONNEE_BARYCENTRIQUE_MAXIMALE
#define PROFONDEUR_MINIMALE_DE_TRONCATION \
GRO1(FRA4(COORDONNEE_BARYCENTRIQUE_MAXIMALE))
#define PROFONDEUR_MAXIMALE_DE_TRONCATION \
COORDONNEE_BARYCENTRIQUE_MAXIMALE
/* Definition du "depth-cueing"... */
#define RAYON_DE_VISUALISATION \
GRO4(FRA1(FU)) \
/* Rayon du disque materialisant une iteration. */
#include xrq/hydrogene.51.I"
/* Rayon minimal du disque materialisant une iteration. ATTENTION : cette constante ne */
/* peut etre exprimee en fonction de 'RAYON_DE_VISUALISATION', car en effet, la premiere */
/* est plutot liee a l'image, alors que la seconde est plutot liee au modele... */
/* ATTENTION, a cause de la derivation formelle d'ordre superieur a 1, cette definition */
/* ne peut etre que dans un fichier d'include... */
#include xrq/Hspherik.15.I"
#include xrs/surfaces.11.I"
dfTRANSFORMAT_31(liste_PARAMETRE_R,fichier_PARAMETRE_R,PARAMETRE_R_IMPLICITE,GRO6(GRO10(PARAMETRE_R)))
/* Definition du fichier du parametre 'r'. */
#define sPARAMETRE_R(numero_de_la_periode) \
FLOT(sTRANSFORMAT_31(numero_de_la_periode,liste_PARAMETRE_R))
/* Formule generale definissant les variations de 'r' au cours du temps. */
#define COLORIER_AVEC_LA_DENSITE_DE_PROBABILITE \
FAUX
DEFV(Local,DEFV(Logical,INIT(colorier_avec_la_densite_de_probabilite,COLORIER_AVEC_LA_DENSITE_DE_PROBABILITE)));
/* Indique si le coloriage, lorsque l'on est en mode 'RVB' se fait a l'aide de la densite */
/* de probabilite ('VRAI') auquel cas le rayon des spheres est constant, ou bien a l'aide */
/* des differentielles sur la sphere ('FAUX') auquel cas le rayon des spheres est alors */
/* proportionnel a la densite de probabilite... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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/sphere.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... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E L ' A T O M E D ' H Y D R O G E N E : */
/* */
/*************************************************************************************************************************************/
#include xrq/hydrogene.25.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S C O U C H E S C O N C E N T R I Q U E S : */
/* */
/*************************************************************************************************************************************/
#define NOMBRE_DE_COUCHES_EXTERNES \
ZERO
DEFV(Local,DEFV(Int,INIT(nombre_de_couches_externes,NOMBRE_DE_COUCHES_EXTERNES)));
#define NOMBRE_DE_COUCHES_INTERNES \
ZERO
DEFV(Local,DEFV(Int,INIT(nombre_de_couches_internes,NOMBRE_DE_COUCHES_INTERNES)));
/* Definition du nombre de couches ; en realite le nombre de couches concentriques qui */
/* vont etre tracees est 'nombre_de_couches_externes + 1 + nombre_de_couches_internes'. */
#define DISTANCE_DES_COUCHES \
FRA10(PARAMETRE_R)
DEFV(Local,DEFV(Float,INIT(distance_des_couches,DISTANCE_DES_COUCHES)));
/* Distance separant deux couches concentriques... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S H A R M O N I Q U E S S P H E R I Q U E S : */
/* */
/*************************************************************************************************************************************/
#include xrq/Laguerre.31.I"
#include xrq/Legendre.31.I"
#include xrq/Hspherik.21.I"
#include xrq/hydrogene.21.I"
#include xrq/Hspherik.12.I"
#include xrq/hydrogene.24.I"
#define ORDRE_n \
ORDRE_n_01
DEFV(Local,DEFV(Int,INIT(ordre_n,ORDRE_n)));
/* Definition de l'ordre 'n'. */
#define ORDRE_l \
ORDRE_l_00
DEFV(Local,DEFV(Int,INIT(ordre_l,ORDRE_l)));
/* Definition de l'ordre 'l'. */
#define ORDRE_m \
ORDRE_m_00
DEFV(Local,DEFV(Int,INIT(ordre_m,ORDRE_m)));
/* Definition de l'ordre 'm'. */
DEFV(Local,DEFV(Float,INIT(rayon_de_visualisation_de_base,FLOT__UNDEF)));
/* Rayon de visualisation de base qui sera module par les fonctions propres... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 xrq/Hspherik.24.I"
/* Pour corriger les initialisations precedentes. ATTENTION, a cause de la derivation */
/* formelle d'ordre superieur a 1, ces 'undef's ne peuvent etre que dans un fichier */
/* d'include... */
#include xrk/attractor.18.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* C H O I X D E C E Q U I S E R A V I S U A L I S E : */
/* */
/*************************************************************************************************************************************/
#include xrv/champs_5.1B.I"
#define CALCULER_AUSSI_LA_PARTIE_IMAGINAIRE \
VRAI
DEFV(Local,DEFV(Logical,INIT(calculer_aussi_la_partie_Imaginaire,CALCULER_AUSSI_LA_PARTIE_IMAGINAIRE)));
/* Faut-il calculer la partie Imaginaire : */
/* */
/* VRAI : c'est le module de la fonction propre qui sera visualise, */
/* FAUX : c'est uniquement la partie Reelle de la fonction propre qui sera visualise. */
/* */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* P R O C E D U R E S U T I L E S : */
/* */
/*************************************************************************************************************************************/
#include xrq/hydrogene.23.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S F O N C T I O N S P R O P R E S D E L ' A T O M E D ' H Y D R O G E N E */
/* E T V I S U A L I S A T I O N S O U S F O R M E D E S P H E R E S T R I D I M E N S I O N N E L L E S : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
#include xrs/surfaces.41.I"
/* Introduit le 20161206103831... */
#include xrs/surfaces.22.I"
/* Introduit le 20221224163321... */
/*..............................................................................................................................*/
INITIALISATIONS_GENERALES;
/* Initialisations generales faites au tout debut... */
iTRANSFORMAT_31(liste_PARAMETRE_R,PARAMETRE_R_IMPLICITE);
/* Initialisation du parametre 'r'. */
#include xrv/champs_5.1A.I"
#include xrs/surfaces.13.I"
GET_ARGUMENTSv(nombre_d_arguments
,BLOC(PROCESS_ARGUMENTS_GEOMETRIQUES;
PROCESS_ARGUMENT_FICHIER("PARAMETRE_R="
,fichier_PARAMETRE_R
,liste_PARAMETRE_R
,PARAMETRE_R_IMPLICITE
,gTRANSFORMAT_31
);
PROCESS_ARGUMENTS_DE_VISUALISATION;
PROCESS_ARGUMENTS_DE_DEFINITION_DU_RAYON_DES_SPHERES_REPRESENTATIVES;
GET_ARGUMENT_I("ordre_n=""n=",ordre_n);
GET_ARGUMENT_I("ordre_l=""l=",ordre_l);
GET_ARGUMENT_I("ordre_m=""m=",ordre_m);
GET_ARGUMENT_L("spheriques=",calculer_la_probabilite_en_coordonnees_spheriques);
GET_ARGUMENT_I("nombre_couches_externes=""ne=""nce=",nombre_de_couches_externes);
GET_ARGUMENT_I("ni=""nci=",nombre_de_couches_internes);
GET_ARGUMENT_F("dc=""dR=",distance_des_couches);
GET_ARGUMENT_L("densite=",colorier_avec_la_densite_de_probabilite);
GET_ARGUMENT_L("imaginaire=",calculer_aussi_la_partie_Imaginaire);
PROCESS_ARGUMENTS_DE_DEFINITION_DES_COORDONNEES_CURVILIGNES;
)
);
VALIDATION_DES_ORDRES_lmn;
/* Validation des ordres demandes. */
EGAL(rayon_de_visualisation_de_base,rayon_de_visualisation);
/* Rayon de visualisation de base qui sera module par les fonctions propres. On notera */
/* que pour se simplifier la vie on prend pour 'rayon_de_visualisation_de_base' la valeur du */
/* parametre interactif 'rayon_de_visualisation'... */
#define NE_PAS_PERMETTRE_LA_DEFINITION_ARBITRAIRE_DES_COORDONNEES_U_ET_V
#include xrs/surfaces.12.I"
/* Visualisation et evolution de la surface au cours du temps... */
VISUALISATION_GENERALE_DE_LA_SURFACE
(BLOC(VIDE;)
,BLOC(Bblock
DEFV(Int,INIT(couches,UNDEF));
/* Pour numeroter les couches... */
DoIn(couches,NEGA(nombre_de_couches_internes),NEUT(nombre_de_couches_externes),I)
Bblock
EGAL(parametre_r
,AXPB(couches
,distance_des_couches
,sPARAMETRE_R(numero_de_la_periode_courante)
)
);
/* Calcul de la valeur du parametre 'r' en fonction de l'instant courant et de la couche */
/* courante... */
GENERATION_D_UNE_IMAGE_DE_LA_SURFACE_AVEC_DIFFERENTIATION
(Fx_Sphere_1,Fy_Sphere_1,Fz_Sphere_1
,BLOC(Bblock
DEFV(Float,INIT(densite_de_probabilite_de_presence,FLOT__UNDEF));
/* Definition de la densite de probabilite de presence... */
CALCUL_D_UNE_DENSITE_DE_PROBABILITE_DE_PRESENCE(ordre_n,ordre_l,ordre_m);
/* Calcul de la densite de probabilite de presence... */
Test(IFET(IL_FAUT(visualiser_en_RVB),IL_FAUT(colorier_avec_la_densite_de_probabilite)))
Bblock
EGAL(densite_de_probabilite_de_presence
,TRON(densite_de_probabilite_de_presence
,COORDONNEE_BARYCENTRIQUE_MINIMALE
,COORDONNEE_BARYCENTRIQUE_MAXIMALE
)
);
/* La densite de probabilite de presence va servir de coordonnee barycentrique... */
EGAL(dcx
,HOMO(BARY(FLOT(NIVR(NOIR))
,FLOT(NIVR(BLANC))
,densite_de_probabilite_de_presence
)
,FLOT(NIVR(NOIR)),FLOT(NIVR(BLANC))
,dXmin_ESPACE,dXmax_ESPACE
)
);
EGAL(dcy
,HOMO(FLOT(NIVR(GRIS_2))
,FLOT(NIVR(NOIR)),FLOT(NIVR(BLANC))
,dYmin_ESPACE,dYmax_ESPACE
)
);
EGAL(dcz
,HOMO(BARY(FLOT(NIVR(BLANC))
,FLOT(NIVR(NOIR))
,densite_de_probabilite_de_presence
)
,FLOT(NIVR(NOIR)),FLOT(NIVR(BLANC))
,dZmin_ESPACE,dZmax_ESPACE
)
);
/* Le choix des couleurs est fait alors de facon a ce que : */
/* */
/* probabilite=0 ==> BLEU, */
/* probabilite=1 ==> ROUGE. */
/* */
EGAL(rayon_de_visualisation
,rayon_de_visualisation_de_base
);
/* Le rayon de visualisation du point courant est constant lorsque le coloriage se fait */
/* avec la densite de probabilite de presence... */
Eblock
ATes
Bblock
EGAL(rayon_de_visualisation
,MUL2(RAYON_DES_SPHERES_REPRESENTATIVES(FU
,densite_de_probabilite_de_presence
)
,rayon_de_visualisation_de_base
)
);
/* Le rayon de visualisation du point courant est proportionnel a la densite de probabilite */
/* de presence elevee a une certaine puissance lorsque le coloriage ne se fait pas avec */
/* cette meme densite de probabilite de presence... */
Eblock
ETes
LIMITATION_DU_RAYON_DE_VISUALISATION;
/* Et on fait en sorte qu'il ne soit pas trop petit... */
Eblock
)
);
Eblock
EDoI
Eblock
)
,Fx_Sphere_1,Fy_Sphere_1,Fz_Sphere_1
/* Fonctions definissant parametriquement la surface... */
,BLOC(VIDE;)
/* Introduit le 20061210113135 car il manquait... */
);
EDITION_DES_DEBORDEMENTS_DU_RAYON_DE_VISUALISATION;
RETU_Commande;
Eblock
ECommande