Annotation of OpenXM/src/k097/debug/Old/toric0.k.19970310, Revision 1.1
1.1 ! maekawa 1:
! 2: /* SSWork/yacc/ip-k/ex3.ccc, 1996, 8/12. This is original. ---> */
! 3: /* debug/toric0.k */
! 4:
! 5: /* load("debug/db.ccc");
! 6: load("help.ccc"); */
! 7:
! 8: def ringd() {
! 9: R = RingD("s,t,a,b,c,d,b1,b2",
! 10: [["s",1,"t",1],["a",100,"b",10,"c",1,"d",150]]);
! 11: /* 定義をやりなおすと別の ring になるよ。
! 12: これは要注意。 */
! 13: SetRingVariables();
! 14: myinput = [ a - s, b - s*t , c - s*t^2, d - s*t^4 ];
! 15: }
! 16: Println("Type in ringd(); to define the ring.");
! 17: Println("Type in toric(myinput) to get the toric ideal.");
! 18:
! 19: sm1(" /mydegree {2 1 roll degree} def ");
! 20: def expvec(f) {
! 21: local vars,n,i,ans;
! 22: vars = [a,b,c,d,Da,Db,Dc,Dd]; /* change here. */
! 23: if (f == Poly("0")) return([ ] );
! 24: sm1(f," /ff.tmp set ", vars ,
! 25: " { ff.tmp mydegree (universalNumber) dc } map /FunctionValue set ");
! 26: }
! 27:
! 28: def in_w(f) {
! 29: local w,top,ans,wtop;
! 30: weight = [-100,-10,-1,0,100,10,1,0]; /* change here. */
! 31: if (f == Poly("0")) return( Poly("0") );
! 32: top = sm1(f, " init ");
! 33: w = expvec(top) * weight ;
! 34: wtop = w;
! 35: ans = top;
! 36: f = f-top;
! 37: while(true) {
! 38: if ( f == Poly("0")) sm1("exit");
! 39: top = sm1(f, " init ");
! 40: w = expvec(top) * weight ;
! 41: if (w < wtop) sm1(" exit ");
! 42: ans = ans + top;
! 43: f = f-top;
! 44: /* Println(ans); */
! 45: }
! 46: return(ans);
! 47: }
! 48:
! 49: def toric(input) {
! 50: local ans ,ans2;
! 51: ans = Replace(Groebner(input),[[h,Poly("1")]]);
! 52: ans2 = sm1( ans, " [(s) (t)] eliminatev ");
! 53: return(ans2);
! 54: }
! 55:
! 56:
! 57:
! 58: def ringd1() {
! 59: local weight,ans;
! 60: weight = ["a",100,"b",10,"c",1,"d",150];
! 61: R = RingD("s,t,a,b,c,d,b1,b2",
! 62: [["s",1,"t",1],weight]);
! 63: SetRingVariables();
! 64: Print("Weight="); Println(weight);
! 65: myinput = [ a - s, b - s*t , c - s*t^2, d - s*t^4 ];
! 66: ans = toric(myinput);
! 67: Println(ans);
! 68: return(ans);
! 69: }
! 70:
! 71: def void ringd2() {
! 72: local weight;
! 73: weight = ["a",1,"b",15,"c",1,"d",4];
! 74: R = RingD("s,t,a,b,c,d,b1,b2",
! 75: [["s",1,"t",1],weight]);
! 76: SetRingVariables();
! 77: Print("Weight="); Println(weight);
! 78: myinput = [ a - s, b - s*t , c - s*t^2, d - s*t^4 ];
! 79: Println(toric(myinput));
! 80: }
! 81:
! 82: /* This is used to find a bug */
! 83: def void ringd3() {
! 84: local weight;
! 85: weight = ["a",1,"b",15,"c",1,"d",4];
! 86: R = RingD("s,t,a,b,c,d,b1,b2",
! 87: [["s",1,"t",1],weight]);
! 88: SetRingVariables();
! 89: /* Print("Weight="); Println(weight);
! 90: myinput = [ a - s, b - s*t , c - s*t^2, d - s*t^4 ]; */
! 91: /* Println(toric(myinput)); */
! 92: }
! 93:
! 94: def MinimalBase(base) {
! 95: sm1( base, " minimalBase /FunctionValue set ");
! 96: }
! 97: HelpAdd(["MinimalBase",
! 98: ["MinimalBase[base] prunes redundant elements in the Grobner basis <<base>>."
! 99: ]]);
! 100:
! 101:
! 102: def IndexedVariables(name,size) {
! 103: local result,i,result2;
! 104: result = [ ];
! 105: for (i=0; i<size-1; i++) {
! 106: result = Append(result,Indexed(name,i));
! 107: result = Append(result,",");
! 108: }
! 109: if (size-1 >= 0) {
! 110: result = Append(result,Indexed(name,size-1));
! 111: }
! 112: result2 = Join(["{"],result);
! 113: result2 = Join(result2,["}"]);
! 114: return(AddString(result2));
! 115: }
! 116: HelpAdd(["IndexedVariables",
! 117: ["IndexedVariables(name,size) returns the string ",
! 118: " {name[0],name[1],...,name[size-1]} which can be used as inputs to ",
! 119: " the function RingD (string name, integer size).",
! 120: " Ex. R = RingD(IndexedVariables(\"a\",3)); ",
! 121: " h = Poly(\"h\");",
! 122: " a = NewArray(3);",
! 123: " for (i=0; i<3; i++) {a[i] = Poly(Indexed(\"a\",i));} ;"]]);
! 124:
! 125:
! 126: def IndexedRingD(vList, size, weightMatrix) {
! 127: local myring,tmp,k00_i,argsize,vListD;
! 128: /* You cannot use these local varialbes as a name of global ring
! 129: variables. Change these names to names that start with k00_ */
! 130: argsize = Length(Arglist);
! 131: if (argsize == 1) {
! 132: Println("Error (IndexedRingD): ");
! 133: return(null);
! 134: }
! 135: if (argsize == 2) {
! 136: vListD = AddString(["D",vList]);
! 137: myring = RingD(IndexedVariables(vList,size));
! 138: SetRingVariables();
! 139: tmp = NewArray(size);
! 140: for (k00_i=0; k00_i<size; k00_i++) {tmp[k00_i]=Poly(Indexed(vList,k00_i));}
! 141: sm1(vList, " (literal) dc ", tmp, " def ");
! 142: tmp = NewArray(size);
! 143: for (k00_i=0; k00_i<size; k00_i++) {tmp[k00_i]=Poly(Indexed(vListD,k00_i));}
! 144: sm1(vListD, " (literal) dc ", tmp, " def ");
! 145: Print("Set the global variables ");
! 146: sm1("[(parse) ",vList," ] extension pop print ");
! 147: sm1("[(parse) ",vListD," ] extension pop print "); Ln();
! 148: return( myring );
! 149: }
! 150: if (argsize == 3) {
! 151: vListD = AddString(["D",vList]);
! 152: myring = RingD(IndexedVariables(vList,size),weightMatrix);
! 153: SetRingVariables();
! 154: tmp = NewArray(size);
! 155: for (k00_i=0;k00_i<size; k00_i++) {tmp[k00_i]=Poly(Indexed(vList,k00_i));}
! 156: sm1(vList, " (literal) dc ", tmp, " def ");
! 157: tmp = NewArray(size);
! 158: for (k00_i=0;k00_i<size; k00_i++) {tmp[k00_i]=Poly(Indexed(vListD,k00_i));}
! 159: sm1(vListD, " (literal) dc ", tmp, " def ");
! 160: Print("Set the global variables ");
! 161: sm1("[(parse) ",vList," ] extension pop print ");
! 162: sm1("[(parse) ",vListD," ] extension pop print "); Ln();
! 163: return( myring );
! 164: }
! 165: return(-1);
! 166: }
! 167:
! 168:
! 169:
! 170: def void foo() {
! 171: Println("IndexedVariables:");
! 172: Println(IndexedVariables("xyz",5));
! 173: }
! 174:
! 175:
! 176:
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to toric0.k.19970310 CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM / src / k097 / debug / Old