我有一个C++函数，它基于一个简单模型计算概率。看起来C++倾向于将非常小的概率四舍五入到0，将非常大的概率四舍五入到1。这在后续计算中（如取log(p)和log(1-p)）会导致问题。

有没有一种方法可以明确表示C++内部可以表示而不进行四舍五入的最大小于1的数字？同样，最小大于0的数字呢？

我可以这样做：

if (probability == 1)     probability = 0.999999999;else if (probability == 0)     probability = 0.000000001;

但这会导致其他数值问题（与累积数字有关）。有没有更原则性的方法，或许使用numeric_limits？

回答：

尽管它的名字可能有点误导，但std::numeric_limits<T>::min正是你所寻找的。对于浮点类型，它会给你大于零的最小值。

对于仍然小于1的最大数字，如果你使用C++11，可以使用std::nexttoward：

代码

#include<cstdio>#include<cmath>#include<limits>int main(){  printf("Near 0: %1.20e\nNear 1: %1.20e\n",    std::numeric_limits<double>::min(),    std::nexttoward(1.0, 0.0)  );  return 0;} 

结果

Near 0: 2.22507385850720138309e-308Near 1: 9.99999999999999888978e-01