diff --git a/src/main.cpp b/src/main.cpp index 4bd3d0bd72dba4c60c223aa54a028f03005d5dc2..28ef2a46e73bff2f49dbb52ed047e89c39d24367 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 94ea4967440a785bc5c93c9a260b17584f6ef763..075fedf94608e6329c8da6fda9ead9d81a9a4e80 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();