#12103: Use MeatAxe as an optional back end for dense matrices over `GF(p^n)`, p
odd, n>1, `p^n<255`
-------------------------------------+-------------------------------------
       Reporter:  SimonKing          |        Owner:  jason, was
           Type:  defect             |       Status:  needs_work
       Priority:  major              |    Milestone:  sage-6.4
      Component:  packages:          |   Resolution:
  experimental                       |    Merged in:
       Keywords:  linear algebra,    |    Reviewers:
  MeatAxe                            |  Work issues:
        Authors:  Simon King         |       Commit:
Report Upstream:  None of the above  |  191477e697d5fb02c0e6bf7f8b80850e1092d4f6
  - read trac for reasoning.         |     Stopgaps:
         Branch:                     |
  u/SimonKing/meataxe                |
   Dependencies:  #19240             |
-------------------------------------+-------------------------------------

Comment (by SimonKing):

 It seems to work. With the commit that I am preparing now and that also
 contains some doctests, I obtain:
 {{{
 sage: M = MatrixSpace(GF(27,'x'),5).random_element(density=0.4)
 sage: M.rank()
 5
 sage: ~M
 [            1 x^2 + 2*x + 2             0             0             0]
 [    x^2 + 2*x     2*x^2 + 2           x^2             0   x^2 + x + 2]
 [      x^2 + 2           x^2         2*x^2     2*x^2 + 1       x^2 + 2]
 [            1       x^2 + x           2*x             0             0]
 [    2*x^2 + x             x             0             0             0]
 sage: M = MatrixSpace(GF(27,'x'),5).random_element(density=0.4)
 sage: M.rank()
 4
 sage: ~M
 ---------------------------------------------------------------------------
 ZeroDivisionError                         Traceback (most recent call
 last)
 <ipython-input-6-1888fd9efc66> in <module>()
 ----> 1 ~M

 /home/king/Sage/git/sage/src/sage/matrix/matrix_gfpn_dense.pyx in
 sage.matrix.matrix_gfpn_dense.Matrix_gfpn_dense.__invert__
 (build/cythonized/sage/matrix/matrix_gfpn_dense.c:12186)()
    1032         OUT._cache = {}
    1033         sig_on()
 -> 1034         OUT.Data = MatInverse(self.Data)
    1035         sig_off()
    1036         if OUT.Data != NULL:

 ZeroDivisionError: Division by zero in file matinv.c (line 50)
 }}}

--
Ticket URL: <http://trac.sagemath.org/ticket/12103#comment:128>
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.
For more options, visit https://groups.google.com/d/optout.

Reply via email to