On Tuesday, January 30, 2018 at 3:50:45 PM UTC, Simon Brandhorst wrote: > > sage: M = Matrix(R,4,[0, 0, 1, 1, 2^19, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1]) > sage: M.inverse() > [ 1048575 2 0 1] > [274878955520 524288 1 1048575] > [ 524289 524287 1 524287] > [274877382656 524289 1048575 524289] > > Works. > Only in the sense that you don't get an error:
sage: min([a.valuation() for a in (M*M^(-1)-M^0).list()]) 1 that's quite far from what it should be. With appropriate pivoting of course: sage: min([a.valuation() for a in (M*(P*(M*P)^(-1))-M^0).list()]) 20 -- You received this message because you are subscribed to the Google Groups "sage-devel" 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 https://groups.google.com/group/sage-devel. For more options, visit https://groups.google.com/d/optout.
