#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.


Reply via email to