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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
|
#
# jython examples for jas.
# $Id$
#
import sys;
from jas import Ring, GF
#from edu.jas.arith import ModIntegerRing
# example from rose (modified) modular computation
#r = Ring( "Mod 19 (U3,U4,A46) L" );
#r = Ring( "Mod 1152921504606846883 (U3,U4,A46) L" ); # 2^60-93
#r = Ring( "Quat(U3,U4,A46) L" );
#r = Ring( "Z(U3,U4,A46) L" );
#r = Ring( "C(U3,U4,A46) L" );
r = Ring( "Rat(A46,U3,U4) 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 = r.ideal( ps );
print "Ideal: " + str(f);
print;
fi = f.toInteger();
print "Ideal: " + str(fi);
print;
# "1152921504606846883"
#mf = ModIntegerRing( str(2**60-93), True );
#mf = ModIntegerRing( str(19), True );
mf = GF( 2**60-93);
fm = fi.toModular( mf );
print "Ideal: " + str(fm);
#sys.exit();
rm = fm.GB();
#print "seq Output:", rm;
#print;
#rg = f.GB();
#print "seq Output:", rg;
#print;
#rgm = rg.toInteger().toModular(mf);
#print "seq Output:", rgm;
#print;
#e = rm == rgm;
#print "e = " + str(e);
sys.exit();
#for i in range(1,9):
for i in range(1,3):
rg = f.parGB(i);
#print "par Output:", rg;
#print;
rg = f.parOldGB(i);
#print "par-old Output:", rg;
#print;
#f.distClient(); # starts in background
#rg = f.distGB(i);
#print "dist Output:", rg;
#print;
#sys.exit();
|
