Linear Algebra and the C Language/a0gd

/* ------------------------------------ */
/*  Save as :   c00b.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void AplsB(int r)
{
double **A  = i_mR(r,r);
double **B  = i_mR(r,r);
double **C  = i_mR(r,r);
double **D  = i_mR(r,r);
double **T  = i_mR(r,r);
double **AB = i_mR(r,r);
double **BA = i_mR(r,r);

  r_commute4_mR(A,B,C,D,9);

  clrscrn();
  printf(" A:");
  p_mR(A, S5,P3,C6); 
                                 
  printf(" B:");
  p_mR(B, S5,P3,C6);
  
  printf(" C:");
  p_mR(C, S5,P3,C6); 
  
  printf(" D:");
  p_mR(D, S5,P3,C6);        
  stop();

  clrscrn();
  printf(" AB:");
  p_mR(mul_mR(A,B,AB), S5,P9,C6);
  printf(" BA:");
  p_mR(mul_mR(B,A,BA), S5,P9,C6);
  stop();

  clrscrn();
  printf(" AC:");
  p_mR(mul_mR(A,C,AB), S5,P9,C6);
  printf(" CA:");
  p_mR(mul_mR(C,A,BA), S5,P9,C6);
  stop();

  clrscrn();
  printf(" BC:");
  p_mR(mul_mR(B,C,AB), S5,P9,C6);
  printf(" CB:");
  p_mR(mul_mR(C,B,BA), S5,P9,C6);   
  stop();

  clrscrn();
  printf(" AD:");
  p_mR(mul_mR(A,D,AB), S5,P9,C6);
  printf(" DA:");
  p_mR(mul_mR(D,A,BA), S5,P9,C6);
  stop();

  clrscrn();
  printf(" BD:");
  p_mR(mul_mR(B,D,AB), S5,P9,C6);
  printf(" DB:");
  p_mR(mul_mR(D,B,BA), S5,P9,C6); 
  stop();

  clrscrn();
  printf(" CD:");
  p_mR(mul_mR(C,D,AB), S5,P9,C6);
  printf(" DC:");
  p_mR(mul_mR(D,C,BA), S5,P9,C6); 
  
  f_mR(A);
  f_mR(B);
  f_mR(C);
  f_mR(D);
  f_mR(T);
  f_mR(AB);
  f_mR(BA);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
  
do
{
  AplsB(RC3);

} while(stop_w());

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

                                                                                       
 A:
-0.455 -0.154 +0.004 
-0.050 -0.619 -0.001 
+0.227 -0.173 -0.652 

 B:
-0.616 +1.046 -0.029 
+0.343 +0.492 +0.007 
-1.542 +1.177 +0.721 

 C:
-0.005 -0.010 +0.000 
-0.003 -0.016 -0.000 
+0.015 -0.011 -0.018 

 D:
+0.005 +0.010 -0.000 
+0.003 +0.016 +0.000 
-0.015 +0.011 +0.018 

 Press return to continue. 


 AB:
+0.221370207 -0.547138127 +0.015412342 
-0.179168471 -0.358519146 -0.003853085 
+0.806258121 -0.615530393 -0.477964794 

 BA:
+0.221370207 -0.547138127 +0.015412342 
-0.179168471 -0.358519146 -0.003853085 
+0.806258121 -0.615530393 -0.477964794 

 Press return to continue. 


 AC:
+0.002957845 +0.006889328 -0.000194066 
+0.002256012 +0.010259563 +0.000048516 
-0.010152056 +0.007750494 +0.011763571 

 CA:
+0.002957845 +0.006889328 -0.000194066 
+0.002256012 +0.010259563 +0.000048516 
-0.010152056 +0.007750494 +0.011763571 

 Press return to continue. 


 BC:
-0.000582702 -0.010052002 +0.000283155 
-0.003291677 -0.011236408 -0.000070789 
+0.014812545 -0.011308502 -0.013430859 

 CB:
-0.000582702 -0.010052002 +0.000283155 
-0.003291677 -0.011236408 -0.000070789 
+0.014812545 -0.011308502 -0.013430859 

 Press return to continue. 
 

 AD:
-0.002957842 -0.006889316 +0.000194065 
-0.002256008 -0.010259546 -0.000048516 
+0.010152037 -0.007750480 -0.011763551 

 DA:
-0.002957842 -0.006889316 +0.000194065 
-0.002256008 -0.010259546 -0.000048516 
+0.010152037 -0.007750480 -0.011763551 

 Press return to continue. 


 BD:
+0.000582698 +0.010051988 -0.000283155 
+0.003291672 +0.011236389 +0.000070789 
-0.014812524 +0.011308486 +0.013430837 

 DB:
+0.000582698 +0.010051988 -0.000283155 
+0.003291672 +0.011236389 +0.000070789 
-0.014812524 +0.011308486 +0.013430837 

 Press return to continue. 


 CD:
-0.000063883 -0.000205158 +0.000005779 
-0.000067182 -0.000281322 -0.000001445 
+0.000302319 -0.000230803 -0.000326110 

 DC:
-0.000063883 -0.000205158 +0.000005779 
-0.000067182 -0.000281322 -0.000001445 
+0.000302319 -0.000230803 -0.000326110 


 Press   return to continue
 Press X return to stop