Bug: det() gives false zero

[Charlie Dyson]

It must be stressed that this matrix is not singular and its det is
not even small, so zero is not expected.

Many thanks for your time and greeting from the isle of wight!


On 24/07/2008, John W. Eaton <jwe at bevo.che.wisc.edu> wrote:
> On 24-Jul-2008, Jaroslav Hajek wrote:
>
> | Looking at src/DLD-FUNCTIONS/det.cc, it seems that this is a feature -
> | if a numerically singular matrix is detected, zero is returned.
> | Perhaps John can comment on why he made this choice.
>
> If I was the one who did that, then I no longer remember.  If it is
> the wrong thing to do, then please submit a patch.  It should not be
> hard to do, as it appears that the Matrix::determinant method already
> returns the value you expect in this case.
>
> Thanks,
>
> jwe
>