HostGrid can provide geomTypes in container not std::vector

Summary

IndexSet is allowed to provide its GeometryTypes list as something else than std::vector<GeometryType>. Therefore the IndexSet::Types alias is available. Most grids just return a std::vector, but, e.g., AlbertaGrid returns a std::array instead. Almost all grids could returns something cheaper than std::vector since there is only a single GeometryType. This MR, fixes the handling of this flexibility in the ParallelGrid wrapper by copying the GeometryTypes into the new storage.

dune/grid/parallelgrid/grid.hh

@@ -315,9 +315,10 @@ namespace Dune
       if( hostGrid.comm().size() > 1 )
         DUNE_THROW( InvalidStateException, "ParallelGrid can only be used with serial host grids." );
 
-      for( int i = 0; i <= HostGrid::dimension; ++i )
-        geomTypes_.push_back( hostGrid.levelGridView( 0 ).indexSet().types( i ) );
-
+      for( int i = 0; i <= HostGrid::dimension; ++i ) {
+        auto types = hostGrid.levelGridView( 0 ).indexSet().types( i );
+        geomTypes_.emplace_back(std::begin(types), std::end(types));
+      }
       sizeCache_.reset();
 
       for (std::size_t i = 0; i < hostGrid.numBoundarySegments(); ++i)
@@ -347,8 +348,10 @@ namespace Dune
       if( hostGrid->comm().size() > 1 )
         DUNE_THROW( InvalidStateException, "ParallelGrid can only be used with serial host grids." );
 
-      for( int i = 0; i <= HostGrid::dimension; ++i )
-        geomTypes_.push_back( hostGrid->levelGridView( 0 ).indexSet().types( i ) );
+      for( int i = 0; i <= HostGrid::dimension; ++i ) {
+        auto types = hostGrid->levelGridView( 0 ).indexSet().types( i );
+        geomTypes_.emplace_back(std::begin(types), std::end(types));
+      }
 
       sizeCache_.reset();