From 741eb7398ae0906c9e1d61272649ed32b7b35106 Mon Sep 17 00:00:00 2001
From: Stephane Del Pino <stephane.delpino44@gmail.com>
Date: Fri, 29 May 2020 16:00:16 +0200
Subject: [PATCH] Fix Cjr sign in 2D
This was buggy since used gmsh meshes were incorrectly oriented ...
Should check cells volume
---
src/mesh/MeshData.hpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/mesh/MeshData.hpp b/src/mesh/MeshData.hpp
index 1647abc0b..5c4a6068a 100644
--- a/src/mesh/MeshData.hpp
+++ b/src/mesh/MeshData.hpp
@@ -106,7 +106,7 @@ class MeshData
int Rp1 = (R + 1) % cell_nodes.size();
int Rm1 = (R + cell_nodes.size() - 1) % cell_nodes.size();
Rd half_xrp_xrm = 0.5 * (xr[cell_nodes[Rp1]] - xr[cell_nodes[Rm1]]);
- Cjr(j, R) = Rd{-half_xrp_xrm[1], half_xrp_xrm[0]};
+ Cjr(j, R) = Rd{half_xrp_xrm[1], -half_xrp_xrm[0]};
}
});
m_Cjr = Cjr;
--
GitLab