===================================================================
RCS file: /home/cvs/OpenXM/src/kan96xx/Kan/kanExport0.c,v
retrieving revision 1.23
retrieving revision 1.45
diff -u -p -r1.23 -r1.45
--- OpenXM/src/kan96xx/Kan/kanExport0.c	2004/07/30 11:21:55	1.23
+++ OpenXM/src/kan96xx/Kan/kanExport0.c	2005/07/03 11:08:53	1.45
@@ -1,5 +1,7 @@
-/* $OpenXM: OpenXM/src/kan96xx/Kan/kanExport0.c,v 1.22 2004/05/13 04:38:28 takayama Exp $  */
+/* $OpenXM: OpenXM/src/kan96xx/Kan/kanExport0.c,v 1.44 2005/06/16 06:21:21 takayama Exp $  */
 #include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
 #include "datatype.h"
 #include "stackm.h"
 #include "extern.h"
@@ -30,7 +32,8 @@ struct object KooAdd(ob1,ob2)
   POLY r;
   int s,i;
   objectp f1,f2,g1,g2;
-  struct object nn,dd;
+  struct object nn = OINIT;
+  struct object dd = OINIT;
   
   switch (Lookup[ob1.tag][ob2.tag]) {
   case SintegerSinteger:
@@ -161,7 +164,8 @@ struct object KooSub(ob1,ob2)
   int s,i;
   objectp f1,f2,g1,g2;
   extern struct coeff *UniversalZero;
-  struct object nn,dd;
+  struct object nn = OINIT;
+  struct object dd = OINIT;
   
   switch (Lookup[ob1.tag][ob2.tag]) {
   case SintegerSinteger:
@@ -292,7 +296,8 @@ struct object KooMult(ob1,ob2)
   POLY r;
   int i,s;
   objectp f1,f2,g1,g2;
-  struct object dd,nn;
+  struct object dd = OINIT;
+  struct object nn = OINIT;
 
   
   switch (Lookup[ob1.tag][ob2.tag]) {
@@ -438,7 +443,7 @@ struct object KoNegate(obj)
 {
   struct object rob = NullObject;
   extern struct ring SmallRing;
-  struct object tob;
+  struct object tob = OINIT;
   switch(obj.tag) {
   case Sinteger:
     rob = obj;
@@ -480,7 +485,7 @@ struct object KoInverse(obj)
   struct object rob = NullObject;
   extern struct coeff *UniversalOne;
   objectp onep;
-  struct object tob;
+  struct object tob = OINIT;
   switch(obj.tag) {
   case Spoly:
     tob.tag = SuniversalNumber;
@@ -546,11 +551,11 @@ struct object KaoMult(aa,bb)
   POLY tmp;
   POLY fik;
   POLY gkj;
-  struct object rob;
+  struct object rob = OINIT;
   int r1,r2;
   int rsize;
-  struct object tob;
-  struct object ob1;
+  struct object tob = OINIT;
+  struct object ob1 = OINIT;
   extern struct ring SmallRing;
 
   m = getoaSize(aa); m2 = getoaSize(bb);
@@ -610,7 +615,9 @@ struct object KaoMult(aa,bb)
   r1 = isMatrix(aa,m,n); r2 = isMatrix(bb,m2,n2);
   if (r1 == -1 || r2 == -1) {
     /* Object multiplication. Elements are not polynomials. */
-    struct object ofik,ogkj,otmp;
+    struct object ofik = OINIT;
+	struct object ogkj = OINIT;
+	struct object otmp = OINIT;
     rob = newObjectArray(m);
     for (i=0; i<m; i++) {
       getoa(rob,i) = newObjectArray(n2);
@@ -685,10 +692,18 @@ KooEqualQ(obj1,obj2)
      struct object obj1;
      struct object obj2;
 {
-  struct object ob;
+  struct object ob = OINIT;
   int i;
+  extern int Verbose;
   if (obj1.tag != obj2.tag) {
     warningKan("KooEqualQ(ob1,ob2): the datatypes of ob1 and ob2  are not same. Returns false (0).\n");
+	if (Verbose & 0x10) {
+	  fprintf(stderr,"obj1(tag:%d)=",obj1.tag);
+	  printObject(obj1,0,stderr);
+	  fprintf(stderr,", obj2(tag:%d)=",obj2.tag);
+	  printObject(obj2,0,stderr);
+	  fprintf(stderr,"\n"); fflush(stderr);
+	}
     return(0);
   }
   switch(obj1.tag) {
@@ -766,7 +781,7 @@ struct object KooGreater(obj1,obj2)
      struct object obj1;
      struct object obj2;
 {
-  struct object ob;
+  struct object ob = OINIT;
   int tt;
   if (obj1.tag != obj2.tag) {
     errorKan1("%s\n","You cannot compare different kinds of objects.");
@@ -797,6 +812,21 @@ struct object KooGreater(obj1,obj2)
     if ( KopDouble(obj1) > KopDouble(obj2) ) return(KpoInteger(1));
     else return(KpoInteger(0));
     break;
+  case Sarray:
+  {
+    int i,m1,m2;
+    struct object rr = OINIT;
+    m1 = getoaSize(obj1); m2 = getoaSize(obj2);
+    for (i=0; i< (m1>m2?m2:m1); i++) {
+      rr=KooGreater(getoa(obj1,i),getoa(obj2,i));
+      if (KopInteger(rr) == 1) return rr;
+      rr=KooGreater(getoa(obj2,i),getoa(obj1,i));
+      if (KopInteger(rr) == 1) return KpoInteger(0);
+    }
+    if (m1 > m2) return KpoInteger(1);
+    else return KpoInteger(0);
+  }
+  break;
   default:
     errorKan1("%s\n","KooGreater() has not supported these objects yet.");
     break;
@@ -838,6 +868,21 @@ struct object KooLess(obj1,obj2)
     if ( KopDouble(obj1) < KopDouble(obj2) ) return(KpoInteger(1));
     else return(KpoInteger(0));
     break;
+  case Sarray:
+  {
+    int i,m1,m2;
+    struct object rr = OINIT;
+    m1 = getoaSize(obj1); m2 = getoaSize(obj2);
+    for (i=0; i< (m1>m2?m2:m1); i++) {
+      rr=KooLess(getoa(obj1,i),getoa(obj2,i));
+      if (KopInteger(rr) == 1) return rr;
+      rr=KooLess(getoa(obj2,i),getoa(obj1,i));
+      if (KopInteger(rr) == 1) return KpoInteger(0);
+    }
+    if (m1 < m2) return KpoInteger(1);
+    else return KpoInteger(0);
+  }
+  break;
   default:
     errorKan1("%s\n","KooLess() has not supported these objects yet.");
     break;
@@ -852,11 +897,12 @@ struct object KdataConversion(obj,key)
 {
   char tmps[128]; /* Assume that double is not more than 128 digits */
   char intstr[100]; /* Assume that int is not more than 100 digits */
-  struct object rob;
+  struct object rob = OINIT;
   extern struct ring *CurrentRingp;
   extern struct ring SmallRing;
   int flag;
-  struct object rob1,rob2;
+  struct object rob1 = OINIT;
+  struct object rob2 = OINIT;
   char *s;
   int i;
   double f;
@@ -889,6 +935,10 @@ struct object KdataConversion(obj,key)
       return(rob);
     }else if (strcmp(key,"poly") == 0) {
       rob = KpoPOLY(ZERO);
+      return rob;
+    }else if (strcmp(key,"array") == 0) {
+      rob = newObjectArray(0);
+      return rob;
     }else{
       warningKan("Sorry. The data conversion from null to this data type has not supported yet.\n");
     }
@@ -914,8 +964,7 @@ struct object KdataConversion(obj,key)
       strcpy(rob.lc.str,intstr);
       return(rob);
     }else if (strcmp(key,"universalNumber")==0) {
-      rob.tag = SuniversalNumber;
-      rob.lc.universalNumber = intToCoeff(obj.lc.ival,&SmallRing);
+      rob = KintToUniversalNumber(obj.lc.ival);
       return(rob);
     }else if (strcmp(key,"double") == 0) {
       rob = KpoDouble((double) (obj.lc.ival));
@@ -986,7 +1035,7 @@ struct object KdataConversion(obj,key)
     if (strcmp(key,"array") == 0) {
       return(rob);
     }else if (strcmp(key,"list") == 0) {
-      rob = *( arrayToList(obj) );
+      rob = KarrayToList(obj);
       return(rob);
     }else if (strcmp(key,"arrayOfPOLY")==0) {
       rob = KpoArrayOfPOLY(arrayToArrayOfPOLY(obj));
@@ -1000,6 +1049,9 @@ struct object KdataConversion(obj,key)
     }else if (strcmp(key,"null") == 0) {
       rob = NullObject;
       return(rob);
+    }else if (strcmp(key,"byteArray") == 0) {
+      rob = newByteArray(getoaSize(obj),obj);
+      return(rob);
     }else {
 	  { /* Automatically maps the elements. */
 		int n,i;
@@ -1079,12 +1131,13 @@ struct object KdataConversion(obj,key)
     break;
   case Slist:
     if (strcmp(key,"array") == 0) {
-      rob = listToArray(&obj);
+      rob = KlistToArray(obj);
       return(rob);
     }
     break;
   case SuniversalNumber:
     if ((strcmp(key,"universalNumber")==0) || (strcmp(key,"numerator")==0)) {
+      rob = obj;
       return(rob);
     }else if (strcmp(key,"integer")==0) {
       rob = KpoInteger(coeffToInt(obj.lc.universalNumber));
@@ -1102,6 +1155,9 @@ struct object KdataConversion(obj,key)
     }else if (strcmp(key,"double") == 0) {
       rob = KpoDouble( toDouble0(obj) );
       return(rob);
+    }else if (strcmp(key,"denominator") == 0) {
+      rob = KintToUniversalNumber(1);
+      return(rob);
     }else{
       warningKan("Sorry. This type of data conversion of universalNumber has not supported yet.\n");
     }
@@ -1173,18 +1229,43 @@ struct object KdataConversion(obj,key)
       warningKan("Sorryl This type of data conversion of ringp has not supported yet.\n");
     }
     break;
+  case SbyteArray:
+    if (strcmp(key,"array") == 0) {
+      rob = byteArrayToArray(obj);
+      return(rob);
+    } else {
+      warningKan("Sorryl This type of data conversion of ringp has not supported yet.\n");
+    }
+    break;
   default:
     warningKan("Sorry. This type of data conversion has not supported yet.\n");
   }
   return(NullObject);
 }
-    
+
+/* cf. macro to_int32 */
+struct object Kto_int32(struct object ob) {
+  int n,i;
+  struct object otmp = OINIT;
+  struct object rob = OINIT;
+  if (ob.tag == SuniversalNumber) return KdataConversion(ob,"integer");
+  if (ob.tag == Sarray) {
+	n = getoaSize(ob);
+	rob = newObjectArray(n);
+	for (i=0; i<n; i++) {
+	  otmp = Kto_int32(getoa(ob,i));
+	  putoa(rob,i,otmp);
+	}
+	return rob;
+  }
+  return ob;
+}    
 /* conversion functions between primitive data and objects.
    If it's not time critical, it is recommended to use these functions */
 struct object KpoInteger(k)
      int k;
 {
-  struct object obj;
+  struct object obj = OINIT;
   obj.tag = Sinteger;
   obj.lc.ival = k; obj.rc.ival = 0;
   return(obj);
@@ -1192,7 +1273,7 @@ struct object KpoInteger(k)
 struct object KpoString(s)
      char *s;
 {
-  struct object obj;
+  struct object obj = OINIT;
   obj.tag = Sdollar;
   obj.lc.str = s; obj.rc.ival = 0;
   return(obj);
@@ -1200,7 +1281,7 @@ struct object KpoString(s)
 struct object KpoPOLY(f)
      POLY f;
 {
-  struct object obj;
+  struct object obj = OINIT;
   obj.tag = Spoly;
   obj.lc.poly = f; obj.rc.ival = 0;
   return(obj);
@@ -1208,7 +1289,7 @@ struct object KpoPOLY(f)
 struct object KpoArrayOfPOLY(ap)
      struct arrayOfPOLY *ap ;
 {
-  struct object obj;
+  struct object obj = OINIT;
   obj.tag = SarrayOfPOLY;
   obj.lc.arrayp = ap; obj.rc.ival = 0;
   return(obj);
@@ -1217,7 +1298,7 @@ struct object KpoArrayOfPOLY(ap)
 struct object KpoMatrixOfPOLY(mp)
      struct matrixOfPOLY *mp ;
 {
-  struct object obj;
+  struct object obj = OINIT;
   obj.tag = SmatrixOfPOLY;
   obj.lc.matrixp = mp; obj.rc.ival = 0;
   return(obj);
@@ -1226,7 +1307,7 @@ struct object KpoMatrixOfPOLY(mp)
 struct object KpoRingp(ringp)
      struct ring *ringp;
 {
-  struct object obj;
+  struct object obj = OINIT;
   obj.tag = Sring;
   obj.lc.ringp = ringp;
   return(obj);
@@ -1235,20 +1316,30 @@ struct object KpoRingp(ringp)
 struct object KpoUniversalNumber(u)
      struct coeff *u;
 {
-  struct object obj;
+  struct object obj = OINIT;
   obj.tag = SuniversalNumber;
   obj.lc.universalNumber = u;
   return(obj);
 }
+struct object KintToUniversalNumber(n)
+	 int n;
+{
+  struct object rob = OINIT;
+  extern struct ring SmallRing;
+  rob.tag = SuniversalNumber;
+  rob.lc.universalNumber = intToCoeff(n,&SmallRing);
+  return(rob);
+}
+
 /*** conversion 2. Data conversions on arrays and matrices. ****/
 struct object arrayOfPOLYToArray(aa)
      struct arrayOfPOLY *aa;
 {
   POLY *a;
   int size;
-  struct object r;
+  struct object r = OINIT;
   int j;
-  struct object tmp;
+  struct object tmp = OINIT;
 
   size = aa->n; a = aa->array;
   r = newObjectArray(size);
@@ -1263,8 +1354,8 @@ struct object arrayOfPOLYToArray(aa)
 struct object matrixOfPOLYToArray(pmat)
      struct matrixOfPOLY *pmat;
 {
-  struct object r;
-  struct object tmp;
+  struct object r = OINIT;
+  struct object tmp = OINIT;
   int i,j;
   int m,n;
   POLY *mat;
@@ -1287,7 +1378,7 @@ struct arrayOfPOLY *arrayToArrayOfPOLY(oa)
   POLY *a;
   int size;
   int i;
-  struct object tmp;
+  struct object tmp = OINIT;
   struct arrayOfPOLY *ap;
   
   if (oa.tag != Sarray) errorKan1("KarrayToArrayOfPOLY(): %s",
@@ -1315,7 +1406,8 @@ struct matrixOfPOLY *arrayToMatrixOfPOLY(oa)
   int i,j;
   struct matrixOfPOLY *ma;
 
-  struct object tmp,tmp2;
+  struct object tmp = OINIT;
+  struct object tmp2 = OINIT;
   if (oa.tag != Sarray) errorKan1("KarrayToMatrixOfPOLY(): %s",
                                   "Argument is not array\n");
   m = getoaSize(oa);
@@ -1353,8 +1445,8 @@ int objArrayToOrderMatrix(oA,order,n,oasize)
 {
   int size;
   int k,j;
-  struct object tmpOa;
-  struct object obj;
+  struct object tmpOa = OINIT;
+  struct object obj = OINIT;
   if (oA.tag != Sarray) {
     warningKan("The argument should be of the form [ [...] [...] ... [...]].");
     return(-1);
@@ -1439,10 +1531,10 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5)
         ob4 = Order matrix
         ob5 = [(keyword) value (keyword) value ....]
      */
-#define RP_LIMIT 500
+#define RP_LIMIT 5000
 {
   int i;
-  struct object ob;
+  struct object ob = OINIT;
   int c,l,m,n;
   int cc,ll,mm,nn;
   int p;
@@ -1465,7 +1557,7 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5)
   extern char *F_mpMult;
   char *fmp_mult_saved;
   char *mpMultName = NULL;
-  struct object rob;
+  struct object rob = OINIT;
   struct ring *savedCurrentRingp;
 
   /* To get the ring structure. */
@@ -1536,7 +1628,8 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5)
       outputVars[i] = i;
     }
   }
-  
+
+  ob4 = Kto_int32(ob4); /* order matrix */
   oasize = getoaSize(ob4);
   order = (int *)sGC_malloc(sizeof(int)*((2*n)*oasize+1));
   if (order == (int *)NULL) errorKan1("%s\n","No memory.");
@@ -1581,6 +1674,8 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5)
   newRingp->degreeShiftSize = 0;
   newRingp->degreeShiftN = 0;
   newRingp->degreeShift = NULL;
+  newRingp->partialEcart = 0;
+  newRingp->partialEcartGlobalVarX = NULL;
 
   if (ob5.tag != Sarray || (getoaSize(ob5) % 2) != 0) {
     errorKan1("%s\n","[(keyword) value (keyword) value ....] should be given.");
@@ -1637,8 +1732,8 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5)
           errorKan1("%s\n","An array of array should be given. (degreeShift)");
         }
         {
-          struct object ods;
-          struct object ods2;
+          struct object ods = OINIT;
+          struct object ods2 = OINIT;
           int dssize,k,j,nn;
           ods=getoa(ob5,i+1);
           if ((getoaSize(ods) < 1) || (getoa(ods,0).tag != Sarray)) {
@@ -1661,6 +1756,36 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5)
             }
           }
         }
+      } else if (strcmp(KopString(getoa(ob5,i)),"partialEcartGlobalVarX") == 0) {
+        if (getoa(ob5,i+1).tag != Sarray) {
+          errorKan1("%s\n","An array of array should be given. (partialEcart)");
+        }
+        {
+          struct object odv = OINIT;
+          struct object ovv = OINIT;
+          int k,j,nn;
+          char *vname;
+          odv=getoa(ob5,i+1);
+          nn = getoaSize(odv);
+          newRingp->partialEcart = nn;
+          newRingp->partialEcartGlobalVarX = (int *) sGC_malloc(sizeof(int)*nn+1);
+          if (newRingp->partialEcartGlobalVarX == NULL) errorKan1("%s\n","No more memory.");
+          for (j=0; j<nn; j++)
+            (newRingp->partialEcartGlobalVarX)[j] = -1;
+          for (j=0; j<nn; j++) {
+            ovv = getoa(odv,j);
+            if (ovv.tag != Sdollar) errorKan1("%s\n","partialEcartGlobalVarX: string is expected.");
+            vname = KopString(ovv);
+            for (k=0; k<n; k++) {
+              if (strcmp(vname,xvars[k]) == 0) {
+                (newRingp->partialEcartGlobalVarX)[j] = k; break;
+              }else{
+                if (k == n-1) errorKan1("%s\n","partialEcartGlobalVarX: no such variable.");
+              }
+            }
+          }
+        }
+
         switch_function("grade","module1v");
         /* Warning: grading is changed to module1v!! */
       } else {
@@ -1715,8 +1840,8 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5)
 struct object KsetVariableNames(struct object ob,struct ring *rp)
 {
   int n,i;
-  struct object ox;
-  struct object otmp;
+  struct object ox = OINIT;
+  struct object otmp = OINIT;
   char **xvars;
   char **dvars;
   if (ob.tag  != Sarray) {
@@ -1761,7 +1886,7 @@ struct object KswitchFunction(ob1,ob2)
      struct object ob1,ob2;
 {
   char *ans ;
-  struct object rob;
+  struct object rob = OINIT;
   int needWarningForAvoidTheSameRing = 0;
   extern int AvoidTheSameRing;
   if ((ob1.tag != Sdollar) || (ob2.tag != Sdollar)) {
@@ -1803,14 +1928,14 @@ struct object KoReplace(of,rule)
      struct object of;
      struct object rule;
 {
-  struct object rob;
+  struct object rob = OINIT;
   POLY f;
   POLY lRule[N0*2];
   POLY rRule[N0*2];
   POLY r;
   int i;
   int n;
-  struct object trule;
+  struct object trule = OINIT;
   
 
   if (rule.tag != Sarray) {
@@ -1866,7 +1991,7 @@ struct object Kparts(f,v)
 {
   POLY ff;
   POLY vv;
-  struct object obj;
+  struct object obj = OINIT;
   struct matrixOfPOLY *co;
   /* check the data type */
   if (f.tag != Spoly || v.tag != Spoly)
@@ -1883,7 +2008,7 @@ struct object Kparts2(f,v)
 {
   POLY ff;
   POLY vv;
-  struct object obj;
+  struct object obj = OINIT;
   struct matrixOfPOLY *co;
   /* check the data type */
   if (f.tag != Spoly || v.tag != Spoly)
@@ -1931,7 +2056,8 @@ struct object Ksp(ob1,ob2)
      struct object ob1,ob2;
 {
   struct spValue sv;
-  struct object rob,cob;
+  struct object rob = OINIT;
+  struct object cob = OINIT;
   POLY f;
   if (ob1.tag != Spoly || ob2.tag != Spoly)
     errorKan1("%s\n","Ksp(): The arguments must be polynomials.");
@@ -1961,7 +2087,7 @@ struct object Keval(obj)
 {
   char *key;
   int size;
-  struct object rob;
+  struct object rob = OINIT;
   rob = NullObject;
   
   if (obj.tag != Sarray)
@@ -2006,8 +2132,8 @@ char *KremoveSpace(str)
 struct object KtoRecords(ob)
      struct object ob;
 {
-  struct object obj;
-  struct object tmp;
+  struct object obj = OINIT;
+  struct object tmp = OINIT;
   int i;
   int size;
   char **argv;
@@ -2103,7 +2229,7 @@ int KtoArgvbyCurryBrace(str,argv,limit)
 }
 
 struct object KstringToArgv(struct object ob) {
-  struct object rob;
+  struct object rob = OINIT;
   char *s;
   int n,wc,i,inblank;
   char **argv;
@@ -2165,7 +2291,7 @@ static void checkDuplicateName(xvars,dvars,n)
 }
 
 struct object KooPower(struct object ob1,struct object ob2) {
-  struct object rob;
+  struct object rob = OINIT;
   /* Bug. It has not yet been implemented. */
   if (QuoteMode) {
     rob = powerTree(ob1,ob2);
@@ -2287,9 +2413,11 @@ struct object KgbExtension(struct object obj)
 {
   char *key;
   int size;
-  struct object keyo;
+  struct object keyo = OINIT;
   struct object rob = NullObject;
-  struct object obj1,obj2,obj3;
+  struct object obj1 = OINIT;
+  struct object obj2 = OINIT;
+  struct object obj3 = OINIT;
   POLY f1;
   POLY f2;
   POLY f3;
@@ -2458,9 +2586,12 @@ struct object KmpzExtension(struct object obj)
 {
   char *key;
   int size;
-  struct object keyo;
+  struct object keyo = OINIT;
   struct object rob = NullObject;
-  struct object obj0,obj1,obj2,obj3;
+  struct object obj0 = OINIT;
+  struct object obj1 = OINIT;
+  struct object obj2 = OINIT;
+  struct object obj3 = OINIT;
   MP_INT *f;
   MP_INT *g;
   MP_INT *h;
@@ -2483,6 +2614,12 @@ struct object KmpzExtension(struct object obj)
     if (size != 3) errorKan1("%s\n","[(gcd)  universalNumber universalNumber] mpzext.");
     obj1 = getoa(obj,1);
     obj2 = getoa(obj,2);
+    if (obj1.tag != SuniversalNumber) {
+      obj1 = KdataConversion(obj1,"universalNumber");
+	}
+    if (obj2.tag != SuniversalNumber) {
+      obj2 = KdataConversion(obj2,"universalNumber");
+	}
     if (obj1.tag != SuniversalNumber || obj2.tag != SuniversalNumber)
       errorKan1("%s\n","[(gcd)  universalNumber universalNumber] mpzext.");
     if (! is_this_coeff_MP_INT(obj1.lc.universalNumber) ||
@@ -2499,6 +2636,12 @@ struct object KmpzExtension(struct object obj)
     if (size != 3) errorKan1("%s\n","[(tdiv_qr)  universalNumber universalNumber] mpzext.");
     obj1 = getoa(obj,1);
     obj2 = getoa(obj,2);
+    if (obj1.tag != SuniversalNumber) {
+      obj1 = KdataConversion(obj1,"universalNumber");
+	}
+    if (obj2.tag != SuniversalNumber) {
+      obj2 = KdataConversion(obj2,"universalNumber");
+	}
     if (obj1.tag != SuniversalNumber || obj2.tag != SuniversalNumber)
       errorKan1("%s\n","[(tdiv_qr)  universalNumber universalNumber] mpzext.");
     if (! is_this_coeff_MP_INT(obj1.lc.universalNumber) ||
@@ -2557,6 +2700,9 @@ struct object KmpzExtension(struct object obj)
     /*  One arg functions  */
     if (size != 2) errorKan1("%s\n","[key num] mpzext");
     obj1 = getoa(obj,1);
+    if (obj1.tag != SuniversalNumber) {
+      obj1 = KdataConversion(obj1,"universalNumber");
+	}
     if (obj1.tag != SuniversalNumber)
       errorKan1("%s\n","[key num] mpzext : num must be a universalNumber.");
     if (! is_this_coeff_MP_INT(obj1.lc.universalNumber))
@@ -2578,6 +2724,12 @@ struct object KmpzExtension(struct object obj)
     if (size != 3) errorKan1("%s\n","[key  num1 num2] mpzext.");
     obj1 = getoa(obj,1);
     obj2 = getoa(obj,2);
+    if (obj1.tag != SuniversalNumber) {
+      obj1 = KdataConversion(obj1,"universalNumber");
+	}
+    if (obj2.tag != SuniversalNumber) {
+      obj2 = KdataConversion(obj2,"universalNumber");
+	}
     if (obj1.tag != SuniversalNumber || obj2.tag != SuniversalNumber)
       errorKan1("%s\n","[key num1 num2] mpzext.");
     if (! is_this_coeff_MP_INT(obj1.lc.universalNumber) ||
@@ -2609,6 +2761,15 @@ struct object KmpzExtension(struct object obj)
     /* three args */
     if (size != 4) errorKan1("%s\n","[key num1 num2 num3] mpzext");
     obj1 = getoa(obj,1); obj2 = getoa(obj,2); obj3 = getoa(obj,3);
+    if (obj1.tag != SuniversalNumber) {
+      obj1 = KdataConversion(obj1,"universalNumber");
+	}
+    if (obj2.tag != SuniversalNumber) {
+      obj2 = KdataConversion(obj2,"universalNumber");
+	}
+    if (obj3.tag != SuniversalNumber) {
+      obj3 = KdataConversion(obj3,"universalNumber");
+	}
     if (obj1.tag != SuniversalNumber ||
         obj2.tag != SuniversalNumber ||
         obj3.tag != SuniversalNumber ) {
@@ -2627,6 +2788,28 @@ struct object KmpzExtension(struct object obj)
     mpz_powm(r1,f,g,h);
     rob.tag = SuniversalNumber;
     rob.lc.universalNumber = mpintToCoeff(r1,SmallRingp);
+  } else if (strcmp(key,"lcm")==0) {
+    if (size != 3) errorKan1("%s\n","[(lcm)  universalNumber universalNumber] mpzext.");
+    obj1 = getoa(obj,1);
+    obj2 = getoa(obj,2);
+    if (obj1.tag != SuniversalNumber) {
+      obj1 = KdataConversion(obj1,"universalNumber");
+	}
+    if (obj2.tag != SuniversalNumber) {
+      obj2 = KdataConversion(obj2,"universalNumber");
+	}
+    if (obj1.tag != SuniversalNumber || obj2.tag != SuniversalNumber)
+      errorKan1("%s\n","[lcm num1 num2] mpzext.");
+    if (! is_this_coeff_MP_INT(obj1.lc.universalNumber) ||
+        ! is_this_coeff_MP_INT(obj2.lc.universalNumber)) {
+      errorKan1("%s\n","[(lcm)  universalNumber universalNumber] mpzext.");
+    }
+    f = coeff_to_MP_INT(obj1.lc.universalNumber);
+    g = coeff_to_MP_INT(obj2.lc.universalNumber);
+    r1 = newMP_INT();
+    mpz_lcm(r1,f,g);
+    rob.tag = SuniversalNumber;
+    rob.lc.universalNumber = mpintToCoeff(r1,SmallRingp); 
   }else {
     errorKan1("%s\n","mpzExtension(): Unknown tag.");
   }
@@ -2637,7 +2820,7 @@ struct object KmpzExtension(struct object obj)
 /** : context   */
 struct object KnewContext(struct object superObj,char *name) {
   struct context *cp;
-  struct object ob;
+  struct object ob = OINIT;
   if (superObj.tag != Sclass) {
     errorKan1("%s\n","The argument of KnewContext must be a Class.Context");
   }
@@ -2656,10 +2839,10 @@ struct object KcreateClassIncetance(struct object ob1,
                                     struct object ob3)
 {
   /* [class-tag super-obj] size [class-tag]  cclass */
-  struct object ob4;
+  struct object ob4 = OINIT;
   int size,size2,i;
-  struct object ob5;
-  struct object rob;
+  struct object ob5 = OINIT;
+  struct object rob = OINIT;
   
   if (ob1.tag != Sarray)
     errorKan1("%s\n","cclass: The first argument must be an array.");
@@ -2712,8 +2895,8 @@ struct object KpoDouble(double a) {
 double toDouble0(struct object ob) {
   double r;
   int r3;
-  struct object ob2;
-  struct object ob3;
+  struct object ob2 = OINIT;
+  struct object ob3 = OINIT;
   switch(ob.tag) {
   case Sinteger:
     return( (double) (KopInteger(ob)) );
@@ -2742,7 +2925,7 @@ double toDouble0(struct object ob) {
 }
 
 struct object KpoGradedPolySet(struct gradedPolySet *grD) {
-  struct object rob;
+  struct object rob = OINIT;
   rob.tag = Sclass;
   rob.lc.ival = CLASSNAME_GradedPolySet;
   rob.rc.voidp = (void *) grD;
@@ -2759,10 +2942,14 @@ static char *getspace0(int a) {
   return(s);
 }
 struct object KdefaultPolyRing(struct object ob) {
-  struct object rob;
+  struct object rob = OINIT;
   int i,j,k,n;
-  struct object ob1,ob2,ob3,ob4,ob5;
-  struct object t1;
+  struct object ob1 = OINIT;
+  struct object ob2 = OINIT;
+  struct object ob3 = OINIT;
+  struct object ob4 = OINIT;
+  struct object ob5 = OINIT;
+  struct object t1 = OINIT;
   char *s1;
   extern struct ring *CurrentRingp;
   static struct ring *a[N0];
@@ -2833,11 +3020,244 @@ struct object KdefaultPolyRing(struct object ob) {
 }
   
   
+struct object Krest(struct object ob) {
+  struct object rob;
+  struct object *op;
+  int n,i;
+  if (ob.tag == Sarray) {
+    n = getoaSize(ob);
+    if (n == 0) return ob;
+    rob = newObjectArray(n-1);
+    for (i=1; i<n; i++) {
+      putoa(rob,i-1,getoa(ob,i));
+    }
+    return rob;
+  }else if ((ob.tag == Slist) || (ob.tag == Snull)) {
+    return Kcdr(ob);
+  }else{
+    errorKan1("%s\n","Krest(ob): ob must be an array or a list.");
+  }
+}
+struct object Kjoin(struct object ob1, struct object ob2) {
+  struct object rob = OINIT;
+  int n1,n2,i;
+  if ((ob1.tag == Sarray) &&  (ob2.tag == Sarray)) {
+    n1 = getoaSize(ob1); n2 = getoaSize(ob2);
+    rob = newObjectArray(n1+n2);
+    for (i=0; i<n1; i++) {
+      putoa(rob,i,getoa(ob1,i));
+    }
+    for (i=n1; i<n1+n2; i++) {
+      putoa(rob,i,getoa(ob2,i-n1));
+    }
+    return rob;
+  }else if ((ob1.tag == Slist) || (ob1.tag == Snull)) {
+	if ((ob2.tag == Slist) || (ob2.tag == Snull)) {
+	  return KvJoin(ob1,ob2);
+	}else{
+	  errorKan1("%s\n","Kjoin: both argument must be a list.");
+	}
+  }else{
+    errorKan1("%s\n","Kjoin: arguments must be arrays.");
+  }
+}
   
-  
+struct object Kget(struct object ob1, struct object ob2) {
+  struct object rob = OINIT;
+  struct object tob = OINIT;
+  int i,j,size,n;
+  if (ob2.tag == Sinteger) {
+    i =ob2.lc.ival;
+  }else if (ob2.tag == SuniversalNumber) {
+    i = KopInteger(KdataConversion(ob2,"integer"));
+  }else if (ob2.tag == Sarray) {
+    n = getoaSize(ob2);
+    if (n == 0) return ob1;
+    rob = ob1;
+    for (i=0; i<n; i++) {
+      rob=Kget(rob,getoa(ob2,i));
+    }
+    return rob;
+  }
+  if (ob1.tag == Sarray) {
+    size = getoaSize(ob1);
+    if ((0 <= i) && (i<size)) {
+      return(getoa(ob1,i));
+    }else{
+      errorKan1("%s\n","Kget: Index is out of bound. (get)\n");
+    }
+  }else if (ob1.tag == Slist) {
+    rob = NullObject;
+    if (i < 0) errorKan1("%s\n","Kget: Index is negative. (get)"); 
+    for (j=0; j<i; j++) {
+      rob = Kcdr(ob1);
+      if ((ob1.tag == Snull) && (rob.tag == Snull)) {
+        errorKan1("%s\n","Kget: Index is out of bound. (get) cdr of null list.\n");
+      }
+      ob1 = rob;
+    }
+    return Kcar(ob1);
+  } else if (ob1.tag == SbyteArray) {
+    size = getByteArraySize(ob1);
+    if ((0 <= i) && (i<size)) {
+      return(KpoInteger(KopByteArray(ob1)[i]));
+    }else{
+      errorKan1("%s\n","Kget: Index is out of bound. (get)\n");
+    }
+  } else if (ob1.tag == Sdollar) {
+    unsigned char *sss;
+    sss = (unsigned char *) KopString(ob1);
+    size = strlen(sss);
+    if ((0 <= i) && (i<size)) {
+      return(KpoInteger(sss[i]));
+    }else{
+      errorKan1("%s\n","Kget: Index is out of bound. (get)\n");
+    }
 
+  }else errorKan1("%s\n","Kget: argument must be an array or a list.");
+}
+
+/* Constructor of byteArray */
+struct object newByteArray(int size,struct object obj) {
+  unsigned char *ba;
+  unsigned char *ba2;
+  struct object rob = OINIT;
+  struct object tob = OINIT;
+  int i,n;
+  ba = NULL;
+  if (size > 0) {
+    ba = (unsigned char *) sGC_malloc(size);
+    if (ba == NULL) errorKan1("%s\n","No more memory.");
+  }
+  rob.tag = SbyteArray; rob.lc.bytes = ba; rob.rc.ival = size;
+  if (obj.tag == SbyteArray) {
+    n = getByteArraySize(obj);
+    ba2 = KopByteArray(obj);
+    for (i=0; i<(n<size?n:size); i++) {
+      ba[i] = ba2[i];
+    }
+    for (i=n; i<size; i++) ba[i] = 0;
+    return rob;
+  }else if (obj.tag == Sarray) {
+    n = getoaSize(obj);
+    for (i=0; i<n; i++) {
+      tob = getoa(obj,i);
+      tob = Kto_int32(tob);
+      if (tob.tag != Sinteger) errorKan1("%s\n","newByteArray: array is not an array of integer or universalNumber.");
+      ba[i] = (unsigned char) KopInteger(tob);
+    }
+    for (i=n; i<size; i++) ba[i] = 0;
+    return rob;
+  }else{
+    for (i=0; i<size; i++) ba[i] = 0;
+    return rob;
+  }
+}
+struct object newByteArrayFromStr(char *s,int size) {
+  unsigned char *ba;
+  struct object rob = OINIT;
+  int i;
+  ba = NULL;
+  if (size > 0) {
+    ba = (unsigned char *) sGC_malloc(size);
+    if (ba == NULL) errorKan1("%s\n","No more memory.");
+  }
+  rob.tag = SbyteArray; rob.lc.bytes = ba; rob.rc.ival = size;
+  for (i=0; i<size; i++) {
+	ba[i] = (char) s[i];
+  }
+  return(rob);
+}
+
+struct object byteArrayToArray(struct object obj) {
+  int n,i; unsigned char *ba;
+  struct object rob = OINIT;
+  if (obj.tag != SbyteArray) errorKan1("%s\n","byteArrayToArray: argument is not an byteArray.");
+  n = getByteArraySize(obj);
+  rob = newObjectArray(n);
+  ba = KopByteArray(obj);
+  for (i=0; i<n; i++) putoa(rob,i,KpoInteger((int) ba[i]));
+  return rob;
+}
+
+struct object KgetAttributeList(struct object ob){
+  struct object rob = OINIT;
+  if (ob.attr != NULL) rob = *(ob.attr);
+  else rob = NullObject;
+  return rob;
+}
+struct object  KsetAttributeList(struct object ob,struct object attr) {
+  ob.attr = newObject();
+  *(ob.attr) = attr;
+  return ob;
+}
+struct object KgetAttribute(struct object ob,struct object key) {
+  struct object rob = OINIT;
+  struct object alist = OINIT;
+  int n,i;
+  struct object tob = OINIT;
+  char *s;
+  rob = NullObject;
+  if (ob.attr == NULL) return rob;
+  alist = *(ob.attr);
+  if (alist.tag != Sarray) return rob;
+  if (key.tag != Sdollar) return rob;
+  s = KopString(key);
+  n = getoaSize(alist);
+  for (i = 0; i < n; i += 2) {
+    tob = getoa(alist,i);
+    if (tob.tag == Sdollar) {
+      if (strcmp(KopString(tob),s) == 0) {
+        if (i+1 < n) rob = getoa(alist,i+1);
+        return rob;
+      }
+    }
+  }
+  return rob;
+}
+/*  ob (key) (value) setAttribute /ob set. They are not destructive. */
+struct object KsetAttribute(struct object ob,struct object key,struct object value) {
+  struct object rob = OINIT;
+  struct object alist = OINIT;
+  int n,i;
+  char *s = "";
+  struct object tob = OINIT;
+  rob = ob;
+  if (ob.attr == NULL) {
+    rob.attr = newObject();
+    *(rob.attr) = newObjectArray(2);
+    putoa((*(rob.attr)),0,key);
+    putoa((*(rob.attr)),1,value);
+    return rob;
+  }
+  alist = *(ob.attr);
+  if (alist.tag != Sarray) return rob;
+  if (key.tag != Sdollar) {
+    s = KopString(key);
+  }
+  n = getoaSize(alist);
+  for (i = 0; i < n; i += 2) {
+    tob = getoa(alist,i);
+    if (tob.tag == Sdollar) {
+      if (strcmp(KopString(tob),s) == 0) {
+        if (i+1 < n) putoa(alist,i+1,value);
+        return rob;
+      }
+    }
+  }
+
+  rob.attr = newObject();
+  *(rob.attr) = newObjectArray(n+2);
+  for (i=0; i<n; i++) {
+    putoa((*(rob.attr)),i,getoa((*(ob.attr)),i));
+  }
+  putoa((*(rob.attr)),n,key);
+  putoa((*(rob.attr)),n+1,value);
+  return rob;
+}
+
 /******************************************************************
-     error handler
+     Error handler
 ******************************************************************/
 
 errorKan1(str,message)
@@ -2847,7 +3267,9 @@ errorKan1(str,message)
   extern char *GotoLabel;
   extern int GotoP;
   extern int ErrorMessageMode;
+  extern int RestrictedMode, RestrictedMode_saved;
   char tmpc[1024];
+  RestrictedMode = RestrictedMode_saved;
   cancelAlarm();
   if (ErrorMessageMode == 1 || ErrorMessageMode == 2) {
     sprintf(tmpc,"\nERROR(kanExport[0|1].c): ");
@@ -2859,6 +3281,7 @@ errorKan1(str,message)
   if (ErrorMessageMode != 1) {
     fprintf(stderr,"\nERROR(kanExport[0|1].c): ");
     fprintf(stderr,str,message);
+    (void) traceShowStack(); traceClearStack();
   }
   /* fprintf(stderr,"Hello "); */
   if (GotoP) {