/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        R E S O L U T I O N   D E   S Y S T E M E S   L I N E A I R E S   B I D I M E N S I O N N E L S  :                         */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Definitions :                                                                                                              */
/*                                                                                                                                   */
/*                    Ce programme resoud des systemes                                                                               */
/*                  lineaires bidimensionnels :                                                                                      */
/*                                                                                                                                   */
/*                                      a11.X + a12.Y = a13                                                                          */
/*                                      a21.X + a22.Y = a23                                                                          */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrv/SysLin_2D.01$K' :                                                                                          */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 20030703092551).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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   SOLUTION_X_INDEFINIE                                                                                                          \
                    FZERO
#define   SOLUTION_Y_INDEFINIE                                                                                                          \
                    FZERO
                                        /* Valeurs a donner aux solutions {X,Y} dans le cas d'un systeme impossible.                 */
#define   PONDERATION_X_IMPLICITE                                                                                                       \
                    FZERO
#define   PONDERATION_Y_IMPLICITE                                                                                                       \
                    FZERO
                                        /* Ponderation de selection des solutions apres resolution.                                  */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D E S   F I C H I E R S  :                                                                           */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xrv/ARITHMET.1d.I"
                                        /* Passage a l'allocation dynamique le 20060214183112...                                     */
#include  xrv/ARITHMET.21.I"
#include  xrv/champs_5.41.I"

#define   A11_IMPLICITE                                                                                                                 \
                    FU
#define   A12_IMPLICITE                                                                                                                 \
                    FZERO
#define   A13_IMPLICITE                                                                                                                 \
                    FU
#define   A21_IMPLICITE                                                                                                                 \
                    FZERO
#define   A22_IMPLICITE                                                                                                                 \
                    FU
#define   A23_IMPLICITE                                                                                                                 \
                    FU

gGENERATION_D_UN_FICHIER(fichier_LISTE_A11,liste_initiale_des_A11);
gGENERATION_D_UN_FICHIER(fichier_LISTE_A12,liste_initiale_des_A12);
gGENERATION_D_UN_FICHIER(fichier_LISTE_A13,liste_initiale_des_A13);
gGENERATION_D_UN_FICHIER(fichier_LISTE_A21,liste_initiale_des_A21);
gGENERATION_D_UN_FICHIER(fichier_LISTE_A22,liste_initiale_des_A22);
gGENERATION_D_UN_FICHIER(fichier_LISTE_A23,liste_initiale_des_A23);
                                        /* Definition en memoire des fichiers des coefficients des systemes lineaires.               */

#define   ELEMENT_DU_FICHIER_LISTE_A11(index)                                                                                           \
                    gELEMENT_DU_FICHIER(liste_initiale_des_A11,index)
#define   ELEMENT_DU_FICHIER_LISTE_A12(index)                                                                                           \
                    gELEMENT_DU_FICHIER(liste_initiale_des_A12,index)
#define   ELEMENT_DU_FICHIER_LISTE_A13(index)                                                                                           \
                    gELEMENT_DU_FICHIER(liste_initiale_des_A13,index)
#define   ELEMENT_DU_FICHIER_LISTE_A21(index)                                                                                           \
                    gELEMENT_DU_FICHIER(liste_initiale_des_A21,index)
#define   ELEMENT_DU_FICHIER_LISTE_A22(index)                                                                                           \
                    gELEMENT_DU_FICHIER(liste_initiale_des_A22,index)
#define   ELEMENT_DU_FICHIER_LISTE_A23(index)                                                                                           \
                    gELEMENT_DU_FICHIER(liste_initiale_des_A23,index)
                                        /* Acces a un element courant des fichiers des coefficients des systemes lineaires.          */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        R E S O L U T I O N   D E   S Y S T E M E S   L I N E A I R E S   B I D I M E N S I O N N E L S  :                         */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
#include  xrv/ARITHMET.22.I"
#include  xci/valeurs.03.I"

     DEFV(Float,INIT(solution_X_indefinie,SOLUTION_X_INDEFINIE));
     DEFV(Float,INIT(solution_Y_indefinie,SOLUTION_Y_INDEFINIE));
                                        /* Valeurs a donner aux solutions {X,Y} dans le cas d'un systeme impossible.                 */
     DEFV(Float,INIT(ponderation_X,PONDERATION_X_IMPLICITE));
     DEFV(Float,INIT(ponderation_Y,PONDERATION_Y_IMPLICITE));
                                        /* Ponderation de selection des solutions apres resolution.                                  */
     /*..............................................................................................................................*/
#include  xrv/champs_5.1A.I"
                                        /* Ceci fut introduit le 20070103172937...                                                   */

     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("LISTE_A11="
                                                 ,fichier_LISTE_A11
                                                 ,liste_initiale_des_A11
                                                 ,A11_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                         PROKESF_ARGUMENT_FICHIER("LISTE_A12="
                                                 ,fichier_LISTE_A12
                                                 ,liste_initiale_des_A12
                                                 ,A12_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                         PROKESF_ARGUMENT_FICHIER("LISTE_A13="
                                                 ,fichier_LISTE_A13
                                                 ,liste_initiale_des_A13
                                                 ,A13_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                         PROKESF_ARGUMENT_FICHIER("LISTE_A21="
                                                 ,fichier_LISTE_A21
                                                 ,liste_initiale_des_A21
                                                 ,A21_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                         PROKESF_ARGUMENT_FICHIER("LISTE_A22="
                                                 ,fichier_LISTE_A22
                                                 ,liste_initiale_des_A22
                                                 ,A22_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );
                         PROKESF_ARGUMENT_FICHIER("LISTE_A23="
                                                 ,fichier_LISTE_A23
                                                 ,liste_initiale_des_A23
                                                 ,A23_IMPLICITE
                                                 ,lTRANSFORMAT_0d
                                                 ,iGENERATION_D_UN_FICHIER
                                                  );

                         GET_ARGUMENT_F("xi=""Xi=",solution_X_indefinie);
                         GET_ARGUMENT_F("yi=""Yi=",solution_Y_indefinie);

                         GET_ARGUMENT_F("px=""pX=""Pond1=",ponderation_X);
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_pasX'...                     */
                         GET_ARGUMENT_F("py=""pY=""Pond2=",ponderation_Y);
                                        /* A ne pas confondre avec 'v $xig/fonct$vv$DEF ENTREE_ARGUMENT_pasY'...                     */

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_3;
                                        /* Cette procedure fut introduite le 20070103172937...                                       */

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_1;
                                        /* Cette procedure fut introduite le 20061226193133...                                       */

                         PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_5;
                                        /* Cette procedure fut introduite le 20211005105203...                                       */
                         )
                    );

     gOPERATION_SUR_LES_FICHIERS(BLOC(
                                      DEFV(Float,INIT(coefficient_A11,ELEMENT_DU_FICHIER_LISTE_A11(index)));
                                      DEFV(Float,INIT(coefficient_A12,ELEMENT_DU_FICHIER_LISTE_A12(index)));
                                      DEFV(Float,INIT(coefficient_A13,ELEMENT_DU_FICHIER_LISTE_A13(index)));
                                      DEFV(Float,INIT(coefficient_A21,ELEMENT_DU_FICHIER_LISTE_A21(index)));
                                      DEFV(Float,INIT(coefficient_A22,ELEMENT_DU_FICHIER_LISTE_A22(index)));
                                      DEFV(Float,INIT(coefficient_A23,ELEMENT_DU_FICHIER_LISTE_A23(index)));
                                        /* Recuperation des coefficients {{a11,a12,a13},{a21,a22,a23}} des systemes lineaires.       */

                                      DEFV(Float,INIT(determinant,FLOT__UNDEF));
                                      DEFV(Float,INIT(determinant_X,FLOT__UNDEF));
                                      DEFV(Float,INIT(determinant_Y,FLOT__UNDEF));
                                        /* Definition des trois determinants necessaires.                                            */

                                      DEFV(Float,INIT(solution_X,solution_X_indefinie));
                                      DEFV(Float,INIT(solution_Y,solution_Y_indefinie));
                                        /* Solution {X,Y} du systeme courant initialisee sur le cas impossible...                    */

                                      EGAL(determinant
                                          ,DET2(coefficient_A11,coefficient_A12
                                               ,coefficient_A21,coefficient_A22
                                                )
                                           );
                                        /* Calcul du determinant "principal".                                                        */

                                      Test(IZNE(determinant))
                                           Bblock
                                        /* Cas ou le systeme peut etre resolu :                                                      */
                                           EGAL(determinant_X
                                               ,DET2(coefficient_A13,coefficient_A12
                                                    ,coefficient_A23,coefficient_A22
                                                     )
                                                );
                                           EGAL(determinant_Y
                                               ,DET2(coefficient_A11,coefficient_A13
                                                    ,coefficient_A21,coefficient_A23
                                                     )
                                                );
                                        /* Calcul des determinants "secondaires".                                                    */

                                           EGAL(solution_X,DIVI(determinant_X,determinant));
                                           EGAL(solution_Y,DIVI(determinant_Y,determinant));
                                        /* Resolution du systeme...                                                                  */
                                           Eblock
                                      ATes
                                           Bblock
                                        /* Cas ou le systeme ne peut etre resolu, les valeurs "indefinie"s sont renvoyees.           */
                                           Eblock
                                      ETes
                                      )
                                ,LIZ2(ponderation_X,solution_X
                                     ,ponderation_Y,solution_Y
                                      )
                                ,EDITER_LA_VALEUR_RESULTANTE_DANS_gOPERATION_SUR_LES_FICHIERS
                                ,nombre_d_exemplaires_du_resultat_de_l_operation_sur_les_valeurs_courantes
                                 );
                                        /* Edition des solutions {X,Y} du systeme lineaire courant.                                  */

     lGENERATION_D_UN_FICHIER(liste_initiale_des_A23,A23_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_initiale_des_A22,A22_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_initiale_des_A21,A21_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_initiale_des_A13,A13_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_initiale_des_A12,A12_IMPLICITE);
     lGENERATION_D_UN_FICHIER(liste_initiale_des_A11,A11_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.