diff --git a/src/main.cpp b/src/main.cpp index c8b6fe135b1a14a035f91d1e09287e205e990d1e..2225d80d2faae3c6da7565c74f5041a0c17d16c0 100644 --- a/src/main.cpp +++ b/src/main.cpp @@ -141,7 +141,7 @@ int main(int argc, char *argv[]) const Kokkos::View<const double*> Vj = mesh_data.Vj(); const Kokkos::View<const Rd**> Cjr = mesh_data.Cjr(); - const double tmax=0.5; + const double tmax=0.2; double t=0.; int itermax=std::numeric_limits<int>::max(); @@ -153,7 +153,7 @@ int main(int argc, char *argv[]) Kokkos::View<double*> gammaj = unknowns.gammaj(); Kokkos::View<double*> cj = unknowns.cj(); Kokkos::View<double*> kj = unknowns.kj(); - + Kokkos::View<Rd*> uj = unknowns.uj(); BlockPerfectGas block_eos(rhoj, ej, pj, gammaj, cj); double c = 0.; @@ -188,7 +188,7 @@ int main(int argc, char *argv[]) t += dt_euler; // ETAPE 2 DU SPLITTING - DIFFUSION - + /* double dt_diff = 0.4*finite_volumes_diffusion.diffusion_dt(rhoj, kj, cj); if (dt_euler <= dt_diff) { @@ -205,7 +205,7 @@ int main(int argc, char *argv[]) t_diff += dt_diff; } } - + */ // DIFFUSION PURE /* @@ -234,7 +234,7 @@ int main(int argc, char *argv[]) //fout << ligne << ' ' << std::fixed << t << '\n'; for (size_t j = 0; j<mesh.numberOfCells(); ++j) { getline(fint, ligne); - fout << ligne << ' ' << std::fixed << rhoj[j] << '\n'; + fout << ligne << ' ' << std::fixed << xj[j][0] << ' ' << std::fixed << rhoj[j] << '\n'; } fint.close(); fout.close(); diff --git a/src/scheme/FiniteVolumesEulerUnknowns.hpp b/src/scheme/FiniteVolumesEulerUnknowns.hpp index 34e9e84a8d7169c808d877d28720c46263786f3a..3579c7d7a82f8a1c7fb3c16c39f229b2ce2a86f0 100644 --- a/src/scheme/FiniteVolumesEulerUnknowns.hpp +++ b/src/scheme/FiniteVolumesEulerUnknowns.hpp @@ -192,7 +192,7 @@ public: Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j){ if (xj[j][0]<0.5) { - m_rhoj[j]=1; + m_rhoj[j]=1.; } else { m_rhoj[j]=0.125; } @@ -242,9 +242,9 @@ public: m_kj[j]=0.; } else { if (xj[j][0]<0.9){ - m_kj[j]=0.5; - } else { m_kj[j]=0.; + } else { + m_kj[j]=0. ; } } });