//
// ------------------------------------------------------------------------
+#include <deal.II/base/mpi_stub.h>
+
#include <deal.II/lac/petsc_block_sparse_matrix.h>
#include <deal.II/lac/petsc_compatibility.h>
+
#ifdef DEAL_II_WITH_PETSC
+# include <petscmat.h>
+
+
+DEAL_II_NAMESPACE_OPEN
+
namespace
{
// A dummy utility routine to create an empty matrix in case we import
}
} // namespace
-DEAL_II_NAMESPACE_OPEN
namespace PETScWrappers
{
std::to_string(c) +
" is completely empty "
"and so it is not possible to determine how many columns it should have."));
- Mat dummy = ::create_dummy_mat(
- comm,
- static_cast<PetscInt>(row_local_sizes[r]),
- static_cast<PetscInt>(row_sizes[r]),
- static_cast<PetscInt>(col_local_sizes[c]),
- static_cast<PetscInt>(col_sizes[c]));
+ Mat dummy =
+ create_dummy_mat(comm,
+ static_cast<PetscInt>(row_local_sizes[r]),
+ static_cast<PetscInt>(row_sizes[r]),
+ static_cast<PetscInt>(col_local_sizes[c]),
+ static_cast<PetscInt>(col_sizes[c]));
this->sub_objects[r][c] = new BlockType(dummy);
// the new object got a reference on dummy, we can safely