From 6aff1078884c820b110f1355c9b0d04a9dd45cc2 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Wed, 4 Feb 2026 16:52:18 +0100 Subject: Semi working state. Broken vtk ordering --- examples/poiseulle_particles_2d_gpu/sim.cpp | 22 ++++++++++++++++++---- 1 file changed, 18 insertions(+), 4 deletions(-) (limited to 'examples/poiseulle_particles_2d_gpu') diff --git a/examples/poiseulle_particles_2d_gpu/sim.cpp b/examples/poiseulle_particles_2d_gpu/sim.cpp index 3084bca..bdd93b0 100644 --- a/examples/poiseulle_particles_2d_gpu/sim.cpp +++ b/examples/poiseulle_particles_2d_gpu/sim.cpp @@ -83,7 +83,7 @@ saw::error_or setup_initial_conditions( {{0u,1u}} ); // - auto& df_f = fields.template get<"dfs">(); + auto& df_f = fields.template get<"dfs_old">(); auto& rho_f = macros.template get<"density">(); auto& vel_f = macros.template get<"velocity">(); @@ -248,6 +248,15 @@ saw::error_or lbm_main(int argc, char** argv){ return eov; } } + { + std::string file_name = "tmp/df_"; + file_name += std::to_string(i.get()); + file_name += ".vtk"; + auto eov = write_vtk_file(file_name, *lbm_data_ptr); + if(eov.is_error()){ + return eov; + } + } { auto eov = step(lsd_view,lsdm_view,i,dev); if(eov.is_error()){ @@ -276,9 +285,14 @@ saw::error_or lbm_main(int argc, char** argv){ if(eov.is_error()){ return eov; } - auto eov2 = write_vtk_file((std::string{"tmp/df_"}+std::to_string(time_steps.get())+std::string{".vtk"}), *lbm_data_ptr); - if(eov2.is_error()){ - return eov2; + } + { + std::string file_name = "tmp/df_"; + file_name += std::to_string(time_steps.get()); + file_name += ".vtk"; + auto eov = write_vtk_file(file_name, *lbm_data_ptr); + if(eov.is_error()){ + return eov; } } -- cgit v1.2.3