From 469e52f6ad1bce83997a1ae0270fbfb6bd0899f2 Mon Sep 17 00:00:00 2001
From: Thomas Toulorge <thomas.toulorge@mines-paristech.fr>
Date: Tue, 18 Mar 2014 09:59:54 +0000
Subject: [PATCH] Fixed OptHOM for HO elements with negative primary Jacobian
---
contrib/HighOrderMeshOptimizer/OptHomMesh.cpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/contrib/HighOrderMeshOptimizer/OptHomMesh.cpp b/contrib/HighOrderMeshOptimizer/OptHomMesh.cpp
index e8cdf9de65..4322850bbe 100644
--- a/contrib/HighOrderMeshOptimizer/OptHomMesh.cpp
+++ b/contrib/HighOrderMeshOptimizer/OptHomMesh.cpp
@@ -108,7 +108,7 @@ Mesh::Mesh(const std::map<MElement*,GEntity*> &element2entity,
_invStraightJac.resize(nEl(),1.);
double dumJac[3][3];
for (int iEl = 0; iEl < nEl(); iEl++)
- _invStraightJac[iEl] = 1. / _el[iEl]->getPrimaryJacobian(0.,0.,0.,dumJac);
+ _invStraightJac[iEl] = 1. / fabs(_el[iEl]->getPrimaryJacobian(0.,0.,0.,dumJac));
}
}
--
GitLab