From e285f88ab4beddbd1044f4cf9b15bc1ca57e2ee4 Mon Sep 17 00:00:00 2001
From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr>
Date: Thu, 26 Apr 2018 16:01:19 +0200
Subject: [PATCH] ajout gnuplot u

---
 src/main.cpp                          | 7 ++++---
 src/scheme/FiniteVolumesDiffusion.hpp | 2 +-
 2 files changed, 5 insertions(+), 4 deletions(-)

diff --git a/src/main.cpp b/src/main.cpp
index e122cb5b9..e4df1a5e7 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -180,10 +180,11 @@ int main(int argc, char *argv[])
     
     { // gnuplot output for density
       const Kokkos::View<const Rd*> xj   = mesh_data.xj();
-      const Kokkos::View<const double*> rhoj = unknowns.rhoj();
-      std::ofstream fout("rho");
+      const Kokkos::View<const Rd*> uj = unknowns.uj();
+      double pi = 4.*std::atan(1.);
+      std::ofstream fout("comparaison u");
       for (size_t j=0; j<mesh.numberOfCells(); ++j) {
-        fout << xj[j][0] << ' ' << rhoj[j] << '\n';
+	fout << xj[j][0] << ' ' << uj[j][0] <<  ' ' << std::sin(pi*xj[j][0])*std::exp(-pi*pi*0.2) <<'\n';
       }
     }
 
diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp
index 93440138c..74e59d14f 100644
--- a/src/scheme/FiniteVolumesDiffusion.hpp
+++ b/src/scheme/FiniteVolumesDiffusion.hpp
@@ -102,7 +102,7 @@ private:
 
         // k = 1
 	 
-	m_Fl(l)= (1./(2*Vl(l)))*sum;
+	m_Fl(l)= (1./Vl(l))*sum;
 
       });
 
-- 
GitLab