#820 Add dev file to compare with Freefem.

Sources/ModelInstances/UnderDevelopment/Poromechanics/ImplicitStepFluid/InnerLoop/VariationalFormulation.hxx

@@ -166,6 +166,19 @@ namespace HappyHeart
                         GetDifferentialDarcyOperator().Assemble(std::make_tuple(std::ref(contribution_with_coeff)));
                     }
                     
+                    const auto filename = parent::GetOutputDirectory(numbering_subset)
+                    + "/"
+                    + DifferentialPreffix(differential::yes)
+                    + "rhs_" + GetIterationTag() + ".m";
+                    
+                    std::cout << "FILE = " << filename << std::endl;
+                    
+                    rhs.View(parent::MpiHappyHeart(),
+                             filename,
+                             __FILE__, __LINE__,
+                             PETSC_VIEWER_ASCII_MATLAB);
+
+                    
                     
                     // \todo #820 At the moment, I accept my norm is not the same as Bruno's; see what Petsc gives when
                     // SNES is applied (system_rhs's is what we got in hyperelastic model).