diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-06-01 17:21:44 +0200 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-06-01 17:21:44 +0200 |
| commit | 7fd9bfd5946472230a3b74c52f88e19c15741faf (patch) | |
| tree | c528bd796f0de524921d06b55d383943eebead08 /lib/core/c++/collision.hpp | |
| parent | 2dd7c95a111a930e8e23140ab3fec074e7de4c8c (diff) | |
| download | libs-lbm-7fd9bfd5946472230a3b74c52f88e19c15741faf.tar.gz | |
I seem to have no clue what I'm doing
Diffstat (limited to 'lib/core/c++/collision.hpp')
| -rw-r--r-- | lib/core/c++/collision.hpp | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/lib/core/c++/collision.hpp b/lib/core/c++/collision.hpp index 9c76c1a..023f61f 100644 --- a/lib/core/c++/collision.hpp +++ b/lib/core/c++/collision.hpp @@ -146,7 +146,8 @@ public: saw::data<sch::Scalar<T>> half; half.at({}).set(0.5); - saw::data<sch::Vector<T,Descriptor::D>> vel = vel_f.at(index) + total_force * ( half / rho ); + auto& vel = vel_f.at(index); + vel = vel + total_force * ( half / rho ); compute_rho_u<T,Descriptor>(dfs_old_f.at(index),rho,vel); auto eq = equilibrium<T,Descriptor>(rho,vel); |