在面向对象的 C++ 中处理大数字变量

Question

我正在使用双精度编译的求解器（编程基于面向对象的 C++）进行数值分析，我的单位是 64 位。我的问题是，当求解器计算一个大数时——比如说 -1.45 的 21 次方，举一个实际的例子——并通过将此值传递给现有变量将其堆叠在分配的内存中，它被转换为 0。所以当然，当我后来在一个部门中使用这个变量时，我得到了一个分段错误。我不明白这个过程是如何工作的，而且因为我使用的是 DP，所以我不知道如何解决这个问题。有人可以帮我解决这个问题吗？

如果有帮助：我刚刚运行了一个测试，我声明 a=-1.45e+21，并“打印”求解器正确返回的值。但是当我不使用“e”指数并输入完整值（有 19 个零）时，我得到 0 作为回报。所以我想问题/限制来自位数，有什么想法吗？谢谢 ！

编辑：我发布了我为计算造成问题的变量之一所经历的步骤的摘要。其他类似定义。 首先我初始化字段指针列表：

PtrList<volScalarField> fInvFluids(fluidRegions.size());

volScalarField 类只是一个双精度数组。然后我填充字段指针列表：

fInvFluids.set
    (
        i,
        new volScalarField
        (
            IOobject
            (
                "fInv",
                runTime.timeName(),
                fluidRegions[i],
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            fluidRegions[i],
            dimensionedScalar
            (
                "fInv",
                dimensionSet(3,1,-9,-1,0,0,0),
                scalar(0)
            )
        )
    );

在此之后，我设置了字段区域：

volScalarField& fInv = fInvFluids[i];

最后我计算了值：

   //   Info<< " ** Calculating fInv            **\n";
    fInv = gT*pow(pow(T/Tlambda, 5.7)*(1 - pow(T/Tlambda, 5.7)), 3);

其中 T 是字段变量，Tlambda 是在运行时定义的标量值。

Answer 1

double 变量可能也不能容纳 19 个零。一个（十进制）数字占用超过 3 位，因此 19 个零将占用至少 57 位。 double 的尾数通常只有 53 位。

但是，这听起来不像您遇到的问题。在 C++ 代码 中，表达式也有类型。 1 与 1.0 不同。第一个是int，第二个是double。虽然您可以将int 转换为double，但它们并不相同。 int 很可能可以容纳高达 20 亿的值，但正式它可能有一个低至 32767 的限制。解决您的问题可能就像添加一个 第二十 个零一样简单：@987654330 @使它成为double。