From 77afc46bf0cb5169edbec626a7c60996c81b2154 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Tue, 3 Mar 2026 16:52:53 +0100 Subject: Confusion about Boundary conditions --- examples/poiseulle_particles_2d_gpu/sim.cpp | 53 ++++++++++++++++------------- 1 file changed, 29 insertions(+), 24 deletions(-) (limited to 'examples') diff --git a/examples/poiseulle_particles_2d_gpu/sim.cpp b/examples/poiseulle_particles_2d_gpu/sim.cpp index 4bc9dfb..e7aec8e 100644 --- a/examples/poiseulle_particles_2d_gpu/sim.cpp +++ b/examples/poiseulle_particles_2d_gpu/sim.cpp @@ -259,22 +259,6 @@ saw::error_or step( }); }).wait(); - q.submit([&](acpp::sycl::handler& h){ - component> stream; - - h.parallel_for(acpp::sycl::range{dim_x,dim_y}, [=](acpp::sycl::id idx){ - saw::data> index; - for(uint64_t i = 0u; i < Desc::D; ++i){ - index.at({{i}}).set(idx[i]); - } - - auto info = info_f.at(index); - - if(info.get() > 0u){ - stream.apply(fields,index,t_i); - } - }); - }).wait(); // Step /* @@ -380,35 +364,56 @@ saw::error_or lbm_main(int argc, char** argv){ } } sycl_q.wait(); - saw::data time_steps{4096ul}; + saw::data time_steps{128ul}; + + auto& info_f = lbm_sycl_data.template get<"info">(); for(saw::data i{0u}; i < time_steps and krun; ++i){ + // BC + Collision { - auto eov = dev.copy_to_host(lbm_sycl_macro_data,*lbm_macro_data_ptr); + auto eov = step(lsd_view,lsdm_view,*lbm_particle_data_ptr,i,dev); if(eov.is_error()){ return eov; } } + sycl_q.wait(); { - { - auto eov = write_vtk_file(out_dir,"t",i.get(), *lbm_macro_data_ptr); - if(eov.is_error()){ - return eov; - } + auto eov = dev.copy_to_host(lbm_sycl_macro_data,*lbm_macro_data_ptr); + if(eov.is_error()){ + return eov; } } { - auto eov = step(lsd_view,lsdm_view,*lbm_particle_data_ptr,i,dev); + auto eov = write_vtk_file(out_dir,"t",time_steps.get(), *lbm_macro_data_ptr); if(eov.is_error()){ return eov; } } + // Stream + sycl_q.submit([&](acpp::sycl::handler& h){ + component> stream; + + h.parallel_for(acpp::sycl::range{dim_x,dim_y}, [=](acpp::sycl::id idx){ + saw::data> index; + for(uint64_t i = 0u; i < Desc::D; ++i){ + index.at({{i}}).set(idx[i]); + } + + auto info = info_f.at(index); + + if(info.get() > 0u){ + stream.apply(fields,index,t_i); + } + }); + }).wait(); wait.poll(); if(print_status){ std::cout<<"Status: "<().get() * 100 / time_steps.get())<<"%"< Date: Tue, 3 Mar 2026 19:03:05 +0100 Subject: Added CSV for simpler script processing --- examples/poiseulle_particles_2d_gpu/sim.cpp | 42 ++++++++++++++++++++++------- 1 file changed, 32 insertions(+), 10 deletions(-) (limited to 'examples') diff --git a/examples/poiseulle_particles_2d_gpu/sim.cpp b/examples/poiseulle_particles_2d_gpu/sim.cpp index e7aec8e..6d012e0 100644 --- a/examples/poiseulle_particles_2d_gpu/sim.cpp +++ b/examples/poiseulle_particles_2d_gpu/sim.cpp @@ -337,6 +337,12 @@ saw::error_or lbm_main(int argc, char** argv){ return eov; } } + { + auto eov = write_vtk_file(out_dir,"initial_state",0u,*lbm_data_ptr); + if(eov.is_error()){ + return eov; + } + } saw::data, encode::Sycl> lbm_sycl_data{sycl_q}; saw::data, encode::Sycl> lbm_sycl_macro_data{sycl_q}; @@ -364,9 +370,9 @@ saw::error_or lbm_main(int argc, char** argv){ } } sycl_q.wait(); - saw::data time_steps{128ul}; + saw::data time_steps{32ul}; - auto& info_f = lbm_sycl_data.template get<"info">(); + auto& info_f = lsd_view.template get<"info">(); for(saw::data i{0u}; i < time_steps and krun; ++i){ // BC + Collision @@ -378,15 +384,31 @@ saw::error_or lbm_main(int argc, char** argv){ } sycl_q.wait(); { - auto eov = dev.copy_to_host(lbm_sycl_macro_data,*lbm_macro_data_ptr); - if(eov.is_error()){ - return eov; + { + auto eov = dev.copy_to_host(lbm_sycl_macro_data,*lbm_macro_data_ptr); + if(eov.is_error()){ + return eov; + } + } + { + auto eov = write_vtk_file(out_dir,"m",i.get(), *lbm_macro_data_ptr); + if(eov.is_error()){ + return eov; + } } } { - auto eov = write_vtk_file(out_dir,"t",time_steps.get(), *lbm_macro_data_ptr); - if(eov.is_error()){ - return eov; + { + auto eov = dev.copy_to_host(lbm_sycl_data,*lbm_data_ptr); + if(eov.is_error()){ + return eov; + } + } + { + auto eov = write_csv_file(out_dir,"lbm",i.get(), *lbm_data_ptr); + if(eov.is_error()){ + return eov; + } } } // Stream @@ -402,7 +424,7 @@ saw::error_or lbm_main(int argc, char** argv){ auto info = info_f.at(index); if(info.get() > 0u){ - stream.apply(fields,index,t_i); + stream.apply(lsd_view,index,i); } }); }).wait(); @@ -416,7 +438,7 @@ saw::error_or lbm_main(int argc, char** argv){ // After Loop sycl_q.wait(); { - auto eov = write_vtk_file(out_dir,"t",time_steps.get(), *lbm_macro_data_ptr); + auto eov = write_vtk_file(out_dir,"m",time_steps.get(), *lbm_macro_data_ptr); if(eov.is_error()){ return eov; } -- cgit v1.2.3 From 16700f013886b8d8cb02e322c2d66c00179d9c89 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Tue, 3 Mar 2026 19:46:31 +0100 Subject: Fixing ZouHe? --- examples/poiseulle_particles_2d_gpu/sim.cpp | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) (limited to 'examples') diff --git a/examples/poiseulle_particles_2d_gpu/sim.cpp b/examples/poiseulle_particles_2d_gpu/sim.cpp index 6d012e0..bb81383 100644 --- a/examples/poiseulle_particles_2d_gpu/sim.cpp +++ b/examples/poiseulle_particles_2d_gpu/sim.cpp @@ -131,7 +131,7 @@ saw::error_or setup_initial_conditions( rho.at({}) = {1}; auto& vel = vel_f.at(index); if(info_f.at(index).get() == 2u){ - vel.at({{0u}}) = 0.01; + vel.at({{0u}}) = 0.0; } auto eq = equilibrium(rho,vel); @@ -433,6 +433,7 @@ saw::error_or lbm_main(int argc, char** argv){ std::cout<<"Status: "<().get() * 100 / time_steps.get())<<"%"<