From 747f1851d5ab36f90e76c24e559bafbe3c259f48 Mon Sep 17 00:00:00 2001
From: Martin Bauer <martin.bauer@fau.de>
Date: Wed, 26 Jul 2017 13:08:20 +0200
Subject: [PATCH] Bugfix in finite-differences: dx forgotten in divergence
---
finitedifferences.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/finitedifferences.py b/finitedifferences.py
index 171260cdc..423b62064 100644
--- a/finitedifferences.py
+++ b/finitedifferences.py
@@ -115,14 +115,14 @@ def discretizeDivergence(vectorTerm, symbolsToFieldDict, dx):
>>> gradX = grad(x, dim=3)
>>> f = Field.createGeneric('f', spatialDimensions=3)
>>> sp.simplify(discretizeDivergence(gradX, {x : f}, dx))
- (f_B - 6*f_C + f_E + f_N + f_S + f_T + f_W)/dx
+ (f_B - 6*f_C + f_E + f_N + f_S + f_T + f_W)/dx**2
"""
dim = len(vectorTerm)
result = 0
for d in range(dim):
for offset in [-1, 1]:
result += offset * discretizeStaggered(vectorTerm[d], symbolsToFieldDict, d, offset, dx, dim)
- return result
+ return result / dx
def __upDownOffsets(d, dim):
--
GitLab