From 694786c1a740a3a9a0144eef8c80e202b2835843 Mon Sep 17 00:00:00 2001 From: guido Date: Thu, 22 Jan 2004 11:33:48 +0000 Subject: [PATCH] orientation git-svn-id: https://svn.dealii.org/trunk@8331 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/base/source/polynomials_bdm.cc | 67 ++++++++++++++------------ 1 file changed, 37 insertions(+), 30 deletions(-) diff --git a/deal.II/base/source/polynomials_bdm.cc b/deal.II/base/source/polynomials_bdm.cc index c4771dc977..2dc82866fb 100644 --- a/deal.II/base/source/polynomials_bdm.cc +++ b/deal.II/base/source/polynomials_bdm.cc @@ -140,42 +140,49 @@ PolynomialsBDM::compute_node_matrix (Table<2,double>& A) const values.resize(n()); for (unsigned int face=0;face<2*dim;++face) - for (unsigned int k=0;k p; - switch (face) - { - case 2: - p(1) = 1.; - case 0: - p(0) = x; - break; - case 1: - p(0) = 1.; - case 3: - p(1) = x; - break; - } + { + double orientation = 1.; + if ((face==0) || (face==3)) + orientation = -1.; + + for (unsigned int k=0;k p; + switch (face) + { + case 2: + p(1) = 1.; + case 0: + p(0) = x; + break; + case 1: + p(0) = 1.; + case 3: + p(1) = x; + break; + } // std::cerr << p // << '\t' << moment_weight[0].value(x) // << '\t' << moment_weight[1].value(x) // ; - - compute (p, values, grads, grad_grads); - - for (unsigned int i=0;i