===================================================================
RCS file: /home/cvs/OpenXM/src/kan96xx/Kan/matrix.c,v
retrieving revision 1.1.1.1
retrieving revision 1.4
diff -u -p -r1.1.1.1 -r1.4
--- OpenXM/src/kan96xx/Kan/matrix.c	1999/10/08 02:12:02	1.1.1.1
+++ OpenXM/src/kan96xx/Kan/matrix.c	2005/07/03 11:08:54	1.4
@@ -1,16 +1,18 @@
+/* $OpenXM: OpenXM/src/kan96xx/Kan/matrix.c,v 1.3 2001/05/04 01:06:24 takayama Exp $ */
 /* matrix.c */
 #include <stdio.h>
+#include <stdlib.h>
 #include "datatype.h"
 #include "stackm.h"
 #include "extern.h"
 #include "extern2.h"
 
 #include "matrix.h"   /* ind() and ind2() are defined. */
-                      /* You must use n and n2. */
+/* You must use n and n2. */
 
 struct arrayOfPOLY *aaAdd(aa,bb)
-struct arrayOfPOLY *aa,*bb;
-/* this function has not been tested yet. */
+     struct arrayOfPOLY *aa,*bb;
+     /* this function has not been tested yet. */
 {
   int i;
   POLY *r;
@@ -32,7 +34,7 @@ struct arrayOfPOLY *aa,*bb;
 }
 
 struct matrixOfPOLY *aaMult(aa,bb)
-struct matrixOfPOLY *aa,*bb;
+     struct matrixOfPOLY *aa,*bb;
 {
   POLY *a;
   POLY *b;
@@ -53,7 +55,7 @@ struct matrixOfPOLY *aa,*bb;
     for (j=0; j<n2; j++) {
       tmp = ZERO;
       for (k=0; k<n; k++) {
-	tmp = ppAddv(tmp, ppMult( a[ind(i,k)], b[ind2(k,j)]));
+        tmp = ppAddv(tmp, ppMult( a[ind(i,k)], b[ind2(k,j)]));
       }
       r[ind2(i,j)] = tmp;
     }
@@ -68,7 +70,7 @@ struct matrixOfPOLY *aa,*bb;
   
 /****************  error handler ************************/
 void errorMatrix(str)
-char *str;
+     char *str;
 {
   fprintf(stderr,"matrix.c: %s\n",str);
   exit(10);
@@ -77,8 +79,8 @@ char *str;
 /********************************************************/
 
 struct arrayOfPOLY *newArrayOfPOLY(size)
-int size;
-/* if size<=0, it returns [0,null]. */
+     int size;
+     /* if size<=0, it returns [0,null]. */
 {
   struct arrayOfPOLY *ap;
   ap = (struct arrayOfPOLY *)sGC_malloc(sizeof(struct arrayOfPOLY));
@@ -94,8 +96,8 @@ int size;
 }
 
 struct matrixOfPOLY *newMatrixOfPOLY(m,n)
-int m,n;
-/* if size<=0, it returns [0,null]. */
+     int m,n;
+     /* if size<=0, it returns [0,null]. */
 {
   struct matrixOfPOLY *ap;
   ap = (struct matrixOfPOLY *)sGC_malloc(sizeof(struct matrixOfPOLY));
@@ -112,9 +114,9 @@ int m,n;
 
 
 struct arrayOfPOLY *carrayToArrayOfPOLY(a,size)
-POLY a[];
-int size;
-/* a[] is read only. */
+     POLY a[];
+     int size;
+     /* a[] is read only. */
 {
   struct arrayOfPOLY *ans;
   int i;