Skip to content
Snippets Groups Projects
Commit 8f2e1c03 authored by Thomas Toulorge's avatar Thomas Toulorge
Browse files

Undid change committed in r19612 that was erroneous

parent fd640cdd
No related branches found
No related tags found
No related merge requests found
...@@ -376,7 +376,9 @@ void GradientBasis::mapFromIdealElement(fullMatrix<double> *dxyzdX, ...@@ -376,7 +376,9 @@ void GradientBasis::mapFromIdealElement(fullMatrix<double> *dxyzdX,
case TYPE_TET: case TYPE_TET:
{ {
static const double cTet[3] = {-1.5/sqrt(6), -0.75/sqrt(2), 3./sqrt(6)}; static const double cTet[3] = {-3./2/std::sqrt(6),
-1./2/std::sqrt(2),
std::sqrt(1.5)};
dxyzdZ->scale(cTet[2]); dxyzdZ->scale(cTet[2]);
dxyzdZ->axpy(*dxyzdX, cTet[0]); dxyzdZ->axpy(*dxyzdX, cTet[0]);
dxyzdZ->axpy(*dxyzdY, cTet[1]); dxyzdZ->axpy(*dxyzdY, cTet[1]);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment