From 1f3e3fdfe841a0fe2e2122d8e6019b974f1e354c Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Thu, 9 Jul 2020 10:47:27 +0200 Subject: interp: fix source index computation --- interp.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/interp.py b/interp.py index c57f417..f0f5693 100644 --- a/interp.py +++ b/interp.py @@ -1,7 +1,7 @@ import numpy as np def interp1d(src_start, src_step, src_val, dst_coords, stencil): - idx_src = (dst_coords / src_step - src_start - (stencil / 2 - 1)).astype(np.int) + idx_src = ((dst_coords - src_start) / src_step - (stencil / 2 - 1)).astype(np.int) src_coord = np.linspace(src_start, src_start + src_step * (src_val.shape[0] - 1), src_val.shape[0], dtype = src_val.dtype) -- cgit v1.2.3