From 3c2edd3ba4db598e46a8be584ac9908d87d7306c Mon Sep 17 00:00:00 2001 From: domfournier Date: Mon, 3 Feb 2025 10:03:16 -0800 Subject: [PATCH] Add small quantity to avoid zero division (cherry picked from commit 925f3e4735060af4e992b1e3711dbad3d20a7f70) --- simpeg/electromagnetics/static/resistivity/simulation.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/simpeg/electromagnetics/static/resistivity/simulation.py b/simpeg/electromagnetics/static/resistivity/simulation.py index 4e65a8d07a..bdd09e1666 100644 --- a/simpeg/electromagnetics/static/resistivity/simulation.py +++ b/simpeg/electromagnetics/static/resistivity/simulation.py @@ -606,7 +606,7 @@ def setBC(self): # TODO: Implement Zhang et al. (1995) r_vec = boundary_faces - source_point - r = np.linalg.norm(r_vec, axis=-1) + r = np.linalg.norm(r_vec, axis=-1) + 1e-12 r_hat = r_vec / r[:, None] r_dot_n = np.einsum("ij,ij->i", r_hat, boundary_normals)