From d0391c45321fd2ed961a7180fcdabb310a084bda Mon Sep 17 00:00:00 2001
From: Amaury Johnan <amjohnen@gmail.com>
Date: Wed, 17 Aug 2016 15:49:11 +0000
Subject: [PATCH] fix bug when optimizing 3D meshes (thanks to Wendy
Merks-Swolfs!)
---
contrib/HighOrderMeshOptimizer/OptHomMesh.cpp | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/contrib/HighOrderMeshOptimizer/OptHomMesh.cpp b/contrib/HighOrderMeshOptimizer/OptHomMesh.cpp
index afaa2d9dd4..1121b97601 100644
--- a/contrib/HighOrderMeshOptimizer/OptHomMesh.cpp
+++ b/contrib/HighOrderMeshOptimizer/OptHomMesh.cpp
@@ -112,6 +112,10 @@ Mesh::Mesh(const std::map<MElement*,GEntity*> &element2entity,
double dumJac[3][3];
for (int iEl = 0; iEl < nEl(); iEl++)
_invStraightJac[iEl] = 1. / fabs(_el[iEl]->getPrimaryJacobian(0.,0.,0.,dumJac));
+
+ // initialize _scaledNormEl with empty fullMatrix for use in function
+ // Mesh::scaledJacAndGradients(..)
+ _scaledNormEl.resize(nEl());
}
}
--
GitLab