/*************************************************************************************************************************************/ /* */ /* 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