From 3598f86bbe0d04e56145ff3be66efdd5a5c3fb0e Mon Sep 17 00:00:00 2001
From: Thomas Toulorge <thomas.toulorge@mines-paristech.fr>
Date: Mon, 9 Jul 2012 15:52:13 +0000
Subject: [PATCH] Fixed bug in 3D Hessian computation
---
Mesh/meshMetric.cpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/Mesh/meshMetric.cpp b/Mesh/meshMetric.cpp
index 5997b9e47a..12e9f7d006 100644
--- a/Mesh/meshMetric.cpp
+++ b/Mesh/meshMetric.cpp
@@ -295,7 +295,7 @@ void meshMetric::computeHessian()
d2udxy = coeffs(1); d2udxz = coeffs(2); d2udyz = coeffs(4);
dudx = d2udx2*x+d2udxy*y+d2udxz*z+coeffs(6);
dudy = d2udxy*x+d2udy2*y+d2udyz*z+coeffs(7);
- dudz = d2udxz*z+d2udyz*z+d2udz2*z+coeffs(8);
+ dudz = d2udxz*x+d2udyz*y+d2udz2*z+coeffs(8);
}
double duNorm = sqrt(dudx*dudx+dudy*dudy+dudz*dudz);
if (duNorm == 0. || _technique == meshMetric::HESSIAN) duNorm = 1.;
--
GitLab