Sorry for replying to myself but after a bit of further digging it appears that numpy didn't have anything to do with this. It was just the route by which I came across the following problem:
sage: C = random_matrix(ZZ, 10, 80, distribution='uniform') sage: C.ncols() - (C.right_kernel().dimension() + C.rank()) -10 On 32bit sage 4.7 on osx the answer is the correct 0. Switching ZZ to QQ or RDF also restores the rank-nullity theorem. Thanks, -- Vegard Lima -- To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
