diff --git a/Solver/linearSystemPETSc.h b/Solver/linearSystemPETSc.h
index e1229bd6ec7323267f7e2631db5ade649d7b5b31..c2ed187e375b31076470efdfdcc1e1d76f55b51a 100644
--- a/Solver/linearSystemPETSc.h
+++ b/Solver/linearSystemPETSc.h
@@ -79,6 +79,7 @@ class linearSystemPETSc : public linearSystem<scalar> {
   std::vector<scalar> getData();
   std::vector<int> getRowPointers();
   std::vector<int> getColumnsIndices();
+  MPI_Comm& getComm() {return _comm;}
 };
 
 class linearSystemPETScBlockDouble : public linearSystem<fullMatrix<double> > {