On 9/8/11 11:00 AM, William Stein wrote:
On Thu, Sep 8, 2011 at 8:49 AM, Jason Grout<jason-s...@creativetrax.com>  wrote:
On 9/8/11 10:30 AM, Dima Pasechnik wrote:

we should also add that both A.solve_right and A.solve_left claim to
solve Ax=b...


Yep.  solve_left should be solve_right, IIRC.  See
http://trac.sagemath.org/sage_trac/ticket/7852.  I think this was an error
from when I did the RDF switch to a numpy backend.  As you can see, Rob
posted a patch to fix it (and enhanced a few other things as well).

Note that A2.solve_left(b) gives the same vector as A2.inverse()*b

So you made solve_left use the numpy backend and actually do
solve_right, and left solve_right to use the generic Gauss elimination
code, so neither is right.   Trac 7852 really needs to get in.
Release manager, merge that sucker.

Yes, that's also my conclusion after trying to unify our two answers.

Jason

--
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to 
sage-devel+unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org

Reply via email to