=================================================================== RCS file: /home/cvs/OpenXM/src/asir-doc/parts/groebner.texi,v retrieving revision 1.12 retrieving revision 1.14 diff -u -p -r1.12 -r1.14 --- OpenXM/src/asir-doc/parts/groebner.texi 2003/12/27 11:52:07 1.12 +++ OpenXM/src/asir-doc/parts/groebner.texi 2004/09/14 01:32:34 1.14 @@ -1,4 +1,4 @@ -@comment $OpenXM: OpenXM/src/asir-doc/parts/groebner.texi,v 1.11 2003/04/28 06:43:10 noro Exp $ +@comment $OpenXM: OpenXM/src/asir-doc/parts/groebner.texi,v 1.13 2004/09/13 09:23:30 noro Exp $ \BJP @node グレブナ基底の計算,,, Top @chapter グレブナ基底の計算 @@ -15,6 +15,7 @@ * 基本的な函数:: * 計算および表示の制御:: * 項順序の設定:: +* Weight:: * 有理式を係数とするグレブナ基底計算:: * 基底変換:: * Weyl 代数:: @@ -26,6 +27,7 @@ * Fundamental functions:: * Controlling Groebner basis computations:: * Setting term orderings:: +* Weight:: * Groebner basis computation with rational function coefficients:: * Change of ordering:: * Weyl algebra:: @@ -1052,6 +1054,139 @@ expressed by variable @code{x}, and the above explanat such a drastic experimental results. In practice, however, optimum ordering for variables may not known beforehand, and some heuristic trial may be inevitable. +\E + +\BJP +@node Weight ,,, グレブナ基底の計算 +@section Weight +\E +\BEG +@node Weight,,, Groebner basis computation +@section Weight +\E +\BJP +前節で紹介した項順序は, 各変数に weight (重み) を設定することで +より一般的なものとなる. +\E +\BEG +Term orderings introduced in the previous section can be generalized +by setting a weight for each variable. +\E +@example +[0] dp_td(<<1,1,1>>); +3 +[1] dp_set_weight([1,2,3])$ +[2] dp_td(<<1,1,1>>); +6 +@end example +\BJP +単項式の全次数を計算する際, デフォルトでは +各変数の指数の和を全次数とする. これは各変数の weight を 1 と +考えていることに相当する. この例では, 第一, 第二, 第三変数の +weight をそれぞれ 1,2,3 と指定している. このため, @code{<<1,1,1>>} +の全次数 (以下ではこれを単項式の weight と呼ぶ) が @code{1*1+1*2+1*3=6} となる. +weight を設定することで, 同じ項順序型のもとで異なる項順序が定義できる. +例えば, weight をうまく設定することで, 多項式を weighted homogeneous +にすることができる場合がある. +\E +\BEG +By default, the total degree of a monomial is equal to +the sum of all exponents. This means that the weight for each variable +is set to 1. +In this example, the weights for the first, the second and the third +variable are set to 1, 2 and 3 respectively. +Therefore the total degree of @code{<<1,1,1>>} under this weight, +which is called the weight of the monomial, is @code{1*1+1*2+1*3=6}. +By setting weights, different term orderings can be set under a type of +term ordeing. In some case a polynomial can +be made weighted homogeneous by setting an appropriate weight. +\E + +\BJP +各変数に対する weight をまとめたものを weight vector と呼ぶ. +すべての成分が正であり, グレブナ基底計算において, 全次数の +代わりに用いられるものを特に sugar weight と呼ぶことにする. +sugar strategy において, 全次数の代わりに使われるからである. +一方で, 各成分が必ずしも正とは限らない weight vector は, +sugar weight として設定することはできないが, 項順序の一般化には +有用である. これらは, 行列による項順序の設定にすでに現れて +いる. すなわち, 項順序を定義する行列の各行が, 一つの weight vector +と見なされる. また, ブロック順序は, 各ブロックの +変数に対応する成分のみ 1 で他は 0 の weight vector による比較を +最初に行ってから, 各ブロック毎の tie breaking を行うことに相当する. +\E + +\BEG +A list of weights for all variables is called a weight vector. +A weight vector is called a sugar weight vector if +its elements are all positive and it is used for computing +a weighted total degree of a monomial, because such a weight +is used instead of total degree in sugar strategy. +On the other hand, a weight vector whose elements are not necessarily +positive cannot be set as a sugar weight, but it is useful for +generalizing term order. In fact, such a weight vector already +appeared in a matrix order. That is, each row of a matrix defining +a term order is regarded as a weight vector. A block order +is also considered as a refinement of comparison by weight vectors. +It compares two terms by using a weight vector whose elements +corresponding to variables in a block is 1 and 0 otherwise, +then it applies a tie breaker. +\E + +\BJP +weight vector の設定は @code{dp_set_weight()} で行うことができる +が, 項順序を指定する際の他のパラメタ (項順序型, 変数順序) と +まとめて設定できることが望ましい. このため, 次のような形でも +項順序が指定できる. +\E +\BEG +A weight vector can be set by using @code{dp_set_weight()}. +However it is more preferable if a weight vector can be set +together with other parapmeters such as a type of term ordering +and a variable order. This is realized as follows. +\E + +@example +[64] B=[x+y+z-6,x*y+y*z+z*x-11,x*y*z-6]$ +[65] dp_gr_main(B|v=[x,y,z],sugarweight=[3,2,1],order=0); +[z^3-6*z^2+11*z-6,x+y+z-6,-y^2+(-z+6)*y-z^2+6*z-11] +[66] dp_gr_main(B|v=[y,z,x],order=[[1,1,0],[0,1,0],[0,0,1]]); +[x^3-6*x^2+11*x-6,x+y+z-6,-x^2+(-y+6)*x-y^2+6*y-11] +[67] dp_gr_main(B|v=[y,z,x],order=[[x,1,y,2,z,3]]); +[x+y+z-6,x^3-6*x^2+11*x-6,-x^2+(-y+6)*x-y^2+6*y-11] +@end example + +\BJP +いずれの例においても, 項順序は option として指定されている. +最初の例では @code{v} により変数順序を, @code{sugarweight} により +sugar weight vector を, @code{order}により項順序型を指定している. +二つ目の例における @code{order} の指定は matrix order と同様である. +すなわち, 指定された weight vector を左から順に使って weight の比較 +を行う. 三つ目の例も同様であるが, ここでは weight vector の要素を +変数毎に指定している. 指定がないものは 0 となる. 三つ目の例では, +@code{order} による指定では項順序が決定しない. この場合には, +tie breaker として全次数逆辞書式順序が自動的に設定される. +この指定方法は, @code{dp_gr_main}, @code{dp_gr_mod_main} など +の組み込み関数でのみ可能であり, @code{gr} などのユーザ定義関数 +では未対応である. +\E +\BEG +In each example, a term ordering is specified as options. +In the first example, a variable order, a sugar weight vector +and a type of term ordering are specified by options @code{v}, +@code{sugarweight} and @code{order} respectively. +In the second example, an option @code{order} is used +to set a matrix ordering. That is, the specified weight vectors +are used from left to right for comparing terms. +The third example shows a variant of specifying a weight vector, +where each component of a weight vector is specified variable by variable, +and unspecified components are set to zero. In this example, +a term order is not determined only by the specified weight vector. +In such a case a tie breaker by the graded reverse lexicographic ordering +is set automatically. +This type of a term ordering specification can be applied only to builtin +functions such as @code{dp_gr_main()}, @code{dp_gr_mod_main()}, not to +user defined functions such as @code{gr()}. \E \BJP