Re: boosting tolerances on gammainc BIST tests

[Mike Miller]

On Tue, Feb 05, 2019 at 09:46:25 -0800, Rik wrote:
> @Mike: No objection from me to relaxing the tolerances if the miss is only
> a further few eps.  If it is larger, maybe 10 eps, then it probably
> deserves a closer look.

The gammainc.m test failures on Christian's system

     ()         0.4906       0.4906      Rel err 1.1315e-15 exceeds tol 
1.1102e-15 by 2e-17
     ()        17.3984      17.3984      Rel err 6.1259e-16 exceeds tol 
4.4409e-16 by 2e-16
     ()        18.0654      18.0654      Rel err 1.5733e-15 exceeds tol 
1.5543e-15 by 2e-17

are from these three tests

%!assert (gammainc (200, 200, "upper"), 0.49059658199276367497217454, -5*eps)
%!assert (gammainc(1000, 1000.5, "scaledlower"), 39.48467539583672271, -2*eps)
%!assert (gammainc (200, 200, "scaledlower"), 18.065406676779221643065, -7*eps)

and look like they will all be fixed by increasing by 1 eps, I'll make
that change on stable soon.

> If you switch back to OpenBLAS from FlexiBLAS do the errors go away?

It should be a pass through, I used the same combination of FlexiBLAS
2.0.0 and OpenBLAS 0.3.5 and was not able to reproduce any of the fails.

-- 
mike

signature.asc
Description: PGP signature