diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-08-25 10:46:46 +0200 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-08-25 10:46:46 +0200 |
| commit | 7af165927a3ebbf70a84605cc3b1ad7fb886e940 (patch) | |
| tree | 2d4b7cf4271dd2e47f98214a218685aa1ddbbcd9 /modules/core/c++ | |
| parent | 8c3107f5a5a603d4f2e2830e38791298ac1ded6c (diff) | |
| parent | 29ab0dd49362add54fae1791e1d0776157ee79b0 (diff) | |
| download | libs-lbm-7af165927a3ebbf70a84605cc3b1ad7fb886e940.tar.gz | |
Merge branch 'dev'
Diffstat (limited to 'modules/core/c++')
| -rw-r--r-- | modules/core/c++/fplbm.hpp | 6 |
1 files changed, 4 insertions, 2 deletions
diff --git a/modules/core/c++/fplbm.hpp b/modules/core/c++/fplbm.hpp index d130c8d..840db99 100644 --- a/modules/core/c++/fplbm.hpp +++ b/modules/core/c++/fplbm.hpp @@ -206,13 +206,15 @@ public: // vel_s is technically time the density of the particle? - force = ( vel_s * rho - vel * rho ) * (two / three) * flip_por; + force = ( vel_s * rho - vel * rho ) * rho * flip_por; + // force = ( vel_s * rho - vel * rho ) * flip_por; + // force = ( vel_s * rho - vel * rho ) /* (two / three) */ * flip_por; force_p = force_p - force; }, aabb.template get<"a">(), aabb.template get<"b">()); auto& pirb_acc = pirb.template get<"acceleration">(); - pirb_acc = force_p / (two * pg.template get<"total_mass">().at({})); + pirb_acc = force_p / (pg.template get<"total_mass">().at({})); for(saw::data<sch::UInt64> i{0u}; i < sub_steps; ++i){ verlet_step_lambda<T,Descriptor::D>(pi,one / sss); |