/*************************************************************************************************************************************/
/* */
/* N O R M A L I S A T I O N D E C H A Q U E E L E M E N T D ' U N F I C H I E R N U M E R I Q U E */
/* P A R R A P P O R T A L A S O M M E D E S E L E M E N T S ( O U D E L E U R C A R R E ) */
/* D E C E L U I - C I : */
/* */
/* */
/* Author of '$xrv/normalise.11$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20041016212956). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 : */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* F I C H I E R S D ' I N C L U D E S : */
/* */
/*************************************************************************************************************************************/
#include INCLUDES_BASE
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* V A L E U R S I M P L I C I T E S D E S P A R A M E T R E S : */
/* */
/*************************************************************************************************************************************/
#define NORMALISER_AVEC_LA_SOMME_DES_CARRES \
FAUX \
/* Faut-il normaliser avec la somme des elements (('FAUX'), ou bien avec la somme des */ \
/* carres ('VRAI') ? Ceci a ete introduit le 20041103164418. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D U F I C H I E R : */
/* */
/*************************************************************************************************************************************/
#include xrv/ARITHMET.1d.I"
/* Passage a l'allocation dynamique le 20060214192645... */
#define NOMBRE_DE_DECIMALES \
SEIZE
#include xrv/ARITHMET.21.I"
#include xrv/champs_5.41.I"
#define VALEUR_IMPLICITE \
FZERO
gGENERATION_D_UN_FICHIER(fichier_des_valeurs,liste_des_valeurs);
/* Definition en memoire du fichier. */
#define ELEMENT_DU_FICHIER(index) \
gELEMENT_DU_FICHIER(liste_des_valeurs,index) \
/* Acces a un element courant du fichier. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S E D I T I O N S : */
/* */
/*************************************************************************************************************************************/
#define FORMAT_EFFECTIF_D_EDITION(message) \
Cara(chain_Aconcaten5(message,INTRODUCTION_FORMAT,valeurs_signees,".*",format_d_edition)) \
/* Format effectif de certaines editions... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* N O R M A L I S A T I O N D E C H A Q U E E L E M E N T D ' U N F I C H I E R N U M E R I Q U E */
/* P A R R A P P O R T A L A S O M M E D E S E L E M E N T S ( O U D E L E U R C A R R E ) */
/* D E C E L U I - C I : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
DEFV(Logical,INIT(normaliser_avec_la_somme_des_carres,NORMALISER_AVEC_LA_SOMME_DES_CARRES));
/* Faut-il normaliser avec la somme des elements (('FAUX'), ou bien avec la somme des */
/* carres ('VRAI') ? Ceci a ete introduit le 20041103164418. */
DEFV(Float,INIT(cumul_de_l_ensemble_des_puissances_1_des_elements_du_fichier,FZERO));
DEFV(Float,INIT(cumul_de_l_ensemble_des_puissances_2_des_elements_du_fichier,FZERO));
/* Somme des elements du fichier et de leurs carres. */
#include xrv/ARITHMET.22.I"
#include xci/valeurs.03.I"
/*..............................................................................................................................*/
#include xrv/champs_5.1A.I"
/* Ceci fut introduit le 20070103174116... */
GET_ARGUMENTS_(nombre_d_arguments
,BLOC(PROCESS_ARGUMENT_I("nombre_elements=""ne=",nombre_d_elements
,BLOC(VIDE;)
,BLOC(Bblock
PRINT_AVERTISSEMENT("'ne=' doit etre defini avant toute entree de fichiers");
Eblock
)
);
PROCESS_ARGUMENTS_DE_DEFINITION_DES_FICHIERS_01;
PROKESF_ARGUMENT_FICHIER("fichier="
,fichier_des_valeurs
,liste_des_valeurs
,VALEUR_IMPLICITE
,lTRANSFORMAT_0d
,iGENERATION_D_UN_FICHIER
);
GET_ARGUMENT_L("p2=""carres=",normaliser_avec_la_somme_des_carres);
GET_ARGUMENT_N("p1=""lineaire=",normaliser_avec_la_somme_des_carres);
PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_3;
/* Cette procedure fut introduite le 20070103174116... */
PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_1;
/* Cette procedure fut introduite le 20061226193412... */
PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_5;
/* Cette procedure fut introduite le 20211005105755... */
)
);
DoIn(index
,PREMIER_ELEMENT_D_UN_FICHIER
,DERNIER_ELEMENT_D_UN_FICHIER
,I
)
Bblock
DEFV(Float,INIT(valeur_courante,ELEMENT_DU_FICHIER(index)));
/* Recuperation de la valeur courante dans le fichier. */
Test(IL_FAUT(normaliser_avec_la_somme_des_carres))
Bblock
INCR(cumul_de_l_ensemble_des_puissances_2_des_elements_du_fichier,EXP2(valeur_courante));
/* Cumul de l'ensemble des carres des elements du fichier... */
Eblock
ATes
Bblock
INCR(cumul_de_l_ensemble_des_puissances_1_des_elements_du_fichier,EXP1(valeur_courante));
/* Cumul de l'ensemble des elements du fichier... */
Eblock
ETes
Eblock
EDoI
gOPERATION_SUR_LES_FICHIERS(BLOC(
DEFV(Float,INIT(valeur_courante,ELEMENT_DU_FICHIER(index)));
/* Recuperation de la valeur courante dans le fichier. */
)
,DIVZ(valeur_courante
,COND(IL_FAUT(normaliser_avec_la_somme_des_carres)
,cumul_de_l_ensemble_des_puissances_2_des_elements_du_fichier
,cumul_de_l_ensemble_des_puissances_1_des_elements_du_fichier
)
)
,EDITER_LA_VALEUR_RESULTANTE_DANS_gOPERATION_SUR_LES_FICHIERS
,nombre_d_exemplaires_du_resultat_de_l_operation_sur_les_valeurs_courantes
);
/* Normalisation du fichier argument. Cela permet de generer un fichier dont la somme */
/* des elements est egale a 1 et donc dont chaque element est dans [0,1] ; le fichier */
/* resultant est donc un ensemble de probabilites. On notera alors que l'on suppose donc */
/* plus ou moins que les elements du fichier sont positifs ou nuls et qu'ainsi on les ramene */
/* effectivement dans [0,1]. On pourrait imaginer que l'on translate prealablement chaque */
/* element par le minimum du fichier ; mais cela n'est pas realise car, en effet, dans le */
/* cas contraire, cela introduirait au moins un element nul (celui -ou ceux- qui */
/* correspondait au minimum) dans le fichier et cela serait ennuyeux dans le cas ou le */
/* fichier resultant serait transmis a 'v $xrv/entropie.01$K' qui prend des logarithmes... */
/* */
/* On notera que l'option 'normaliser_avec_la_somme_des_carres' introduite le 20041103164418 */
/* permet de generer deux types de vecteur : */
/* */
/* p1=VRAI : ensemble de probabilites (de somme egale a 1), */
/* p2=VRAI : composantes d'un vecteur normalise (norme egale a 1). */
/* */
lGENERATION_D_UN_FICHIER(liste_des_valeurs,VALEUR_IMPLICITE);
RETU_Commande;
Eblock
ECommande
Copyright © Jean-François COLONNA, 2019-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2019-2024.