From 9e7ad10e0ab676813f7f0822e842df5ae1af1992 Mon Sep 17 00:00:00 2001
From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr>
Date: Tue, 22 May 2018 14:43:44 +0200
Subject: [PATCH] correction coefficient tau
---
src/main.cpp | 4 ++--
src/scheme/FiniteVolumesEulerUnknowns.hpp | 1 -
2 files changed, 2 insertions(+), 3 deletions(-)
diff --git a/src/main.cpp b/src/main.cpp
index 4bd3d0bd7..28ef2a46e 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -260,7 +260,7 @@ int main(int argc, char *argv[])
{ // gnuplot output for density
const Kokkos::View<const Rd*> xj = mesh_data.xj();
const Kokkos::View<const Rd*> uj = unknowns.uj();
- double h = std::sqrt(1. - (0.2*0.2)/((50./6.)*(50./6.)));
+ double h = std::sqrt(1. - (0.2*0.2)/(50./9.));
std::ofstream fout("rho");
fout.precision(15);
for (size_t j=0; j<mesh.numberOfCells(); ++j) {
@@ -287,7 +287,7 @@ int main(int argc, char *argv[])
const Kokkos::View<const Rd*> xj = mesh_data.xj();
const Kokkos::View<const double*> ej = unknowns.ej();
//double pi = 4.*std::atan(1.);
- double h = std::sqrt(1. - (0.2*0.2)/((50./6.)*(50./6.)));
+ double h = std::sqrt(1. - (0.2*0.2)/(50./9.));
std::ofstream fout("e");
fout.precision(15);
for (size_t j=0; j<mesh.numberOfCells(); ++j) {
diff --git a/src/scheme/FiniteVolumesEulerUnknowns.hpp b/src/scheme/FiniteVolumesEulerUnknowns.hpp
index 94ea49674..075fedf94 100644
--- a/src/scheme/FiniteVolumesEulerUnknowns.hpp
+++ b/src/scheme/FiniteVolumesEulerUnknowns.hpp
@@ -224,7 +224,6 @@ public:
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] = 2.;
});
const Kokkos::View<const double*> Vj = m_mesh_data.Vj();
--
GitLab