From efc1be484234ee621763ae590488ee62ce61bbec Mon Sep 17 00:00:00 2001
From: Brian Helenbrook <bhelenbr@clarkson.edu>
Date: Wed, 19 Nov 2008 08:44:39 +0000
Subject: [PATCH] Fixed bug in getfacevertices for prism case.
---
Mesh/HighOrder.cpp | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/Mesh/HighOrder.cpp b/Mesh/HighOrder.cpp
index 6608363e63..06dbd39b1b 100644
--- a/Mesh/HighOrder.cpp
+++ b/Mesh/HighOrder.cpp
@@ -17,6 +17,7 @@
#include "GFaceCompound.h"
#include "GmshMatrix.h"
#include "FunctionSpace.h"
+#include <iostream>
#define SQU(a) ((a)*(a))
@@ -865,6 +866,7 @@ void getFaceVertices(GRegion *gr, MElement *ele,
Double_Matrix points;
+
int start = 0;
switch (nPts){
@@ -911,7 +913,7 @@ void getFaceVertices(GRegion *gr, MElement *ele,
else{
std::vector<MVertex*>& vtcs = faceVertices[face];
- if(face.getNumVertices() == 3){ // triangles
+ if(face.getNumVertices() == 3 && nPts > 1){ // triangles
// construct incomplete element to take into account curved edges on surface boundaries
--
GitLab