be more friendly for floating point rounding issue (#7298)

jangorecki · web-flow · commit ed2f36f68975 · 2025-09-09T11:00:59.000+02:00
diff --git a/src/froll.c b/src/froll.c
@@ -185,8 +185,8 @@ void frollmeanFast(double *x, uint64_t nx, ans_t *ans, int k, double fill, bool
     SUM_WINDOW_STEP_FRONT                                       // i==k-1
     MEAN_WINDOW_STEP_VALUE
     for (uint64_t i=k; i<nx; i++) {                             // loop over obs, complete window, all remaining after partial window
-      SUM_WINDOW_STEP_FRONT
       SUM_WINDOW_STEP_BACK
+      SUM_WINDOW_STEP_FRONT
       MEAN_WINDOW_STEP_VALUE
     }
   }
@@ -375,8 +375,8 @@ if (nc == 0) {                                                 \
     SUM_WINDOW_STEP_FRONT                                       // i==k-1
     SUM_WINDOW_STEP_VALUE
     for (uint64_t i=k; i<nx; i++) {                             // loop over obs, complete window, all remaining after partial window
-      SUM_WINDOW_STEP_FRONT
       SUM_WINDOW_STEP_BACK
+      SUM_WINDOW_STEP_FRONT
       SUM_WINDOW_STEP_VALUE
     }
   }

Original file line number	Diff line number	Diff line change
`@@ -185,8 +185,8 @@ void frollmeanFast(double x, uint64_t nx, ans_t ans, int k, double fill, bool`
`185`	`185`	`SUM_WINDOW_STEP_FRONT // i==k-1`
`186`	`186`	`MEAN_WINDOW_STEP_VALUE`
`187`	`187`	`for (uint64_t i=k; i<nx; i++) { // loop over obs, complete window, all remaining after partial window`
`188`		`- SUM_WINDOW_STEP_FRONT`
`189`	`188`	`SUM_WINDOW_STEP_BACK`
	`189`	`+ SUM_WINDOW_STEP_FRONT`
`190`	`190`	`MEAN_WINDOW_STEP_VALUE`
`191`	`191`	`}`
`192`	`192`	`}`
`@@ -375,8 +375,8 @@ if (nc == 0) { \`
`375`	`375`	`SUM_WINDOW_STEP_FRONT // i==k-1`
`376`	`376`	`SUM_WINDOW_STEP_VALUE`
`377`	`377`	`for (uint64_t i=k; i<nx; i++) { // loop over obs, complete window, all remaining after partial window`
`378`		`- SUM_WINDOW_STEP_FRONT`
`379`	`378`	`SUM_WINDOW_STEP_BACK`
	`379`	`+ SUM_WINDOW_STEP_FRONT`
`380`	`380`	`SUM_WINDOW_STEP_VALUE`
`381`	`381`	`}`
`382`	`382`	`}`