Comparación de punto flotante en STL, BOOST

Question

¿Hay en el STL o en Boost un conjunto de funciones simples de comparación genéricas?

El que encontré siempre requieren parámetros de plantilla y/o instanciación de una plantilla de estructura .

estoy buscando algo con una sintaxis como:

if (is_equal(x,y)) 
{ 
    ... 
} 

que podría ser implementado como:

template <typename T> 
bool is_equal(const T& x, const T& y) 
{ 
    return (fabs(x - y) < Precision<T>::eps); 
} 

EDIT: He cambiado el operador de igualdad. (Ver comentarios abajo)

Answer 1

no sé de cualquier biblioteca que lo hace, tal vez porque es tan simple como una sola línea o tal vez porque se olvida ...

Como generalidad va sin embargo, está seguro de que' Me gusta configurar el épsilon para un tipo dado a un valor dado ... ¿a lo largo de la aplicación? Personalmente, me gustaría personalizarlo en función de las operaciones que estoy realizando (aunque un valor predeterminado sería bueno).

En cuanto a sus operadores, ¿por qué no idearlos usted mismo?

template <class T> 
bool rough_eq(T lhs, T rhs, T epsilon = Precision<T>::epsilon) // operator== 
{ 
    return fabs(lhs - rhs) < epsilon; 
} 

template <class T> 
bool rough_lt(T lhs, T rhs, T epsilon = Precision<T>::epsilon) // operator< 
{ 
    return rhs - lhs >= epsilon; 
     // tricky >= because if the difference is equal to epsilon 
     // then they are not equal per the rough_eq method 
} 

template <class T> 
bool rough_lte(T lhs, T rhs, T epsilon = Precision<T>::epsilon) // operator<= 
{ 
    return rhs - lhs > -epsilon; 
} 

La desigualdad y más que los métodos se pueden derivar trivialmente de esto.

El parámetro adicional significa que puede desear especificar otro valor para un conjunto dado de cálculos ... una configuración de toda la aplicación es demasiado estricta.

Answer 2

Desde el comentario de Marcelo Cantos:

... entonces también necesita una prueba que < vuelve falsa, incluso si el primer número es siempre tan ligeramente más pequeña que el segundo.

Me imagino que la implementación sería:

return !roughly_equal(a, b) && a < b; 

Answer 3

se puede encontrar algunas de las razones de la lógica comparacion "complejo" , en la documentación de la biblioteca de prueba de impulso.