/* ------------------------------------ */
/*  Save as :   cq3.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **U   = r_mR(i_mR(r,r),9.);
double **V   = r_mR(i_mR(r,r),9.);
double **A   =      i_mR(r,r);
double **B   =      i_mR(r,r);
double **AB  =      i_mR(r,r);
double **ABT =      i_mR(r,r);
double **T   =      i_mR(r,r);

  clrscrn();
  
  r_Quv_mR(U,A);    
  printf(" A :                an orthogonal matrix ");
  p_mR(A,3,3,6);
  
  r_Quv_mR(V,B);    
  printf(" B :                an orthogonal matrix ");
  p_mR(B,3,3,6);
   
  printf("   A * B  = AB  :   AB is an orthogonal matrix");
  mul_mR(A,B,AB);
  p_mR(AB,3,3,6);    

  printf(" AB * ABT = id :    inv(AB) = ABT ");
  transpose_mR(AB,ABT);
  mul_mR(AB,ABT,T);
  p_mR(T,3,3,6);    
  
  f_mR(U);
  f_mR(V);
  f_mR(A);
  f_mR(B);
  f_mR(AB);
  f_mR(ABT);
  f_mR(T);

}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun(rp_I(3)+R2);

} while(stop_w());

  return 0;
}

/* ------------------------------------ */
/* ------------------------------------ */