]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fixed issue with periodicity.
authorLuca Heltai <luca.heltai@sissa.it>
Sat, 26 Mar 2016 09:47:59 +0000 (10:47 +0100)
committerLuca Heltai <luca.heltai@sissa.it>
Wed, 6 Apr 2016 11:00:14 +0000 (13:00 +0200)
source/grid/manifold.cc

index 4fc83bf99dbce8e0adcace97ea6bea90e0543f90..361057fbf9f87ba520b7bc94d50daa028cc7f850 100644 (file)
@@ -423,7 +423,8 @@ get_new_point (const Quadrature<spacedim> &quad) const
     if (periodicity[d] > 0)
       for (unsigned int i=0; i<surrounding_points.size(); ++i)
         {
-          minP[d] = std::min(minP[d], surrounding_points[i][d]);
+          if (weights[i] > tolerance)
+            minP[d] = std::min(minP[d], surrounding_points[i][d]);
           Assert( (surrounding_points[i][d] < periodicity[d]+tolerance*periodicity.norm()) ||
                   (surrounding_points[i][d] >= -tolerance*periodicity.norm()),
                   ExcPeriodicBox(d, surrounding_points[i], periodicity, tolerance*periodicity.norm()));

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.