[cig-commits] commit 1899 by bangerth to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 1899

Only call add_periodicity() if there is actually something to be added.

U   trunk/aspect/source/geometry_model/box.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=1899&peg=1899

Diff:
Modified: trunk/aspect/source/geometry_model/box.cc
===================================================================
--- trunk/aspect/source/geometry_model/box.cc	2013-09-17 18:43:09 UTC (rev 1898)
+++ trunk/aspect/source/geometry_model/box.cc	2013-09-17 18:53:25 UTC (rev 1899)
@@ -48,11 +48,12 @@
       std::vector<std_cxx1x::tuple< typename parallel::distributed::Triangulation<dim>::cell_iterator, unsigned int,
                                     typename parallel::distributed::Triangulation<dim>::cell_iterator, unsigned int> >
                                    periodicity_vector;
-      for( unsigned int i=0; i<dim; ++i)
+      for (unsigned int i=0; i<dim; ++i)
         if (periodic[i])
           GridTools::identify_periodic_face_pairs(coarse_grid, 2*i, 2*i+1, i, periodicity_vector);
 
-      coarse_grid.add_periodicity(periodicity_vector);
+      if (periodicity_vector.size() > 0)
+        coarse_grid.add_periodicity (periodicity_vector);
 #else
       for( unsigned int i=0; i<dim; ++i)
         AssertThrow(!periodic[i],