From d68684b70ae65cb3a2a5f707a2e2b9e8d4fe4fb3 Mon Sep 17 00:00:00 2001
From: Amaury Johnan <amjohnen@gmail.com>
Date: Wed, 2 Oct 2013 13:50:55 +0000
Subject: [PATCH] fix bug tet p3 serendip
---
Numeric/pointsGenerators.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/Numeric/pointsGenerators.cpp b/Numeric/pointsGenerators.cpp
index 324ace50ca..77f97d94e1 100644
--- a/Numeric/pointsGenerators.cpp
+++ b/Numeric/pointsGenerators.cpp
@@ -286,7 +286,7 @@ fullMatrix<double> gmshGenerateMonomialsTetrahedron(int order, bool serendip)
}
}
- if (order > 2) {
+ if (!serendip && order > 2) {
fullMatrix<double> dudv = gmshGenerateMonomialsTriangle(order - 3);
dudv.add(1);
@@ -311,7 +311,7 @@ fullMatrix<double> gmshGenerateMonomialsTetrahedron(int order, bool serendip)
}
}
- if (!serendip && order > 3) {
+ if (order > 3) {
fullMatrix<double> inner = gmshGenerateMonomialsTetrahedron(order - 4);
inner.add(1);
monomials.copy(inner, 0, nbMonomials - index, 0, 3, index, 0);
--
GitLab