[TrilinosApplication] Improve TrilinosSpace, its python bindings and MPI test coverage

applications/TrilinosApplication/custom_python/add_trilinos_space_to_python.cpp

@@ -250,6 +250,61 @@ Vector GatherValues(TrilinosSparseSpaceType& dummy,
     dummy.sublist(sublist_name).set(name, value);
 }
 
+bool WriteMatrixMarketMatrixHelper(TrilinosSparseSpaceType& dummy, const char* filename, const TrilinosSparseSpaceType::MatrixType& rM, const bool symmetric)
+{
+    return TrilinosSparseSpaceType::WriteMatrixMarketMatrix(filename, rM, symmetric);
+}
+
+bool WriteMatrixMarketVectorHelper(TrilinosSparseSpaceType& dummy, const char* filename, const TrilinosSparseSpaceType::VectorType& rV)
+{
+    return TrilinosSparseSpaceType::WriteMatrixMarketVector(filename, rV);
+}
+
+void GlobalAssembleMatrix(TrilinosSparseSpaceType& dummy, TrilinosSparseSpaceType::MatrixType& rA)
+{
+    TrilinosSparseSpaceType::GlobalAssemble(rA);
+}
+
+void GlobalAssembleVector(TrilinosSparseSpaceType& dummy, TrilinosSparseSpaceType::VectorType& rV)
+{
+    TrilinosSparseSpaceType::GlobalAssemble(rV);
+}
+
+void SetVector(TrilinosSparseSpaceType& dummy, TrilinosSparseSpaceType::VectorType& rX, double value)
+{
+    TrilinosSparseSpaceType::Set(rX, value);
+}
+
+double GetDiagonalNormHelper(TrilinosSparseSpaceType& dummy, const TrilinosSparseSpaceType::MatrixType& rA)
+{
+    return TrilinosSparseSpaceType::GetDiagonalNorm(rA);
+}
+
+double GetAveragevalueDiagonalHelper(TrilinosSparseSpaceType& dummy, const TrilinosSparseSpaceType::MatrixType& rA)
+{
+    return TrilinosSparseSpaceType::GetAveragevalueDiagonal(rA);
+}
+
+double GetMaxDiagonalHelper(TrilinosSparseSpaceType& dummy, const TrilinosSparseSpaceType::MatrixType& rA)
+{
+    return TrilinosSparseSpaceType::GetMaxDiagonal(rA);
+}
+
+double GetMinDiagonalHelper(TrilinosSparseSpaceType& dummy, const TrilinosSparseSpaceType::MatrixType& rA)
+{
+    return TrilinosSparseSpaceType::GetMinDiagonal(rA);
+}
+
+void CopyMatrixValuesHelper(TrilinosSparseSpaceType& dummy, TrilinosSparseSpaceType::MatrixType& rA, const TrilinosSparseSpaceType::MatrixType& rB)
+{
+    TrilinosSparseSpaceType::CopyMatrixValues(rA, rB);
+}
+
+double GetValue(TrilinosSparseSpaceType& dummy, const TrilinosSparseSpaceType::VectorType& rX, std::size_t I)
+{
+    return TrilinosSparseSpaceType::GetValue(rX, I);
+}
+
 void  AddBasicOperations(pybind11::module& m)
 {
     py::class_< Epetra_MpiComm > (m,"Epetra_MpiComm")
@@ -310,10 +365,22 @@ void  AddBasicOperations(pybind11::module& m)
     .def("CreateEmptyVectorPointer", CreateEmptyVectorPointer)
     .def("ReadMatrixMarketMatrix", ReadMatrixMarketMatrix)
     .def("ReadMatrixMarketVector", ReadMatrixMarketVector)
+    .def("WriteMatrixMarketMatrix", WriteMatrixMarketMatrixHelper)
+    .def("WriteMatrixMarketVector", WriteMatrixMarketVectorHelper)
     .def("SetValue", SetValue)
     //.def("GetValue", GetValue) //deliberately commented out. Only works for local Ids
     .def("GatherValues", GatherValues)
+    .def("GlobalAssembleMatrix", GlobalAssembleMatrix)
+    .def("GlobalAssembleVector", GlobalAssembleVector)
+    .def("Set", SetVector)
+    .def("TransposeMult", TransposeMult)
+    .def("GetDiagonalNorm", GetDiagonalNormHelper)
+    .def("GetAveragevalueDiagonal", GetAveragevalueDiagonalHelper)
+    .def("GetMaxDiagonal", GetMaxDiagonalHelper)
+    .def("GetMinDiagonal", GetMinDiagonalHelper)
+    .def("CopyMatrixValues", CopyMatrixValuesHelper)
     .def_static("IsDistributed", &TrilinosSparseSpaceType::IsDistributed)
+    .def_static("IsDistributedSpace", &TrilinosSparseSpaceType::IsDistributedSpace)
     .def_static("FastestDirectSolverList", &TrilinosSparseSpaceType::FastestDirectSolverList)
     ;