diff options
Diffstat (limited to 'src/PJ_horner.c')
| -rw-r--r-- | src/PJ_horner.c | 8 |
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; } |
