On Tue, Jul 31, 2012 at 12:31 PM, Daniel J Sebald
<address@hidden> wrote:
On 07/31/2012 02:22 PM, Ed Meyer wrote:
As I pointed out in the "make check" thread the problem is that we
should not be
using absolute tolerances because that does not account for the size of
the matrix
elements. What we should do is to estimate if the result lies in the
interval we would
get with roundoff perturbations in the matrix. If it does, that's as
good as can be expected.
Eispack used to do something like that with their "performance index"
but I think you can
get better estimates.
Give us an idea of what code you are suggesting Ed. Something like
eps*size(s,1)*size(s,2)
? Something concerning the rank of the input matrix?
Dan
PS: Did you mean to CC to the maintainers list?