From 85beea0a8c566bb075638843f6472b7b9885ea5e Mon Sep 17 00:00:00 2001 From: Luca Heltai Date: Tue, 12 Jul 2016 13:29:00 +0200 Subject: [PATCH] Fixed bug in mapping_manifold --- source/fe/mapping_manifold.cc | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/source/fe/mapping_manifold.cc b/source/fe/mapping_manifold.cc index e21a20c1a9..e86c10ea2b 100644 --- a/source/fe/mapping_manifold.cc +++ b/source/fe/mapping_manifold.cc @@ -457,11 +457,12 @@ namespace internal Assert(pi >=0 && pi <= 1.0, ExcInternalError("Was expecting a quadrature point " "inside the unit reference element.")); - const Point np(pi > .5 ? p-pi *ei : p+(1-pi)*ei); // In the length L, we store also the direction sign, // which is positive, if the coordinate is < .5, const double L = pi > .5 ? -pi: 1-pi; + + const Point np(p + L*ei); // Get the weights to compute the np point in real space for (unsigned int j=0; j::vertices_per_cell; ++j) -- 2.39.5