diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2025-09-03 18:04:36 +0200 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2025-09-03 18:04:36 +0200 |
| commit | 616b0011c3638eceaae8890c3a66037625587f1d (patch) | |
| tree | 0487d2ddddcaedaa3006fd83ddfaf6520c41f419 /examples | |
| parent | dd90f1efc84a5f0d92e97fce6085763697092eae (diff) | |
Error in guo forcing seems to persist
Some debug output left there, but there's no force. Might be that the
higher precision is impacting this
Diffstat (limited to 'examples')
| -rw-r--r-- | examples/cavity_2d_gpu.cpp | 1 | ||||
| -rw-r--r-- | examples/poiseulle_particles_channel_2d.cpp | 6 |
2 files changed, 4 insertions, 3 deletions
diff --git a/examples/cavity_2d_gpu.cpp b/examples/cavity_2d_gpu.cpp index b24ca38..19563e2 100644 --- a/examples/cavity_2d_gpu.cpp +++ b/examples/cavity_2d_gpu.cpp @@ -199,7 +199,6 @@ void lbm_step( */ component<sch::T, sch::D2Q9, cmpt::BGK> coll{0.59}; component<sch::T, sch::D2Q9, cmpt::BounceBack> bb; - component<sch::T, sch::D2Q9, cmpt::MovingWall> bb_lid; bb_lid.lid_vel = {0.1,0.0}; diff --git a/examples/poiseulle_particles_channel_2d.cpp b/examples/poiseulle_particles_channel_2d.cpp index e495043..fed78b2 100644 --- a/examples/poiseulle_particles_channel_2d.cpp +++ b/examples/poiseulle_particles_channel_2d.cpp @@ -643,7 +643,9 @@ int main(int argc, char** argv){ compute_rho_u<sch::T,sch::D2Q9>(dfs,rho,vel); rho = rho * saw::data<sch::T>{dfi::cs2}; - force = + for(uint64_t d = 0u; d < sch::D2Q9::D; ++d){ + force.at({d}) = cell.template get<"force">()({d}); + } part_mask.set(0u); }, {{0u,0u}}, meta); @@ -683,7 +685,7 @@ int main(int argc, char** argv){ // Cast down to get lower corner. // Before casting shift by 0.5 for closest pick saw::data<sch::FixedArray<sch::UInt64,2u>> p_cell_pos {{ - static_cast<uint64_t>(p_pos_lie.at(saw::data<T>{{0u}}).get()+0.5), + static_cast<uint64_t>(p_pos_lie.at({{0u}}).get()+0.5), static_cast<uint64_t>(p_pos_lie.at({{1u}}).get()+0.5) }}; |