Linear Algebra and the C Language/a0gl

/* ------------------------------------ */
/*  Save as :  c00e.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define  ARRAY  A2 
#define    RCA  R5  
/* ------------------------------------ */
void fun(void)
{
double **A[ARRAY];
double **InvA[ARRAY];

double **P            =   r_mR(   i_mR(RCA,RCA),99); 
double **InvP         = inv_mR(P, i_mR(RCA,RCA)   );
double **T            =           i_mR(RCA,RCA);
double **InvP_InvA0_P =           i_mR(RCA,RCA);
int i;

  for(i=A0; i<ARRAY; i++)
     {   
          A[i] = i_mR(RCA,RCA); 
       InvA[i] = i_mR(RCA,RCA);   
	   }
	    
  rsymmetric_mR(A[R0],99);
  
  mul_mR(InvP,A[A0],T); 
  mul_mR(T,P,A[A0+A1]);

  clrscrn();
  printf(" The two similar matrices :\n\n"
         " A[%d] = InvP A[%d] P   \n\n",A1,A0);
  for(i=A0; i<ARRAY; i++)
     {
	  printf(" A[%d] : ",i);      
      p_mR(A[i], S12,P4,C6); 
     }   
  stop();    

   clrscrn();
   
   for(i=A0; i<ARRAY; i++)
      {
       inv_mR(A[i],InvA[i]); 
       }       
          
   mul_mR(InvP,InvA[A0],T); 
   mul_mR(T,P,InvP_InvA0_P);  
     
   printf(" If   A[%d] = InvP A[%d] P    Then   "
          " InVA[%d] = InvP InVA[%d] P \n\n",A1,A0 ,A1,A0); 
   printf(" InVA[%d] : ",A1);      
   pE_mR(InvA[1], S12,P4,C6); 
   
   printf(" InvP InVA[%d] P : ",A0);      
   pE_mR(InvP_InvA0_P, S12,P4,C6);  

  for(i=A0; i<ARRAY; i++)
     { 
	   f_mR(A[i]);
	   f_mR(InvA[i]);
	   }
	   
  f_mR(P);		 
  f_mR(InvP); 	   
  f_mR(T); 
  f_mR(InvP_InvA0_P);       
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
  
  do{
        fun();
        
  }while(stop_w());

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

                                                                                       
 The two similar matrices :

 A[1] = InvP A[0] P   

 A[0] : 
    +82.0000     +65.0000     +86.0000     -99.0000     +15.0000 
    +65.0000     -75.0000     -85.0000     -42.0000     -22.0000 
    +86.0000     -85.0000      -1.0000     +75.0000      +8.0000 
    -99.0000     -42.0000     +75.0000     +31.0000     -97.0000 
    +15.0000     -22.0000      +8.0000     -97.0000     -40.0000 

 A[1] : 
     +5.4781     +82.9444    -180.5620    +212.9966    +218.8392 
   +230.4664    -359.8709    -788.2735    +917.4393    +239.7937 
    +55.2520    -231.7830     -58.2085     +50.3652    -178.5631 
   +234.1154    -518.1748    -458.9903    +584.3454    -141.5204 
   -200.7368    +430.8753    +760.6541    -949.9923    -174.7442 

 Press return to continue. 


 If   A[1] = InvP A[0] P    Then    InVA[1] = InvP InVA[0] P 

 InVA[1] : 
 +2.0339e-02  +4.3383e-04  -3.8811e-03  +1.7726e-02  +1.5677e-02 
 -1.4125e-01  +4.3604e-02  -7.0015e-02  -3.6664e-02  -1.5825e-02 
 +7.1233e-02  -3.2594e-02  +1.6698e-02  +2.9501e-02  +3.5259e-03 
 -3.9873e-02  +1.9200e-03  -3.0451e-02  -4.5130e-03  -1.2529e-02 
 +1.5518e-01  -4.5299e-02  +7.0049e-02  +4.2186e-02  +2.0709e-02 

 InvP InVA[0] P : 
 +2.0339e-02  +4.3383e-04  -3.8811e-03  +1.7726e-02  +1.5677e-02 
 -1.4125e-01  +4.3604e-02  -7.0015e-02  -3.6664e-02  -1.5825e-02 
 +7.1233e-02  -3.2594e-02  +1.6698e-02  +2.9501e-02  +3.5259e-03 
 -3.9873e-02  +1.9200e-03  -3.0451e-02  -4.5130e-03  -1.2529e-02 
 +1.5518e-01  -4.5299e-02  +7.0049e-02  +4.2186e-02  +2.0709e-02 


 Press   return to continue
 Press X return to stop