Mercurial > hg > octave-nkf
diff scripts/plot/draw/fplot.m @ 18775:37c300acfcfd stable
don't truncate when casting char to uintN values (bug #42054)
* oct-inttypes.h (octave_int<T>::octave_int (char)):
New special case constructor.
* oct-inttypes.cc: New tests.
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Sat, 05 Apr 2014 15:09:09 -0400 |
| parents | d63878346099 |
| children | 6fdd3ab55b78 |
line wrap: on
line diff
--- a/scripts/plot/draw/fplot.m +++ b/scripts/plot/draw/fplot.m @@ -146,10 +146,10 @@ ## problems with the current solution. while (n < 2^18) # Something is wrong if we need more than 250K points - yi = interp1 (x0, y0, x, "linear"); + yi = interp1 (x0, y0, x, "linear") ## relative error calculation using average of [yi,y] as reference ## since neither estimate is known a priori to be better than the other. - err = 0.5 * max (abs ((yi - y) ./ (yi + y))(:)); + err = 0.5 * max (abs ((yi - y) ./ (yi + y))(:)) if (err < tol || abs (err - err0) < tol/2) ## Either relative tolerance has been met OR ## algorithm has stopped making any reasonable progress per iteration.