On 28/04/2008, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> The multiplication of sigma and
>  sigma_inv should be 1 on diagonal and zero elsewhere. I am using
>  double and am expecting the epsilon of the result I get should be
>  small (< 1e-10) but actually some of them are very large (~0.001).

What is the condition number of your matrix? Can we see your matrix, please?

Thanks,
- Jordi G. H.


_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl

Reply via email to