Curve Fitting: correct initial handle sign

In practice the initial values are almost never used.

diff --git a/extern/curve_fit_nd/intern/curve_fit_cubic_refit.c b/extern/curve_fit_nd/intern/curve_fit_cubic_refit.c
index b51535bab106798c2cb21e2e8af4f7fb3d935039..d22d042bff50a546c6e35fb4088f6235d8565e87 100644 (file)
--- a/extern/curve_fit_nd/intern/curve_fit_cubic_refit.c
+++ b/extern/curve_fit_nd/intern/curve_fit_cubic_refit.c
@@ -1048,7 +1048,6 @@ int curve_fit_cubic_to_points_refit_db(
 {
        const uint knots_len = points_len;
        struct Knot *knots = malloc(sizeof(Knot) * knots_len);
-       knots[0].next = NULL;
 
 #ifndef USE_CORNER_DETECT
        (void)r_corner_index_array;
@@ -1085,7 +1084,6 @@ int curve_fit_cubic_to_points_refit_db(
 
                        knots[i].heap_node = NULL;
                        knots[i].index = i;
-                       knots[i].index = i;
                        knots[i].can_remove = true;
                        knots[i].is_removed = false;
                        knots[i].is_corner = false;
@@ -1155,8 +1153,8 @@ int curve_fit_cubic_to_points_refit_db(
                        add_vn_vnvn(k->tan[0], tan_prev, tan_next, dims);
                        normalize_vn(k->tan[0], dims);
                        copy_vnvn(k->tan[1], k->tan[0], dims);
-                       k->handles[0] = len_prev / 3;
-                       k->handles[1] = len_next / 3;
+                       k->handles[0] = len_prev /  3;
+                       k->handles[1] = len_next / -3;
                }
 #else
                if (knots_len < 2) {
@@ -1185,8 +1183,8 @@ int curve_fit_cubic_to_points_refit_db(
                                add_vn_vnvn(k->tan[0], tan_prev, tan_next, dims);
                                normalize_vn(k->tan[0], dims);
                                copy_vnvn(k->tan[1], k->tan[0], dims);
-                               k->handles[0] = len_prev / 3;
-                               k->handles[1] = len_next / 3;
+                               k->handles[0] = len_prev /  3;
+                               k->handles[1] = len_next / -3;
 
                                copy_vnvn(tan_prev, tan_next, dims);
                                len_prev = len_next;
@@ -1201,8 +1199,8 @@ int curve_fit_cubic_to_points_refit_db(
                                tan_prev, &points[0 * dims], &points[1 * dims], dims);
                        copy_vnvn(knots[0].tan[0], tan_prev, dims);
                        copy_vnvn(knots[0].tan[1], tan_prev, dims);
-                       knots[0].handles[0] = len_prev / 3;
-                       knots[0].handles[1] = len_prev / 3;
+                       knots[0].handles[0] = len_prev /  3;
+                       knots[0].handles[1] = len_prev / -3;
 
                        for (uint i_curr = 1, i_next = 2; i_next < knots_len; i_curr = i_next++) {
                                struct Knot *k = &knots[i_curr];
@@ -1215,8 +1213,8 @@ int curve_fit_cubic_to_points_refit_db(
                                add_vn_vnvn(k->tan[0], tan_prev, tan_next, dims);
                                normalize_vn(k->tan[0], dims);
                                copy_vnvn(k->tan[1], k->tan[0], dims);
-                               k->handles[0] = len_prev / 3;
-                               k->handles[1] = len_next / 3;
+                               k->handles[0] = len_prev /  3;
+                               k->handles[1] = len_next / -3;
 
                                copy_vnvn(tan_prev, tan_next, dims);
                                len_prev = len_next;
@@ -1224,8 +1222,8 @@ int curve_fit_cubic_to_points_refit_db(
                        copy_vnvn(knots[knots_len - 1].tan[0], tan_next, dims);
                        copy_vnvn(knots[knots_len - 1].tan[1], tan_next, dims);
 
-                       knots[knots_len - 1].handles[0] = len_next / 3;
-                       knots[knots_len - 1].handles[1] = len_next / 3;
+                       knots[knots_len - 1].handles[0] = len_next /  3;
+                       knots[knots_len - 1].handles[1] = len_next / -3;
                }
 #endif
        }