From: guido Date: Thu, 22 Jan 2004 11:33:48 +0000 (+0000) Subject: orientation X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=694786c1a740a3a9a0144eef8c80e202b2835843;p=dealii-svn.git orientation git-svn-id: https://svn.dealii.org/trunk@8331 0785d39b-7218-0410-832d-ea1e28bc413d --- 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