8
8
The basic command syntax is as follows::
9
9
10
10
sage: import sage.matrix.benchmark as b
11
- sage: print "starting"; import sys; sys.stdout.flush(); b.report([b.det_ZZ], 'Test', systems=['sage'])
11
+ sage: print( "starting") ; import sys; sys.stdout.flush(); b.report([b.det_ZZ], 'Test', systems=['sage'])
12
12
starting...
13
13
======================================================================
14
14
Test
15
15
======================================================================
16
16
...
17
17
======================================================================
18
18
"""
19
+ from __future__ import print_function
20
+
19
21
from constructor import random_matrix , Matrix
20
22
from sage .rings .all import ZZ , QQ , GF
21
23
from sage .misc .misc import cputime
@@ -41,7 +43,7 @@ def report(F, title, systems = ['sage', 'magma'], **kwds):
41
43
EXAMPLES::
42
44
43
45
sage: import sage.matrix.benchmark as b
44
- sage: print "starting"; import sys; sys.stdout.flush(); b.report([b.det_ZZ], 'Test', systems=['sage'])
46
+ sage: print( "starting") ; import sys; sys.stdout.flush(); b.report([b.det_ZZ], 'Test', systems=['sage'])
45
47
starting...
46
48
======================================================================
47
49
Test
@@ -52,15 +54,15 @@ def report(F, title, systems = ['sage', 'magma'], **kwds):
52
54
import os
53
55
if len (systems ) > 2 :
54
56
raise NotImplementedError ("at most two systems ('sage' or 'magma')" )
55
- print '=' * 70
56
- print ' ' * 10 + title
57
- print '=' * 70
57
+ print ( '=' * 70 )
58
+ print ( ' ' * 10 + title )
59
+ print ( '=' * 70 )
58
60
os .system ('uname -a' )
59
- print '\n '
61
+ print ( '\n ' )
60
62
for f in F :
61
- print "-" * 70
62
- print f .__doc__ .strip ()
63
- print ( '%15s' * len (systems ))% tuple (systems )
63
+ print ( "-" * 70 )
64
+ print ( f .__doc__ .strip () )
65
+ print (( '%15s' * len (systems )) % tuple (systems ) )
64
66
w = []
65
67
for s in systems :
66
68
alarm (timeout )
@@ -77,8 +79,8 @@ def report(F, title, systems = ['sage', 'magma'], **kwds):
77
79
w .append (w [0 ]/ w [1 ])
78
80
79
81
w = tuple (w )
80
- print ( '%15.3f' * len (w ))% w
81
- print '=' * 70
82
+ print (( '%15.3f' * len (w )) % w )
83
+ print ( '=' * 70 )
82
84
83
85
84
86
#######################################################################
@@ -97,7 +99,7 @@ def report_ZZ(**kwds):
97
99
EXAMPLES::
98
100
99
101
sage: import sage.matrix.benchmark as b
100
- sage: print "starting"; import sys; sys.stdout.flush(); b.report_ZZ(systems=['sage']) # long time (15s on sage.math, 2012)
102
+ sage: print( "starting") ; import sys; sys.stdout.flush(); b.report_ZZ(systems=['sage']) # long time (15s on sage.math, 2012)
101
103
starting...
102
104
======================================================================
103
105
Dense benchmarks over ZZ
@@ -147,7 +149,7 @@ def nullspace_ZZ(n=200, min=0, max=2**32, system='sage'):
147
149
K := Kernel(A);
148
150
s := Cputime(t);
149
151
""" % (n ,min ,max )
150
- if verbose : print code
152
+ if verbose : print ( code )
151
153
magma .eval (code )
152
154
return float (magma .eval ('s' ))
153
155
else :
@@ -186,7 +188,7 @@ def charpoly_ZZ(n=100, min=0, max=9, system='sage'):
186
188
K := CharacteristicPolynomial(A);
187
189
s := Cputime(t);
188
190
""" % (n ,min ,max )
189
- if verbose : print code
191
+ if verbose : print ( code )
190
192
magma .eval (code )
191
193
return float (magma .eval ('s' ))
192
194
else :
@@ -225,7 +227,7 @@ def rank_ZZ(n=700, min=0, max=9, system='sage'):
225
227
K := Rank(A);
226
228
s := Cputime(t);
227
229
""" % (n ,min ,max )
228
- if verbose : print code
230
+ if verbose : print ( code )
229
231
magma .eval (code )
230
232
return float (magma .eval ('s' ))
231
233
else :
@@ -263,7 +265,7 @@ def rank2_ZZ(n=400, min=0, max=2**64, system='sage'):
263
265
K := Rank(A);
264
266
s := Cputime(t);
265
267
""" % (n ,min ,max )
266
- if verbose : print code
268
+ if verbose : print ( code )
267
269
magma .eval (code )
268
270
return float (magma .eval ('s' ))
269
271
else :
@@ -303,7 +305,7 @@ def smithform_ZZ(n=128, min=0, max=9, system='sage'):
303
305
K := ElementaryDivisors(A);
304
306
s := Cputime(t);
305
307
""" % (n ,min ,max )
306
- if verbose : print code
308
+ if verbose : print ( code )
307
309
magma .eval (code )
308
310
return float (magma .eval ('s' ))
309
311
else :
@@ -348,7 +350,7 @@ def matrix_multiply_ZZ(n=300, min=-9, max=9, system='sage', times=1):
348
350
end for;
349
351
s := Cputime(t);
350
352
""" % (n ,min ,max ,times )
351
- if verbose : print code
353
+ if verbose : print ( code )
352
354
magma .eval (code )
353
355
return float (magma .eval ('s' ))/ times
354
356
else :
@@ -394,7 +396,7 @@ def matrix_add_ZZ(n=200, min=-9, max=9, system='sage', times=50):
394
396
end for;
395
397
s := Cputime(t);
396
398
""" % (n ,min ,max ,times )
397
- if verbose : print code
399
+ if verbose : print ( code )
398
400
magma .eval (code )
399
401
return float (magma .eval ('s' ))/ times
400
402
else :
@@ -454,7 +456,7 @@ def det_ZZ(n=200, min=1, max=100, system='sage'):
454
456
d := Determinant(A);
455
457
s := Cputime(t);
456
458
""" % (n ,min ,max )
457
- if verbose : print code
459
+ if verbose : print ( code )
458
460
magma .eval (code )
459
461
return float (magma .eval ('s' ))
460
462
else :
@@ -494,7 +496,7 @@ def det_QQ(n=300, num_bound=10, den_bound=10, system='sage'):
494
496
d := Determinant(A);
495
497
s := Cputime(t);
496
498
""" % (n ,- num_bound , num_bound , den_bound )
497
- if verbose : print code
499
+ if verbose : print ( code )
498
500
magma .eval (code )
499
501
return float (magma .eval ('s' ))
500
502
else :
@@ -541,7 +543,7 @@ def vecmat_ZZ(n=300, min=-9, max=9, system='sage', times=200):
541
543
end for;
542
544
s := Cputime(t);
543
545
""" % (n ,min ,max ,times )
544
- if verbose : print code
546
+ if verbose : print ( code )
545
547
magma .eval (code )
546
548
return float (magma .eval ('s' ))/ times
547
549
else :
@@ -572,7 +574,7 @@ def report_GF(p=16411, **kwds):
572
574
EXAMPLES::
573
575
574
576
sage: import sage.matrix.benchmark as b
575
- sage: print "starting"; import sys; sys.stdout.flush(); b.report_GF(systems=['sage'])
577
+ sage: print( "starting") ; import sys; sys.stdout.flush(); b.report_GF(systems=['sage'])
576
578
starting...
577
579
======================================================================
578
580
Dense benchmarks over GF with prime 16411
@@ -617,7 +619,7 @@ def nullspace_GF(n=300, p=16411, system='sage'):
617
619
K := Kernel(A);
618
620
s := Cputime(t);
619
621
""" % (n ,p )
620
- if verbose : print code
622
+ if verbose : print ( code )
621
623
magma .eval (code )
622
624
return magma .eval ('s' )
623
625
else :
@@ -656,7 +658,7 @@ def charpoly_GF(n=100, p=16411, system='sage'):
656
658
K := CharacteristicPolynomial(A);
657
659
s := Cputime(t);
658
660
""" % (n ,p )
659
- if verbose : print code
661
+ if verbose : print ( code )
660
662
magma .eval (code )
661
663
return magma .eval ('s' )
662
664
else :
@@ -697,7 +699,7 @@ def matrix_add_GF(n=1000, p=16411, system='sage',times=100):
697
699
end for;
698
700
s := Cputime(t);
699
701
""" % (n ,p ,p ,times )
700
- if verbose : print code
702
+ if verbose : print ( code )
701
703
magma .eval (code )
702
704
return magma .eval ('s' )
703
705
else :
@@ -743,7 +745,7 @@ def matrix_multiply_GF(n=100, p=16411, system='sage', times=3):
743
745
end for;
744
746
s := Cputime(t);
745
747
""" % (n ,p ,times )
746
- if verbose : print code
748
+ if verbose : print ( code )
747
749
magma .eval (code )
748
750
return float (magma .eval ('s' ))/ times
749
751
else :
@@ -780,7 +782,7 @@ def rank_GF(n=500, p=16411, system='sage'):
780
782
K := Rank(A);
781
783
s := Cputime(t);
782
784
""" % (n ,p )
783
- if verbose : print code
785
+ if verbose : print ( code )
784
786
magma .eval (code )
785
787
return float (magma .eval ('s' ))
786
788
else :
@@ -816,7 +818,7 @@ def rank2_GF(n=500, p=16411, system='sage'):
816
818
K := Rank(A);
817
819
s := Cputime(t);
818
820
""" % (n ,p )
819
- if verbose : print code
821
+ if verbose : print ( code )
820
822
magma .eval (code )
821
823
return float (magma .eval ('s' ))
822
824
else :
@@ -853,7 +855,7 @@ def det_GF(n=400, p=16411 , system='sage'):
853
855
d := Determinant(A);
854
856
s := Cputime(t);
855
857
""" % (n ,p )
856
- if verbose : print code
858
+ if verbose : print ( code )
857
859
magma .eval (code )
858
860
return float (magma .eval ('s' ))
859
861
else :
@@ -915,7 +917,7 @@ def echelon_QQ(n=100, min=0, max=9, system='sage'):
915
917
K := EchelonForm(A);
916
918
s := Cputime(t);
917
919
""" % (n ,min ,max )
918
- if verbose : print code
920
+ if verbose : print ( code )
919
921
magma .eval (code )
920
922
return float (magma .eval ('s' ))
921
923
else :
@@ -954,7 +956,7 @@ def inverse_QQ(n=100, min=0, max=9, system='sage'):
954
956
K := A^(-1);
955
957
s := Cputime(t);
956
958
""" % (n ,min ,max )
957
- if verbose : print code
959
+ if verbose : print ( code )
958
960
magma .eval (code )
959
961
return float (magma .eval ('s' ))
960
962
else :
@@ -1000,7 +1002,7 @@ def matrix_multiply_QQ(n=100, bnd=2, system='sage', times=1):
1000
1002
end for;
1001
1003
s := Cputime(t);
1002
1004
""" % (n , A .name (), times )
1003
- if verbose : print code
1005
+ if verbose : print ( code )
1004
1006
magma .eval (code )
1005
1007
return float (magma .eval ('s' ))/ times
1006
1008
else :
@@ -1037,7 +1039,7 @@ def det_hilbert_QQ(n=80, system='sage'):
1037
1039
s := Cputime(tinit);
1038
1040
delete h;
1039
1041
""" % n
1040
- if verbose : print code
1042
+ if verbose : print ( code )
1041
1043
magma .eval (code )
1042
1044
return float (magma .eval ('s' ))
1043
1045
@@ -1071,7 +1073,7 @@ def invert_hilbert_QQ(n=40, system='sage'):
1071
1073
s := Cputime(tinit);
1072
1074
delete h;
1073
1075
""" % n
1074
- if verbose : print code
1076
+ if verbose : print ( code )
1075
1077
magma .eval (code )
1076
1078
return float (magma .eval ('s' ))
1077
1079
@@ -1114,7 +1116,7 @@ def MatrixVector_QQ(n=1000,h=100,system='sage',times=1):
1114
1116
end for;
1115
1117
s := Cputime(t);
1116
1118
""" % (n ,h ,times )
1117
- if verbose : print code
1119
+ if verbose : print ( code )
1118
1120
magma .eval (code )
1119
1121
return float (magma .eval ('s' ))
1120
1122
else :
@@ -1162,7 +1164,7 @@ def nullspace_RR(n=300, min=0, max=10, system='sage'):
1162
1164
K := Kernel(A);
1163
1165
s := Cputime(t);
1164
1166
""" % (n ,min ,max )
1165
- if verbose : print code
1167
+ if verbose : print ( code )
1166
1168
magma .eval (code )
1167
1169
return float (magma .eval ('s' ))
1168
1170
else :
@@ -1202,7 +1204,7 @@ def nullspace_RDF(n=300, min=0, max=10, system='sage'):
1202
1204
K := Kernel(A);
1203
1205
s := Cputime(t);
1204
1206
""" % (n ,min ,max )
1205
- if verbose : print code
1207
+ if verbose : print ( code )
1206
1208
magma .eval (code )
1207
1209
return float (magma .eval ('s' ))
1208
1210
else :
0 commit comments