Aller au contenu

Mathc matrices/e1608

Un livre de Wikilivres.


Étude du code


Installer et compiler ces fichiers dans votre répertoire de travail.


gj08.c
/* ------------------------------------ */
/*  Save as :   gj08.c                 */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
void xx_gj3_T_mR(
double **Ab,
int below
)
{
int r= R1;
int c= C1;

   while( (r<Ab[R_SIZE][C0]) &&
           (c<Ab[C_SIZE_A][C0]) )
       {
        if(pivotbest_gj3Ab_mR(Ab,r,c)>ERROR_E) 

           zero_under_pivot_gj3Ab_mR(Ab,r,c);
 
        else r--;
        
        r++;
        c++;
       }

 r = Ab[R_SIZE][C0];
 
 if(below)
 
 while(r>R2)
  {
   r--;
   c=C1;
   
   while(fabs(Ab[r][c])<ERROR_E  &&  c<(Ab[C_SIZE_A][C0]-C1) ) c++;
         
         
   if(fabs(Ab[r][c])>ERROR_E)
   {
	  printf(" Ab : pivot[%d][[%d] ",r,c); 
      p_mR(Ab,S8,P4,C8);	 
       	   
      zero_below_pivot_gj3Ab_mR(Ab,r,c);
      
      printf(" Ab : pivot[%d][[%d]"
             " The work of  zero_below_pivot_gj1Ab_mR();",r,c); 
      p_mR(Ab,S8,P4,C8);	
      stop();
      clrscrn();          
    }
  }
}
/* ------------------------------------ */
double **XX_gj_PP_mR(
double **Ab,
int below
)
{
   xx_gj3_T_mR(Ab,below);
  
return(Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */
void fun(int r)
{
double **A  = i_mR(r,r);
double **b  = i_mR(r,C1);
double **Ab = i_Abr_Ac_bc_mR(r,r,C1);
/*            i_Abr_Ac_bc_mR(R_Ab,C_A,C_b));  */

 do
  {
   r_mR(A,999.);
   printf(".");
  }while(!det_R(A));

  r_mR(b,999.);
  c_A_b_Ab_mR(A,b,Ab);
  
  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mR(Ab,"Ab",P0);
  printf("\n rref(Ab,.00000000001)\n\n");
  stop();

  clrscrn();
  XX_gj_PP_mR(Ab,YES);
  printf("  gj_PP_mR(Ab,YES) :");
  p_mR(Ab,S10,P4,C7);
  
  
  f_mR(Ab);
  f_mR(b);
  f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
do
{
  fun(R5);


} while(stop_w());

  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */


Cette fois ci on met des zéros au dessus des pivots.


Exemple de sortie écran :
 ------------------------------------ 
 Copy/Past into the octave window.

 Ab=[
-807,-885,+368,+990,-643,+400;
-801,-229,-989,-441,-287,-243;
+52,-129,-983,-405,-223,-969;
-785,-561,+380,-969,-653,-633;
+996,+38,-905,+56,-815,+406]


 rref(Ab,.00000000001)

 Press return to continue. 



 Ab : pivot[5][[5] 
 +1.0000  +0.0382  -0.9086  +0.0562  -0.8183  +0.4076 
 +0.0000  +1.0000  +0.4276  -1.2121  +1.5258  -0.8534 
 +0.0000  -0.0000  +1.0000  +0.3900  +0.3920  +0.0526 
 +0.0000  -0.0000  -0.0000  +1.0000  +0.2904  +0.4981 
 +0.0000  +0.0000  +0.0000  -0.0000  +1.0000  -2.2005 

 Ab : The work of  zero_below_pivot_gj1Ab_mR();
 +1.0000  +0.0382  -0.9086  +0.0562  -0.0000  -1.3930 
 -0.0000  +1.0000  +0.4276  -1.2121  +0.0000  +2.5042 
 -0.0000  -0.0000  +1.0000  +0.3900  +0.0000  +0.9151 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  +1.1371 
 +0.0000  +0.0000  +0.0000  -0.0000  +1.0000  -2.2005 

 Press return to continue. 



 Ab : pivot[4][[4] 
 +1.0000  +0.0382  -0.9086  +0.0562  -0.0000  -1.3930 
 -0.0000  +1.0000  +0.4276  -1.2121  +0.0000  +2.5042 
 -0.0000  -0.0000  +1.0000  +0.3900  +0.0000  +0.9151 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  +1.1371 
 +0.0000  +0.0000  +0.0000  -0.0000  +1.0000  -2.2005 

 Ab : The work of  zero_below_pivot_gj1Ab_mR();
 +1.0000  +0.0382  -0.9086  +0.0000  -0.0000  -1.4569 
 -0.0000  +1.0000  +0.4276  -0.0000  +0.0000  +3.8824 
 -0.0000  +0.0000  +1.0000  +0.0000  +0.0000  +0.4716 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  +1.1371 
 +0.0000  +0.0000  +0.0000  -0.0000  +1.0000  -2.2005 

 Press return to continue. 



 Ab : pivot[3][[3] 
 +1.0000  +0.0382  -0.9086  +0.0000  -0.0000  -1.4569 
 -0.0000  +1.0000  +0.4276  -0.0000  +0.0000  +3.8824 
 -0.0000  +0.0000  +1.0000  +0.0000  +0.0000  +0.4716 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  +1.1371 
 +0.0000  +0.0000  +0.0000  -0.0000  +1.0000  -2.2005 

 Ab : The work of  zero_below_pivot_gj1Ab_mR();
 +1.0000  +0.0382  +0.0000  +0.0000  -0.0000  -1.0284 
 -0.0000  +1.0000  +0.0000  -0.0000  +0.0000  +3.6807 
 -0.0000  +0.0000  +1.0000  +0.0000  +0.0000  +0.4716 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  +1.1371 
 +0.0000  +0.0000  +0.0000  -0.0000  +1.0000  -2.2005 

 Press return to continue. 



 Ab : pivot[2][[2] 
 +1.0000  +0.0382  +0.0000  +0.0000  -0.0000  -1.0284 
 -0.0000  +1.0000  +0.0000  -0.0000  +0.0000  +3.6807 
 -0.0000  +0.0000  +1.0000  +0.0000  +0.0000  +0.4716 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  +1.1371 
 +0.0000  +0.0000  +0.0000  -0.0000  +1.0000  -2.2005 

 Ab : The work of  zero_below_pivot_gj1Ab_mR();
 +1.0000  +0.0000  +0.0000  +0.0000  -0.0000  -1.1688 
 -0.0000  +1.0000  +0.0000  -0.0000  +0.0000  +3.6807 
 -0.0000  +0.0000  +1.0000  +0.0000  +0.0000  +0.4716 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  +1.1371 
 +0.0000  +0.0000  +0.0000  -0.0000  +1.0000  -2.2005 

 Press return to continue. 



  gj_PP_mR(Ab,YES) :
  +1.00000   +0.00000   +0.00000   +0.00000   -0.00000   -1.16882 
  -0.00000   +1.00000   +0.00000   -0.00000   +0.00000   +3.68071 
  -0.00000   +0.00000   +1.00000   +0.00000   +0.00000   +0.47164 
  -0.00000   -0.00000   -0.00000   +1.00000   +0.00000   +1.13706 
  +0.00000   +0.00000   +0.00000   -0.00000   +1.00000   -2.20053 


 Press return to continue
 Press X      to stop