[Octave-bug-tracker] [bug #55587] splinefit.m BIST may fail, slightly out of tolerance result

[Mike Miller]

Follow-up Comment #2, bug #55587 (project octave):

Here is a specific example that fails for me on both stable and default. I
inserted a 'save' into the BIST to get the exact yb vector that causes a
failure


>> xb = 0:2:10;
>> yb = [-0.85659524457527947 -0.43282809207909179 -0.82873770262843049
0.85540706303396452 -0.93702483313540752 -0.46131091926622292];
>> x = 0:0.1:10;
>> y = interp1 (xb, yb, x, "spline");
>> assert (ppval (splinefit (x, y, xb, "order", 3), x), y, 15 * eps ());
error: ASSERT errors for:  assert (ppval (splinefit (x, y, xb, "order", 3),
x),y,15 * eps ())

  Location  |  Observed  |  Expected  |  Reason
   (101)       -0.46131     -0.46131     Abs err 3.4417e-15 exceeds tol
3.3307e-15 by 1e-16


This was found on the default branch using the following


i = 0;
n = 3;
while (n == 3)
  i++;
  [n, ~] = test ('splinefit.m');
endwhile
i

output:
i =  5361


    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?55587>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/