Annotation of OpenXM/src/k097/lib/minimal/minimal-note-ja.txt, Revision 1.12
1.12 ! takayama 1: $OpenXM: OpenXM/src/k097/lib/minimal/minimal-note-ja.txt,v 1.11 2000/08/02 05:14:30 takayama Exp $
1.1 takayama 2:
3: SpairAndReduction() :
4: $BM?$($i$l$?(B pair $B$r(B reduction $B$9$k(B.
5: V-minimal $B$KI,MW$+$I$&$+$NH=Dj$b$9$k(B.
6:
7: SpairAndReduction2():
8: tower2 = StowerOf(tower,level-1);
9: SsetTower(tower2);
10: /** sm1(" show_ring "); */
11:
12: $BM?$($i$l$?(B pair $B$r(B reduction $B$9$k$?$a$N(B schreyer order
13: $B$r@_Dj$9$k(B. Resolution $B$N?<$5$K1~$8$F(B, tower $B$b?<$/$9$kI,MW$,$"$k(B.
14:
15:
16: if (IsConstant(t_syz[i])){
17:
18: Syzygy $B$r$_$F(B, $BDj?t@.J,$,$J$$$+(B check.
19: t_syz[i] $B$,Dj?t@.J,$G$"$l$P(B, $B0l$DA0$N(B GB $B$N9=@.MWAG$G$"$k(B
20: g_i $B$,M>J,$J(B GB $B$G$"$k2DG=@-$,$?$+$$(B.
21: SpairAndReduction() ( LaScala-Stillman $B$NJ}K!(B) $B$H$N@09g@-$r$H$k$?$a(B
22: g_i $B$r(B tmp[0] $B$KBeF~$7(B ( reduction $B$G$-$J$+$C$?$U$j$r$9$k(B )
23: g_i $B$N(B V-degree $B$r$7$i$Y$k(B.
24:
25:
26: Sannfs2_laScala2()
27: Sannfs3_laScala2() $B$r:n$k(B.
28:
29: $BFs$D$N%"%k%4%j%:%`$NHf3S(B.
30: In(11)=sm1_pmat(a1[1]); $B$N=gHV$r$+$($k(B.
31: [
32: [ 3*Dx^2*h , 0 , Dy , -Dz ]
33: [ 6*x*Dy*Dz^2-9*y*Dx^2*h , -2*x*Dy*Dz^2+3*y*Dx^2*h , -2*x*Dx-3*y*Dy , 0]
34: [ 2*x*Dx+3*z*Dz-h^2 , y*Dy-z*Dz , 0 , 0 ]
35: [ 2*x*Dy*Dz , 0 , z , -y ]
36:
37: [ 0 , 2*x*Dy^2*Dz-3*z*Dx^2*h , 0 , 2*x*Dx+3*z*Dz ]
38: ]
39: In(12)=sm1_pmat(a2[1]);
40: [
41: [ 3*Dx^2*h , 0 , Dy , -Dz ]
42: [ 6*x*Dy*Dz^2-9*y*Dx^2*h , -2*x*Dy*Dz^2+3*y*Dx^2*h , -2*x*Dx-3*y*Dy, 0 ]
43: [ 2*x*Dx+3*z*Dz-h^2 , y*Dy-z*Dz , 0 , 0 ]
44: [ 2*x*Dy*Dz , 0 , z , -y ]
45:
46: [ 9*z*Dx^2*h , 2*x*Dy^2*Dz-3*z*Dx^2*h , 3*z*Dy , 2*x*Dx ]
47: [ 2*x*Dx*Dz^2+3*z*Dz^3+5*Dz^2*h^2 , y*Dy*Dz^2-z*Dz^3-2*Dz^2*h^2 , 0 , 0 ]
48: ]
49: In(13)=
50:
51: ----------------------
52: In(16)=sm1_pmat(a1[2]);
53: [
54: [ -2*x*Dx-3*y*Dy-3*z*Dz-6*h^2 , -Dy , -Dz , 3*Dx^2*h , 3*Dy^2 , 3*Dy*Dz , -2*x*Dy , 2*x*Dz , 0 ]
55: [ 3*y*z , z , y , -2*x*Dy*Dz , -3*z*Dy , 2*x*Dx , 2*x*z , -2*x*y , 0 ]
56: ]
57: In(17)=sm1_pmat(a2[2]);
58: [
59: [ -y , 2*x*Dy*Dz , z , 0 , 2*x*Dx , 0 ]
60: [ -Dz , 3*Dx^2*h , Dy , -2*x*Dx-3*y*Dy-3*h^2 , -3*Dy*Dz , 0 ]
61: ]
62: In(18)=
63:
1.2 takayama 64: ---------------------------
65:
66: May 22, (Tue), 5:50 (Spain local time, 12:50 JST)
67:
68: kan96xx/Kan/resol.c $B$G(B,
69: RemoveRedundantInSchreyerSkelton = 0
70: $B$KJQ$($F(B ($B$3$N(B option $B$b$"$?$i$7$/2C$($k(B), schreyer $B$,@5$7$/F0$/$+(B
71: $BD4$Y$k$3$H$K$9$k(B.
72: ( commit $B$O(B kan96xx $B$H(B k097 $BN>J}$9$Y$7(B.)
73:
74: test8() $B$G(B sm1 $B$G=q$$$?J}$N(B Schreyer $B$r8+$k$H(B,
75: RemoveRedundantInSchreyerSkelton = 1
76: $B$G$b(B,
77: kernel = image
78: $B$H$J$C$F$$$k$N$G0J8e$3$N(B option $B$O(B 1 $B$N$^$^;H$&$3$H$H$9$k(B.
79: $BMW$9$k$K(B k0 $B$N%3!<%I$,$I$&$d$i$*$+$7$$$i$7$$(B.
1.4 takayama 80: ==>
81: 6/8 $B$N%N!<%H$h$j(B.
82: syzygy $B$r(B homogenization $B$r2p$7$F7W;;$9$k$N$OLdBj$"$j(B.
83: --> usage of isExact
84: $BMW$9$k$K(B kernel = image $B$N%3!<%I$bJQ(B. Homogenized $B$N$^$^$d$kI,MW$"$j(B.
1.3 takayama 85:
86: -----------------------------------
87: June 8, 2000 (Thu), 9:10 (Spain local time)
88: hol.sm1 : gb_h, syz_h, isSameIdeal, isSameIdeal_h
89: complex.sm1 : isExact, isExact_h
90:
91: syzygy $B$r(B homogenization $B$r2p$7$F7W;;$9$k$N$OLdBj$"$j(B.
92: --> usage of isExact
93:
94: [(Homogenize_vec) 0] system_variable : vector $B$N(B homogenize $B$r$7$J$$(B.
95: (grade) (module1v) switch_function : vector $BJQ?t$O(B, total
96: degree $B$K?t$($J$$(B.
97: ==> $BL58B%k!<%W$KCm0U(B ---> gb_h, syz_h $B$N(B usage.
98:
99: minimal-test.k $B$N(B ann(x^3-y^2*z^2) $B$N(B laplace $BJQ49$N(B
100: betti $B?t$,JQ(B, exact $B$G$J$$(B, $B$r(B isExact_h $B$G(B check
101: $B$7$h$&(B.
102:
103: minimal-test.k
104: test10();
105: LaScala-Stillman $B$NJ}K!$G$D$/$C$?(B, schreyer resol $B$,(B exact $B$+(B
106: $BD4$Y$k(B.
107: $BNcBj$O(B, ann(1/(x^3-y^2 z^2)) $B$N(B Laplace $BJQ49(B.
1.4 takayama 108: ==> OK. IsExact_h $B$G$7$i$Y$k(B. (IsExact $B$O$@$a$h(B)
109:
110: June 8, 2000 (Thu), 19:35
111: load["minimal-test.k"];;
112: test11();
113: LaScala-Stillman $B$NJ}K!$G$D$/$C$?(B, minimal resol $B$,(B exact $B$+(B
114: $BD4$Y$k(B.
115: $BNcBj$O(B, ann(1/(x^3-y^2 z^2)) $B$N(B Laplace $BJQ49(B.
116:
117: SwhereInTower $B$r;H$&$H$-$O(B,
118: SsetTower() $B$G(B gbList $B$rJQ99$7$J$$$H$$$1$J$$(B.
119: $B$b$A$m$s;HMQ$7$?$i(B, $B$=$l$rLa$9$3$H(B.
120: SpairAndReduction, SpairAndReduction2 $B$G(B,
121: SsetTower(StowerOf(tower,level));
122: pos = SwhereInTower(syzHead,tower[level]);
123:
124: SsetTower(StowerOf(tower,level-1));
125: pos2 = SwhereInTower(tmp[0],tower[level-1]);
126: $B$H(B, SwhereInTower $B$NA0$K(B setTower $B$r$/$o$($?(B.
127: ( $B0c$&%l%Y%k$G$NHf3S$N$?$a(B.)
128:
129: IsExact_h $B$O(B, 0 $B%Y%/%H%k$r4^$`>l9g(B, $B$?$@$7$/F0:n$7$J$$$h$&$@(B.
130: test11().
131: test11a() $B$G(B, 0 $B%Y%/%H%k$r<j$G=|$$$?9TNs$N(B exactness $B$r%A%'%C%/(B. ==> OK.
132:
133:
134: ---------------------------------
135: June 9, 6:20
136: SpairAndReduction
137: $B$H(B
138: SpairAndReduction2
139: $B$N0c$$(B.
140: SpairAndReduction : SlaScala (LaScala-Stillman's algorithm $B$G;H$&(B)
141: SpairAndReduction2 : Sschreyer (schreyer algorithm $B$G;H$&(B, laScala $B$O$J$7(B.)
142:
143: 0 $B$r<+F0$G=|$/%3!<%I$r=q$3$&(B.
144:
145: SpruneZeroRow() $B$r(B Sminimal() $B$K2C$($?(B.
146: test11() $B$b@5$7$/F0:n$9$k$O$:(B.
147: IsExact_h $B$O(B schreyer $B$r(B off $B$7$F(B, ReParse $B$7$F$+$i(B,
148: $B8F$S=P$9$3$H(B.
149:
150:
151: #ifdef TOTAL_STRATEGY
152: return(StotalDegree(f)+Sdegree(tower[level-2,i],tower,level-1));
153: #endif
154: /* Strategy must be compatible with ordering. */
155: /* Weight vector must be non-negative, too. */
156: /* See Sdegree, SgenerateTable, reductionTable. */
157: wd = Sord_w(f,ww);
158: return(wd+Sdegree(tower[level-2,i],tower,level-1));
159: TOTAL_STRATEGY $B$rMQ$$$kI,MW$,$"$k$N$G$O(B??
160: Example 1: Sweyl("x,y",[["x",-1,"y",-1,"Dx",1,"Dy",1]]);
161: v=[[2*x*Dx + 3*y*Dy+6, 0],
162: [3*x^2*Dy + 2*y*Dx, 0],
163: [0, x^2+y^2],
164: [0, x*y]];
165: a=Sminimal(v);
166: strategy $B$,$*$+$7$$$H$$$C$F$H$^$k(B. $BM}M3$O(B?
167:
168: a=test_ann3("x^3+y^3+z^3); $B$O;~4V$,$+$+$j$=$&(B.
169: a=test_ann3("x^3+y^3"); OK.
170: a=test_ann3("x^2+y^2+z"); OK.
171:
172:
173: $B>e$N(B example 1 $B$N%(%i!<(B $B$N8+J}(B:
174: Processing [ 1 , 3 ] Strategy = 2
175: 1 $B$N(B 3 $BHVL\$N(B spair $B$N(B reduction $B$r=hM}Cf(B.
176: In(7)=reductionTable:
177: [[ 1 , 1 , 1 , 2 , 2 , 3 ] , [ 3 , 2 , 1 , 2 , 3 ] , [ 2 ] ]
178: -- $B$3$l(B.
179: SpairAndReduction:
180: [ p and bases , [ [ 0 , 3 ] , [ y*h , -x ] ] , [ 2*x*Dx+3*y*Dy+6*h^2 , e_*x^2+e_*y^2 , e_*x*y , 2*y*Dx*h+3*x^2*Dy , e_*y^3 , %[null] ] ]
181: 0 $B$N(B 0 $BHVL\$H(B 3 $BHVL\(B $B$N(B spair $B$r7W;;$7$F(B, 0 $B%l%Y%k$N(B gb $B$G(B reduction.
182: [ 1 , 1 , 1 , 2 , 2 , 3 ] $B$K$"$k$h$&$K(B, strategy 3 $B0J30$O7W;;$:$_(B.
183: ( $B7W;;$7$F$J$$$b$N$O(B %[null] $B$H$J$C$F$k(B. )
184: [ level= , 1 ]
185: [ tower2= , [ [ ] ] ] ( $B0lHV2<$J$N$G(B, tower $B$O$J$7$h(B. )
186: [ y*h , -es^3*x ]
187: [gi, gj] = [ 2*x*Dx+3*y*Dy+6*h^2 , 2*y*Dx*h+3*x^2*Dy ]
188: 1
189: Reduce the element 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy
190: by [ 2*x*Dx+3*y*Dy+6*h^2 , e_*x^2+e_*y^2 , e_*x*y , 2*y*Dx*h+3*x^2*Dy , e_*y^3 , %[null] ]
191: result is [ 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy , 1 , [ 0 , 0 , 0 , 0 , 0 , 0 ] ]
192: vdegree of the original = -1
193: vdegree of the remainder = -1
194: [ 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy , [ y*h , 0 , 0 , -x , 0 , 0 ] , 3 , 5 , -1 , -1 ]
195:
196: In(11)=freeRes:
197: [ [ 2*x*Dx+3*y*Dy+6*h^2 , e_*x^2+e_*y^2 , e_*x*y , 2*y*Dx*h+3*x^2*Dy , e_*y^3 , 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy ] , [ %[null] , [ 0 , 0 , y^2 , 0 , -x , 0 ] , [ 0 , -y , x , 0 , 1 , 0 ] , [ -y*h , 0 , 0 , x , 0 , 1 ] , %[null] ] , [ %[null] ] ]
198: $B$r$_$l$P$o$+$k$h$&$K(B, SlaScala $B$G(B, freeRes $B$K$3$N85$,(B [0,5] $B$K2C$((B
199: $B$i$l$?(B.
200:
201: $B<!$K(B SnextI $B$,(B SlaScala $B$h$j8F$P$l$F$3$N%(%i!<(B.
202: i = SnextI(reductionTable_tmp,strategy,redundantTable,
203: skel,level,freeRes);
204: In(22)=reductionTable:
205: [ [ 1 , 1 , 1 , 2 , 2 , 3 ] , [ 3 , 2 , 1 , 2 , 3 ] , [ 2 ] ]
206: $B$J$N$G(B, $B:G8e(B $B$N(B 2 $B$,=hM}$5$l$k$O$:$@$,(B,
207: In(25)=skel[2]:
208: [ [ [ 0 , 2 ] , [ 1 , -y^2 ] ] ]
209: $B$N$h$&$K(B, 0 $BHVL\$H(B, 2 $BHVL\$N(B spair.
210: $B$7$+$7(B,
211: In(26)=bases:
212: [ %[null] , [ 0 , 0 , y^2 , 0 , -x , 0 ] , [ 0 , -y , x , 0 , 1 , 0 ] , [ -y*h , 0 , 0 , x , 0 , 1 ] , %[null] ]
213: $B$N$h$&$K(B, 0 $BHVL\$O(B strategy 3 $B$J$N$G(B, $B$^$@$b$H$^$C$F$$$J$$(B.
214:
215: reductionTable_tmp=[ 2 ]
216: See also reductionTable, strategy, level,i
217: ERROR(sm): error operator : SnextI: bases[i] or bases[j] is null for all combinations.
218: --- Engine error or interrupt : In function : Error of class PrimitiveObject
219:
220: Type in Cleards() to exit the debug mode and Where() to see the stack trace.
221: In(7)=reductionTable:
222: [ [ 1 , 1 , 1 , 2 , 2 , 3 ] , [ 3 , 2 , 1 , 2 , 3 ] , [ 2 ] ]
223: In(8)=strategy:
224: 2
225: In(9)=level:
226: 2
227:
228: RemoveRedundantInSchreyerSkelton = 0
229: $B$H$7$F$bF1$8%(%i!<(B.
230:
231: -------------------------------------------------
232: test_ann3("x*y+y*z+z*x"); OK.
233:
234: 6/9 (Fri)
235: Sminimal $B$N<BAu$KAjJQ$o$i$:6lO+$7$F$^$9(B.
236: Sevilla $B$G$$$m$$$m$HD>$7$?7k2L(B,
237: Sminimal $B$O$&$^$/$&$4$1$P@5$7$$Ez$($r$@$7$F$k$_$?$$$G$9$,(B
238: (D<h> : homogenized Weyl $B$G(B ker = im $B$r(B check $B$7$F$k(B,
239: V-adapted (strict) $B$+$I$&$+$N(B check routing $B$O$^$@=q$$$F$J$$(B),
240: strategy $B$,$&$^$/$&$4$+$J$/$F$H$^$k>l9g$b$"$j$^$9(B
241: ( strategy = 2 $B$N(B sp $B$r7W;;$9$k$N$K(B, strategy 3 $B$N(B $B85$rI,MW$H(B
242: $B$7$?$j$9$k>l9g$"$j(B).
243:
244:
245: strategy $B$O(B
246: def Sdegree(f,tower,level) {
247: local i,ww, wd;
248: /* extern WeightOfSweyl; */
249: ww = WeightOfSweyl;
250: f = Init(f);
251: if (level <= 1) return(StotalDegree(f));
252: i = Degree(f,es);
253: return(StotalDegree(f)+Sdegree(tower[level-2,i],tower,level-1));
254: }
255: $B$rMQ$$$F(B,
256: ans_at_each_floor[j] = Sdegree(tower[i,j],tower,i+1)-(i+1)
257: $B$G7W;;$7$F$^$9(B.
258:
259: $B$$$/$D$+=PNO$r$D$1$F$*$-$^$9$N$G(B, $B8!F$(B!!!
260:
261: $BNc(B 1:
262: load["minimal-test.k"];;
263: a=test_ann3("x^3-y^2*z^2"); $B0z?t$N(B annihilating ideal $B$N(B laplace $BJQ49$N(B
264: homogenization $B$N(B resolution.
265: weight vector $B$O(B (-1,-1,-1,1,1,1)
266:
267: In(4)=sm1_pmat(a[1]);
268: [
269: [ 0 $B<!(B
270: [ y*Dy-z*Dz ]
271: [ -2*x*Dx-3*z*Dz+h^2 ]
272: [ 2*x*Dy*Dz^2-3*y*Dx^2*h ]
273: [ 2*x*Dy^2*Dz-3*z*Dx^2*h ]
274: ]
275: [ 1 $B<!(B
276: [ 3*Dx^2*h , 0 , Dy , -Dz ]
277: [ 6*x*Dy*Dz^2-9*y*Dx^2*h , -2*x*Dy*Dz^2+3*y*Dx^2*h , -2*x*Dx-3*y*Dy , 0 ]
278: [ 0 , 2*x*Dy^2*Dz-3*z*Dx^2*h , 0 , 2*x*Dx+3*z*Dz ]
279: [ 2*x*Dx+3*z*Dz-h^2 , y*Dy-z*Dz , 0 , 0 ]
280: [ 2*x*Dy*Dz , 0 , z , -y ]
281: ]
282: [ 2 $B<!(B
283: [ -2*x*Dx-3*y*Dy-3*z*Dz-6*h^2 , -Dy , -Dz , 3*Dx^2*h , 3*Dy*Dz ]
284: [ 3*y*z , z , y , -2*x*Dy*Dz , 2*x*Dx ]
285: ]
286: ]
287: In(5)=
288:
289: $BNc(B 2:
290: load["minimal-test.k"];;
291: a=test_ann3("x*y+y*z+z*x");
292: In(6)=sm1_pmat(a[1]);
293: [
294: [ 0 $B<!(B
295: [ 2*x*Dx+x*Dz-y*Dz+z*Dz+h^2 ]
296: [ -2*y*Dy+x*Dz-y*Dz-z*Dz-h^2 ]
297: [ -2*x*Dy+2*z*Dy+x*Dz-y*Dz+3*z*Dz+h^2 ]
298: [ -2*y*Dx+2*z*Dx-x*Dz+y*Dz+3*z*Dz+h^2 ]
299: ]
300: [ 1 $B<!(B
301: [ y-z , x-z , -y , x ]
302: [ 2*Dy-2*Dz , 2*Dx-2*Dz , 2*Dx+2*Dz , -2*Dy-2*Dz ]
303: [ 2*y*Dx-2*z*Dx+x*Dz-y*Dz-3*z*Dz-2*h^2 , 0 , 0 , 2*x*Dx+x*Dz-y*Dz+z*Dz+2*h^2 ]
304: [ 2*y*Dy-2*z*Dy+y*Dz-z*Dz+h^2 , 2*x*Dz-y*Dz+2*z*Dz+h^2 , -x*Dz+z*Dz , 2*x*Dy+x*Dz ]
305: [ -2*y*Dy+2*z*Dy+y*Dz-z*Dz , y*Dz-4*z*Dz , -2*y*Dx+2*z*Dx-y*Dz+2*z*Dz , -2*z*Dy+y*Dz-3*z*Dz ]
306: ]
307: [ 2 $B<!(B
308: [ -2*y*Dx+2*z*Dx-y*Dz+2*z*Dz , x*y-x*z-y*z+z^2 , y-z , y , x+y-z ]
309: [ -6*Dx*Dz-2*Dz^2 , x*Dz+y*Dz-5*z*Dz-4*h^2 , -2*Dy+2*Dz , 2*Dx+2*Dz , 4*Dz ]
310: ]
311: ]
312: In(7)=
313:
314: $BNc(B 3: $B$&$^$/9T$+$J$$Nc(B:
315:
316: Example 1: Sweyl("x,y",[["x",-1,"y",-1,"Dx",1,"Dy",1]]);
317: v=[[2*x*Dx + 3*y*Dy+6, 0],
318: [3*x^2*Dy + 2*y*Dx, 0],
319: [0, x^2+y^2],
320: [0, x*y]];
321: a=Sminimal(v);
322: strategy $B$,$*$+$7$$$H$$$C$F$H$^$k(B. $BM}M3$O(B?
323: Negative weight vector $B$r;H$o$J$$$H$-$A$s$HF0$-$^$9(B.
324:
325:
326: DEBUG $B=PNO(B:
327: rf= [
328: [
329: [ Schreyer frame.
330: [ 0 , y^3 , 0 , 0 , -x^2 , 0 ]
331: [ 0 , 0 , y^2 , 0 , -x , 0 ]
332: [ 0 , y , -x , 0 , 0 , 0 ]
333: [ y*h , 0 , 0 , -x , 0 , 0 ]
334: [ 0 , 0 , 0 , 3*y*Dy , 0 , -2*Dx ]
335: ]
336: [
337: [ 1 , 0 , -y^2 , 0 , 0 ]
338: ]
339: [ ]
340: ]
341: [
342: [ 2*x*Dx , e_*x^2 , e_*x*y , 2*y*Dx*h , e_*y^3 , 3*y^2*Dy*h ]
343: [ es*y^3 , es^2*y^2 , es*y , y*h , 3*es^3*y*Dy ]
344: [ 1 ]
345: ]
346: [
347: [ ]
348: [
349: [
350: [ 1 , 4 ]
351: [ y^3 , -x^2 ]
352: ]
353: [
354: [ 2 , 4 ]
355: [ y^2 , -x ]
356: ]
357: [
358: [ 1 , 2 ]
359: [ y , -x ]
360: ]
361: [
362: [ 0 , 3 ]
363: [ y*h , -x ]
364: ]
365: [
366: [ 3 , 5 ]
367: [ 3*y*Dy , -2*Dx ]
368: ]
369: ]
370: [
371: [
372: [ 0 , 2 ]
373: [ 1 , -y^2 ]
374: ]
375: ]
376: [ ]
377: ]
378: [ resolution $B$9$Y$-(B $BItJ,2C72(B e_ $B$O(B $B%Y%/%H%k@.J,$N%^!<%/(B.
379: [ 2*x*Dx+3*y*Dy+6*h^2 , e_*x^2+e_*y^2 , e_*x*y , 2*y*Dx*h+3*x^2*Dy , e_*y^3 , 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy ]
380: ]
381: ]
382:
383: $BN,(B
384: Processing [ 1 , 3 ] Strategy = 2
385: 1 $B$N(B 3 $BHVL\$N(B spair $B$N(B reduction $B$r=hM}Cf(B.
386: In(7)=reductionTable:
387: [[ 1 , 1 , 1 , 2 , 2 , 3 ] , [ 3 , 2 , 1 , 2 , 3 ] , [ 2 ] ]
388: -- $B$3$l(B.
389: SpairAndReduction:
390: [ p and bases , [ [ 0 , 3 ] , [ y*h , -x ] ] , [ 2*x*Dx+3*y*Dy+6*h^2 , e_*x^2+e_*y^2 , e_*x*y , 2*y*Dx*h+3*x^2*Dy , e_*y^3 , %[null] ] ]
391: 0 $B$N(B 0 $BHVL\$H(B 3 $BHVL\(B $B$N(B spair $B$r7W;;$7$F(B, 0 $B%l%Y%k$N(B gb $B$G(B reduction.
392: [ 1 , 1 , 1 , 2 , 2 , 3 ] $B$K$"$k$h$&$K(B, strategy 3 $B0J30$O7W;;$:$_(B.
393: ( $B7W;;$7$F$J$$$b$N$O(B %[null] $B$H$J$C$F$k(B. )
394: [ level= , 1 ]
395: [ tower2= , [ [ ] ] ] ( $B0lHV2<$J$N$G(B, tower $B$O$J$7$h(B. )
396: [ y*h , -es^3*x ]
397: [gi, gj] = [ 2*x*Dx+3*y*Dy+6*h^2 , 2*y*Dx*h+3*x^2*Dy ]
398: 1
399: Reduce the element 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy
400: by [ 2*x*Dx+3*y*Dy+6*h^2 , e_*x^2+e_*y^2 , e_*x*y , 2*y*Dx*h+3*x^2*Dy , e_*y^3 , %[null] ]
401: result is [ 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy , 1 , [ 0 , 0 , 0 , 0 , 0 , 0 ] ]
402: vdegree of the original = -1
403: vdegree of the remainder = -1
404: [ 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy , [ y*h , 0 , 0 , -x , 0 , 0 ] , 3 , 5 , -1 , -1 ]
405:
406: In(11)=freeRes:
407: [ [ 2*x*Dx+3*y*Dy+6*h^2 , e_*x^2+e_*y^2 , e_*x*y , 2*y*Dx*h+3*x^2*Dy , e_*y^3 , 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy ] , [ %[null] , [ 0 , 0 , y^2 , 0 , -x , 0 ] , [ 0 , -y , x , 0 , 1 , 0 ] , [ -y*h , 0 , 0 , x , 0 , 1 ] , %[null] ] , [ %[null] ] ]
408: $B$r$_$l$P$o$+$k$h$&$K(B, SlaScala $B$G(B, freeRes $B$K$3$N85$,(B [0,5] $B$K2C$((B
409: $B$i$l$?(B.
410:
411: $B<!$K(B SnextI $B$,(B SlaScala $B$h$j8F$P$l$F$3$N%(%i!<(B.
412: i = SnextI(reductionTable_tmp,strategy,redundantTable,
413: skel,level,freeRes);
414: In(22)=reductionTable:
415: [ [ 1 , 1 , 1 , 2 , 2 , 3 ] , [ 3 , 2 , 1 , 2 , 3 ] , [ 2 ] ]
416: $B$J$N$G(B, $B:G8e(B $B$N(B 2 $B$,=hM}$5$l$k$O$:$@$,(B,
417: In(25)=skel[2]:
418: [ [ [ 0 , 2 ] , [ 1 , -y^2 ] ] ]
419: $B$N$h$&$K(B, 0 $BHVL\$H(B, 2 $BHVL\$N(B spair.
420: $B$7$+$7(B,
421: In(26)=bases:
422: [ %[null] , [ 0 , 0 , y^2 , 0 , -x , 0 ] , [ 0 , -y , x , 0 , 1 , 0 ] , [ -y*h , 0 , 0 , x , 0 , 1 ] , %[null] ]
423: $B$N$h$&$K(B, 0 $BHVL\$O(B strategy 3 $B$J$N$G(B, $B$^$@$b$H$^$C$F$$$J$$(B.
424:
425: reductionTable_tmp=[ 2 ]
426: See also reductionTable, strategy, level,i
427: ERROR(sm): error operator : SnextI: bases[i] or bases[j] is null for all combinations.
428: --- Engine error or interrupt : In function : Error of class PrimitiveObject
429:
430: Type in Cleards() to exit the debug mode and Where() to see the stack trace.
431: In(7)=reductionTable:
432: [ [ 1 , 1 , 1 , 2 , 2 , 3 ] , [ 3 , 2 , 1 , 2 , 3 ] , [ 2 ] ]
433: In(8)=strategy:
434: 2
435: In(9)=level:
436: 2
437: $B$3$N;~E@$^$G$G$b$H$^$C$?(B basis
438: [
439: [ 2*x*Dx+3*y*Dy+6*h^2 , e_*x^2+e_*y^2 , e_*x*y , 2*y*Dx*h+3*x^2*Dy , e_*y^3 , 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy ]
440: [ %[null] , [ 0 , 0 , y^2 , 0 , -x , 0 ] , [ 0 , -y , x , 0 , 1 , 0 ] , [ -y*h , 0 , 0 , x , 0 , 1 ] , %[null] ]
441: [ %[null] ]
442: ]
443:
444: -------------------------------------
445:
446: Sweyl("x,y",[["x",-1,"y",-1,"Dx",1,"Dy",1]]);
447: a=Sminimal([x^2+y^2,x*y]);
448: $B$3$l$G$b;w$?$h$&$J%(%i!<$r$@$;$k(B.
449: $B$3$NJ}$,(B debug $B$7$d$9$$(B:
450: Sweyl("x,y",[["x",-1,"y",-1,"Dx",1,"Dy",1]]);
451: a=Sminimal([x*y,x^2+y^2]);
452: $B$G$O%(%i!<$,$G$J$$$N$,IT;W5D(B.
453: pruneZero $B$,F0$$$F$J$$$N$,JQ(B.
454:
455: rf= [
456: [
457: [
458: [ y^3 , 0 , -x^2 ]
459: [ 0 , y^2 , -x ]
460: [ y , -x , 0 ]
461: ]
462: [
463: [ 1 , 0 , -y^2 ]
464: ]
465: [ ]
466: ]
467: [
468: [ x^2 , x*y , y^3 ]
469: [ y^3 , es*y^2 , y ]
470: [ 1 ]
471: ]
472: [
473: [ ]
474: [
475: [
476: [ 0 , 2 ]
477: [ y^3 , -x^2 ]
478: ]
479: [
480: [ 1 , 2 ]
481: [ y^2 , -x ]
482: ]
483: [
484: [ 0 , 1 ]
485: [ y , -x ]
486: ]
487: ]
488: [
489: [
490: [ 0 , 2 ]
491: [ 1 , -y^2 ]
492: ]
493: ]
494: [ ]
495: ]
496: [
497: [ x^2+y^2 , x*y , y^3 ]
498: ]
499: ]
500: [ 0 , 0 ]
501: Processing [ 0 , 0 ] Strategy = 1
502: [ 0 , 1 ]
503: Processing [ 0 , 1 ] Strategy = 1
504: [ 1 , 2 ]
505: Processing [ 1 , 2 ] Strategy = 1
506: SpairAndReduction:
507: [ p and bases , [ [ 0 , 1 ] , [ y , -x ] ] , [ x^2+y^2 , x*y , %[null] ] ]
508: [ level= , 1 ]
509: [ tower2= , [ [ ] ] ]
510: [ y , -es*x ]
511: [gi, gj] = [ x^2+y^2 , x*y ]
512: 1
513: Reduce the element y^3
514: by [ x^2+y^2 , x*y , %[null] ]
515: result is [ y^3 , 1 , [ 0 , 0 , 0 ] ]
516: vdegree of the original = -3
517: vdegree of the remainder = -3
518: [ y^3 , [ y , -x , 0 ] , 2 , 2 , -3 , -3 ]
519: [ 0 , 2 ]
520: Processing [ 0 , 2 ] Strategy = 2
521: [ 1 , 1 ]
522: Processing [ 1 , 1 ] Strategy = 2
523: SpairAndReduction:
524: [ p and bases , [ [ 1 , 2 ] , [ y^2 , -x ] ] , [ x^2+y^2 , x*y , y^3 ] ]
525: [ level= , 1 ]
526: [ tower2= , [ [ ] ] ]
527: [ es*y^2 , -es^2*x ]
528: [gi, gj] = [ x*y , y^3 ]
529: 1
530: Reduce the element 0
531: by [ x^2+y^2 , x*y , y^3 ]
532: result is [ 0 , 1 , [ 0 , 0 , 0 ] ]
533: vdegree of the original = -4
534: vdegree of the remainder = %[null]
535: [ 0 , [ 0 , y^2 , -x ] , 1 , -1 , -4 , %[null] ]
536: reductionTable_tmp=[ 2 ]
537: See also reductionTable, strategy, level,i
538: ERROR(sm): error operator : SnextI: bases[i] or bases[j] is null for all combinations.
539: --- Engine error or interrupt : In function : Error of class PrimitiveObject
540:
541: Type in Cleards() to exit the debug mode and Where() to see the stack trace.
542: In(10)=reductionTable :
543: [ [ 1 , 1 , 2 ] , [ 3 , 2 , 1 ] , [ 2 ] ]
544: In(11)=bases:
545: [ %[null] , [ 0 , y^2 , -x ] , [ -y , x , 1 ] ]
546: In(12)= $B$3$l$O(B, [3, 2, 1] $B$N85$N$&$A(B, 2,1 $B$,$b$H$^$C$F$$$k(B.
1.6 takayama 547: $B:G8e$N(B [ 2 ] $B$N7W;;$K(B 0 $BHVL\$,I,MW$G$3$l$,$^$@$J$$(B.
548: $BMW$9$k$K(B 1 $BHVL\$H(B 3 $BHVL\$r>C$9(B operator [1, 0, -y^2]
549: [ y^3 , 0 , -x^2 ]
550: [ 0 , y^2 , -x ]
551: [ y , -x , 0 ]
552: $B$N(B reduction $B$,I,MW(B.
1.4 takayama 553:
1.5 takayama 554: -----------------------------------------
555: June 11, 2000 (Tue), 20:05
556: V-strict $B$+$I$&$+$r%A%'%C%/$9$k4X?t$r=q$-$?$$(B.
557: $B0BA4$K(B ring (schreyer order) $B$rDj5A$9$k4X?t$,M_$7$$(B.
558: $B0BA4$K(B parse $B$9$k4X?t$bM_$7$$(B.
559: $B%Y%/%H%k$H(B es $BI=8=$NJQ494X?t$b$$$k(B.
560:
561: AvoidTheSameRing == 1 $B$J$i(B, schreyer $B$N(B gbList $B$bJQ99$G$-$J$$$h$&$K(B
562: $B$9$Y$-$+!)(B
563: $B4XO"JQ?t(B:
564: needWarningForAvoidTheSameRing
565: isTheSameRing() : ring $B$,F1$8$+(B check. pointer $B$G$J$/Cf?H$^$G$_$k(B.
566: see poly4.c. $B$3$3$N%3%a%s%H$O;29M$K$J$k(B.
567: 3.If Schreyer = 1, then the system always generates a new ring.
568:
569: define_ring $B$K(B gbList $B$bEO$;$k$N(B?
570: ==> set_up_ring@ $B$r8+$k(B. grep set_up_ring ==>
571: primitive.c KsetUpRing() grep KsetUpRing ==>
572: keyword gbListTower $B$,;H$($k$,(B, list $B$GM?$($J$$$H$$$1$J$$(B.
573: list $B$KJQ49$9$k$N$O(B, (list) dc.
574:
575: tparse $B$NI,MW$J$o$1(B?
576: ?? $B$*$b$$$@$;$J$$(B.
577:
578: ring_def $B$G(B ring (schreyer order) $B$rDj5A$9$k$H(B, $B7W;;$N$H$-$N(B
579: order $B$b(B tower $B$G$d$C$F$/$l$k$N(B?
580: $BB?J,(B NO.
581: grep ppAdd *.c ==>
582: poly2.c
583: checkRing(f,g);
584:
585: while (f != POLYNULL && g != POLYNULL) {
586: /*printf("%s + %s\n",POLYToString(f,'*',1),POLYToString(g,'*',1));*/
587: checkRing2(f,g); /* for debug */
588: gt = (*mmLarger)(f,g);
589:
590: mmLarger $B$OJQ$($F$J$$$h$&$K8+$($k(B. checkRing $B$O%^%/%m(B.
591:
592: mmLarger_tower $B$O(B
593: if (!(f->m->ringp->schreyer) || !(g->m->ringp->schreyer))
594: return(mmLarger_matrix(f,g));
595: $B$H$J$C$F$k$N$G(B mmLarger_tower $B$r(B default $B$K$7$F$*$1$P?4G[$J$$$h$&$K8+$($k(B.
596:
597: ring_def $B$O@5$7$/F0$/(B?
598:
1.6 takayama 599: TODO:
600: $B4X?t$N;EMM(B: ( new.sm1 $B$^$?$O(B complex.sm1 $B$K$*$$$H$/(B )
1.5 takayama 601: mmLarger $B$O(B tower $B$KJQ$($F$7$^$&(B.
602: $BJQ?tL>(B, weight vector, $B%7%U%H%Y%/%H%k(B m $B$rM?$($k$H(B ring (with schreyer order)
603: $B$r:n$k(B. ==> weyl<m>, weyl
1.6 takayama 604: parser $B$O$H$/$K:n$kI,MW$,$J$$$h$&$K8+$($k$,(B...(tparse) ==> name
1.5 takayama 605: $B%Y%/%H%k(B <---> es $BI=8=(B cf. toVectors, [(toe_) f] gbext ==> name
606: $BE,@Z$J(B homogenization $B4X?t(B ==> homogenize<m>
607: ord_w $B$N(B schreyer $BHG(B ==> ord_w<m>
608: init $B$N(B schreyer $BHG(B ==> init<m>
609: gb_h, syz_h $B$NBP1~HG(B ==> [ ii vv ww m] syz_h
610: resolution $B$+$i(B shift vector $B$r7W;;$9$k4X?t(B.
611:
1.6 takayama 612: $B7k2L$N(B check $B$r$9$k(B assert $B4X?t$bI,MW(B.
613:
1.5 takayama 614: $B>e$N(B $B%7%U%H%Y%/%H%kBP1~HG$N4X?t$OEvJ,(B new.sm1 $B$X(B. $B$=$N$"$H(B complex.sm1 $B$X(B.
615:
616: cohom.sm1 $B$N(B interface $B4X?t$O(B cohom.k $B$X(B.
617: Help key word $B$O(B (Cohom.deRham) $B$_$?$$$K(B, . $B$G$/$.$C$F=q$/(B.
1.6 takayama 618:
619: ----------------------
620: $B%(%i!<$N860x$,$h$&$d$/$o$+$k(B: June 14, 19:00
621: Schreyer frame $B$NCJ3,$G(B syz $B$K(B 1 $B$,$"$k$H(B strategy $B$,(B
622: $B$O$?$i$+$J$$(B.
623:
624: test13() GKZ $B$N(B minimal free resolution. 2 $BEY<B9T$9$k$HJQ(B.
625: grade $B$,JQ99$5$l$k$H(B, $BJQ$J$3$H$,$*$-$k$N$G(B,
626: ScheckIfSchreyer() $B4X?t$G(B, $B$3$l$r(B scheck $B$9$k$3$H$K$7$?(B.
627: sm1(" (report) (mmLarger) switch_function /ss set ");
628: $B$O$^$@$d$a$H$/(B. matrix $B$K$J$C$F$k$N$G(B.
629:
630: ------------------------------------------
631: June 15, 2000
632: TODO:
633: 1.if (IdenfityIntegerAndUniversalNumber) $B$N$H$-(B --- default
634: lt, gt, eq $B$G(B integer $B$H(B universalNumber $B$NHf3S$,$G$-$k$h$&$K$9$k(B.
635: rational $B$H$NHf3S$b2DG=$K$9$k(B.
636:
637: 2. sm1_push_int0 $B$KBP1~$9$k$3$H$r(B, sm1 $B$NB&$G$d$k(B.
638: $B%^%/%mL>(B obj to_int --> Done.
639: weight_vector $B$N(B universalNumber ==> $B$^$@(B. $B%(%i!<$r$@$5$J$$$N$,$3$o$$(B.
640: s_weight_vector
641: weightv
642: ord_w
643: toVectors
644: define_ring
645: init
646: gkz
647:
648: -------------
649: Schreyer skelton $B$,$I$&$7$F(B 1 $B$rMWAG$K$b$D$+$7$i$Y$k(B.
1.7 takayama 650:
651: June 24 (Sat), 22:30 at Posthouse (Heathrow) www.posthouse-hotels.com
652: Sevilla $BBZ:_(B, Mega $B$b$h$&$d$/$*$o$j(B minimal resolution $B$N(B check $B$KLa$k(B.
653: resol1.c $B$K<!$N(B line $B$r2C$($?(B.
654: /* If isConstant(sv.a) is added, (x^3 - y^2 z^2) deRham stops
655: with an error. I've not yet understood the reason.
656: At Posthouse at Heathrow. June 24, 2000 */
657: if (isConstant(sv.b)) {
658: s->deleted = 1;
659: }
1.8 takayama 660: ===> $B$*$+$7$$$N$G:o=|(B.
1.7 takayama 661:
662: isConstant(sv.a) $B$,$J$$$H(B, $B$3$s$I$O(B,
663: Sminimal([x^2+y^2,x*y]); $B$,%(%i!<$G$H$^$k(B.
664: (x,y $B$N(B weight $B$O(B -1).
665: LaScala-Stillman $B$NO@J8$r$b$&0lEY$J$,$a$h$&(B.
666:
1.8 takayama 667: commit $B$9$Y$-(B: misc/mega2000 (cvs-misc add) Done.
668: OpenXM/src/kan96xx Done.
669: OpenXM/src/k097/lib/minimal Done.
670:
671: July 26.
672: resol.c $B$N(B schreyerSkelton0 $B$G(B, skelton $B$,(B minimal $B$K$J$k$h$&$K(B
673: $B%3!<%I$rA^F~(B.
674: $B%F%9%H$O(B
675: cd src/k097/lib/minimal
676: k0
677: load["minimal.k"];;
678: Sweyl("x,y",[["x",-1,"y",-1,"Dx",1,"Dy",1]]);
679: Sminimal([x^2+y^2,x*y]);
680: $B$G(B.
681:
682: LaScala-Stillman $B$NO@J85U$G(B i<j $B$J$i(B e_i > e_j $B$H$J$k(B.
683: (order.c mmLarger_tower())
684:
685: $B%F%9%H(B 2.
686: cd src/k097/lib/minimal
687: k0
688: load["minimal-test.k"];;
689: v:
690: Sminimal(v);
691:
692: test11(); /* a = test_ann3("x^3-y^2*z^2"); */
693: test14(); /* gkz (1,2,3) */
694:
1.9 takayama 695: July 30. Removed unnecessary code.
696: $BNc(B:
697: Sminimal("x^3-y^2");
698: test12() ( x^3-y^2 z^2)
699: test15() GKZ 1,2,3 with a check.
700: test15b() toric
1.10 takayama 701: test15c() (u,v) = (-1,1)
702:
703: August 1.
704: (u,v)-minimal $B$N%F%9%H%3!<%I$r$$$l$?(B.
705: IsExact_h $B$G(B $BJQ?t(B c $B$NCM$,$+$o$k(B. $B860xITL@(B.
706: c=Sinit_w(b,w);
707: Println("Resolution (b)----");
708: sm1_pmat(b);
709: Println("Initial (c)----");
710: sm1_pmat(c); cc=c;
711: Println("Exactness of the resolution ---");
712: Println(IsExact_h(b,v)); /* IsExact_h breaks the variable c.
713: THIS BUG SHOULD BE FIXED. */
714: $B$3$N$"$H$J$<$+(B, c $B$,(B b $B$NCM$K$+$o$C$F$7$^$&(B.
715: $B$J$*(B def IsExact(c,...) $B$HDj5A$5$l$F$*$j(B, $B$3$N(B c $B$rJL$NJQ?tL>$K(B
716: $BJQ$($l$P$3$NLdBj$O$*$-$J$$(B.
717: Println("Why is the initial c rewritten by b? (buggy) ");sm1_pmat(c[0]);
718:
719: ===> complex.sm1 $B$N(B isExact_h (isExact) $B$G(B popVariables $B$rK:$l$F$?$@$1(B.
720:
721: betti $B?t$O(B, $B9TNs$N>C5n$r$d$k$^$G$o$+$i$J$$$N(B?
722: SbettiTable().
723:
724: Sminimal $B$O(B [(Homogenize_vec) 0] system_variable $B$K$9$k$h$&$G(B,
725: $B$3$l$,(B cohomology $B$N7W;;$K$O<YKb(B.
1.11 takayama 726:
727: August 2, 2000.
728:
729: Sminimal $B$O(B [(Homogenize_vec) 0] system_variable $B$K$9$k$h$&$G(B,
730: $B$3$l$,(B cohomology $B$N7W;;$K$O<YKb(B.
731: ( cf. $BBg0$5W;a$N%9%/%j%W%H(B. $B8=:_?@8M$KBZ:_Cf(B. )
732:
733: /restoreEnvAfterResolution {
734: [(AvoidTheSameRing)] pushEnv
735: [ [(AvoidTheSameRing) 0] system_variable
736: [(gbListTower) [[ ]] (list) dc] system_variable
737: ] pop popEnv
738: setupEnvForResolution.opts restoreOptions <=== $BJQ99(B. opts $B$O$$$m$s$J$H$3$m$G;H$C$F$k(B.
739: } def
740:
741: $B$3$N%^%/%m$r$h$Y$P$$$$$N$+!)(B
742: sm1(" restoreEnvAfterResolution ");
743: $B$r(B Sminimal $B$N$*$o$j$K8F$V$h$&$KJQ$($?(B.
744: test17b(), test18() $B$O@5>oF0:n(B.
745:
1.10 takayama 746:
1.12 ! takayama 747: August 7, Mon 13:00JST ( 5:00 St.Andrews, Scotland, 4039 $B9f<<(B)
! 748: example-ja.tex $B$r=q$/$?$a$N=PNO(B.
! 749:
! 750: % k0
! 751: sm1>macro package : dr.sm1, 9/26,1995 --- Version 6/15, 2000.
! 752: sm1>macro package : module1.sm1, 1994 -- Nov 8, 1998
! 753: This is kan/k0 Version 1998,12/15
! 754: WARNING: This is an EXPERIMENTAL version
! 755: sm1>var.sm1 : Version 3/7, 1997
! 756:
! 757:
! 758: In(1)=Loading startup files (startup.k) 1997, 3/11.
! 759: sm1 version = 3.000726
! 760: Default ring is Z[x,h].
! 761: WARNING(sm): You rewrited the protected symbol pushVariables.
! 762: WARNING(sm): You rewrited the protected symbol popVariables.
! 763: In(2)=load["minimal-test.k"];;
! 764: cpp: -lang-c++: linker input file unused since linking not done
! 765: cpp: -lang-c++: linker input file unused since linking not done
! 766: cohom.sm1 is the top of an experimental package to compute restrictions
! 767: of all degrees based on restall.sm1 and restall_s.sm1
! 768: See, http://www.math.kobe-u.ac.jp to get these files of the latest version.
! 769: Note that the package b-function.sm1 cannot be used with this package.
! 770: r-interface.sm1 (C) N.Takayama, restriction, deRham
! 771:
! 772: oxasir.sm1, --- open asir protocol module 3/1 1998, 6/5 1999
! 773: asirconnect, asir, fctr, primadec, (C) M.Noro, N.Takayama
! 774: ox.sm1, --- open sm1 protocol module 11/11,1999 (C) N.Takayama. oxhelp for help
! 775: hol.sm1, basic package for holonomic systems (C) N.Takayama, 2000, 06/08
! 776: rank characteristic ch rrank gb pgb syz genericAnn annfs gb_h syz_h isSameIdeal isSameIdeal_h
! 777: sm1>gkz.sm1 generates gkz systems (C) N.Takayama, 1998, 11/8, cf. rrank in hol.sm1
! 778: gkz
! 779: sm1>appell.sm1 generates Appell hypergeometric differential equations (C) N.Takayama, 1998, 11/8, cf. rank in hol.sm1
! 780: appell1 appell4
! 781: sm1>resol0.sm1, package to construct schreyer resolutions -- not minimal
! 782: (C) N.Takayama, 1999, 5/18. resol0, resol1
! 783: complex.sm1 : 1999, 9/28, res-div, res-solv, res-kernel-image, res-dual
! 784: 2000, 6/8, isExact_h, isExact
! 785: In this package, complex is expressed in terms of matrices.
! 786: restall.sm1 ... compute all the cohomology groups of the restriction
! 787: of a D-module to tt = (t_1,...,t_d) = (0,...,0).
! 788: non-Schreyer Version: 19980415 by T.Oaku
! 789: usage: [(P1)...] [(t1)...] bfm --> the b-function
! 790: [(P1)...] [(t1)...] k0 k1 deg restall --> cohomologies of restriction
! 791: [(P1)...] [(t1)...] intbfm --> the b-function for integration
! 792: [(P1)...] [(t1)...] k0 k1 deg intall --> cohomologies of integration
! 793: restall_s.sm1...compute all the cohomology groups of the restriction
! 794: of a D-module to tt = (t_1,...,t_d) = (0,...,0).
! 795: Schreyer Version: 19990521 by N.Takayama & T.Oaku
! 796: usage: [(P1)...] [(t1)...] k0 k1 deg restall_s -> cohomologies of restriction
! 797: [(P1)...] [(t1)...] k0 k1 deg intall_s --> cohomologies of integration
! 798: No truncation from below in restall
! 799: The variable Schreyer is set to 2.
! 800: Loading tower.sm1 in the standard context. You cannot use Schyrer 1. It is controlled from cohom.sm1
! 801:
! 802: oxpath.oxlog.xterm is set to /home/nobuki/OpenXM/lib/sm1/bin/oxlog
! 803: In(3)=a=Sannfs2("x^3-y^2");
! 804: Starting ox_asir server.
! 805: Hello from open. serverName is localhost and portnumber is 0
! 806: Done the initialization. port =1024
! 807: Hello from open. serverName is localhost and portnumber is 0
! 808: Done the initialization. port =1025
! 809: [ 7 , 1025 , 6 , 1024 ]
! 810: [1] 250
! 811: Trying to accept from localhost... len= 16
! 812: 4 2 7f 0 0 1 0 0 0 0 0 0 0 0 8 0
! 813: Authentification: localhost is allowed to be accepted.
! 814: Accepted.
! 815: Trying to accept from localhost... len= 16
! 816: 4 3 7f 0 0 1 0 0 0 0 0 0 0 0 6 0
! 817: Authentification: localhost is allowed to be accepted.
! 818: Accepted.
! 819:
! 820: Control port 1024 : Connected.
! 821:
! 822: Stream port 1025 : Connected.
! 823: Byte order for control process is network byte order.
! 824: Byte order for engine process is network byte order.
! 825: WeightOfSweyl=[ x , -1 , y , -1 , Dx , 1 , Dy , 1 ]
! 826: Automatic homogenization.
! 827: [ -2*x*Dx , -3*y*Dx^2 , -9*y^2*Dx*Dy , -27*y^3*Dy^2 ]
! 828: Warning: Homogenization and ReduceLowerTerms options are automatically turned off.
! 829: ....Done. betti=4
! 830: Warning: Homogenization and ReduceLowerTerms options are automatically turned ON.
! 831: Warning: Homogenization and ReduceLowerTerms options are automatically turned off.
! 832: .Done. betti=1
! 833: Warning: Homogenization and ReduceLowerTerms options are automatically turned ON.
! 834: Warning: Homogenization and ReduceLowerTerms options are automatically turned off.
! 835: Done. betti=0
! 836: Warning: Homogenization and ReduceLowerTerms options are automatically turned ON.
! 837: rf= [
! 838: [
! 839: [
! 840: [ -9*y^2*Dy , 0 , 2*x , 0 ]
! 841: [ 0 , 0 , -3*y*Dy , Dx ]
! 842: [ 0 , -3*y*Dy , Dx , 0 ]
! 843: [ -3*y*Dx , 2*x , 0 , 0 ]
! 844: ]
! 845: [
! 846: [ -Dx , 0 , 0 , 3*y*Dy ]
! 847: ]
! 848: [ ]
! 849: ]
! 850: [
! 851: [ -2*x*Dx , -3*y*Dx^2 , -9*y^2*Dx*Dy , -27*y^3*Dy^2 ]
! 852: [ -9*y^2*Dy , -3*es^2*y*Dy , -3*es*y*Dy , -3*y*Dx ]
! 853: [ -Dx ]
! 854: ]
! 855: [
! 856: [ ]
! 857: [
! 858: [
! 859: [ 0 , 2 ]
! 860: [ -9*y^2*Dy , 2*x ]
! 861: ]
! 862: [
! 863: [ 2 , 3 ]
! 864: [ -3*y*Dy , Dx ]
! 865: ]
! 866: [
! 867: [ 1 , 2 ]
! 868: [ -3*y*Dy , Dx ]
! 869: ]
! 870: [
! 871: [ 0 , 1 ]
! 872: [ -3*y*Dx , 2*x ]
! 873: ]
! 874: ]
! 875: [
! 876: [
! 877: [ 0 , 3 ]
! 878: [ -Dx , 3*y*Dy ]
! 879: ]
! 880: ]
! 881: [ ]
! 882: ]
! 883: [
! 884: [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h , -9*y^2*Dx*Dy-3*y*Dx*h^2-4*x^2*Dy*h , -27*y^3*Dy^2-27*y^2*Dy*h^2+3*y*h^4+8*x^3*Dy*h ]
! 885: ]
! 886: ]
! 887: Generating reduction table which gives an order of reduction.
! 888: WeghtOfSweyl=[ x , -1 , y , -1 , Dx , 1 , Dy , 1 ]
! 889: tower[ [ -2*x*Dx , -3*y*Dx^2 , -9*y^2*Dx*Dy , -27*y^3*Dy^2 ] , [ -9*y^2*Dy , -3*es^2*y*Dy , -3*es*y*Dy , -3*y*Dx ] , [ -Dx ] ]
! 890: reductionTable= [
! 891: [ 1 , 2 , 3 , 4 ]
! 892: [ 3 , 4 , 3 , 2 ]
! 893: [ 3 ]
! 894: ]
! 895: [ 0 , 0 ]
! 896: Processing [level,i]= [ 0 , 0 ] Strategy = 1
! 897: [ 0 , 1 ]
! 898: Processing [level,i]= [ 0 , 1 ] Strategy = 2
! 899: [ 1 , 3 ]
! 900: Processing [level,i]= [ 1 , 3 ] Strategy = 2
! 901: SpairAndReduction:
! 902: [ p and bases , [ [ 0 , 1 ] , [ -3*y*Dx , 2*x ] ] , [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h , %[null] , %[null] ] ]
! 903: [ level= , 1 ]
! 904: [ tower2= , [ [ ] ] ]
! 905: [ -3*y*Dx , 2*es*x ]
! 906: [gi, gj] = [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h ]
! 907: 1
! 908: Reduce the element 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h
! 909: by [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h , %[null] , %[null] ]
! 910: result is [ 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h , 1 , [ 0 , 0 , 0 , 0 ] ]
! 911: vdegree of the original = 0
! 912: vdegree of the remainder = 0
! 913: [ 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h , [ -3*y*Dx , 2*x , 0 , 0 ] , 3 , 2 , 0 , 0 ]
! 914: [ 0 , 2 ]
! 915: Processing [level,i]= [ 0 , 2 ] Strategy = 3
! 916: [ 1 , 0 ]
! 917: Processing [level,i]= [ 1 , 0 ] Strategy = 3
! 918: SpairAndReduction:
! 919: [ p and bases , [ [ 0 , 2 ] , [ -9*y^2*Dy , 2*x ] ] , [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h , 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h , %[null] ] ]
! 920: [ level= , 1 ]
! 921: [ tower2= , [ [ ] ] ]
! 922: [ 9*y^2*Dy , 2*es^2*x ]
! 923: [gi, gj] = [ -2*x*Dx-3*y*Dy+h^2 , 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h ]
! 924: 1
! 925: Reduce the element -27*y^3*Dy^2+6*x*y*Dx*h^2-18*y^2*Dy*h^2+8*x^3*Dy*h
! 926: by [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h , 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h , %[null] ]
! 927: result is [ 27*y^3*Dy^2+27*y^2*Dy*h^2-3*y*h^4-8*x^3*Dy*h , -1 , [ -3*y*h^2 , 0 , 0 , 0 ] ]
! 928: vdegree of the original = -1
! 929: vdegree of the remainder = -1
! 930: [ 27*y^3*Dy^2+27*y^2*Dy*h^2-3*y*h^4-8*x^3*Dy*h , [ -9*y^2*Dy-3*y*h^2 , 0 , -2*x , 0 ] , 0 , 3 , -1 , -1 ]
! 931: [ 1 , 2 ]
! 932: Processing [level,i]= [ 1 , 2 ] Strategy = 3
! 933: SpairAndReduction:
! 934: [ p and bases , [ [ 1 , 2 ] , [ -3*y*Dy , Dx ] ] , [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h , 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h , 27*y^3*Dy^2+27*y^2*Dy*h^2-3*y*h^4-8*x^3*Dy*h ] ]
! 935: [ level= , 1 ]
! 936: [ tower2= , [ [ ] ] ]
! 937: [ 3*es*y*Dy , es^2*Dx ]
! 938: [gi, gj] = [ -3*y*Dx^2+2*x*Dy*h , 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h ]
! 939: 1
! 940: Reduce the element -6*y*Dx^2*h^2+4*x^2*Dx*Dy*h+6*x*y*Dy^2*h+8*x*Dy*h^3
! 941: by [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h , 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h , 27*y^3*Dy^2+27*y^2*Dy*h^2-3*y*h^4-8*x^3*Dy*h ]
! 942: result is [ 0 , 1 , [ 2*x*Dy*h , -2*h^2 , 0 , 0 ] ]
! 943: vdegree of the original = 1
! 944: vdegree of the remainder = %[null]
! 945: [ 0 , [ 2*x*Dy*h , 3*y*Dy-2*h^2 , Dx , 0 ] , 2 , -1 , 1 , %[null] ]
! 946: [ 2 , 0 ]
! 947: Processing [level,i]= [ 2 , 0 ] Strategy = 3
! 948: SpairAndReduction:
! 949: [ p and bases , [ [ 0 , 3 ] , [ -Dx , 3*y*Dy ] ] , [ [ 9*y^2*Dy+3*y*h^2 , 0 , 2*x , 1 ] , %[null] , [ 2*x*Dy*h , 3*y*Dy-2*h^2 , Dx , 0 ] , [ 3*y*Dx , -2*x , 1 , 0 ] ] ]
! 950: [ level= , 2 ]
! 951: [ tower2= , [ [ -2*x*Dx , -3*y*Dx^2 , -9*y^2*Dx*Dy , -27*y^3*Dy^2 ] ] ]
! 952: [ Dx , -3*es^3*y*Dy ]
! 953: [gi, gj] = [ 9*y^2*Dy+2*es^2*x+es^3+3*y*h^2 , 3*y*Dx-2*es*x+es^2 ]
! 954: 1
! 955: Reduce the element 6*es*x*y*Dy+2*es^2*x*Dx-3*es^2*y*Dy+es^3*Dx-6*y*Dx*h^2+2*es^2*h^2
! 956: by [ [ 9*y^2*Dy+3*y*h^2 , 0 , 2*x , 1 ] , %[null] , [ 2*x*Dy*h , 3*y*Dy-2*h^2 , Dx , 0 ] , [ 3*y*Dx , -2*x , 1 , 0 ] ]
! 957: result is [ -3*es^2*y*Dy+es^3*Dx+4*es^2*h^2-4*x^2*Dy*h , 1 , [ 0 , 0 , -2*x , 2*h^2 ] ]
! 958: vdegree of the original = 0
! 959: vdegree of the remainder = 0
! 960: [ -3*es^2*y*Dy+es^3*Dx+4*es^2*h^2-4*x^2*Dy*h , [ Dx , 0 , -2*x , -3*y*Dy+2*h^2 ] , 0 , 1 , 0 , 0 ]
! 961: [ 0 , 3 ]
! 962: Processing [level,i]= [ 0 , 3 ] Strategy = 4
! 963: [ 1 , 1 ]
! 964: Processing [level,i]= [ 1 , 1 ] Strategy = 4
! 965: Betti numbers are ------
! 966: [ 2 , 1 , 0 ]
! 967: [seq,level,q]=[ 3 , 1 , 1 ]
! 968: [ level, q = , 1 , 1 ]
! 969: bases=
! 970: [
! 971: [ -Dx , 1 , 2*x , 3*y*Dy-2*h^2 ]
! 972: ]
! 973: dr=
! 974: [ Dx , -1 , -2*x , -3*y*Dy+2*h^2 ]
! 975: newbases=
! 976: [
! 977: [ 0 , 0 , 0 , 0 ]
! 978: ]
! 979: [seq,level,q]=[ 2 , 0 , 3 ]
! 980: [ level, q = , 0 , 3 ]
! 981: bases=
! 982: [
! 983: [ 9*y^2*Dy+3*y*h^2 , 0 , 2*x , 1 ]
! 984: [ -4*x^2*Dy*h , 0 , -3*y*Dy+4*h^2 , Dx ]
! 985: [ 2*x*Dy*h , 3*y*Dy-2*h^2 , Dx , 0 ]
! 986: [ 3*y*Dx , -2*x , 1 , 0 ]
! 987: ]
! 988: dr=
! 989: [ -9*y^2*Dy-3*y*h^2 , 0 , -2*x , -1 ]
! 990: newbases=
! 991: [
! 992: [ 0 , 0 , 0 , 0 ]
! 993: [ -9*y^2*Dx*Dy-3*y*Dx*h^2-4*x^2*Dy*h , 0 , -2*x*Dx-3*y*Dy+2*h^2 , 0 ]
! 994: [ 2*x*Dy*h , 3*y*Dy-2*h^2 , Dx , 0 ]
! 995: [ 3*y*Dx , -2*x , 1 , 0 ]
! 996: ]
! 997: [seq,level,q]=[ 1 , 0 , 2 ]
! 998: [ level, q = , 0 , 2 ]
! 999: bases=
! 1000: [
! 1001: [ 0 , 0 , 0 , 0 ]
! 1002: [ -9*y^2*Dx*Dy-3*y*Dx*h^2-4*x^2*Dy*h , 0 , -2*x*Dx-3*y*Dy+2*h^2 , 0 ]
! 1003: [ 2*x*Dy*h , 3*y*Dy-2*h^2 , Dx , 0 ]
! 1004: [ 3*y*Dx , -2*x , 1 , 0 ]
! 1005: ]
! 1006: dr=
! 1007: [ -3*y*Dx , 2*x , -1 , 0 ]
! 1008: newbases=
! 1009: [
! 1010: [ 0 , 0 , 0 , 0 ]
! 1011: [ 6*x*y*Dx^2-4*x^2*Dy*h , -4*x^2*Dx-6*x*y*Dy , 0 , 0 ]
! 1012: [ -3*y*Dx^2+2*x*Dy*h , 2*x*Dx+3*y*Dy , 0 , 0 ]
! 1013: [ 0 , 0 , 0 , 0 ]
! 1014: ]
! 1015: [ level= , 0 ]
! 1016: [
! 1017: [ -2*x*Dx-3*y*Dy+h^2 ]
! 1018: [ -3*y*Dx^2+2*x*Dy*h ]
! 1019: ]
! 1020: [
! 1021: [ -2*x*Dx-3*y*Dy+h^2 ]
! 1022: [ -3*y*Dx^2+2*x*Dy*h ]
! 1023: ]
! 1024: [ level= , 1 ]
! 1025: [
! 1026: [ 0 , 0 , 0 , 0 ]
! 1027: [ -3*y*Dx^2+2*x*Dy*h , 2*x*Dx+3*y*Dy , 0 , 0 ]
! 1028: [ 0 , 0 , 0 , 0 ]
! 1029: ]
! 1030: [
! 1031: [ 0 , 0 ]
! 1032: [ -3*y*Dx^2+2*x*Dy*h , 2*x*Dx+3*y*Dy ]
! 1033: [ 0 , 0 ]
! 1034: ]
! 1035: [ level= , 2 ]
! 1036: [
! 1037: [ 0 , 0 , 0 , 0 ]
! 1038: ]
! 1039: [
! 1040: [ 0 , 0 , 0 ]
! 1041: ]
! 1042: ------------ Note -----------------------------
! 1043: To get shift vectors, use Reparse and SgetShifts(resmat,w)
! 1044: To get initial of the complex, use Reparse and Sinit_w(resmat,w)
! 1045: 0: minimal resolution, 3: Schreyer resolution
! 1046: ------------ Resolution Summary --------------
! 1047: Betti numbers : [ 2 , 1 ]
! 1048: Betti numbers of the Schreyer frame: [ 4 , 4 , 1 ]
! 1049: -----------------------------------------------
! 1050: In(4)=sm1_pmat(a);
! 1051: [
! 1052: [
! 1053: [
! 1054: [ -2*x*Dx-3*y*Dy+h^2 ]
! 1055: [ -3*y*Dx^2+2*x*Dy*h ]
! 1056: ]
! 1057: [
! 1058: [ -3*y*Dx^2+2*x*Dy*h , 2*x*Dx+3*y*Dy ]
! 1059: ]
! 1060: ]
! 1061: [
! 1062: [
! 1063: [ -2*x*Dx-3*y*Dy+h^2 ]
! 1064: [ -3*y*Dx^2+2*x*Dy*h ]
! 1065: ]
! 1066: [
! 1067: [ 0 , 0 ]
! 1068: [ -3*y*Dx^2+2*x*Dy*h , 2*x*Dx+3*y*Dy ]
! 1069: [ 0 , 0 ]
! 1070: ]
! 1071: [
! 1072: [ 0 , 0 , 0 ]
! 1073: ]
! 1074: ]
! 1075: [
! 1076: [
! 1077: [
! 1078: [ -2*x*Dx-3*y*Dy+h^2 ]
! 1079: [ -3*y*Dx^2+2*x*Dy*h ]
! 1080: [ 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h ]
! 1081: [ 27*y^3*Dy^2+27*y^2*Dy*h^2-3*y*h^4-8*x^3*Dy*h ]
! 1082: ]
! 1083: [
! 1084: [ 0 , 0 , 0 , 0 ]
! 1085: [ 6*x*y*Dx^2-4*x^2*Dy*h , -4*x^2*Dx-6*x*y*Dy , 0 , 0 ]
! 1086: [ -3*y*Dx^2+2*x*Dy*h , 2*x*Dx+3*y*Dy , 0 , 0 ]
! 1087: [ 0 , 0 , 0 , 0 ]
! 1088: ]
! 1089: [
! 1090: [ 0 , 0 , 0 , 0 ]
! 1091: ]
! 1092: ]
! 1093: [
! 1094: [ 0 , 0 , 1 , 2 ]
! 1095: [ 0 , 3 , 0 , 0 ]
! 1096: [ 0 ]
! 1097: ]
! 1098: [
! 1099: [ %[null] , %[null] , [ -3*y*Dx , 2*x , -1 , 0 ] , [ -9*y^2*Dy-3*y*h^2 , 0 , -2*x , -1 ] ]
! 1100: [ %[null] , [ Dx , -1 , -2*x , -3*y*Dy+2*h^2 ] , %[null] , %[null] ]
! 1101: [ %[null] ]
! 1102: ]
! 1103: [ 1 , 4 , 4 , 1 ]
! 1104: [
! 1105: [ 0 , 0 , 1 , 2 ]
! 1106: [ 0 , 3 , %[null] , 0 ]
! 1107: [ 0 ]
! 1108: ]
! 1109: ]
! 1110: [
! 1111: [
! 1112: [ -2*x*Dx-3*y*Dy+h^2 ]
! 1113: [ -3*y*Dx^2+2*x*Dy*h ]
! 1114: [ 9*y^2*Dx*Dy+3*y*Dx*h^2+4*x^2*Dy*h ]
! 1115: [ 27*y^3*Dy^2+27*y^2*Dy*h^2-3*y*h^4-8*x^3*Dy*h ]
! 1116: ]
! 1117: [
! 1118: [ 9*y^2*Dy+3*y*h^2 , 0 , 2*x , 1 ]
! 1119: [ -4*x^2*Dy*h , 0 , -3*y*Dy+4*h^2 , Dx ]
! 1120: [ 2*x*Dy*h , 3*y*Dy-2*h^2 , Dx , 0 ]
! 1121: [ 3*y*Dx , -2*x , 1 , 0 ]
! 1122: ]
! 1123: [
! 1124: [ -Dx , 1 , 2*x , 3*y*Dy-2*h^2 ]
! 1125: ]
! 1126: ]
! 1127: [
! 1128: [
! 1129: [
! 1130: [ -9*y^2*Dy , 0 , 2*x , 0 ]
! 1131: [ 0 , 0 , -3*y*Dy , Dx ]
! 1132: [ 0 , -3*y*Dy , Dx , 0 ]
! 1133: [ -3*y*Dx , 2*x , 0 , 0 ]
! 1134: ]
! 1135: [
! 1136: [ -Dx , 0 , 0 , 3*y*Dy ]
! 1137: ]
! 1138: [ ]
! 1139: ]
! 1140: [
! 1141: [ -2*x*Dx , -3*y*Dx^2 , -9*y^2*Dx*Dy , -27*y^3*Dy^2 ]
! 1142: [ -9*y^2*Dy , -3*es^2*y*Dy , -3*es*y*Dy , -3*y*Dx ]
! 1143: [ -Dx ]
! 1144: ]
! 1145: [
! 1146: [ ]
! 1147: [
! 1148: [
! 1149: [ 0 , 2 ]
! 1150: [ -9*y^2*Dy , 2*x ]
! 1151: ]
! 1152: [
! 1153: [ 2 , 3 ]
! 1154: [ -3*y*Dy , Dx ]
! 1155: ]
! 1156: [
! 1157: [ 1 , 2 ]
! 1158: [ -3*y*Dy , Dx ]
! 1159: ]
! 1160: [
! 1161: [ 0 , 1 ]
! 1162: [ -3*y*Dx , 2*x ]
! 1163: ]
! 1164: ]
! 1165: [
! 1166: [
! 1167: [ 0 , 3 ]
! 1168: [ -Dx , 3*y*Dy ]
! 1169: ]
! 1170: ]
! 1171: [ ]
! 1172: ]
! 1173: [
! 1174: [ -2*x*Dx-3*y*Dy+h^2 , -3*y*Dx^2+2*x*Dy*h , -9*y^2*Dx*Dy-3*y*Dx*h^2-4*x^2*Dy*h , -27*y^3*Dy^2-27*y^2*Dy*h^2+3*y*h^4+8*x^3*Dy*h ]
! 1175: ]
! 1176: ]
! 1177: ]
! 1178: In(5)=
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to minimal-note-ja.txt CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM / src / k097 / lib / minimal