Question

的Boost实施中

enter image description here

在标题<boost/math/special_functions/erf.hpp>中

result = z * 1.125 + z * 0.003379167095512573896158903121545171688L;

在下面附带的代码段中。我的问题是：为什么没有

→

此添加的第一个组件中的长后缀，而第二个中有一个？它背后的理由是什么？它引入了什么含义？

template <class T, class Policy>
T erf_imp(T z, bool invert, const Policy& pol, const mpl::int_<113>& t)
{
   BOOST_MATH_STD_USING
   BOOST_MATH_INSTRUMENT_CODE("113-bit precision erf_imp called");

   if(z < 0){
      if(!invert)
         return -erf_imp(-z, invert, pol, t);
      else if(z < -0.5)
         return 2 - erf_imp(-z, invert, pol, t);
      else
         return 1 + erf_imp(-z, false, pol, t);
   }

   T result;

   // Big bunch of selection statements now
   if(z < 0.5){
      // We're going to calculate erf:
      if(z == 0){
         result = 0;
      }else if(z < 1e-20){
         result = z * 1.125 + z * 0.003379167095512573896158903121545171688L;
                          ^^^ no L?                                       ^^^

Answer 1

如果z已经是long double，那么乘法z * 1.125已经是long double次乘法。常量1.125可以完全表示为double，它的类型。事实上，它可以完全代表float。因此，无需为此常量指示类型long double。如果z具有更宽的浮点类型，则无论如何都会在乘法之前提升常量。

相比之下，以十进制写的实数0.003379167095512573896158903121545171688在任何二进制浮点精度中都不能完全表示。如果0.003379167095512573896158903121545171688L在程序中写入0.003379167095512573896158903121545171688，它将表示与指示的十进制序列最接近的double，这与最近的long double不同，并且远离真正的价值。

您可以观察下面C程序的不同之处：

#include <stdio.h>

int main(int c, char **v)
{
  printf("%s\n%.24Lf\n%.24Lf\n\n%La\n%La\n", 
     "0.003379167095512573896158903121545171688",
     (long double) 0.003379167095512573896158903121545171688,
     0.003379167095512573896158903121545171688L,
     (long double) 0.003379167095512573896158903121545171688,
     0.003379167095512573896158903121545171688L);
}

结果：

0.003379167095512573896158903121545171688
0.003379167095512573739530
0.003379167095512573896231

0xd.d750429b6d118p-12
0xd.d750429b6d11ae4p-12

boost erf（）中的L（长）后缀：何时需要，何时不需要？

1 个答案: