#14007: When computing determinant over GF(p), don't lift to ZZ
----------------------------------+-----------------------------------------
Reporter: jdemeyer | Owner: jason, was
Type: enhancement | Status: needs_work
Priority: blocker | Milestone: sage-5.7
Component: linear algebra | Resolution:
Keywords: | Work issues:
Report Upstream: N/A | Reviewers: Charles Bouillaguet
Authors: Jeroen Demeyer | Merged in:
Dependencies: | Stopgaps:
----------------------------------+-----------------------------------------
Comment (by slabbe):
with OSX 10.5.8, 32-bit, sage-5.6.rc0 + 14007_matrix_modp.patch :
{{{
sage: p = next_prime(2^32-60)
sage: p < 2^32
True
sage: M = random_matrix(Integers(p), 75*5)
sage: M2 = matrix(ZZ, M)
sage: time d = M.det()
Time: CPU 7.49 s, Wall: 7.53 s
sage: time d = M2.det()
Time: CPU 9.55 s, Wall: 9.49 s
sage: p = next_prime(2^32)
sage: M = random_matrix(Integers(p), 75*5)
sage: M2 = matrix(ZZ, M)
sage: time d = M.det()
Time: CPU 14.15 s, Wall: 14.59 s
sage: time d = M2.det()
Time: CPU 9.81 s, Wall: 9.72 s
sage: p = next_prime(2^64-60)
sage: M = random_matrix(Integers(p), 75*5)
sage: M2 = matrix(ZZ, M)
sage: time d = M.det()
Time: CPU 15.82 s, Wall: 16.08 s
sage: time d = M2.det()
Time: CPU 14.91 s, Wall: 15.72 s
sage: p = next_prime(2^64)
sage: M = random_matrix(Integers(p), 75*5)
sage: M2 = matrix(ZZ, M)
sage: time d = M.det()
Time: CPU 18.14 s, Wall: 18.23 s
sage: time d = M2.det()
Time: CPU 14.59 s, Wall: 14.63 s
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/14007#comment:10>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.