From 7e218e56945b8369109978009e062665e62f2d5c Mon Sep 17 00:00:00 2001
From: Jonathan Lambrechts <jonathan.lambrechts@uclouvain.be>
Date: Fri, 22 Oct 2010 08:42:33 +0000
Subject: [PATCH] fake 0D gauss rule
---
Numeric/Gauss.cpp | 5 +++++
Numeric/polynomialBasis.cpp | 2 +-
2 files changed, 6 insertions(+), 1 deletion(-)
diff --git a/Numeric/Gauss.cpp b/Numeric/Gauss.cpp
index 1f6e54b962..50fbb050db 100644
--- a/Numeric/Gauss.cpp
+++ b/Numeric/Gauss.cpp
@@ -61,6 +61,11 @@ void gaussIntegration::get(int elementType, int order, fullMatrix<double> &pts,
case TYPE_TET : pts2fullMatrix(getNGQTetPts(order), getGQTetPts(order), pts, weights); break;
case TYPE_HEX : pts2fullMatrix(getNGQHPts(order), getGQHPts(order), pts, weights); break;
case TYPE_PRI : pts2fullMatrix(getNGQPriPts(order), getGQPriPts(order), pts, weights); break;
+ case TYPE_PNT :
+ weights.resize(1,1);
+ weights(0,0) = 1.;
+ pts.resize(1,3);
+ break;
default : Msg::Error("No integration rules defined for type %i", elementType);
}
}
diff --git a/Numeric/polynomialBasis.cpp b/Numeric/polynomialBasis.cpp
index 97d67876e4..79d31aca83 100644
--- a/Numeric/polynomialBasis.cpp
+++ b/Numeric/polynomialBasis.cpp
@@ -1032,7 +1032,7 @@ const polynomialBasis *polynomialBases::find(int tag)
F.numFaces = 1;
F.monomials = generate1DMonomials(0);
F.points = generate1DPoints(0);
- F.parentType = MSH_PNT;
+ F.parentType = TYPE_PNT;
break;
case MSH_LIN_2 :
F.numFaces = 2;
--
GitLab