From 40f0f478154a7e5558f5acab4d8387d070c08e35 Mon Sep 17 00:00:00 2001 From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr> Date: Mon, 7 May 2018 11:33:50 +0200 Subject: [PATCH] Correction solution exacte et condition initiale de E --- src/main.cpp | 2 +- src/scheme/FiniteVolumesEulerUnknowns.hpp | 1 - 2 files changed, 1 insertion(+), 2 deletions(-) diff --git a/src/main.cpp b/src/main.cpp index 55a0362a6..5ee448520 100644 --- a/src/main.cpp +++ b/src/main.cpp @@ -216,7 +216,7 @@ int main(int argc, char *argv[]) double pi = 4.*std::atan(1.); std::ofstream fout("comparaison E"); for (size_t j=0; j<mesh.numberOfCells(); ++j) { - fout << xj[j][0] << ' ' << Ej[j] << ' ' << (-(std::cos(pi*xj[j][0]))*(std::cos(pi*xj[j][0]))+(std::sin(pi*xj[j][0])*std::sin(pi*xj[j][0])))*(pi*0.5)*std::exp(-4.*pi*0.2) <<'\n'; + fout << xj[j][0] << ' ' << Ej[j] << ' ' << (-(std::cos(pi*xj[j][0])*std::cos(pi*xj[j][0]))+(std::sin(pi*xj[j][0])*std::sin(pi*xj[j][0])))*(pi*0.5)*(std::exp(-4.*pi*0.2)-1.) + 2. <<'\n'; // cas k constant } } diff --git a/src/scheme/FiniteVolumesEulerUnknowns.hpp b/src/scheme/FiniteVolumesEulerUnknowns.hpp index 1db575470..a39dd4157 100644 --- a/src/scheme/FiniteVolumesEulerUnknowns.hpp +++ b/src/scheme/FiniteVolumesEulerUnknowns.hpp @@ -227,7 +227,6 @@ void initializeSod() Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j){ //m_Ej[j] = m_ej[j]+0.5*(m_uj[j],m_uj[j]); - //m_Ej[j] = (std::cos(pi*xj[j][0]))*(std::cos(pi*xj[j][0]))*std::sin(pi*xj[j][0])*((pi*pi*pi)*0.5); m_Ej[j] = 2.; }); -- GitLab