From 4aeb7380203f6b2c53feff9dd7e09624f1a21f38 Mon Sep 17 00:00:00 2001
From: Jonathan Lambrechts <jonathan.lambrechts@uclouvain.be>
Date: Thu, 5 May 2011 13:11:34 +0000
Subject: [PATCH] fix prisms shape functions
---
Numeric/polynomialBasis.cpp | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/Numeric/polynomialBasis.cpp b/Numeric/polynomialBasis.cpp
index 67c16b6f9d..d77b0644ee 100644
--- a/Numeric/polynomialBasis.cpp
+++ b/Numeric/polynomialBasis.cpp
@@ -1605,14 +1605,14 @@ const polynomialBasis *polynomialBases::find(int tag)
case TYPE_PRI :
F.numFaces = 5;
F.dimension = 3;
- F.monomials = generatePascalPrism(0);
- F.points = gmshGeneratePointsPrism(0, F.serendip);
+ F.monomials = generatePascalPrism(F.order);
+ F.points = gmshGeneratePointsPrism(F.order, F.serendip);
if (F.order == 0) {
generateClosureOrder0(F.closures,48);
generateClosureOrder0(F.fullClosures,48);
F.closureRef.resize(48, 0);
} else {
- generateFaceClosurePrism(F.closures, 1);
+ generateFaceClosurePrism(F.closures, F.order);
generateFaceClosurePrismFull(F.fullClosures, F.closureRef, F.order);
}
break;
--
GitLab