[PL] Enable local matrix debug output

ProcessLib/VectorMatrixAssembler.cpp

@@ -100,6 +100,9 @@ void VectorMatrixAssembler::assemble(
         assert(_local_b_data.size() == num_r_c);
         b.add(indices, _local_b_data);
     }
+
+    _local_output(t, process_id, mesh_item_id, _local_M_data, _local_K_data,
+                  _local_b_data);
 }
 
 void VectorMatrixAssembler::assembleWithJacobian(
@@ -180,6 +183,9 @@ void VectorMatrixAssembler::assembleWithJacobian(
             "No Jacobian has been assembled! This might be due to programming "
             "errors in the local assembler of the current process.");
     }
+
+    _local_output(t, process_id, mesh_item_id, _local_M_data, _local_K_data,
+                  _local_b_data, &_local_Jac_data);
 }
 
 }  // namespace ProcessLib

ProcessLib/VectorMatrixAssembler.h

@@ -11,14 +11,14 @@
 #pragma once
 
 #include <vector>
-#include "NumLib/NumericsConfig.h"
+
 #include "AbstractJacobianAssembler.h"
-#include "CoupledSolutionsForStaggeredScheme.h"
+#include "Assembly/MatrixOutput.h"
 
 namespace NumLib
 {
 class LocalToGlobalIndexMap;
-}  // NumLib
+}  // namespace NumLib
 
 namespace ProcessLib
 {
@@ -74,6 +74,8 @@ private:
 
     //! Used to assemble the Jacobian.
     std::unique_ptr<AbstractJacobianAssembler> _jacobian_assembler;
+
+    Assembly::LocalMatrixOutput _local_output;
 };
 
 }  // namespace ProcessLib