楼层函数返回 EXC_BAD_ACCESS

Question

我使用的 cod 包含这些 sn-ps 代码。我用参数 2455343.50000 调用 ThetaG_JD，这只是一个儒略日期示例。每次我运行程序时，我都会在指定的行上收到一个 EXC_BAD_ACCESS。当使用 gdb 并打印出中间值并通过 floor 函数传递它们时，我没有收到任何错误，但是当使用 Frac() 时它总是返回错误。

double Frac(double arg)
{
    /* Returns fractional part of double argument */
    return arg - floor(arg);
}

double ThetaG_JD(double jd)
{
    /* Reference:  The 1992 Astronomical Almanac, page B6. */

    double UT=0, TU=0, GMST=0;
    //THIS LINE
    UT=Frac(jd+0.5);
    // THAT ONE ^^
    jd=jd-UT;
    TU=(jd-2451545.0)/36525;
    GMST=24110.54841+TU*(8640184.812866+TU*(0.093104-TU*6.2E-6));
    GMST=Modulus(GMST+secday*omega_E*UT,secday);

    return (twopi*GMST/secday);
}

Answer 1

EXC_BAD_ACCESS 让我有些费解，但这听起来有点像浮点异常。已经有一段时间了，但我记得在 x87 硬件上，您可以生成溢出/下溢/NaN，并且处理器不会让您知道异常，直到 next FP 操作可能在代码的完全不同的部分。你可以试试jd += 0.5 之类的东西，而不是调用 Frac，看看它是否仍然死机。

此外，x87 状态寄存器将能够向您显示是否存在错误状态，您应该能够在 gbd 中看到它。