Differences between revisions 4 and 5
⇤ ← Revision 4 as of 2008-02-08 13:37:45 →
Size: 990
Editor: was
Comment:
← Revision 5 as of 2008-02-08 13:38:37 → ⇥
Size: 946
Editor: was
Comment:
Deletions are marked like this. Additions are marked like this.
Line 46: Line 46:

---- /!\ '''End of edit conflict''' ----

Fast Hermite Normal Form over ZZ

Benchmark 1: Random 200x200 single-digit matrix

Mathematica code: 

sage: mathematica.eval('a = Table[RandomInteger[{0,9}], {i,200}, {j,200}];')

sage: mathematica.eval('Timing[HermiteDecomposition[a];]')
        {98.9791, Null}

Maple code: 

> with(LinearAlgebra); n := 200: k := 9: A := RandomMatrix(n,n,generator=rand(-k..k)): time( HermiteForm(A,output=['H', 'U']) );

60.431 seconds

Sage (via PARI): 

sage: a = random_matrix(ZZ,200)
sage: time v = a.echelon_form()
CPU times: user 43.72 s, sys: 0.27 s, total: 43.98 s
Wall time: 44.36

sage: a = random_matrix(ZZ,200)
sage: z = magma(a)
sage: t = magma.cputime()
sage: time w = z.HermiteForm()
CPU times: user 0.00 s, sys: 0.00 s, total: 0.00 s
Wall time: 0.82
sage: magma.cputime(t)
0.68999999999999995

Benchmark 1: Random 200x200 matrix with 32-bit entries

days7/sprints/linalg (last edited 2008-11-14 13:42:08 by anonymous)