===================================================================
RCS file: /home/cvs/OpenXM/src/kan96xx/Kan/dr.sm1,v
retrieving revision 1.3
retrieving revision 1.31
diff -u -p -r1.3 -r1.31
--- OpenXM/src/kan96xx/Kan/dr.sm1 2000/02/02 09:07:59 1.3
+++ OpenXM/src/kan96xx/Kan/dr.sm1 2004/08/27 08:33:27 1.31
@@ -1,9 +1,9 @@
-% $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.2 1999/11/08 00:16:05 takayama Exp $
+% $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.30 2004/08/23 08:33:55 takayama Exp $
%% dr.sm1 (Define Ring) 1994/9/25, 26
%% This file is error clean.
@@@.quiet { }
-{ (macro package : dr.sm1, 9/26,1995 --- Version 2/2, 2000. ) message } ifelse
+{ (macro package : dr.sm1, 9/26,1995 --- Version 12/10, 2000. ) message } ifelse
/ctrlC-hook {
%%% define your own routing in case of error.
@@ -50,7 +50,7 @@
} for
] /dList set
- [(H)] xList join [@@@.esymbol] join /xList set
+ [@@@.Hsymbol] xList join [@@@.esymbol] join /xList set
[(h)] dList join [@@@.Esymbol] join /dList set
[0 %% dummy characteristic
<< xList length >> << xList length >> << xList length >>
@@ -167,7 +167,7 @@
vars reverse /xList set %[z y x]
vars {@@@.Dsymbol 2 1 roll 2 cat_n} map
reverse /dList set %[Dz Dy Dx]
- [(H)] xList join [@@@.esymbol] join /xList set
+ [@@@.Hsymbol] xList join [@@@.esymbol] join /xList set
[(h)] dList join [@@@.Esymbol] join /dList set
[0 1 1 1 << xList length >>
1 1 1 << xList length 1 sub >> ] /param set
@@ -257,6 +257,8 @@
/arg1 set
[/vars /n /i /xList /dList /param] pushVariables
[
+ (This is an obsolete macro. Use ring_of_differential_difference_operators)
+ error
(mmLarger) (matrix) switch_function
(mpMult) (difference) switch_function
(red@) (module1) switch_function
@@ -267,7 +269,7 @@
vars reverse /xList set %[z y x]
vars {@@@.diffEsymbol 2 1 roll 2 cat_n} map
reverse /dList set %[Dz Dy Dx]
- [(H)] xList join [@@@.esymbol] join /xList set
+ [@@@.Hsymbol] xList join [@@@.esymbol] join /xList set
[(h)] dList join [@@@.Esymbol] join /dList set
[0 1 1 << xList length >> << xList length >>
1 1 << xList length 1 sub >> << xList length >> ] /param set
@@ -278,7 +280,44 @@
} def
+/ring_of_differential_difference_operators {
+ /arg1 set
+ [/vars /n /i /xList /dList /param /dvar /evar /vars2 ] pushVariables
+ [
+ /vars arg1 def
+ vars tag 6 eq not {
+ ( List is expected as the argument for ring_of_differential_difference_operators ) error
+ } { } ifelse
+ vars 0 get /dvar set
+ vars 1 get /evar set
+ (mmLarger) (matrix) switch_function
+ (mpMult) (difference) switch_function
+ (red@) (module1) switch_function
+ (groebner) (standard) switch_function
+ (isSameComponent) (x) switch_function
+ [dvar to_records pop] /vars set %[x y z]
+ vars reverse /xList set %[z y x]
+
+ [evar to_records pop] /vars2 set %[s1 s2]
+
+ vars2 reverse {@@@.Esymbol 2 1 roll 2 cat_n} map
+ xList
+ join /xList set %[Es2 Es1 z y x]
+
+ vars2 reverse
+ vars {@@@.Dsymbol 2 1 roll 2 cat_n} map
+ reverse join /dList set %[s2 s1 Dz Dy Dx]
+ [@@@.Hsymbol] xList join [@@@.esymbol] join /xList set
+ [(h)] dList join [@@@.Esymbol] join /dList set
+ [0 1 1 << vars2 length 1 add >> << xList length >>
+ 1 1 << vars2 length 1 add >> << xList length 1 sub >> ] /param set
+ [ xList dList param ] /arg1 set
+ ] pop
+ popVariables
+ arg1
+} def
+
/reverse {
/arg1 set
arg1 length 1 lt
@@ -1474,6 +1513,13 @@
(type?) data_conversion RationalFunctionP eq
} def
+[(isRing)
+ [(obj isRing bool) ]
+] putUsages
+/isRing {
+ (type?) data_conversion RingP eq
+} def
+
/toString.tmp {
/arg1 set
[/obj /fname] pushVariables
@@ -1720,7 +1766,11 @@
/arg1 set
[/f /k /i /ans] pushVariables
[
- /ans (1).. def
+ /ans (1).. def
+ [(QuoteMode)] system_variable {
+ /f arg1 def /k arg2 def
+ [(ooPower) f k] extension /ans set
+ } {
/f arg1 def /k arg2 ..int def
k 0 lt {
1 1 << 0 k sub >> {
@@ -1733,7 +1783,8 @@
/ans f ans {mul} sendmsg2 def
} for
} ifelse
- /arg1 ans def
+ } ifelse
+ /arg1 ans def
] pop
popVariables
arg1
@@ -2276,7 +2327,7 @@ newline
/tolower.aux {
/arg1 set
- arg1 64 gt arg1 96 lt and
+ arg1 64 gt arg1 91 lt and
{ arg1 32 add }
{ arg1 } ifelse
} def
@@ -2460,7 +2511,7 @@ newline
[(ord_w)
[(ff [v1 w1 v2 w2 ... vm wm] ord_w d)
(poly ff; string v1; integer w1; ...)
- (order of ff by the weight vector [w1 w2 ...])
+ (order of the initial of ff by the weight vector [w1 w2 ...])
(Example: [(x,y) ring_of_polynomials 0] define_ring )
( (x^2 y^3-x). [(x) 2 (y) 1] ord_w ::)
]
@@ -2488,6 +2539,46 @@ newline
arg1
} def
+[(ord_w_all)
+ [(ff [v1 w1 v2 w2 ... vm wm] ord_w d)
+ (poly ff; string v1; integer w1; ...)
+ (order of ff by the weight vector [w1 w2 ...])
+ (Example: [(x,y,t) ring_of_polynomials 0] define_ring )
+ ( (x^2 y^3-x-t). [(t) 1 ] ord_w_all ::)
+ ]
+] putUsages
+/ord_w_all {
+ /arg2 set /arg1 set
+ [/ord_w_all-in /fff /fff-in /www /rrr /iii /ddd /zzz /ddd-tmp] pushVariables
+ [
+ /fff arg1 def
+ /www arg2 def
+ fff (0). eq { /ddd -intInfinity def /ord_w_all.LLL goto} { } ifelse
+ /ddd -intInfinity def
+ fff (ring) dc /rrr set
+ /zzz (0) rrr ,, def
+ fff init /fff-in set
+ fff fff-in sub /fff set
+ {
+ /ddd-tmp 0 def
+ 0 2 www length 1 sub {
+ /iii set
+ fff-in << www iii get rrr ,, >> degree
+ << www iii 1 add get >> mul
+ ddd-tmp add /ddd-tmp set
+ } for
+ ddd-tmp ddd gt { /ddd ddd-tmp def } { } ifelse
+ fff zzz eq { exit } { } ifelse
+ fff init /fff-in set
+ fff fff-in sub /fff set
+ } loop
+ /ord_w_all.LLL
+ /arg1 ddd def
+ ] pop
+ popVariables
+ arg1
+} def
+
[(laplace0)
[
(f [v1 ... vn] laplace0 g)
@@ -2676,6 +2767,7 @@ newline
[/in-ngcd /nlist /g.ngcd /ans] pushVariables
[
/nlist arg1 def
+ nlist to_univNum /nlist set
nlist length 2 lt
{ /ans nlist 0 get def /L.ngcd goto }
{
@@ -3236,6 +3328,10 @@ newline
]] putUsages
/cancelCoeff {
+ /arg1 set
+ [(reduceContent) arg1] gbext 0 get
+} def
+/cancelCoeff_org {
/arg1 set
[/in-cancelCoeff /ff /gg /dd /dd2] pushVariables
[ /ff arg1 def
@@ -3532,6 +3628,577 @@ $ [ff ff] fromVectors :: $
[(nl)
[(nl is the newline character.)
$Example: [(You can break line) nl (here.)] cat message$
+]] putUsages
+
+/to_int {
+ /arg1 set
+ [/to-int /ob /ans] pushVariables
+ [
+ /ob arg1 def
+ /ans ob def
+ ob isArray {
+ ob {to_int} map /ans set
+ /LLL.to_int goto
+ } { } ifelse
+ ob isInteger {
+ ob (universalNumber) dc /ans set
+ /LLL.to_int goto
+ } { } ifelse
+ /LLL.to_int
+ /arg1 ans def
+ ] pop
+ popVariables
+ arg1
+} def
+[(to_int)
+[(obj to_int obj2)
+ (All integers in obj are changed to universalNumber.)
+ (Example: /ff [1 2 [(hello) (0).]] def ff { tag } map ::)
+ ( ff to_int { tag } map :: )
+]] putUsages
+
+/define_ring_variables {
+ [/in-define_ring_variables /drv._v /drv._p /drv._v0] pushVariables
+%% You cannot use these names for names for polynomials.
+ [
+ /drv._v getVariableNames def
+ /drv._v0 drv._v def
+ drv._v { dup /drv._p set (/) 2 1 roll ( $) drv._p ($. def ) } map cat
+ /drv._v set
+% drv._v message
+ [(parse) drv._v] extension
+ ] pop
+ popVariables
+} def
+[(define_ring_variables)
+[(It binds a variable <> in the current ring to the sm1 variable <>.)
+ (For example, if x is a variable in the current ring, it defines the sm1)
+ (variable x by /x (x) def)
+]] putUsages
+
+/boundp {
+ /arg1 set
+ [/a /ans] pushVariables
+ [
+ /a arg1 def
+ [(parse) [(/) a ( load tag 0 eq { /ans 0 def } )
+ ( { /ans 1 def } ifelse )] cat ] extension
+ /arg1 ans def
+ ] pop
+ popVariables
+ arg1
+} def
+[(boundp)
+ [( a boundp b)
+ (string a, b is 0 or 1.)
+ (If the variable named << a >> is bounded to a value,)
+ (it returns 1 else it returns 0.)
+ $Example: (hoge) boundp ::$
+]] putUsages
+[(isSubstr)
+ [
+ (s1 s2 isSubstr pos)
+ (If s1 is a substring of s2, isSubstr returns the position in s2 from which)
+ (s1 is contained in s2.)
+ (If s1 is not a substring of s2, then isSubstr returns -1.)
+ ]
+] putUsages
+/isSubstr {
+ /arg2 set /arg1 set
+ [/in-isSubstr /s1 /s2 /i1 /i2 /n1 /n2
+ /ans /flg
+ ] pushVariables
+ [
+ /s1 arg1 def
+ /s2 arg2 def
+ s1 (array) dc /s1 set
+ s2 (array) dc /s2 set
+ /n1 s1 length def
+ /n2 s2 length def
+ /ans -1 def
+ 0 1 n2 n1 sub {
+ /i2 set
+ /flg 1 def
+ 0 1 n1 1 sub {
+ /i1 set
+ s1 i1 get s2 i2 i1 add get eq {
+ } {
+ /flg 0 def exit
+ } ifelse
+ } for
+ flg {
+ /ans i2 def
+ /isSubstr.L2 goto
+ } { /ans -1 def } ifelse
+ } for
+ /isSubstr.L2
+ /arg1 ans def
+ ] pop
+ popVariables
+ arg1
+} def
+
+[(execve)
+ [
+ (command execve)
+ ([arg0 arg1 arg2 ...] execve )
+ (It executes the command by the system call execve.)
+ (cf. system, forkExec)
+ ]
+] putUsages
+
+/execve {
+ /execve.arg set
+ [(forkExec) execve.arg [ ] 1] extension
+} def
+
+[(beginEcart)
+ [
+ (beginEcart)
+ (Set the environments for the ecart division algorithm.)
+ ]
+] putUsages
+
+/ecart.debug_reduction1 0 def
+/beginEcart {
+ (red@) (ecart) switch_function
+ [(Ecart) 1] system_variable
+ [(CheckHomogenization) 0] system_variable
+ [(ReduceLowerTerms) 0] system_variable
+ [(AutoReduce) 0] system_variable
+ [(EcartAutomaticHomogenization) 0] system_variable
+ ecart.debug_reduction1 {
+ (red@) (debug) switch_function
+ } { } ifelse
+} def
+
+[(endEcart)
+ [
+ (endEcart)
+ (End of using the ecart division algorithm.)
+ ]
+] putUsages
+
+/endEcart {
+ (red@) (standard) switch_function
+ [(Ecart) 0] system_variable
+ [(degreeShift) (reset)] homogenize pop
+} def
+
+/ord_ws_all {
+ /arg2 set /arg1 set
+ [(ord_ws_all) arg1 arg2] gbext
+} def
+[(ord_ws_all)
+ [
+ (fv wv ord_ws_all degree)
+ ( ord_ws_all returns the ord with respect to the weight vector wv.)
+ $Example: [(x,y) ring_of_differential_operators 0] define_ring $
+ $ (Dx^2+x*Dx*Dy+2). [(Dx) 1 (Dy) 1] weightv ord_ws_all :: $
+ ( )
+ (fv [wv shiftv] ord_ws_all degree)
+ ( ord_ws_all returns the ord with respect to the weight vector wv and)
+ ( the shift vector shiftv.)
+ $Example: [(x,y) ring_of_differential_operators 0] define_ring $
+ $ [(Dx^2+x*Dx*Dy+2). (Dx).] [[(Dx) 1 (Dy) 1] weightv [0 2]] ord_ws_all ::$
+ ( )
+ (cf: init, gbext. Obsolete: ord_w, ord_w_all)
+ ]
+] putUsages
+
+[(newVector)
+ [( n newVector vec)
+]] putUsages
+/newVector {
+ /arg1 set
+ [/in-newVector /n] pushVariables
+ [
+ /n arg1 def
+ [(newVector) n] extension /arg1 set
+ ] pop
+ popVariables
+ arg1
+} def
+
+[(newMatrix)
+ [( [m n] newMatrix mat)
+]] putUsages
+/newMatrix {
+ /arg1 set
+ [/in-newMatrix /n] pushVariables
+ [
+ /n arg1 def
+ [(newMatrix) n 0 get n 1 get] extension /arg1 set
+ ] pop
+ popVariables
+ arg1
+} def
+
+/addStdoutStderr {
+ [(>) (stringOut://@@@stdout) (2>) (stringOut://@@@stderr)] join
+} def
+
+[(,,,)
+[(reparse a polynomial or polynomials)]
+] putUsages
+/,,, {
+ /arg1 set
+ [/in-reparse /ff] pushVariables
+ [
+ /ff arg1 def
+ ff tag 6 eq {
+ ff { ,,, } map /arg1 set
+ } {
+ ff toString . /arg1 set
+ } ifelse
+ ] pop
+ popVariables
+ arg1
+} def
+
+/to_univNum {
+ /arg1 set
+ [/rr ] pushVariables
+ [
+ /rr arg1 def
+ rr isArray {
+ rr { to_univNum } map /rr set
+ } {
+ } ifelse
+ rr isInteger {
+ rr (universalNumber) dc /rr set
+ } {
+ } ifelse
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+[(to_univNum)
+[(obj to_univNum obj2)
+ (Example. [ 2 (3).. ] to_univNum)
+ (cf. to_int)
+]] putUsages
+
+[(lcm)
+ [ ([a b c ...] lcm r)
+ (cf. polylcm, mpzext)
+ ]
+] putUsages
+/lcm {
+ /arg1 set
+ [/aa /bb /rr /pp /i] pushVariables
+ [
+ /aa arg1 def
+ /rr (1).. def
+ /pp 0 def % isPolynomial array?
+ 0 1 aa length 1 sub {
+ /i set
+ aa i get isPolynomial {
+ /pp 1 def
+ exit
+ } { } ifelse
+ } for
+
+ 0 1 aa length 1 sub {
+ /i set
+ pp {
+ [rr aa i get] polylcm /rr set
+ } {
+ [(lcm) rr aa i get ] mpzext /rr set
+ } ifelse
+ } for
+
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+[(gcd)
+ [ ([a b c ...] gcd r)
+ (cf. polygcd, mpzext)
+ ]
+] putUsages
+/gcd {
+ /arg1 set
+ [/aa /bb /rr /pp /i] pushVariables
+ [
+ /aa arg1 def
+ /rr (1).. def
+ /pp 0 def % isPolynomial array?
+ 0 1 aa length 1 sub {
+ /i set
+ aa i get isPolynomial {
+ /pp 1 def
+ /rr aa i get def
+ exit
+ } { } ifelse
+ } for
+
+ pp {
+ 0 1 aa length 1 sub {
+ /i set
+ [rr aa i get] polygcd /rr set
+ } for
+ } {
+ aa ngcd /rr set
+ } ifelse
+
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+
+[(denominator)
+ [ ([a b c ...] denominator r)
+ ( a denominator r )
+ (cf. dc, numerator)
+ (Output is Z or a polynomial.)
+ ]
+] putUsages
+% test data.
+% [(1).. (2).. div (1).. (3).. div ] denominator
+% [(2).. (3).. (4).. ] denominator
+/denominator {
+ /arg1 set
+ [/pp /dd /ii /rr] pushVariables
+ [
+ /pp arg1 def
+ pp to_univNum /pp set
+ {
+ pp isArray {
+ pp { denominator } map /dd set
+ /rr dd lcm def % rr = lcm(dd[0], dd[1], ... )
+ rr /dd set
+ exit
+ } { } ifelse
+
+ pp (denominator) dc /dd set
+ exit
+
+ } loop
+ /arg1 dd def
+ ] pop
+ popVariables
+ arg1
+} def
+
+[(numerator)
+ [ ([a b c ...] numerator r)
+ ( a numerator r )
+ (cf. dc, denominator)
+ (Output is a list of Z or polynomials.)
+ ]
+] putUsages
+% test data.
+/numerator {
+ /arg1 set
+ [/pp /dd /ii /rr] pushVariables
+ [
+ /pp arg1 def
+ pp to_univNum /pp set
+ {
+ pp isArray {
+ pp denominator /dd set
+ pp dd mul /rr set
+ rr cancel /rr set
+ exit
+ } { } ifelse
+
+ pp (numerator) dc /rr set
+ exit
+
+ } loop
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+
+/cancel.Q {
+ /arg1 set
+ [/aa /rr /nn /dd /gg] pushVariables
+ [
+ /aa arg1 def
+ {
+ aa isRational {
+ [(cancel) aa] mpzext /rr set
+ rr (denominator) dc (1).. eq {
+ /rr rr (numerator) dc def
+ exit
+ } { } ifelse
+ rr (denominator) dc (-1).. eq {
+ /rr rr (numerator) dc (-1).. mul def
+ } { } ifelse
+ exit
+ } { } ifelse
+
+ /rr aa def
+ exit
+ } loop
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+
+/cancel.one {
+ /arg1 set
+ [/aa /rr /nn /dd /gg] pushVariables
+ [
+ /aa arg1 def
+ {
+ aa isRational {
+ aa (numerator) dc /nn set
+ aa (denominator) dc /dd set
+ nn isUniversalNumber dd isUniversalNumber and {
+ /rr aa cancel.Q def
+ exit
+ } { (cancel: not implemented) error } ifelse
+ } { } ifelse
+
+ /rr aa def
+ exit
+ } loop
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+
+[(cancel)
+ [ (obj cancel r)
+ (Cancel numerators and denominators)
+ (The implementation has not yet been completed. It works only for Q.)
+]] putUsages
+/cancel {
+ /arg1 set
+ [/aa /rr] pushVariables
+ [
+ /aa arg1 def
+ aa isArray {
+ aa {cancel} map /rr set
+ } {
+ aa cancel.one /rr set
+ } ifelse
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+
+/nnormalize_vec {
+ /arg1 set
+ [/pp /rr /dd ] pushVariables
+ [
+ /pp arg1 def
+ pp denominator /dd set
+ dd (0).. lt { (nnormalize_vec: internal error) error } { } ifelse
+ pp numerator dd mul cancel /pp set
+ /@@@.nnormalize_vec_c dd def
+ pp gcd /dd set
+ dd (0).. lt { (nnormalize_vec: internal error) error } { } ifelse
+ pp (1).. dd div mul cancel /rr set
+ @@@.nnormalize_vec_c dd div cancel /@@@.nnormalize_vec_c set
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+[(nnormalize_vec)
+[(pp nnormalize_vec npp)
+ (It normalizes a given vector of Q into a vector of Z with relatively prime)
+ (entries by multiplying a postive number.)
+]] putUsages
+
+/getNode {
+ /arg2 set
+ /arg1 set
+ [/in-getNode /ob /key /rr /rr /ii] pushVariables
+ [
+ /ob arg1 def
+ /key arg2 def
+ /rr null def
+ {
+ ob isClass {
+ ob (array) dc /ob set
+ } { exit } ifelse
+ ob 0 get key eq {
+ /rr ob def
+ exit
+ } { } ifelse
+ ob 2 get /ob set
+ 0 1 ob length 1 sub {
+ /ii set
+ ob ii get key getNode /rr set
+ rr tag 0 eq { } { exit } ifelse
+ } for
+ exit
+ } loop
+ /arg1 rr def
+ ] pop
+ popVariables
+ arg1
+} def
+[(getNode)
+[(ob key getNode)
+ (ob is a class object.)
+ (The operator getNode returns the node with the key in ob.)
+ (The node is an array of the format [key attr-list node-list])
+ (Example:)
+ ( /dog [(dog) [[(legs) 4] ] [ ]] [(class) (tree)] dc def)
+ ( /man [(man) [[(legs) 2] ] [ ]] [(class) (tree)] dc def)
+ ( /ma [(mammal) [ ] [man dog]] [(class) (tree)] dc def)
+ ( ma (dog) getNode )
+]] putUsages
+
+
+/usages {
+ /arg1 set
+ [/name /flag /n /k /slist /m /i /sss] pushVariables
+ [
+ /name arg1 def
+ /flag true def
+
+ %BUG: should use regular expression in a future.
+ name tag 1 eq {
+ @.usages { 0 get } map shell { (, ) nl } map /sss set
+ } {
+
+ /sss [ ] def
+ @.usages length /n set
+ 0 1 << n 1 sub >>
+ {
+ /k set
+ name << @.usages k get 0 get >> eq
+ {
+ /slist @.usages k get 1 get def
+ /m slist length def
+ 0 1 << m 1 sub >> {
+ /i set
+ sss slist i get append nl append /sss set
+ } for
+ /flag false def
+ }
+ { }
+ ifelse
+ } for
+
+ %BUG: cannot get usages of primitives.
+ flag
+ {name Usage /sss [(Usage of ) name ( could not obtained.) nl ] def}
+ { }
+ ifelse
+ } ifelse
+ /arg1 sss cat def
+ ] pop
+ popVariables
+ arg1
+} def
+[(usages)
+ [(key usages usages-as-a-string)
+ (num usages list-of-key-words)
]] putUsages
;