=================================================================== RCS file: /home/cvs/OpenXM/src/kan96xx/Kan/kanExport0.c,v retrieving revision 1.39 retrieving revision 1.49 diff -u -p -r1.39 -r1.49 --- OpenXM/src/kan96xx/Kan/kanExport0.c 2004/09/20 02:26:56 1.39 +++ OpenXM/src/kan96xx/Kan/kanExport0.c 2015/10/08 11:49:37 1.49 @@ -1,5 +1,7 @@ -/* $OpenXM: OpenXM/src/kan96xx/Kan/kanExport0.c,v 1.38 2004/09/20 02:11:22 takayama Exp $ */ +/* $OpenXM: OpenXM/src/kan96xx/Kan/kanExport0.c,v 1.48 2012/09/16 01:53:08 takayama Exp $ */ #include +#include +#include #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; im2?m2:m1); i++) { rr=KooGreater(getoa(obj1,i),getoa(obj2,i)); @@ -864,7 +871,7 @@ struct object KooLess(obj1,obj2) case Sarray: { int i,m1,m2; - struct object rr; + 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)); @@ -890,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; @@ -932,6 +940,7 @@ struct object KdataConversion(obj,key) rob = newObjectArray(0); return rob; }else{ + /* fprintf(stderr,"key=%s\n",key); */ warningKan("Sorry. The data conversion from null to this data type has not supported yet.\n"); } break; @@ -1238,8 +1247,8 @@ struct object KdataConversion(obj,key) /* cf. macro to_int32 */ struct object Kto_int32(struct object ob) { int n,i; - struct object otmp; - struct object rob; + struct object otmp = OINIT; + struct object rob = OINIT; if (ob.tag == SuniversalNumber) return KdataConversion(ob,"integer"); if (ob.tag == Sarray) { n = getoaSize(ob); @@ -1257,7 +1266,7 @@ struct object Kto_int32(struct object ob) { 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); @@ -1265,7 +1274,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); @@ -1273,7 +1282,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); @@ -1281,7 +1290,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); @@ -1290,7 +1299,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); @@ -1299,7 +1308,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); @@ -1308,7 +1317,7 @@ 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); @@ -1316,7 +1325,7 @@ struct object KpoUniversalNumber(u) struct object KintToUniversalNumber(n) int n; { - struct object rob; + struct object rob = OINIT; extern struct ring SmallRing; rob.tag = SuniversalNumber; rob.lc.universalNumber = intToCoeff(n,&SmallRing); @@ -1329,9 +1338,9 @@ struct object arrayOfPOLYToArray(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); @@ -1346,8 +1355,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; @@ -1370,7 +1379,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", @@ -1398,7 +1407,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); @@ -1436,8 +1446,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); @@ -1522,10 +1532,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; @@ -1548,7 +1558,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. */ @@ -1650,6 +1660,7 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5) newRingp->cc = cc; newRingp->x = xvars; newRingp->D = dvars; + newRingp->Dsmall = makeDsmall(dvars,n); /* You don't need to set order and orderMatrixSize here. It was set by setOrder(). */ setFromTo(newRingp); @@ -1723,8 +1734,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)) { @@ -1752,8 +1763,8 @@ int KsetUpRing(ob1,ob2,ob3,ob4,ob5) errorKan1("%s\n","An array of array should be given. (partialEcart)"); } { - struct object odv; - struct object ovv; + struct object odv = OINIT; + struct object ovv = OINIT; int k,j,nn; char *vname; odv=getoa(ob5,i+1); @@ -1831,8 +1842,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) { @@ -1877,7 +1888,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)) { @@ -1919,14 +1930,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) { @@ -1982,7 +1993,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) @@ -1999,7 +2010,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) @@ -2047,7 +2058,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."); @@ -2077,7 +2089,7 @@ struct object Keval(obj) { char *key; int size; - struct object rob; + struct object rob = OINIT; rob = NullObject; if (obj.tag != Sarray) @@ -2122,8 +2134,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; @@ -2219,7 +2231,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; @@ -2259,6 +2271,56 @@ struct object KstringToArgv(struct object ob) { return(rob); } +struct object KstringToArgv2(struct object ob,struct object oseparator) { + struct object rob = OINIT; + char *s; + int n,wc,i,inblank; + char **argv; + int separator; + if (ob.tag != Sdollar) + errorKan1("%s\n","KstringToArgv2(): the argument must be a string."); + if (oseparator.tag == Sinteger) { + separator = KopInteger(oseparator); + }else if (oseparator.tag == Sdollar) { + s = KopString(oseparator); + separator=s[0]; + }else { + errorKan1("%s\n","KstringToArgv2(ob,separator):the argument must be strings."); + } + n = strlen(KopString(ob)); + s = (char *) sGC_malloc(sizeof(char)*(n+2)); + if (s == NULL) errorKan1("%s\n","KstringToArgv(): No memory."); + strcpy(s,KopString(ob)); + inblank = 1; wc = 0; + for (i=0; i 0) { @@ -3133,9 +3205,25 @@ struct object newByteArray(int size,struct object obj) 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