version 1.29, 2004/08/22 12:52:34 |
version 1.30, 2004/08/23 08:33:55 |
|
|
% $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.28 2004/05/13 05:33:10 takayama Exp $ |
% $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.29 2004/08/22 12:52:34 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. |
|
|
Line 3954 $ [ff ff] fromVectors :: $ |
|
Line 3954 $ [ff ff] fromVectors :: $ |
|
[ ([a b c ...] denominator r) |
[ ([a b c ...] denominator r) |
( a denominator r ) |
( a denominator r ) |
(cf. dc, numerator) |
(cf. dc, numerator) |
|
(Output is Z or a polynomial.) |
] |
] |
] putUsages |
] putUsages |
% test data. |
% test data. |
Line 3964 $ [ff ff] fromVectors :: $ |
|
Line 3965 $ [ff ff] fromVectors :: $ |
|
[/pp /dd /ii /rr] pushVariables |
[/pp /dd /ii /rr] pushVariables |
[ |
[ |
/pp arg1 def |
/pp arg1 def |
|
pp to_univNum /pp set |
{ |
{ |
pp isArray { |
pp isArray { |
pp { denominator } map /dd set |
pp { denominator } map /dd set |
Line 3986 $ [ff ff] fromVectors :: $ |
|
Line 3988 $ [ff ff] fromVectors :: $ |
|
[ ([a b c ...] numerator r) |
[ ([a b c ...] numerator r) |
( a numerator r ) |
( a numerator r ) |
(cf. dc, denominator) |
(cf. dc, denominator) |
|
(Output is a list of Z or polynomials.) |
] |
] |
] putUsages |
] putUsages |
% test data. |
% test data. |
Line 3994 $ [ff ff] fromVectors :: $ |
|
Line 3997 $ [ff ff] fromVectors :: $ |
|
[/pp /dd /ii /rr] pushVariables |
[/pp /dd /ii /rr] pushVariables |
[ |
[ |
/pp arg1 def |
/pp arg1 def |
|
pp to_univNum /pp set |
{ |
{ |
pp isArray { |
pp isArray { |
pp denominator /dd set |
pp denominator /dd set |
pp dd mul /rr set |
pp dd mul /rr set |
rr reduce /rr set |
rr cancel /rr set |
exit |
exit |
} { } ifelse |
} { } ifelse |
|
|
Line 4012 $ [ff ff] fromVectors :: $ |
|
Line 4016 $ [ff ff] fromVectors :: $ |
|
arg1 |
arg1 |
} def |
} def |
|
|
/reduce.Q { |
/cancel.Q { |
/arg1 set |
/arg1 set |
[/aa /rr /nn /dd /gg] pushVariables |
[/aa /rr /nn /dd /gg] pushVariables |
[ |
[ |
Line 4039 $ [ff ff] fromVectors :: $ |
|
Line 4043 $ [ff ff] fromVectors :: $ |
|
arg1 |
arg1 |
} def |
} def |
|
|
/reduce.one { |
/cancel.one { |
/arg1 set |
/arg1 set |
[/aa /rr /nn /dd /gg] pushVariables |
[/aa /rr /nn /dd /gg] pushVariables |
[ |
[ |
Line 4049 $ [ff ff] fromVectors :: $ |
|
Line 4053 $ [ff ff] fromVectors :: $ |
|
aa (numerator) dc /nn set |
aa (numerator) dc /nn set |
aa (denominator) dc /dd set |
aa (denominator) dc /dd set |
nn isUniversalNumber dd isUniversalNumber and { |
nn isUniversalNumber dd isUniversalNumber and { |
/rr aa reduce.Q def |
/rr aa cancel.Q def |
exit |
exit |
} { (reduce: not implemented) error } ifelse |
} { (cancel: not implemented) error } ifelse |
} { } ifelse |
} { } ifelse |
|
|
/rr aa def |
/rr aa def |
Line 4063 $ [ff ff] fromVectors :: $ |
|
Line 4067 $ [ff ff] fromVectors :: $ |
|
arg1 |
arg1 |
} def |
} def |
|
|
[(reduce) |
[(cancel) |
[ (obj reduce r) |
[ (obj cancel r) |
(Cancel numerators and denominators) |
(Cancel numerators and denominators) |
(The implementation has not yet been completed. It works only for Q.) |
(The implementation has not yet been completed. It works only for Q.) |
]] putUsages |
]] putUsages |
/reduce { |
/cancel { |
/arg1 set |
/arg1 set |
[/aa /rr] pushVariables |
[/aa /rr] pushVariables |
[ |
[ |
/aa arg1 def |
/aa arg1 def |
aa isArray { |
aa isArray { |
aa {reduce} map /rr set |
aa {cancel} map /rr set |
} { |
} { |
aa reduce.one /rr set |
aa cancel.one /rr set |
} ifelse |
} ifelse |
/arg1 rr def |
/arg1 rr def |
] pop |
] pop |
popVariables |
popVariables |
arg1 |
arg1 |
} def |
} 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 |
|
|
/usages { |
/usages { |
/arg1 set |
/arg1 set |