MueLu RefMaxwell: Dump BC info for Kokkos

trilinos · Jun 16, 2018 · e6f4bd9 · e6f4bd9
1 parent e89a132
commit e6f4bd9
Showing 1 changed file with 12 additions and 0 deletions.
diff --git 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<LO>(outBCrows, "\n"));
       std::ofstream outBCcols("BCcols.mat");
       std::copy(BCcols_.begin(), BCcols_.end(), std::ostream_iterator<LO>(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<SC, LO, GlobalOrdinal, Node>::Write(std::string("nullspace.mat"), *Nullspace_);
       if (Coords_ != Teuchos::null)