From ca6625ba0966cb22c60050b9ce158ef0db5b02df Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Thu, 31 May 2018 09:20:11 +0200 Subject: diff: add the missing power for higher order derivatives --- diff.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/diff.py b/diff.py index bd9a135..553ee69 100644 --- a/diff.py +++ b/diff.py @@ -113,7 +113,7 @@ class FiniteDifferenceUniform(object): for s, c in zip(slices_src, coeff): ret[tuple(slicelist_ret)] += c * arr[tuple(s)] - ret[tuple(slicelist_ret)] /= (div * dx) + ret[tuple(slicelist_ret)] /= (div * (dx ** self.diff_order)) if bound_lower: if self.acc_order > 2: -- cgit v1.2.3