correction coefficient tau

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) {

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();