Skip to content
Snippets Groups Projects
Commit 3e5d06ea authored by Fanny CHOPOT's avatar Fanny CHOPOT
Browse files

test avec k = 2 et rho = 1, OK

parent e285f88a
No related branches found
No related tags found
No related merge requests found
...@@ -184,10 +184,12 @@ int main(int argc, char *argv[]) ...@@ -184,10 +184,12 @@ int main(int argc, char *argv[])
double pi = 4.*std::atan(1.); double pi = 4.*std::atan(1.);
std::ofstream fout("comparaison u"); std::ofstream fout("comparaison u");
for (size_t j=0; j<mesh.numberOfCells(); ++j) { for (size_t j=0; j<mesh.numberOfCells(); ++j) {
fout << xj[j][0] << ' ' << uj[j][0] << ' ' << std::sin(pi*xj[j][0])*std::exp(-pi*pi*0.2) <<'\n'; fout << xj[j][0] << ' ' << uj[j][0] << ' ' << std::sin(pi*xj[j][0])*std::exp(-2.*pi*pi*0.2) <<'\n';
} }
} }
} }
Kokkos::finalize(); Kokkos::finalize();
......
...@@ -100,9 +100,9 @@ private: ...@@ -100,9 +100,9 @@ private:
//Fl(l) = ((kj(cell_nodes(j,r)) + kj(cell_nodes(j-1,r)))/(2*Vl(j))) //Fl(l) = ((kj(cell_nodes(j,r)) + kj(cell_nodes(j-1,r)))/(2*Vl(j)))
//* ((uj(j,r),Cjr(j,r))+ (uj(j-1,r),Cjr(j-1,r))) ; //tensorProduct(uj(j,r),Cjr(j,r)) ? //* ((uj(j,r),Cjr(j,r))+ (uj(j-1,r),Cjr(j-1,r))) ; //tensorProduct(uj(j,r),Cjr(j,r)) ?
// k = 1 // k = 2
m_Fl(l)= (1./Vl(l))*sum; m_Fl(l)= (2./Vl(l))*sum;
}); });
......
...@@ -222,7 +222,7 @@ void initializeSod() ...@@ -222,7 +222,7 @@ void initializeSod()
}); });
Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j){ Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j){
m_kj[j] = 1; m_kj[j] = 2;
}); });
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment