Mercurial > hg > octave-nkf

comparison libinterp/corefcn/__contourc__.cc @ 15220:61822c866ba1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
use std::numeric_limits<T>::epsilon in C++ code * __contourc__.cc, __qp__.cc, quadcc.cc, qz.cc, eigs.cc, mex.cc, graphics.cc, CMatrix.cc, CollocWt.cc, DASPK-opts.in, DASRT-opts.in, DASSL-opts.in, LSODE-opts.in, Quad-opts.in, Range.cc, dMatrix.cc, eigs-base.cc, fCMatrix.cc, fMatrix.cc: Replace DBL_EPSILON and FLT_EPSILON with std::numeric_limits<T>::epsilon.
author John W. Eaton <jwe@octave.org>
date Thu, 23 Aug 2012 12:13:59 -0400
parents 2fc554ffbc28
children 704e15f8fecd
comparison
equal deleted inserted replaced
15219:aeba1adfd76b 15220:61822c866ba1
227 f[1] = Z(r, c+1) - lvl; 227 f[1] = Z(r, c+1) - lvl;
228 f[3] = Z(r+1, c) - lvl; 228 f[3] = Z(r+1, c) - lvl;
229 f[2] = Z(r+1, c+1) - lvl; 229 f[2] = Z(r+1, c+1) - lvl;
230 230
231 for (unsigned int i = 0; i < 4; i++) 231 for (unsigned int i = 0; i < 4; i++)
232 if (fabs(f[i]) < DBL_EPSILON) 232 if (fabs(f[i]) < std::numeric_limits<double>::epsilon ())
233 f[i] = DBL_EPSILON; 233 f[i] = std::numeric_limits<double>::epsilon ();
234 234
235 if (f[1] * f[2] < 0) 235 if (f[1] * f[2] < 0)
236 mark(r, c) += 2; 236 mark(r, c) += 2;
237 237
238 if (f[0] * f[3] < 0) 238 if (f[0] * f[3] < 0)
246 f[1] = Z(r, c+1) - lvl; 246 f[1] = Z(r, c+1) - lvl;
247 f[3] = Z(r+1, c) - lvl; 247 f[3] = Z(r+1, c) - lvl;
248 f[2] = Z(r+1, c+1) - lvl; 248 f[2] = Z(r+1, c+1) - lvl;
249 249
250 for (unsigned int i = 0; i < 4; i++) 250 for (unsigned int i = 0; i < 4; i++)
251 if (fabs(f[i]) < DBL_EPSILON) 251 if (fabs(f[i]) < std::numeric_limits<double>::epsilon ())
252 f[i] = DBL_EPSILON; 252 f[i] = std::numeric_limits<double>::epsilon ();
253 253
254 if (f[0] * f[1] < 0) 254 if (f[0] * f[1] < 0)
255 mark(r, c) += 1; 255 mark(r, c) += 1;
256 256
257 if (f[2] * f[3] < 0) 257 if (f[2] * f[3] < 0)