From fd640cdda23b9cdff09bb22ced64140421bb9f02 Mon Sep 17 00:00:00 2001
From: Thomas Toulorge <thomas.toulorge@mines-paristech.fr>
Date: Thu, 2 Oct 2014 14:01:28 +0000
Subject: [PATCH] Fixed mapFromIdealElement for tets in JacobianBasis
---
Numeric/JacobianBasis.cpp | 4 +---
1 file changed, 1 insertion(+), 3 deletions(-)
diff --git a/Numeric/JacobianBasis.cpp b/Numeric/JacobianBasis.cpp
index 828b8973c9..a35231fd4a 100644
--- a/Numeric/JacobianBasis.cpp
+++ b/Numeric/JacobianBasis.cpp
@@ -376,9 +376,7 @@ void GradientBasis::mapFromIdealElement(fullMatrix<double> *dxyzdX,
case TYPE_TET:
{
- static const double cTet[3] = {-3./2/std::sqrt(6),
- -1./2/std::sqrt(2),
- std::sqrt(1.5)};
+ static const double cTet[3] = {-1.5/sqrt(6), -0.75/sqrt(2), 3./sqrt(6)};
dxyzdZ->scale(cTet[2]);
dxyzdZ->axpy(*dxyzdX, cTet[0]);
dxyzdZ->axpy(*dxyzdY, cTet[1]);
--
GitLab