#include <stdio.h>
#include "proc/set_mtrx.proc"
FILE *fout;

main() {
  void set_p_mtrx(), set_q_mtrx(), set_ca_mtrx(), set_cu_mtrx(),
       set_cx_mtrx(), init_to_zero(), fpr_mtrx();
  int i,j;
  float p[3][8][8], q[2][8][8], ca[3][8], cu[8], cx[2][8];

  fout=fopen("/mit/smmadana/mtrx.dat","w");
  for(i=0;i<6;i++) {
    init_to_zero(p,3);
      set_p_mtrx(p,i);
    init_to_zero(q,2);
      set_q_mtrx(q,i);
    fprintf(fout,"                                    P=%d\n\n",i);
    for(j=0;j<3;j++) {
      fprintf(fout,"p[%d]\n",j);
      fpr_mtrx(p[j],8,8);
    }
    fprintf(fout,"\n\n\n                                    Q=%d\n\n",i);
    for(j=0;j<2;j++) {
      fprintf(fout,"q[%d]\n",j);
      fpr_mtrx(q[j],8,8);
    }
    fprintf(fout,"\f\n");
  }
  fprintf(fout,"\n\n                                      CA\n\n",i);
  set_ca_mtrx(ca);
    fpr_mtrx(ca,3,8);
  fprintf(fout,"\n\n                                      CU\n\n",i);
  set_cu_mtrx(cu);
    fpr_mtrx(cu,1,8);
  fprintf(fout,"\n\n                                      CX\n\n",i);
  set_cx_mtrx(cx);
    fpr_mtrx(cx,2,8);
  fclose(fout);
}  
/*****************************************************************************/
void init_to_zero(a,i)
short i;
float a[][8][8];{
  short loop1, loop2, loop3;

  for(loop1=0;loop1<i;loop1++)  
    for(loop2=0;loop2<8;loop2++)
	for(loop3=0;loop3<8;loop3++)
	  a[loop1][loop2][loop3]=0.0;
}
/*****************************************************************************/
void fpr_mtrx (mtrx,x,y)
float mtrx [8][8]; 
int x,y; {
  int i,j;
 
  for(i=x-1;i>(-1);i--) {
    for(j=y-1;j>(-1);j--)
      fprintf(fout,"%10.3f",mtrx[i][j]);
   fprintf(fout,"\n");
  }
}