From 75fda9a23cd30efa72bd212d5642c3d30ee34913 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" 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