感觉 double 精度不够用啊

[tool2d]

这几天程序有一个 bug ，查出来是计算一个三角形面积，理论上应该是负数，但是函数算出来是正数，百思不得其解。

后来把计算面积公式里的 double 换成了 doubledouble 类，把精度加倍，结果就正确了。

寻思坐标也不是什么很大的值，怎么 double 就当场挂掉了呢？

=========
测试代码如下:

double buf[6];
buf[0*2+0] = 830776.16442871094; buf[0*2+1] = 430600.56744384766;
buf[1*2+0] = 830776.37707519531; buf[1*2+1] = 430600.68286132813;
buf[2*2+0] = 830776.16430664063; buf[2*2+1] = 430600.56732177734;

int numvertex = 3;

double area = 0.0;
for (int p=numvertex-1,q=0; q<numvertex; p=q++)
{
double a1 = buf[p*2+0] * buf[q*2+1];
double a2 = buf[q*2+0] * buf[p*2+1];
area += a1 - a2;
}

最终 area 结果
如果是 double ，那么值是 0.0000610 （错误）
如果是 doubledouble ，那么值是-0.0000118 （正确）

看起来是计算误差问题，但结果一个正，另一个负，足以把程序逻辑搞得天翻地覆。

[zmxnv123]

为什么三角形的面积可以是负数

[MeteorCat]

用 bigdecimal 相关库绝对没错，其他技巧都是假的

[realpg]

不是科班出身就容易踩这种坑
能不用 float 就不用 float
必须得用 float ，每次任意运算之前都得先问问自己精度的问题