I think this might need some work: S=RealBallField(100) M=Matrix(S,2,1,[1,1]) M.solve_right(vector([1,2]))
There's enough information here to conclude there is no solution; or in a rather deranged way, perhaps it should give a rather large ball back so that the multiplication results in a vector consisting of balls that actually contain the values 1 and 2. Currently, it just returns a wrong answer. For systems that are even generically overdetermined, I don't think it makes sense to default to a least squares solution in a routine that is supposed to detect inconsistent systems. There the stable numerical behaviour actually IS being inconsistent. -- 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 sage-devel+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sage-devel/0be0053d-501b-442c-bd3c-23e662b4e49c%40googlegroups.com.
