version 1.34, 2004/09/01 05:06:09 |
version 1.37, 2004/09/10 13:20:23 |
|
|
% $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.33 2004/08/31 05:30:20 takayama Exp $ |
% $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.36 2004/09/09 11:42:22 takayama Exp $ |
%% dr.sm1 (Define Ring) 1994/9/25, 26 |
%% dr.sm1 (Define Ring) 1994/9/25, 26 |
%% This file is error clean. |
%% This file is error clean. |
|
|
|
|
( Pointer to the ring. ) |
( Pointer to the ring. ) |
(Example: [$x,y$ ring_of_q_difference_operators $Qx,Qy$ elimination_order) |
(Example: [$x,y$ ring_of_q_difference_operators $Qx,Qy$ elimination_order) |
( 0] define_qring ) |
( 0] define_qring ) |
(cf. define_ring, set_up_ring@ <coefficient ring>, ring_def, << ,, >>) |
(cf. define_ring, set_up_ring@ <coefficient ring>, ring_def, << __ >>) |
] |
] |
] putUsages |
] putUsages |
/define_qring { |
/define_qring { |
|
|
(one may use the command ) |
(one may use the command ) |
( f (ring) data_conversion /R set) |
( f (ring) data_conversion /R set) |
(cf. define_ring, define_qring, system_variable, poly (ring) data_conversion) |
(cf. define_ring, define_qring, system_variable, poly (ring) data_conversion) |
(cf. << ,, >>) |
(cf. << __ >>) |
] |
] |
] putUsages |
] putUsages |
|
|
|
|
gg (0). eq |
gg (0). eq |
{ 0 } |
{ 0 } |
{ gg (ring) data_conversion /rr set |
{ gg (ring) data_conversion /rr set |
gg << var rr ,, >> degree |
gg << var rr __ >> degree |
} ifelse |
} ifelse |
} map def |
} map def |
%%degs message |
%%degs message |
|
|
[(ClassP) @.datatypeConstant.usage ] putUsages |
[(ClassP) @.datatypeConstant.usage ] putUsages |
[(DoubleP) @.datatypeConstant.usage ] putUsages |
[(DoubleP) @.datatypeConstant.usage ] putUsages |
|
|
[(,,) |
[(__) |
[( string ring ,, polynomial) |
[( string ring __ polynomial) |
(Parse the <<string>> as an element in the <<ring>> and returns) |
(Parse the <<string>> as an element in the <<ring>> and returns) |
(the polynomial.) |
(the polynomial.) |
(cf. define_ring, define_qring, ring_def) |
(cf. define_ring, define_qring, ring_def) |
(Example: [(x,y) ring_of_polynomials [[(x) 1]] weight_vector 7]define_ring) |
(Example: [(x,y) ring_of_polynomials [[(x) 1]] weight_vector 7]define_ring) |
( /myring set) |
( /myring set) |
( ((x+y)^4) myring ,, /f set) |
( ((x+y)^4) myring __ /f set) |
]] putUsages |
]] putUsages |
|
|
/,, { |
/__ { |
/arg2 set /arg1 set |
/arg2 set /arg1 set |
[/rrr] pushVariables |
[/rrr] pushVariables |
[ arg1 tag StringP eq |
[ arg1 tag StringP eq |
|
|
/arg1 arg1 expand def |
/arg1 arg1 expand def |
[(CurrentRingp) rrr] system_variable |
[(CurrentRingp) rrr] system_variable |
} |
} |
{(Argument Error for ,, ) error } |
{(Argument Error for __ ) error } |
ifelse |
ifelse |
] pop |
] pop |
popVariables |
popVariables |
|
|
} def |
} def |
%%end of function |
%%end of function |
|
|
/rest { % returns remainder of a given list |
/rest { |
[ 2 1 roll aload length -1 roll pop ] |
/arg1 set [(Krest) arg1] extension |
} def |
} def |
[(rest) |
[(rest) |
[(array rest the-rest-of-the-array) |
[(array rest the-rest-of-the-array) |
|
|
error |
error |
} ifelse |
} ifelse |
/myring base 0 get (ring) dc def |
/myring base 0 get (ring) dc def |
/zero (0) myring ,, def |
/zero (0) myring __ def |
base length 1 sub /n set |
base length 1 sub /n set |
/minbase [ 0 1 n { /i set base i get } for ] def |
/minbase [ 0 1 n { /i set base i get } for ] def |
0 1 n { |
0 1 n { |
|
|
[(variableNames) M2N 0 get] system_variable |
[(variableNames) M2N 0 get] system_variable |
f toString /f2 set |
f toString /f2 set |
[(variableNames) M2N 3 get] system_variable |
[(variableNames) M2N 3 get] system_variable |
f2 M2N 2 get ,, /f2 set |
f2 M2N 2 get __ /f2 set |
} ifelse |
} ifelse |
[(CurrentRingp) cp] system_variable |
[(CurrentRingp) cp] system_variable |
/arg1 f2 def |
/arg1 f2 def |
|
|
$/R2 set$ |
$/R2 set$ |
$[[(x) (Dx)] [((t-1) Dt) (z)]] /r0 set$ |
$[[(x) (Dx)] [((t-1) Dt) (z)]] /r0 set$ |
$r0 R1 R2 makeRingMap /maptable set$ |
$r0 R1 R2 makeRingMap /maptable set$ |
$(Dx-1) R1 ,, /ff set$ |
$(Dx-1) R1 __ /ff set$ |
$ ff maptable ringmap :: $ |
$ ff maptable ringmap :: $ |
] |
] |
] putUsages |
] putUsages |
|
|
/rrr set |
/rrr set |
base { . } map /base set |
base { . } map /base set |
} ifelse |
} ifelse |
vlist { dup isPolynomial { } { rrr ,, } ifelse } map /vlist set |
vlist { dup isPolynomial { } { rrr __ } ifelse } map /vlist set |
|
|
[(hilbert) base vlist] extension /ff set |
[(hilbert) base vlist] extension /ff set |
[(CurrentRingp) rrrorg] system_variable |
[(CurrentRingp) rrrorg] system_variable |
|
|
/vv arg2 def |
/vv arg2 def |
/f hhh 1 get def |
/f hhh 1 get def |
f (0). eq { /ans [0] def /hilbReduce.label goto } { } ifelse |
f (0). eq { /ans [0] def /hilbReduce.label goto } { } ifelse |
f vv << f (ring) dc >> ,, degree /vv set |
f vv << f (ring) dc >> __ degree /vv set |
hhh 0 get /d set |
hhh 0 get /d set |
d d (integer) dc factorial /d set |
d d (integer) dc factorial /d set |
d << vv (universalNumber) dc vv factorial >> idiv /d set |
d << vv (universalNumber) dc vv factorial >> idiv /d set |
|
|
{ /fn (0). def } |
{ /fn (0). def } |
{ |
{ |
f (ring) dc /rrr set |
f (ring) dc /rrr set |
v toString (^) n toString 3 cat_n rrr ,, |
v toString (^) n toString 3 cat_n rrr __ |
f mul |
f mul |
[[v (0).] [(h) rrr ,, (1) rrr ,,]] replace /fn set |
[[v (0).] [(h) rrr __ (1) rrr __]] replace /fn set |
} ifelse |
} ifelse |
fn /arg1 set |
fn /arg1 set |
] pop |
] pop |
|
|
/ddd 0 def |
/ddd 0 def |
0 2 www length 1 sub { |
0 2 www length 1 sub { |
/iii set |
/iii set |
fff << www iii get rrr ,, >> degree |
fff << www iii get rrr __ >> degree |
<< www iii 1 add get >> mul |
<< www iii 1 add get >> mul |
ddd add /ddd set |
ddd add /ddd set |
} for |
} for |
|
|
fff (0). eq { /ddd -intInfinity def /ord_w_all.LLL goto} { } ifelse |
fff (0). eq { /ddd -intInfinity def /ord_w_all.LLL goto} { } ifelse |
/ddd -intInfinity def |
/ddd -intInfinity def |
fff (ring) dc /rrr set |
fff (ring) dc /rrr set |
/zzz (0) rrr ,, def |
/zzz (0) rrr __ def |
fff init /fff-in set |
fff init /fff-in set |
fff fff-in sub /fff set |
fff fff-in sub /fff set |
{ |
{ |
/ddd-tmp 0 def |
/ddd-tmp 0 def |
0 2 www length 1 sub { |
0 2 www length 1 sub { |
/iii set |
/iii set |
fff-in << www iii get rrr ,, >> degree |
fff-in << www iii get rrr __ >> degree |
<< www iii 1 add get >> mul |
<< www iii 1 add get >> mul |
ddd-tmp add /ddd-tmp set |
ddd-tmp add /ddd-tmp set |
} for |
} for |
|
|
/v0 vv ii get (string) dc def |
/v0 vv ii get (string) dc def |
v0 (array) dc 0 get Dascii eq %% If the first character is D? |
v0 (array) dc 0 get Dascii eq %% If the first character is D? |
{ rule %% Dx-->x |
{ rule %% Dx-->x |
[v0 rr ,, |
[v0 rr __ |
v0 (array) dc rest { (string) dc} map aload length cat_n rr ,,] |
v0 (array) dc rest { (string) dc} map aload length cat_n rr __] |
append /rule set |
append /rule set |
} |
} |
{ rule %% x --> -Dx |
{ rule %% x --> -Dx |
[v0 rr ,, |
[v0 rr __ |
(0). |
(0). |
[Dascii] v0 (array) dc join { (string) dc } map aload length |
[Dascii] v0 (array) dc join { (string) dc } map aload length |
cat_n rr ,, sub |
cat_n rr __ sub |
] |
] |
append /rule set |
append /rule set |
} ifelse |
} ifelse |
} ifelse |
} ifelse |
} for |
} for |
% rule message |
% rule message |
ff rule replace [[(h) rr ,, (1) rr ,,]] replace /ans1 set |
ff rule replace [[(h) rr __ (1) rr __]] replace /ans1 set |
} ifelse |
} ifelse |
} |
} |
{ |
{ |
|
|
{ |
{ |
/xx xx (string) dc def |
/xx xx (string) dc def |
/dxx [@@@.Dsymbol xx] cat def |
/dxx [@@@.Dsymbol xx] cat def |
/xx xx f (ring) dc ,, def |
/xx xx f (ring) dc __ def |
/dxx dxx f (ring) dc ,, def |
/dxx dxx f (ring) dc __ def |
/one (1) f (ring) dc ,, def |
/one (1) f (ring) dc __ def |
|
|
{ |
{ |
/g f init def |
/g f init def |
|
|
f (0). eq { } |
f (0). eq { } |
{ |
{ |
/rr f (ring) dc def |
/rr f (ring) dc def |
xx {toString rr ,, } map /xx set |
xx {toString rr __ } map /xx set |
dx {toString rr ,, } map /dx set |
dx {toString rr __ } map /dx set |
ss {toString rr ,, } map /ss set |
ss {toString rr __ } map /ss set |
/n xx length def |
/n xx length def |
0 1 n 1 sub { |
0 1 n 1 sub { |
/i set |
/i set |
|
|
[ |
[ |
/f arg1 def /xx arg2 def /dx arg3 def /ss arg4 def |
/f arg1 def /xx arg2 def /dx arg3 def /ss arg4 def |
f (ring) dc /rr set |
f (ring) dc /rr set |
/one (1) rr ,, def %% |
/one (1) rr __ def %% |
/ww [ xx toString -1 dx toString 1 ] weightv def |
/ww [ xx toString -1 dx toString 1 ] weightv def |
f ww init f sub (0). eq { } |
f ww init f sub (0). eq { } |
{ [(destraction2.1 : inhomogeneous with respect to ) |
{ [(destraction2.1 : inhomogeneous with respect to ) |
Line 3874 $ [ff ff] fromVectors :: $ |
|
Line 3874 $ [ff ff] fromVectors :: $ |
|
[(>) (stringOut://@@@stdout) (2>) (stringOut://@@@stderr)] join |
[(>) (stringOut://@@@stdout) (2>) (stringOut://@@@stderr)] join |
} def |
} def |
|
|
[(,,,) |
[(___) |
[(reparse a polynomial or polynomials)] |
[(reparse a polynomial or polynomials)] |
] putUsages |
] putUsages |
/,,, { |
/___ { |
/arg1 set |
/arg1 set |
[/in-reparse /ff] pushVariables |
[/in-reparse /ff] pushVariables |
[ |
[ |
/ff arg1 def |
/ff arg1 def |
ff tag 6 eq { |
ff tag 6 eq { |
ff { ,,, } map /arg1 set |
ff { ___ } map /arg1 set |
} { |
} { |
ff toString . /arg1 set |
ff toString . /arg1 set |
} ifelse |
} ifelse |
Line 4189 $ [ff ff] fromVectors :: $ |
|
Line 4189 $ [ff ff] fromVectors :: $ |
|
( ma (dog) getNode ) |
( ma (dog) getNode ) |
]] putUsages |
]] putUsages |
|
|
|
/cons { |
|
/arg2 set /arg1 set |
|
[/aa /bb] pushVariables |
|
[ |
|
/aa arg1 def /bb arg2 def |
|
[aa] (list) dc bb join /arg1 set |
|
] pop |
|
popVariables |
|
arg1 |
|
} def |
|
[(cons) |
|
[(obj list cons list) |
|
]] putUsages |
|
|
/usages { |
/usages { |
/arg1 set |
/arg1 set |