在 Java 中使用 double 计算视锥体的侧面时如何处理浮点精度

问：

我正在尝试使用下面的代码计算视锥体的侧表面。

输入：2 个节点，包含 x、y、z 值（绝对位置）和半径（位置半径）均为双精度

到目前为止我正在做什么：
1. 计算视锥体的长度
2.计算视锥体的侧表面

我的问题：
假设减去 2 个具有相似大小的浮点数是有问题的，因为精度损失。（我没有遇到这个问题）

我的问题：
我能做些什么来改善最终结果？

我发现/想到的可能性：
- 使用 BigDecimal（我不想这样做，因为运行时间更长） - 将 （r1-r2）*（r1-r2）
替换为 r1^2 - 2*r1*r2 + r2^2
- 检查两个双精度值的接近程度，如果非常接近，则假设它们的差值为 0。这甚至会提高我最终结果的准确性吗？减法结果的绝对误差不是比假设为 0 的结果小吗？

public static double calculateLateralSurface(Node node1, Node node2) {
    double x, y, z, length, r1, r2, lateralSurface;

    // calculate length of frustum cone
    x = node1.getPosX() - node2.getPosX();
    y = node1.getPosY() - node2.getPosY();
    z = node1.getPosZ() - node2.getPosZ();
    length = Math.sqrt(x * x + y * y + z * z);
    r1 = node1.getRadius();
    r2 = node2.getRadius();

    // calculate lateral surface of frustum cone
    lateralSurface = (r1+r2)*Math.PI*Math.sqrt((r1-r2)*(r1-r2)+length*length);
    return lateralSurface;
}

我希望有人能帮助我:)