什么时候会: std::abs(x -y) < std::numeric_limits:<double>:min()
When will: std::abs(x -y) < std::numeric_limits<double>::min()
在此链接中,有一个函数可以检查2浮点值的平等:
template<class T>
typename std::enable_if<!std::numeric_limits<T>::is_integer, bool>::type
almost_equal(T x, T y, int ulp)
{
// the machine epsilon has to be scaled to the magnitude of the values used
// and multiplied by the desired precision in ULPs (units in the last place)
return std::abs(x-y) <= std::numeric_limits<T>::epsilon() * std::abs(x+y) * ulp
// unless the result is subnormal
|| std::abs(x-y) < std::numeric_limits<T>::min();
}
但是我不完全了解std::abs(x-y) < std::numeric_limits<T>::min()
实际上会发生什么?有例子吗?谢谢。
std::numeric_limits<T>::min()
返回可以表示的最小的"正常"浮点值。IEEE754可以将值表示为0和std::numeric_limits<T>::min()
之间,为"亚正常"数字。看看这个问题,该问题有几个解释的答案。
您可以轻松地生成一个亚正态值:
// Example program
#include <iostream>
#include <limits>
#include <cmath>
int main()
{
std::cout << "double min: " << std::numeric_limits<double>::min() << " subnormal min: ";
std::cout << std::numeric_limits<double>::denorm_min() << 'n';
double superSmall = std::numeric_limits<double>::min();
std::cout << std::boolalpha << "superSmall is normal: " << std::isnormal(superSmall) << 'n';
double evenSmaller = superSmall/2.0;
std::cout << "evenSmaller = " << evenSmaller << 'n';
std::cout << std::boolalpha << "evenSmaller is normal: " << std::isnormal(evenSmaller);
std::cout << std::boolalpha << "std::abs(superSmall - evenSmaller) < std::numeric_limits<double>::min(): " << (std::abs(superSmall - evenSmaller) < std::numeric_limits<double>::min());
}
在我的计算机上生成以下内容:
double min: 2.22507e-308 subnormal min: 4.94066e-324
superSmall is normal: true
evenSmaller = 1.11254e-308
evenSmaller is normal: false
std::abs(superSmall - evenSmaller) < std::numeric_limits<double>::min(): true
相关文章:
- 使用std::multimap迭代器创建std::list
- C++中std::resize(n)和std::shrink_to_fit之间的区别
- 来自 std::list 的迭代器 .end() 按预期返回"0xcdcdcdcdcdcdcdcd"但 .begin()
- 错误:调用"std::vector<:vector<int>>::p ush_back(std::vector<std::__cxx11::basic_string<
- C 建造者Clang STD :: Sill,找不到超载的操作员&lt;
- 为什么STD :: MAP需要操作员&lt;以及我如何写一个
- std::vector::reserve(未知m),我知道m<<;N(通常)并且知道N
- std::vector<;uint8_t>;当C++11/14启用时,手动复制而不是调用memcpy
- C++运算符<<调用::ostream而不是std::osttream
- 是std :: set&lt; std :: future&gt;不可能存在
- 在修改列表后,std :: list&lt; t&gt; :: end()的值是否会更改
- BOOST ::变体无法解决运算符&lt;&lt;对于STD :: Ostream
- C++重载<<具有typedef'd std::vector
- 以x的倍数填充前导零,使用std::cout<<std::十六进制
- 错误:没有匹配'运算符<<"在'std::cout
- std::pair的默认构造函数<>将基本类型(int等)设置为零
- std::ostream&operator< & lt; (std:: ostream&压力,压力& &;val)
- 将std::endl传递给std::operator<<
- std::映射<>或std::vector<>在处理大型标志集时
- 重载& lt; & lt;使用命名空间std