PARMG cleanup, write out deformationfunction to a file, Continue the trustregionsolver with a smaller trustregion radius in case of an exception

• Clean up PARMG compiler-flags
• Write out the deformation to a file; this can then be read by dune-gfe
• Continue the trustregionsolver with a smaller trustregion radius in case of an exception @oliver.sander_at_tu-dresden.de: Can you have a look?

src/finite-strain-elasticity.cc

@@ -42,7 +42,9 @@
 #include <dune/elasticity/materials/neumannenergy.hh>
 #include <dune/elasticity/materials/sumenergy.hh>
 
+#if HAVE_DUNE_VTK
 #include <dune/vtk/writers/vtkunstructuredgridwriter.hh>
+#endif
 
 #include <iostream>
 #include <fstream>
@@ -418,21 +420,25 @@ int main (int argc, char *argv[]) try
     vtkWriter.addVertexData(displacementFunction, VTK::FieldInfo("displacement", VTK::FieldInfo::Type::vector, dim));
     vtkWriter.write(resultPath + "finite-strain_homotopy_" + std::to_string(i+1));
   }
-
-  displacement = x;
-  displacement -= identity;
-  auto finalDisplacementFunction = Dune::Functions::makeDiscreteGlobalBasisFunction<FieldVector<double,dim>>(feBasis, displacement);
-  //Write out the grid using an unstructured GridWriter
-  VtkUnstructuredGridWriter<GridView> vtkUnstructuredGridWriter(gridView, Vtk::ASCII);
-  vtkUnstructuredGridWriter.addPointData(finalDisplacementFunction, "displacement");
-  vtkUnstructuredGridWriter.write(resultPath + "finite-strain_homotopy_original");
-  std::ofstream file;
-  file.open("deformation");
-  for (auto& d : displacement){
-    file << d << "\n";
+#if HAVE_DUNE_VTK
+  if (grid->leafGridView().comm().size() == 1) {
+    displacement = x;
+    displacement -= identity;
+    auto finalDisplacementFunction = Dune::Functions::makeDiscreteGlobalBasisFunction<FieldVector<double,dim>>(feBasis, displacement);
+    //Write out the grid using an unstructured GridWriter
+    VtkUnstructuredGridWriter<GridView> vtkUnstructuredGridWriter(gridView, Vtk::ASCII);
+    vtkUnstructuredGridWriter.addPointData(finalDisplacementFunction, "displacement");
+    vtkUnstructuredGridWriter.write(resultPath + "finite-strain_homotopy_original");
+    std::ofstream file;
+    file.open("deformation");
+    for (auto& d : displacement){
+      file << d << "\n";
+    }
+    file.close();
+  } else if (mpiHelper.rank() == 0) {
+    std::cout << "Not writing out the deformation to a file, this currently only works for sequential calculations" << std::endl;
   }
-  file.close();
-
+#endif
   if (mpiHelper.rank()==0)
     std::cout << "Complete duration: " << homotopyTimer.elapsed() << " sec." << std::endl;