Application


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

c00a.c
/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int rc)
{
double **A      =         r_mR(    i_mR(rc,rc),9);
double **Q      =                  i_mR(rc,rc);
double **R      =                  i_mR(rc,rc);	

double **A_T    = transpose_mR(A,  i_mR(rc,rc));
double **Q_T    =                  i_mR(rc,rc);
double **R_T    =                  i_mR(rc,rc);	

  clrscrn();
  printf(" A : ");
  p_mR(A,S4,P0,C4);
  QR_mR(A,Q,R);    
  printf(" Q :");
  p_mR(Q,S8,P3,C4);  
  printf(" R :");
  p_mR(R,S8,P3,C4);
  stop();
  
  clrscrn();
  printf(" A_T : ");
  p_mR(A_T,S4,P0,C4);
  QR_mR(A_T,Q_T,R_T);    
  printf(" Q_T :");
  p_mR(Q_T,S8,P3,C4);  
  printf(" R_T :");
  p_mR(R_T,S8,P3,C4);
  stop();

  clrscrn();
  printf(" Q :");
  p_mR(Q,S8,P3,C4);  
  printf(" R :");
  p_mR(R,S8,P3,C4);   
  printf(" Q_T :");
  p_mR(Q_T,S8,P3,C4);  
  printf(" R_T :");
  p_mR(R_T,S8,P3,C4);
  stop();
    
  clrscrn();
  printf(" A :");
  p_mR(A,S3,P0,C4); 
  printf(" A = Q * R  :");  
  mul_mR(Q,R, A);
  p_mR(A,S8,P3,C4);  

  printf(" A_T :");
  p_mR(A_T,S3,P0,C4); 
  printf(" A_T = Q_T * R_T  :");  
  mul_mR(Q_T,R_T, A_T);
  p_mR(A_T,S8,P3,C4); 
       
  f_mR(A);
  f_mR(Q);
  f_mR(R);
  f_mR(A_T);
  f_mR(Q_T);
  f_mR(R_T);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
  
   do 
         fun(rp_I(RC3)+RC1);
        
    while(stop_w());
        
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */


Exemple de sortie écran :

 A : 
  +8   +5   -9 
  +1   -8   -4 
  +5   +3   +8 

 Q :
  +0.843   +0.096   -0.529 
  +0.105   -0.994   -0.012 
  +0.527   +0.045   +0.849 

 R :
  +9.487   +4.954   -3.795 
  -0.000   +8.571   +3.477 
  +0.000   +0.000  +11.598 

 Press return to continue. 


 A_T : 
  +8   +1   +5 
  +5   -8   +3 
  -9   -4   +8 

 Q_T :
  +0.614   +0.090   +0.784 
  +0.383   -0.902   -0.196 
  -0.690   -0.421   +0.588 

 R_T :
 +13.038   +0.307   -1.304 
  +0.000   +8.995   -5.625 
  +0.000   +0.000   +8.041 

 Press return to continue. 


 Q :
  +0.843   +0.096   -0.529 
  +0.105   -0.994   -0.012 
  +0.527   +0.045   +0.849 

 R :
  +9.487   +4.954   -3.795 
  -0.000   +8.571   +3.477 
  +0.000   +0.000  +11.598 

 Q_T :
  +0.614   +0.090   +0.784 
  +0.383   -0.902   -0.196 
  -0.690   -0.421   +0.588 

 R_T :
 +13.038   +0.307   -1.304 
  +0.000   +8.995   -5.625 
  +0.000   +0.000   +8.041 

 Press return to continue. 


 A :
 +8  +5  -9 
 +1  -8  -4 
 +5  +3  +8 

 A = Q * R  :
  +8.000   +5.000   -9.000 
  +1.000   -8.000   -4.000 
  +5.000   +3.000   +8.000 

 A_T :
 +8  +1  +5 
 +5  -8  +3 
 -9  -4  +8 

 A_T = Q_T * R_T  :
  +8.000   +1.000   +5.000 
  +5.000   -8.000   +3.000 
  -9.000   -4.000   +8.000 


 Press   return to continue
 Press X return to stop