diff options
Diffstat (limited to 'examples')
| -rw-r--r-- | examples/poiseulle_particles_2d_gpu/sim.cpp | 27 |
1 files changed, 7 insertions, 20 deletions
diff --git a/examples/poiseulle_particles_2d_gpu/sim.cpp b/examples/poiseulle_particles_2d_gpu/sim.cpp index 691ec93..6d5452e 100644 --- a/examples/poiseulle_particles_2d_gpu/sim.cpp +++ b/examples/poiseulle_particles_2d_gpu/sim.cpp @@ -26,11 +26,14 @@ template<typename T, typename Desc> using ScalarChunk = Chunk<Scalar<T>, 0u, dim_x, dim_y>; template<typename T, typename Desc> +using VectorChunk = Chunk<Vector<T,Desc::D>, 0u, dim_x, dim_y>; + +template<typename T, typename Desc> using ChunkStruct = Struct< Member<InfoChunk, "info">, Member<DfChunk<T,Desc>, "dfs">, Member<DfChunk<T,Desc>, "dfs_old">, - Member<ScalarChunk<T,Desc>, "particle_N">, + Member<VectorChunk<T,Desc>, "particle_N">, Member<ScalarChunk<T,Desc>, "particle_D"> >; @@ -102,11 +105,13 @@ saw::error_or<void> setup_initial_conditions( auto& df_f = fields.template get<"dfs_old">(); auto& rho_f = macros.template get<"density">(); auto& vel_f = macros.template get<"velocity">(); + auto& por_f = macros.template get<"porosity">(); iterator<Desc::D>::apply( [&](auto& index){ auto& df = df_f.at(index); auto& rho = rho_f.at(index); + por_f.at(index).at({}) = {1}; rho.at({}) = {1}; auto& vel = vel_f.at(index); auto eq = equilibrium<T,Desc>(rho,vel); @@ -324,7 +329,7 @@ saw::error_or<void> lbm_main(int argc, char** argv){ } } sycl_q.wait(); - saw::data<sch::UInt64> time_steps{particle_size}; + saw::data<sch::UInt64> time_steps{256ul}; for(saw::data<sch::UInt64> i{0u}; i < time_steps and krun; ++i){ { @@ -348,15 +353,6 @@ saw::error_or<void> lbm_main(int argc, char** argv){ }*/ } { - std::string file_name = "df_"; - file_name += std::to_string(i.get()); - file_name += ".vtk"; - auto eov = write_vtk_file(out_dir/file_name, *lbm_data_ptr); - if(eov.is_error()){ - return eov; - } - } - { auto eov = step<T,Desc>(lsd_view,lsdm_view,i,dev); if(eov.is_error()){ return eov; @@ -391,15 +387,6 @@ saw::error_or<void> lbm_main(int argc, char** argv){ return eov; } } - { - std::string file_name = "df_"; - file_name += std::to_string(time_steps.get()); - file_name += ".vtk"; - auto eov = write_vtk_file(out_dir/file_name, *lbm_data_ptr); - if(eov.is_error()){ - return eov; - } - } /* iterator<Desc::D>::apply( |