[Octave-bug-tracker] [bug #48365] __run-test_suite_ FAIL for assert.m (tolerances)

[Rik]

Follow-up Comment #4, bug #48365 (project octave):

Is it possible to write a pure C++ function to test this Octave code?


x = [-40:0];
 y1 = (10.^x).*(10.^x);
 y2 = 10.^(2*x);
 assert (y1, y2, eps (y1));


It seems possible that we are also running in to issues with BLAS for the
multiplication.  Looking at y1


 y1 = (10.^x).*(10.^x);


I could imagine that the first '^' operator introduces 1 eps, the second power
operator introduces 1 eps, and the multiplication of the two produces a
further 1 eps for a total of 3 eps.


    _______________________________________________________

Reply to this item at:

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

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