Aller au contenu

Mathc matrices/e1508

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"
/* ------------------------------------ */
/* ------------------------------------ */
double **xx_gj1_T_mR(
double **Ab
)
{
int r= R1;

    while(r<Ab[R_SIZE][C0])
       {        	   
        if(fabs(pivotbest_gj1Ab_mR(Ab,r))>ERROR_E)  
 
          zero_under_pivot_gj1Ab_mR(Ab,r);
                    
        r++;
       } 
                
    while(r>R2)
    {
		printf(" Ab : pivot[%d][[%d] \n",r-R1,r-R1); /* See : --r */
        p_mR(Ab,S8,P4,C8);	 
           
        zero_below_pivot_gj1Ab_mR(Ab,--r);
        
        printf(" Ab : pivot[%d][[%d]"
               " The work of  zero_below_pivot_gj1Ab_mR();\n",
                 r,r);
                 
        p_mR(Ab,S8,P4,C8);	
        stop();
        clrscrn();    
     }
        
return(Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */
double **XX_gj_TP_mR(
double **Ab
)
{
   xx_gj1_T_mR(Ab); 
    
   sort_c_mR(Ab);  
   sort_r_mR(Ab);
  
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_TP_mR(Ab);
  printf("  gj_TP_mR(Ab) :\n");
  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=[
-943,-389,+384,-897,+264,+938;
+938,+794,+536,-605,-917,+342;
+538,-905,-235,-921,-971,+296;
+20,+400,+758,-543,-883,+820;
-735,+514,-669,+966,-539,-499]


 rref(Ab,.00000000001)

 Press return to continue. 


 Ab : pivot[5][[5] 
 +1.0000  +0.9320  +0.9485  -0.5541  +0.2420  -0.3048 
 +0.0000  +1.0000  +0.1606  +0.2608  +0.4597  +0.0379 
 +0.0000  +0.0000  +1.0000  -0.9883  -0.7654  -0.5341 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.1130  -0.2909 
 +0.0000  +0.0000  +0.0000  +0.0000  +1.0000  +0.6323 

 Ab : The work of  zero_below_pivot_gj1Ab_mR();
 +1.0000  +0.9320  +0.9485  -0.5541  +0.0000  -0.4579 
 +0.0000  +1.0000  +0.1606  +0.2608  +0.0000  -0.2528 
 +0.0000  +0.0000  +1.0000  -0.9883  +0.0000  -0.0501 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  -0.3624 
 +0.0000  +0.0000  +0.0000  +0.0000  +1.0000  +0.6323 

 Press return to continue. 



 Ab : pivot[4][[4] 
 +1.0000  +0.9320  +0.9485  -0.5541  +0.0000  -0.4579 
 +0.0000  +1.0000  +0.1606  +0.2608  +0.0000  -0.2528 
 +0.0000  +0.0000  +1.0000  -0.9883  +0.0000  -0.0501 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  -0.3624 
 +0.0000  +0.0000  +0.0000  +0.0000  +1.0000  +0.6323 

 Ab : The work of  zero_below_pivot_gj1Ab_mR();
 +1.0000  +0.9320  +0.9485  +0.0000  +0.0000  -0.6587 
 +0.0000  +1.0000  +0.1606  +0.0000  +0.0000  -0.1583 
 +0.0000  +0.0000  +1.0000  +0.0000  +0.0000  -0.4083 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  -0.3624 
 +0.0000  +0.0000  +0.0000  +0.0000  +1.0000  +0.6323 

 Press return to continue. 



 Ab : pivot[3][[3] 
 +1.0000  +0.9320  +0.9485  +0.0000  +0.0000  -0.6587 
 +0.0000  +1.0000  +0.1606  +0.0000  +0.0000  -0.1583 
 +0.0000  +0.0000  +1.0000  +0.0000  +0.0000  -0.4083 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  -0.3624 
 +0.0000  +0.0000  +0.0000  +0.0000  +1.0000  +0.6323 

 Ab : The work of  zero_below_pivot_gj1Ab_mR();
 +1.0000  +0.9320  +0.0000  +0.0000  +0.0000  -0.2714 
 +0.0000  +1.0000  +0.0000  +0.0000  +0.0000  -0.0927 
 +0.0000  +0.0000  +1.0000  +0.0000  +0.0000  -0.4083 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  -0.3624 
 +0.0000  +0.0000  +0.0000  +0.0000  +1.0000  +0.6323 

 Press return to continue. 


 Ab : pivot[2][[2] 
 +1.0000  +0.9320  +0.0000  +0.0000  +0.0000  -0.2714 
 +0.0000  +1.0000  +0.0000  +0.0000  +0.0000  -0.0927 
 +0.0000  +0.0000  +1.0000  +0.0000  +0.0000  -0.4083 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  -0.3624 
 +0.0000  +0.0000  +0.0000  +0.0000  +1.0000  +0.6323 

 Ab : The work of  zero_below_pivot_gj1Ab_mR();
 +1.0000  +0.0000  +0.0000  +0.0000  +0.0000  -0.1850 
 +0.0000  +1.0000  +0.0000  +0.0000  +0.0000  -0.0927 
 +0.0000  +0.0000  +1.0000  +0.0000  +0.0000  -0.4083 
 -0.0000  -0.0000  -0.0000  +1.0000  +0.0000  -0.3624 
 +0.0000  +0.0000  +0.0000  +0.0000  +1.0000  +0.6323 

 Press return to continue. 


  gj_TP_mR(Ab) :
  +1.00000   -0.00000   +0.00000   -0.00000   -0.00000   -0.36240 
  +0.00000   +1.00000   +0.00000   +0.00000   +0.00000   -0.09275 
  +0.00000   +0.00000   +1.00000   +0.00000   +0.00000   +0.63234 
  +0.00000   +0.00000   +0.00000   +1.00000   +0.00000   -0.40825 
  +0.00000   +0.00000   +0.00000   +0.00000   +1.00000   -0.18500 


 Press return to continue
 Press X      to stop