=================================================================== RCS file: /home/cvs/OpenXM/src/k097/lib/minimal/minimal-note-ja.txt,v retrieving revision 1.3 retrieving revision 1.4 diff -u -p -r1.3 -r1.4 --- OpenXM/src/k097/lib/minimal/minimal-note-ja.txt 2000/06/08 08:37:53 1.3 +++ OpenXM/src/k097/lib/minimal/minimal-note-ja.txt 2000/06/09 08:04:54 1.4 @@ -1,4 +1,4 @@ -$OpenXM: OpenXM/src/k097/lib/minimal/minimal-note-ja.txt,v 1.2 2000/05/24 15:24:54 takayama Exp $ +$OpenXM: OpenXM/src/k097/lib/minimal/minimal-note-ja.txt,v 1.3 2000/06/08 08:37:53 takayama Exp $ SpairAndReduction() : $BM?$($i$l$?(B pair $B$r(B reduction $B$9$k(B. @@ -77,6 +77,11 @@ test8() $B$G(B sm1 $B$G=q$$$?J}$N(B Schreyer $B$r kernel = image $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. $BMW$9$k$K(B k0 $B$N%3!<%I$,$I$&$d$i$*$+$7$$$i$7$$(B. +==> +6/8 $B$N%N!<%H$h$j(B. +syzygy $B$r(B homogenization $B$r2p$7$F7W;;$9$k$N$OLdBj$"$j(B. +--> usage of isExact +$BMW$9$k$K(B kernel = image $B$N%3!<%I$bJQ(B. Homogenized $B$N$^$^$d$kI,MW$"$j(B. ----------------------------------- June 8, 2000 (Thu), 9:10 (Spain local time) @@ -100,5 +105,447 @@ test10(); LaScala-Stillman $B$NJ}K!$G$D$/$C$?(B, schreyer resol $B$,(B exact $B$+(B $BD4$Y$k(B. $BNcBj$O(B, ann(1/(x^3-y^2 z^2)) $B$N(B Laplace $BJQ49(B. + ==> OK. IsExact_h $B$G$7$i$Y$k(B. (IsExact $B$O$@$a$h(B) - \ No newline at end of file +June 8, 2000 (Thu), 19:35 +load["minimal-test.k"];; +test11(); + LaScala-Stillman $B$NJ}K!$G$D$/$C$?(B, minimal resol $B$,(B exact $B$+(B + $BD4$Y$k(B. + $BNcBj$O(B, ann(1/(x^3-y^2 z^2)) $B$N(B Laplace $BJQ49(B. + +SwhereInTower $B$r;H$&$H$-$O(B, +SsetTower() $B$G(B gbList $B$rJQ99$7$J$$$H$$$1$J$$(B. +$B$b$A$m$s;HMQ$7$?$i(B, $B$=$l$rLa$9$3$H(B. +SpairAndReduction, SpairAndReduction2 $B$G(B, + SsetTower(StowerOf(tower,level)); + pos = SwhereInTower(syzHead,tower[level]); + + SsetTower(StowerOf(tower,level-1)); + pos2 = SwhereInTower(tmp[0],tower[level-1]); +$B$H(B, SwhereInTower $B$NA0$K(B setTower $B$r$/$o$($?(B. +( $B0c$&%l%Y%k$G$NHf3S$N$?$a(B.) + +IsExact_h $B$O(B, 0 $B%Y%/%H%k$r4^$`>l9g(B, $B$?$@$7$/F0:n$7$J$$$h$&$@(B. +test11(). +test11a() $B$G(B, 0 $B%Y%/%H%k$r OK. + + +--------------------------------- +June 9, 6:20 +SpairAndReduction +$B$H(B +SpairAndReduction2 +$B$N0c$$(B. +SpairAndReduction : SlaScala (LaScala-Stillman's algorithm $B$G;H$&(B) +SpairAndReduction2 : Sschreyer (schreyer algorithm $B$G;H$&(B, laScala $B$O$J$7(B.) + +0 $B$r<+F0$G=|$/%3!<%I$r=q$3$&(B. + +SpruneZeroRow() $B$r(B Sminimal() $B$K2C$($?(B. +test11() $B$b@5$7$/F0:n$9$k$O$:(B. +IsExact_h $B$O(B schreyer $B$r(B off $B$7$F(B, ReParse $B$7$F$+$i(B, +$B8F$S=P$9$3$H(B. + + +#ifdef TOTAL_STRATEGY + return(StotalDegree(f)+Sdegree(tower[level-2,i],tower,level-1)); +#endif + /* Strategy must be compatible with ordering. */ + /* Weight vector must be non-negative, too. */ + /* See Sdegree, SgenerateTable, reductionTable. */ + wd = Sord_w(f,ww); + return(wd+Sdegree(tower[level-2,i],tower,level-1)); +TOTAL_STRATEGY $B$rMQ$$$kI,MW$,$"$k$N$G$O(B?? +Example 1: Sweyl("x,y",[["x",-1,"y",-1,"Dx",1,"Dy",1]]); + v=[[2*x*Dx + 3*y*Dy+6, 0], + [3*x^2*Dy + 2*y*Dx, 0], + [0, x^2+y^2], + [0, x*y]]; + a=Sminimal(v); +strategy $B$,$*$+$7$$$H$$$C$F$H$^$k(B. $BM}M3$O(B? + +a=test_ann3("x^3+y^3+z^3); $B$O;~4V$,$+$+$j$=$&(B. +a=test_ann3("x^3+y^3"); OK. +a=test_ann3("x^2+y^2+z"); OK. + + +$B>e$N(B example 1 $B$N%(%i!<(B $B$N8+J}(B: +Processing [ 1 , 3 ] Strategy = 2 + 1 $B$N(B 3 $BHVL\$N(B spair $B$N(B reduction $B$r=hM}Cf(B. + In(7)=reductionTable: + [[ 1 , 1 , 1 , 2 , 2 , 3 ] , [ 3 , 2 , 1 , 2 , 3 ] , [ 2 ] ] + -- $B$3$l(B. +SpairAndReduction: +[ 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] ] ] +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. +[ 1 , 1 , 1 , 2 , 2 , 3 ] $B$K$"$k$h$&$K(B, strategy 3 $B0J30$O7W;;$:$_(B. +( $B7W;;$7$F$J$$$b$N$O(B %[null] $B$H$J$C$F$k(B. ) +[ level= , 1 ] +[ tower2= , [ [ ] ] ] ( $B0lHV2<$J$N$G(B, tower $B$O$J$7$h(B. ) +[ y*h , -es^3*x ] +[gi, gj] = [ 2*x*Dx+3*y*Dy+6*h^2 , 2*y*Dx*h+3*x^2*Dy ] +1 +Reduce the element 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy +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] ] +result is [ 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy , 1 , [ 0 , 0 , 0 , 0 , 0 , 0 ] ] +vdegree of the original = -1 +vdegree of the remainder = -1 +[ 3*y^2*Dy*h+6*y*h^3-3*x^3*Dy , [ y*h , 0 , 0 , -x , 0 , 0 ] , 3 , 5 , -1 , -1 ] + +In(11)=freeRes: +[ [ 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] ] ] +$B$r$_$l$P$o$+$k$h$&$K(B, SlaScala $B$G(B, freeRes $B$K$3$N85$,(B [0,5] $B$K2C$((B +$B$i$l$?(B. + +$B$7$?7k2L(B, +Sminimal $B$O$&$^$/$&$4$1$P@5$7$$Ez$($r$@$7$F$k$_$?$$$G$9$,(B +(D : homogenized Weyl $B$G(B ker = im $B$r(B check $B$7$F$k(B, + V-adapted (strict) $B$+$I$&$+$N(B check routing $B$O$^$@=q$$$F$J$$(B), +strategy $B$,$&$^$/$&$4$+$J$/$F$H$^$k>l9g$b$"$j$^$9(B +( strategy = 2 $B$N(B sp $B$r7W;;$9$k$N$K(B, strategy 3 $B$N(B $B85$rI,MW$H(B + $B$7$?$j$9$k>l9g$"$j(B). + + +strategy $B$O(B +def Sdegree(f,tower,level) { + local i,ww, wd; + /* extern WeightOfSweyl; */ + ww = WeightOfSweyl; + f = Init(f); + if (level <= 1) return(StotalDegree(f)); + i = Degree(f,es); + return(StotalDegree(f)+Sdegree(tower[level-2,i],tower,level-1)); +} +$B$rMQ$$$F(B, + ans_at_each_floor[j] = Sdegree(tower[i,j],tower,i+1)-(i+1) +$B$G7W;;$7$F$^$9(B. + +$B$$$/$D$+=PNO$r$D$1$F$*$-$^$9$N$G(B, $B8!F$(B!!! + +$BNc(B 1: +load["minimal-test.k"];; +a=test_ann3("x^3-y^2*z^2"); $B0z?t$N(B annihilating ideal $B$N(B laplace $BJQ49$N(B + homogenization $B$N(B resolution. + weight vector $B$O(B (-1,-1,-1,1,1,1) + +In(4)=sm1_pmat(a[1]); + [ + [ 0 $B