e^x 的无限和收敛到不正确的值

我正在尝试实现一种不使用阶乘来计算 e^x 的方法。我通过获取每两个连续项之间的比率来做到这一点，并通过将该比率与上一项相乘来计算下一项。所有这些都会添加到结果总和中，直到该项足够小而无关紧要。

对于较小的 x 值，这是有效的。为了“极端”x 的值，这会由于某种原因而中断。我尝试了 -50，它应该打印 1.9287498e-22 （以固定形式）。 我得到的是 GNU GCC 中的 -56676.4235303065 和 VC++ 中的 2041.8329628977 。两者都大错特错。这是我的代码：

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;

int main()
{
    double x;
    cout << fixed << setprecision(10) ;
    cout << "Enter x: " << endl;
    cin >> x;
    double sum = 0;
    long long i = 0;
    double term = 1;
    do
    {
        sum += term;
        term *= (x / (++i));
    } while (fabs(term) > 1e-10);
    cout << sum << endl;
}

这是我的问题：导致出现此错误的问题是什么only对于 x 的大值，我该如何修复这个错误？为什么 GNU GCC 和 VC++2017 得到完全不同的答案？

该求和中的最大项比精度高约 10 个数量级double: ~ 2.92e+20相比于double13 平方英尺或者。这些术语中的误差幅度本身比求和结果大约 30 个数量级。

因此，您的级数无法正确收敛也就不足为奇了，因为连续项（相反符号）可能不会按理论量抵消。即使使用一些数值技巧，例如 Kahan-Neumaier 求和以及在相加之前对项进行排序，结果仍然只能简化为~6000。请注意，对于积极的情况，不会发生这种情况x因为连续的任期不需要取消。

克服这个问题的一种方法是对x，并使用乘方求幂来放大到正确的值x.

更新：上述方法的实现：

// integer exponentiation by squaring (won't explain here)
double pow_square(double x, unsigned a)
{
   double r = 1.0;
   while (a > 0)
   {
      if (a % 2 == 1)
      {
         a--;
         r *= x;
      }
      a /= 2;
      x *= x;
   }
   return r;
}

// original method
double exp_original(double x, double e)
{
   double sum = 0.0;
   unsigned i = 0;
   double term = 1.0;
   do
   {
      sum += term;
      term *= (x / (++i));
   } while (fabs(term) > e);
   return sum;
}

// new adaptive method
double exp_new(double x, double e)
{
   static const double min_X = -3;

   // if within limit, simply use original function
   if (x >= min_X)
      return exp_original(x, e);

   // compute smallest possible scaling coefficient
   unsigned s = (unsigned)(x / (-min_X) + 0.5);
   double p = exp_original(x / s, e);
   return pow_square(p, s);
}

测试范围大x值证实了新方法处理极值（negative）情况要好得多：

x    | exp (C-library)    exp_original       exp_new
-------------------------------------------------------------
-10  | 4.53999297625e-05  4.53998989141e-05  4.53999299001e-05
-30  | 9.35762296884e-14  6.10299992426e-06  9.35762292245e-14
-50  | 1.92874984796e-22  2041.8329629       1.92874983803e-22
-60  | 8.7565107627e-27   722745700.93       8.75651067587e-27
-80  | 1.80485138785e-35  2.45082011705e+17  1.80485137011e-35
-100 | 3.72007597602e-44  8.1446527451e+25   3.72007589785e-44
-150 | 7.17509597316e-66  -9.14622659954e+47 7.1750957953e-66
-200 | 1.38389652674e-87  7.69097143891e+69  1.38389648613e-87