From d9bbdfd8316642aaafc309c86e46b560317decfe Mon Sep 17 00:00:00 2001
From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr>
Date: Thu, 17 May 2018 16:03:34 +0200
Subject: [PATCH] ajout gnuplot densite rho
---
src/main.cpp | 10 ++++++++++
1 file changed, 10 insertions(+)
diff --git a/src/main.cpp b/src/main.cpp
index a0acdd5f7..e939a34cd 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -225,6 +225,16 @@ int main(int argc, char *argv[])
//method_cost_map["AcousticSolverWithMesh"] = timer.seconds();
method_cost_map["FiniteVolumesDiffusionWithMesh"] = timer.seconds();
+ { // gnuplot output for density
+ const Kokkos::View<const Rd*> xj = mesh_data.xj();
+ const Kokkos::View<const Rd*> uj = unknowns.uj();
+ std::ofstream fout("resultat rho");
+ fout.precision(15);
+ for (size_t j=0; j<mesh.numberOfCells(); ++j) {
+ fout << xj[j][0] << ' ' << rhoj[j] << '\n';
+ }
+ }
+
{ // gnuplot output for vitesse
const Kokkos::View<const Rd*> xj = mesh_data.xj();
const Kokkos::View<const Rd*> uj = unknowns.uj();
--
GitLab