]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fixed bug in mapping_manifold
authorLuca Heltai <luca.heltai@sissa.it>
Tue, 12 Jul 2016 11:29:00 +0000 (13:29 +0200)
committerLuca Heltai <luca.heltai@sissa.it>
Thu, 4 Aug 2016 08:44:10 +0000 (10:44 +0200)
source/fe/mapping_manifold.cc

index e21a20c1a90c98ee10596a1db3cd20b7629daf99..e86c10ea2b13729a06e7a30d9a18b090aaf1bec4 100644 (file)
@@ -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<dim> 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<dim> np(p + L*ei);
 
                   // Get the weights to compute the np point in real space
                   for (unsigned int j=0; j<GeometryInfo<dim>::vertices_per_cell; ++j)

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.