#1533 Add a Petsc wrapper around MatEqual.

Sources/ThirdParty/Wrappers/Petsc/Matrix/Matrix.cpp

@@ -666,6 +666,29 @@ namespace MoReFEM
                 mpi.Barrier();
             }
 
+
+            bool AreStrictlyEqual(const Matrix & lhs,
+                                  const Matrix& rhs,
+                                  const char* invoking_file, int invoking_line)
+            {
+                PetscBool value;
+                int error_code = MatEqual(lhs.Internal(), rhs.Internal(), &value);
+
+                if (error_code)
+                    throw ExceptionNS::Exception(error_code, "MatEqual", invoking_file, invoking_line);
+
+                switch (value)
+                {
+                    case PETSC_TRUE:
+                        return true;
+                    case PETSC_FALSE:
+                        return false;
+                }
+
+                assert(false);
+                exit(EXIT_FAILURE);
+            }
+
             
         } // namespace Petsc
         

Sources/ThirdParty/Wrappers/Petsc/Matrix/Matrix.hpp

@@ -621,6 +621,23 @@ namespace MoReFEM
             void Swap(Matrix& , Matrix& );
 
 
+            /*!
+             * \brief Checks whether Petsc matrices are equal.
+             *
+             * This is a thin wrapper over PETSc's \a MatEqual, it probably won't work if there are rounding errors (see the namesake for PETSc vectors -
+             * I had to resort to something else than \a VecEqual for this reason).
+             *
+             * \copydoc doxygen_hide_lhs_rhs_arg
+             * \copydoc doxygen_hide_invoking_file_and_line
+             *
+             * \return True if lhs and rhs are identical (same pointer or bitwise equality).
+             */
+            bool AreStrictlyEqual(const Matrix & lhs,
+                                  const Matrix& rhs,
+                                  const char* invoking_file, int invoking_line);
+
+
+
         } // namespace Petsc