Linear Algebra and the C Language/a0cg

/* ------------------------------------ */
/*  Save as :   c00b.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define RCA      RC4
/* ------------------------------------ */
void fun(void)
{
double **A = rEsymmetric_mR(i_mR(RCA,RCA),999,+1.E-3);  
     
double **LogA_EValue = i_mR(RCA,RCA);
double **LogA        = i_mR(RCA,RCA);
     
double **Exp_LogA_EValue = i_mR(RCA,RCA);
double **Exp_LogA        = i_mR(RCA,RCA);

double **V       = i_mR(RCA,RCA);
double **V_T     = i_mR(RCA,RCA);
double **EValue  = i_mR(RCA,RCA); 
double **T       = i_mR(RCA,RCA);

/* ================================================================== */           
   clrscrn();
   printf(" Copy/Paste into the octave window.       \n\n");
   p_Octave_mR(A,"A", P4);
   printf(" logm (A)\n\n\n");
  
   eigs_V_mR(A,V);
   transpose_mR(V,V_T);
    
/* EValue = V_T * A * V */   
   mul_mR(V_T,A,T);
   mul_mR(T,V,EValue); 
    
   printf(" LogA");
   f_eigs_mR(log,EValue, LogA_EValue);     
      mul_mR(V,LogA_EValue,T);
      mul_mR(T,V_T,LogA); 
        p_mR(LogA, S9,P4, C6);
   stop(); 
   
/* ================================================================== */   
   clrscrn();     
   printf(" Copy/Paste into the octave window.       \n\n");
   p_Octave_mR(LogA, "LogA", P4);
   printf(" expm (LogA)\n\n\n");
  
   eigs_V_mR(LogA,V);
   transpose_mR(V,V_T);   
   
/* LogA_EValue = V_T * LogA * V */   
   mul_mR(V_T,LogA,T);
   mul_mR(T,V,LogA_EValue); 
   
   printf(" Exp_LogA = A");
   f_eigs_mR(exp,LogA_EValue, Exp_LogA_EValue);    
      mul_mR(V,Exp_LogA_EValue,T);
      mul_mR(T,V_T,Exp_LogA);   
   p_mR(Exp_LogA, S9,P4, C6);

   printf(" A :");
   p_mR(A, S9,P4, C6);  
          
   f_mR(A);
   
   f_mR(LogA_EValue);
   f_mR(LogA);
      
   f_mR(V);
   f_mR(V_T);
   f_mR(EValue); 
   f_mR(T);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
    fun();
    
} while(stop_w());

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

                                                                                       
 Copy/Paste into the octave window.       

 A=[
+0.9278,+0.0302,-0.1843,+0.3124;
+0.0302,+0.8712,+0.4961,+0.4018;
-0.1843,+0.4961,+0.9266,+0.5917;
+0.3124,+0.4018,+0.5917,+0.7404]

 logm (A)


 LogA
  -0.3164   +0.0339   -0.5592   +0.7361 
  +0.0339   -0.3705   +0.5283   +0.3694 
  -0.5592   +0.5283   -0.7673   +1.1333 
  +0.7361   +0.3694   +1.1333   -1.1409 

 Press return to continue. 


 Copy/Paste into the octave window.       

 LogA=[
-0.3164,+0.0339,-0.5592,+0.7361;
+0.0339,-0.3705,+0.5283,+0.3694;
-0.5592,+0.5283,-0.7673,+1.1333;
+0.7361,+0.3694,+1.1333,-1.1409]

 expm (LogA)


 Exp_LogA = A
  +0.9278   +0.0302   -0.1843   +0.3124 
  +0.0302   +0.8712   +0.4961   +0.4018 
  -0.1843   +0.4961   +0.9266   +0.5917 
  +0.3124   +0.4018   +0.5917   +0.7404 

 A :
  +0.9278   +0.0302   -0.1843   +0.3124 
  +0.0302   +0.8712   +0.4961   +0.4018 
  -0.1843   +0.4961   +0.9266   +0.5917 
  +0.3124   +0.4018   +0.5917   +0.7404 


 Press   return to continue
 Press X return to stop