aboutsummaryrefslogtreecommitdiff
path: root/src/PJ_horner.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/PJ_horner.c')
-rw-r--r--src/PJ_horner.c8
1 files changed, 6 insertions, 2 deletions
diff --git a/src/PJ_horner.c b/src/PJ_horner.c
index 09554a7f..a6a26e52 100644
--- a/src/PJ_horner.c
+++ b/src/PJ_horner.c
@@ -270,8 +270,12 @@ summing the tiny high order elements first.
double u, v, N, E;
/* Double Horner's scheme: N = n*Cy*e -> yout, E = e*Cx*n -> xout */
- for (N = *--tcy, E = *--tcx; r > 0; r--) {
- for (c = g, u = *--tcy, v = *--tcx; c >= r; c--) {
+ N = *--tcy;
+ E = *--tcx;
+ for (; r > 0; r--) {
+ u = *--tcy;
+ v = *--tcx;
+ for (c = g; c >= r; c--) {
u = n*u + *--tcy;
v = e*v + *--tcx;
}