Parallel grid wip

dune/elasticity/assemblers/feassembler.hh

@@ -71,7 +71,7 @@ getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const
 
     nb.resize(n, n);
 
-    for (const auto& element : elements(basis_.getGridView(), Dune::Partitions::interior))
+    for (const auto& element : elements(basis_.getGridView()))
     {
         const typename Basis::LocalFiniteElement& lfe = basis_.getLocalFiniteElement(element);
 
@@ -112,7 +112,7 @@ assembleGradientAndHessian(const VectorType& sol,
     gradient.resize(sol.size());
     gradient = 0;
 
-    for (const auto& element : elements(basis_.getGridView(), Dune::Partitions::interior))
+    for (const auto& element : elements(basis_.getGridView()))
     {
         const int numOfBaseFct = basis_.getLocalFiniteElement(element).localBasis().size();
 
@@ -174,7 +174,10 @@ computeEnergy(const VectorType& sol) const
 
     }
 
-    return energy;
-
+    auto&& grid = basis_.getGridView().grid();
+    double totalEnergy = grid.comm().sum(energy);
+    std::cout << "rank " << grid.comm().rank() << " gives energy of " << energy << " and total = " << totalEnergy << std::endl;
+    grid.comm().broadcast(&totalEnergy, 1, 0);
+    return totalEnergy;
 }
 #endif