Implemented multiple_output backend

Author: mdelorme

IOManager.h

@@ -20,15 +20,80 @@ class IOManager {
   ~IOManager() = default;
 
   void saveSolution(const Array &Q, int iteration, real_t t, real_t dt) {
+    if (params.multiple_outputs)
+      saveSolutionMultiple(Q, iteration, t, dt);
+    else
+      saveSolutionUnique(Q, iteration, t, dt);
+  }
+
+  void saveSolutionMultiple(const Array &Q, int iteration, real_t t, real_t dt) {
+    std::ostringstream oss;
+    
+    oss << params.filename_out << "_" << std::setw(4) << std::setfill('0') << iteration << ".h5";
+
+    File file(oss.str(), File::Truncate);
+
+    file.createAttribute("Ntx", params.Ntx);
+    file.createAttribute("Nty", params.Nty);
+    file.createAttribute("Nx", params.Nx);
+    file.createAttribute("Ny", params.Ny);
+    file.createAttribute("ibeg", params.ibeg);
+    file.createAttribute("iend", params.iend);
+    file.createAttribute("jbeg", params.jbeg);
+    file.createAttribute("jend", params.jend);
+    file.createAttribute("problem", params.problem);
+
+    std::vector<real_t> x;
+    for (int i=params.ibeg; i < params.iend; ++i)
+      x.push_back(((i-params.ibeg)+0.5) * params.dx);
+    file.createDataSet("x", x);
+    std::vector<real_t> y;
+    for (int j=params.jbeg; j < params.jend; ++j)
+      y.push_back(((j-params.jbeg)+0.5) * params.dy);
+    file.createDataSet("y", y);
+
+    using Table = std::vector<std::vector<real_t>>;
+
+    auto Qhost = Kokkos::create_mirror(Q);
+    Kokkos::deep_copy(Qhost, Q);
+
+    Table trho, tu, tv, tprs;
+    for (int j=params.jbeg; j<params.jend; ++j) {
+      std::vector<real_t> rrho, ru, rv, rprs;
+
+      for (int i=params.ibeg; i<params.iend; ++i) {
+        real_t rho = Qhost(j, i, IR);
+        real_t u   = Qhost(j, i, IU);
+        real_t v   = Qhost(j, i, IV);
+        real_t p   = Qhost(j, i, IP);
+
+        rrho.push_back(rho);
+        ru.push_back(u);
+        rv.push_back(v);
+        rprs.push_back(p);
+      }
+
+      trho.push_back(rrho);
+      tu.push_back(ru);
+      tv.push_back(rv);
+      tprs.push_back(rprs);
+    }
+
+    file.createDataSet("rho", trho);
+    file.createDataSet("u", tu);
+    file.createDataSet("v", tv);
+    file.createDataSet("prs", tprs);
+    file.createAttribute("time", t);
+  }
+
+  void saveSolutionUnique(const Array &Q, int iteration, real_t t, real_t dt) {
     std::ostringstream oss;
 
-    std::setw(4);
-    std::setfill('0');
-    oss << "ite_" << iteration;
+    oss << "ite_" << std::setw(4) << std::setfill('0') << iteration;
     std::string path = oss.str();
     auto flag = (iteration == 0 ? File::Truncate : File::ReadWrite);
 
-    File file(params.filename_out, flag);
+    File file(params.filename_out + ".h5", flag);
 
     if (iteration == 0) {
       file.createAttribute("Ntx", params.Ntx);

SimInfo.h

@@ -69,14 +69,16 @@ enum ViscosityMode {
 struct Params {
   real_t save_freq;
   real_t tend;
-  std::string filename_out = "run.h5";
+  std::string filename_out = "run";
   BoundaryType boundary_x = BC_REFLECTING;
   BoundaryType boundary_y = BC_REFLECTING;
   ReconstructionType reconstruction = PCM; 
   RiemannSolver riemann_solver = HLL;
   TimeStepping time_stepping = TS_EULER;
   real_t CFL = 0.1;
 
+  bool multiple_outputs = false;
+
   // Parallel stuff
   ParallelRange range_tot;
   ParallelRange range_dom;
@@ -179,6 +181,7 @@ Params readInifile(std::string filename) {
 
   // Run
   res.tend = reader.GetFloat("run", "tend", 1.0);
+  res.multiple_outputs = reader.GetBoolean("run", "multiple_outputs", false);
   res.save_freq = reader.GetFloat("run", "save_freq", 1.0e-1);
   res.filename_out = reader.Get("run", "output_filename", "run.h5");
 

settings/blast.ini

@@ -9,7 +9,7 @@ ymax=1.0
 [run]
 tend=1.0
 save_freq=0.01
-
+multiple_outputs=true
 boundaries_x=periodic
 boundaries_y=periodic