diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-02-09 11:21:55 +0100 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-02-09 11:21:55 +0100 |
| commit | 6b10cae1ffa2521ce3c7e6aec85513a24d25934b (patch) | |
| tree | 97f9b3662df3d689da2a21d68723a60cceeb6b54 /examples/poiseulle_particles_2d_gpu/sim.cpp | |
| parent | 22659a4a4af8be137bcd702b89f1f44cd0fe2455 (diff) | |
| download | libs-lbm-6b10cae1ffa2521ce3c7e6aec85513a24d25934b.tar.gz | |
Set equilibrium boundary for sim
Diffstat (limited to 'examples/poiseulle_particles_2d_gpu/sim.cpp')
| -rw-r--r-- | examples/poiseulle_particles_2d_gpu/sim.cpp | 14 |
1 files changed, 12 insertions, 2 deletions
diff --git a/examples/poiseulle_particles_2d_gpu/sim.cpp b/examples/poiseulle_particles_2d_gpu/sim.cpp index d53a00c..6fdff7b 100644 --- a/examples/poiseulle_particles_2d_gpu/sim.cpp +++ b/examples/poiseulle_particles_2d_gpu/sim.cpp @@ -119,6 +119,14 @@ saw::error_or<void> step( // Need nicer things to handle the flow. I see improvement here component<T,Desc,cmpt::BGK,encode::Sycl<saw::encode::Native>> collision{0.6}; component<T,Desc,cmpt::BounceBack,encode::Sycl<saw::encode::Native>> bb; + + saw::data<sch::Scalar<T>> rho_b; + rho_b.at({}) = 1.0; + saw::data<sch::Vector<T,Desc::D>> vel_b; + vel_b.at({{0u}}) = 0.01; + + component<T,Desc,cmpt::Equilibrium,encode::Sycl<saw::encode::Native>> equi{rho_b,vel_b}; + component<T,Desc,cmpt::ZouHeHorizontal<true>,encode::Sycl<saw::encode::Native>> flow_in{ [&](){ uint64_t target_t_i = 256u; @@ -148,11 +156,13 @@ saw::error_or<void> step( collision.apply(fields,macros,index,t_i); break; case 3u: - flow_in.apply(fields,index,t_i); + equi.apply(fields,index,t_i); + // flow_in.apply(fields,index,t_i); collision.apply(fields,macros,index,t_i); break; case 4u: - flow_out.apply(fields,index,t_i); + equi.apply(fields,index,t_i); + // flow_out.apply(fields,index,t_i); collision.apply(fields,macros,index,t_i); break; default: |