minor changes.

examples/Makefile.am

@@ -3,7 +3,7 @@
 GRIDTYPE = ALUGRID_CUBE
 #GRIDTYPE = ALUGRID_SIMPLEX
 #GRIDTYPE = ALUGRID_CONFORM
-GRIDDIM = 3
+GRIDDIM = 2
 
 SUBDIRS = callback loadbalance dgf backuprestore communication testEfficiency
 

examples/dgf/unitcube2d.dgf

@@ -3,7 +3,7 @@ DGF
 Interval
 0  0                                    % first corner
 1  1                                    % second corner
-8  8                                    % 4 cells in each direction
+16 16                                    % 4 cells in each direction
 #
 
 BoundaryDomain

examples/main.cc

@@ -24,6 +24,18 @@
 #include "diagnostics.hh"
 #include "paralleldgf.hh"
 
+  template <class Grid>
+  inline void printSpaceFillingCurve( const Grid& grid, std::ostream& out )
+  {
+    const auto& gridView = grid.leafGridView();
+    for( auto it = gridView.template begin<0>(), end = gridView.template end<0>();
+         it != end; ++it )
+    {
+      out << (*it).geometry().center() << std::endl;
+    }
+  }
+
+
 // method
 // ------
 void method ( int problem, int startLvl, int maxLvl,
@@ -51,6 +63,12 @@ void method ( int problem, int startLvl, int maxLvl,
 
   Grid &grid = *gridPtr;
 
+  if( grid.comm().rank() == 0 )
+  {
+    std::ofstream file( "curve.gnu" );
+    printSpaceFillingCurve( grid, file );
+  }
+
 #ifndef BALL
   if ( grid.comm().size() > 1 &&
        (grid.overlapSize(0)==0 && grid.ghostSize(0)==0)