From 915509778e8b1a9f1ececbe2d8a0ce33c9d9f14d Mon Sep 17 00:00:00 2001
From: Nicolas Marsic <nicolas.marsic@gmail.com>
Date: Tue, 21 Jan 2014 17:00:56 +0000
Subject: [PATCH] Now TetNodeBasis is coherent with TriNodeBasis -- 3D Free
space case seems ok now \o/
---
FunctionSpace/TetNodeBasis.cpp | 11 ++++-------
1 file changed, 4 insertions(+), 7 deletions(-)
diff --git a/FunctionSpace/TetNodeBasis.cpp b/FunctionSpace/TetNodeBasis.cpp
index bd66ecdff0..71fa5d7923 100644
--- a/FunctionSpace/TetNodeBasis.cpp
+++ b/FunctionSpace/TetNodeBasis.cpp
@@ -79,8 +79,8 @@ TetNodeBasis::TetNodeBasis(size_t order){
for(size_t l = 1; l < order; l++){
basis[s][i] =
new Polynomial(intLegendre[l].compose
- (lagrange[edgeIdx[s][e][0]] -
- lagrange[edgeIdx[s][e][1]]
+ (lagrange[edgeIdx[s][e][1]] -
+ lagrange[edgeIdx[s][e][0]]
,
lagrange[edgeIdx[s][e][0]] +
lagrange[edgeIdx[s][e][1]]));
@@ -89,9 +89,6 @@ TetNodeBasis::TetNodeBasis(size_t order){
}
}
- // Face Based //
- // TO CHECK: Are Triangles face matching tets ?
-
for(size_t s = 0; s < nRefSpace; s++){
size_t i = nVertex + nEdge;
@@ -105,8 +102,8 @@ TetNodeBasis::TetNodeBasis(size_t order){
basis[s][i] =
new Polynomial(intLegendre[l1].compose
- (lagrange[faceIdx[s][f][0]] -
- lagrange[faceIdx[s][f][1]]
+ (lagrange[faceIdx[s][f][1]] -
+ lagrange[faceIdx[s][f][0]]
,
lagrange[faceIdx[s][f][0]] +
lagrange[faceIdx[s][f][1]])
--
GitLab