[Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps

[Nir Krakauer]

Follow-up Comment #23, bug #47800 (project octave):

I get 3 tolerances not met with the new D on Octave 4.0.2:

***** test
 assert(D(10,10),0.12511003572113329898476497894772544708420990097708600,-7*eps)
!!!!! test failed
ASSERT errors for:  assert (D (10,
10),0.12511003572113329898476497894772544708420990097708600,-7 * eps)

  Location  |  Observed  |  Expected  |  Reason
     ()        0.12511      0.12511      Rel err 3.3277e-15 exceeds tol
1.5543e-15

***** test
 assert(D(10,-10),6.069903492836947893782490257188118487200593203967e7,-8*eps)
!!!!! test failed
ASSERT errors for:  assert (D (10,
-10),6.069903492836947893782490257188118487200593203967e7,-8 * eps)

  Location  |  Observed  |  Expected  |  Reason
     ()      60699034.928 60699034.928   Rel err 5.1553e-15 exceeds tol
1.7764e-15

***** test
 assert(D(9,-10),-6.069903492836947893782490257188118487200593203967e7,-8*eps)
!!!!! test failed
ASSERT errors for:  assert (D (9,
-10),-6.069903492836947893782490257188118487200593203967e7,-8 * eps)

  Location  |  Observed  |  Expected  |  Reason
     ()      -60699034.928 -60699034.928   Rel err 5.1553e-15 exceeds tol
1.7764e-15

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?47800>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/