diff --git a/packages/muelu/adapters/xpetra/MueLu_RefMaxwell_def.hpp b/packages/muelu/adapters/xpetra/MueLu_RefMaxwell_def.hpp index 38985e1b7934..2e741012d965 100644 --- a/packages/muelu/adapters/xpetra/MueLu_RefMaxwell_def.hpp +++ b/packages/muelu/adapters/xpetra/MueLu_RefMaxwell_def.hpp @@ -454,6 +454,18 @@ namespace MueLu { std::copy(BCrows_.begin(), BCrows_.end(), std::ostream_iterator(outBCrows, "\n")); std::ofstream outBCcols("BCcols.mat"); std::copy(BCcols_.begin(), BCcols_.end(), std::ostream_iterator(outBCcols, "\n")); +#else + std::ofstream outBCrows("BCrows.mat"); + auto BCrows = Kokkos::create_mirror_view (BCrows_); + Kokkos::deep_copy(BCrows , BCrows_); + for (size_t i = 0; i < BCrows.size(); i++) + outBCrows << BCrows[i] << "\n"; + + std::ofstream outBCcols("BCcols.mat"); + auto BCcols = Kokkos::create_mirror_view (BCcols_); + Kokkos::deep_copy(BCcols , BCcols_); + for (size_t i = 0; i < BCcols.size(); i++) + outBCcols << BCcols[i] << "\n"; #endif Xpetra::IO::Write(std::string("nullspace.mat"), *Nullspace_); if (Coords_ != Teuchos::null)