2

我有一个在 PC 上运行的用 Visual C++ 编写的由来已久的例程，它根据纬度/经度坐标计算地球表面上两点之间的距离。

我已将其转换为 Java 作为 Android 移植的一部分，但我得到了不同的结果。为了调查其原因，我将一个复杂的三角方程拆分为多个位，以跟踪每一步发生的情况。

但是，在 Visual C++ 中，逐步计算的结果与复杂语句得到的结果不同，尽管它们显然应该是相同的。我想知道为什么！

正如预期的那样，Java 版本在这两种情况下都会产生相同的结果。

这是 Visual C++ 代码：

double CCoord::GetDistanceRadians(double dLatRadians, double dLongRadians)
{
    double dResult;

    // This is the step by step code
    double r1, r2, r3, r4, r5;
    r1 = sin(m_dLatRadians);
    r2 = sin(dLatRadians);
    r3 = cos(dLatRadians);
    r4 = m_dLongRadians - dLongRadians;
    r5 = cos(r4);

    r1 *= r2;
    r2 = r3 * r3;
    r2 *= r5;
    r4 = r1 + r2;

    // This is the original statement: r4 and dResult should be equal
    // (but they aren't)
    dResult = sin(m_dLatRadians) * sin(dLatRadians) +
        cos(m_dLatRadians) * cos(dLatRadians) * cos(m_dLongRadians - dLongRadians);

    // Protect against maths processor inaccuracies
    if ((dResult > 1.0) || (dResult < -1.0))
        dResult = 0.0;
    else
        dResult = acos(dResult);

    return dResult;
}

这是Java代码：

    // Return angular distance in radians
private double getDistanceRadians(double dLatRadians, double dLongRadians)
{
    double result;

    // Do JIT test
    checkRadianValues();

            // This is the step by step code
    double r1, r2, r3, r4, r5;
    r1 = Math.sin(m_dLatRadians);
    r2 = Math.sin(dLatRadians);
    r3 = Math.cos(dLatRadians);
    r4 = m_dLongRadians - dLongRadians;
    r5 = Math.cos(r4);

    r1 *= r2;
    r2 = r3 * r3;
    r2 *= r5;
    r4 = r1 + r2;

    // This is the original statement: r4 and result should be equal
    // (and they are)
    result = Math.sin(m_dLatRadians) * Math.sin(dLatRadians) +
                Math.cos(dLatRadians) * Math.cos(dLatRadians) *
                Math.cos(m_dLongRadians - dLongRadians);

    // Protect against maths processor inaccuracies
    if ((result > 1.0d) || (result < -1.0d))
        result = 0.0d;
    else
        result = Math.acos(result);

    return result;
}

这些是输入和结果值：

Variable         VC++ value         Java value
m_dLatRadians    0.896808347        0.896808378
m_dLongRadians  -0.047414778(0)    -0.047414778(6)
dLatRadians      0.896192633(9)     0.896192588(9)
dLongRadians    -0.026897463       -0.026897463
r4               1.000218106        1.0002181433
dResult/result   0.999917766        1.0002181433

请任何人都可以建议为什么 r4 和 dResult 在 Visual C++ 版本中应该不同？

PS 我从 VC++ 版本中获得了良好的最终距离结果，而从 Java 中获得了糟糕的最终结果 - 几乎所有距离都为零（调用这些例程的更高级别的代码没有显示，因为它有很多！但是这是差异我有兴趣了解）。

4

1 回答 1

1

在 C++ 代码中

cos(m_dLatRadians) * cos(dLatRadians)

应该

cos(dLatRadians) * cos(dLatRadians)

r2 = r3 * r3;与您在 Java 中所做的事情保持一致并匹配。

此外，正如@CaseyB指出的那样，dLatRadians参数是不同的（就像 一样m_dLatRadians）。

于 2012-09-19T21:07:21.327 回答