From a818ffd82d9fbe13ba5e74b01fff4f69870f763e Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@ulg.ac.be>
Date: Thu, 20 Nov 2008 07:01:20 +0000
Subject: [PATCH] *** empty log message ***
---
Mesh/HighOrder.cpp | 10 +++-------
1 file changed, 3 insertions(+), 7 deletions(-)
diff --git a/Mesh/HighOrder.cpp b/Mesh/HighOrder.cpp
index 67263e0ce8..5e950238b3 100644
--- a/Mesh/HighOrder.cpp
+++ b/Mesh/HighOrder.cpp
@@ -7,7 +7,6 @@
// Koen Hillewaert
//
-#include <iostream>
#include "HighOrder.h"
#include "gmshSmoothHighOrder.h"
#include "MElement.h"
@@ -389,8 +388,8 @@ static void getEdgeVertices(GRegion *gr, MElement *ele, std::vector<MVertex*> &v
static void getFaceVertices(GFace *gf, MElement *incomplete, MElement *ele,
std::vector<MVertex*> &vf, faceContainer &faceVertices,
bool linear, int nPts = 1, gmshHighOrderSmoother *displ2D = 0,
- gmshHighOrderSmoother *displ3D = 0) {
-
+ gmshHighOrderSmoother *displ3D = 0)
+{
Double_Matrix points;
int start = 0;
@@ -641,7 +640,6 @@ static void getFaceVertices(GRegion *gr, MElement *ele, std::vector<MVertex*> &v
edgeContainer &edgeVertices, bool linear, int nPts = 1)
{
Double_Matrix points;
-
int start = 0;
switch (nPts){
@@ -1134,7 +1132,6 @@ void SetOrderN(GModel *m, int order, bool linear, bool incomplete)
setHighOrder(*it, edgeVertices, faceVertices, linear, incomplete, nPts,
displ2D, displ3D);
-
// now we smooth mesh the internal vertices of the faces
// we do that model face by model face
std::vector<MElement*> bad;
@@ -1156,13 +1153,12 @@ void SetOrderN(GModel *m, int order, bool linear, bool incomplete)
displ3D->smooth(*it);
}
- printJacobians(m, "detjIni.pos");
-
if(displ2D){
delete displ2D;
delete displ3D;
}
+ printJacobians(m, "detjIni.pos");
checkHighOrderTriangles(m, bad, worst);
double t2 = Cpu();
--
GitLab