Cómo encontrar el valor doble siguiente/anterior más próximo (numeric_limits :: épsilon para el número dado)

Question

El título es bastante autoexplicativo, la entrada tiene doble valor y quiero agregar/sustraer la menor cantidad posible.

Answer 1

Si su compilador implementa funciones matemáticas de C99/C++ 11, se puede utilizar el nextafter:

#include <cfloat> // DBL_MAX 
#include <cmath> // std::nextafter 

double x = 0.1; 

// next representable number after x in the direction of DBL_MAX 
double xPlusSmallest = std::nextafter(x, DBL_MAX); 

Incluso si su compilador no lo soporta, probablemente tiene un intrínseco para eso. (MSVC ha tenido _nextafter desde el año 2005, por ejemplo CCG probablemente lo implementa como estándar..)

Si el compilador no lo soporta, pero Boost está disponible para usted, usted puede hacer esto:

#include <boost/math/special_functions/next.hpp> // boost::float_next 

double x = 0.1; 

// next representable number after x 
double xPlusSmallest = boost::math::float_next(x); 

lo que equivale a esto (C99 emulando):

#include <boost/math/special_functions/next.hpp> // boost::nextafter 
#include <cfloat> // DBL_MAX 

double x = 0.1; 

// next representable number after x in the direction of DBL_MAX 
double xPlusSmallest = boost::math::nextafter(x, DBL_MAX); 

y si ninguno de los trabajos para usted, sólo tendrá que abrir una grieta en la cabecera Boost y copiarlo.

Answer 2

#define FLT_MIN   1.175494351e-38F  /* min positive value */ 
#define FLT_MAX   3.402823466e+38F  /* max value */ 
#define DBL_MIN   2.2250738585072014e-308 /* min positive value */ 
#define DBL_MAX   1.7976931348623158e+308 /* max value */ 

http://xona.com/2006/07/26.html

Answer 3

Aquí hay un truco muy sucio que no es legal y solo funciona si su plataforma usa flotadores IEEE754: La representación binaria del flotador se ordena de la misma manera que el valor flotante, por lo que puede incrementar la representación binaria:

double x = 1.25; 

uint64_t * const p = reinterpret_cast<uint64_t*>(&x); 

++*p; // undefined behaviour! but it gets the next value 

// now x has the next value 

se puede lograr el mismo efecto totalmente contrario por hacer las habituales gimnasia copia binaria para obtener un valor adecuado uint64_t. Asegúrate también de verificar cero, infinito y NaN correctamente.

Answer 4

¿Qué tal:

x += fabs(x) * std::numeric_limits<double>::epsilon();