1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56
|
#
# jython examples for jas.
# $Id$
#
import sys;
from jas import Ring
from jas import Ideal
from jas import startLog
# example from rose (modified)
# optimal is
# Rat(A46, U3, U4)
#
r = Ring( "Rat(U3,U4,A46) G" );
print "Ring: " + str(r);
print;
ps = """
(
( U4^4 - 20/7 A46^2 ),
( A46^2 U3^4 + 7/10 A46 U3^4 + 7/48 U3^4 - 50/27 A46^2 - 35/27 A46 - 49/216 ),
( A46^5 U4^3 + 7/5 A46^4 U4^3 + 609/1000 A46^3
U4^3 + 49/1250 A46^2 U4^3 - 27391/800000 A46 U4^3
- 1029/160000 U4^3 + 3/7 A46^5 U3 U4^2 + 3/5 A46^6
U3 U4^2 + 63/200 A46^3 U3 U4^2 + 147/2000 A46^2
U3 U4^2 + 4137/800000 A46 U3 U4^2 - 7/20 A46^4
U3^2 U4 - 77/125 A46^3 U3^2 U4 - 23863/60000 A46^2
U3^2 U4 - 1078/9375 A46 U3^2 U4 - 24353/1920000
U3^2 U4 - 3/20 A46^4 U3^3 - 21/100 A46^3 U3^3
- 91/800 A46^2 U3^3 - 5887/200000 A46 U3^3
- 343/128000 U3^3 )
)
""";
f = Ideal( r, ps );
print "Ideal: " + str(f);
print;
startLog();
o = f.optimize();
print "optimized Ideal: " + str(o);
print;
#rg = f.GB();
#print "Output:", rg;
#print;
org = o.GB();
print "opt Output:", org;
print;
|