迭代时浮点数不精确

Question

curves_error curves_bezier(curves_PointList* list, curves_Point* dest, curves_float t) {
    /* ... */
    if (t < 0 || t > 1)
        return curves_invalid_args;
    /* ... */
    return curves_no_error;
}

for (t=0; t <= 1.0; t += 0.1) {
    curves_error error = curves_bezier(points, point, t);
    if (error != curves_no_error)
        printf("Error with t = %f.\n", t);
    else
        printf("t = %f is ok.\n", t);
}

Answer 1

for (t=0; t <= 1.0; t += 0.1) {

您的问题是二进制浮点数不能精确地表示0.1。

最接近的32位单精度IEEE754浮点数是0.100000001490116119384765625和最接近的64位双精度浮点数0.10000000000000000000055511151231257827021181583404541015625。如果严格以32位精度执行算术，则将0.1f加十次到0的结果为

1.00000011920928955078125

如果以比float更高的精度执行中间计算，则可以得到精确的1.0或甚至更小的数字。

若要修复问题，在本例中，您可以使用

for(k = 0; k <= 10; ++k) {
    t = k*0.1;

因为10*0.1f正好是1.0。

另一种选择是在curves_bezier函数中使用小容差,

if (t > 1 && t < 1 + epsilon) {
    t = 1;
}