diff options
| -rw-r--r-- | src/interpolate.c | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/interpolate.c b/src/interpolate.c index 4457f70..d990cf5 100644 --- a/src/interpolate.c +++ b/src/interpolate.c @@ -231,14 +231,14 @@ CCTK_REAL interpolate_eno(CCTK_REAL x0, CCTK_REAL dx, CCTK_REAL xr; /* Find seed index */ - while (x > x0 + dx * ((float)seed+0.5)) seed++; + while (x > x0 + dx * ((CCTK_REAL)seed+0.5)) seed++; if (seed!=2) { /* Not enough stencil, only perform linear interpolation */ seed = 0; - while (x > x0 + dx * ((float)(seed+1))) seed++; + while (x > x0 + dx * ((CCTK_REAL)(seed+1))) seed++; if (seed==4) seed=3; - return y[seed] + (y[seed+1]-y[seed])/dx * (x-x0-(float)seed*dx); + return y[seed] + (y[seed+1]-y[seed])/dx * (x-x0-(CCTK_REAL)seed*dx); } /* Calculate the undivided differences */ |