diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-08-25 16:49:40 +0200 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-08-25 16:49:40 +0200 |
| commit | 9a1ee9e0b348f75f5f4f8fb41ea901a0dea02f13 (patch) | |
| tree | eebcd9403ac04a0e99b2f57fbfc895b4344a356b /modules/core | |
| parent | 29ab0dd49362add54fae1791e1d0776157ee79b0 (diff) | |
| download | libs-lbm-9a1ee9e0b348f75f5f4f8fb41ea901a0dea02f13.tar.gz | |
Weird things are happening?
Diffstat (limited to 'modules/core')
| -rw-r--r-- | modules/core/c++/fplbm.hpp | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/modules/core/c++/fplbm.hpp b/modules/core/c++/fplbm.hpp index 840db99..9514c93 100644 --- a/modules/core/c++/fplbm.hpp +++ b/modules/core/c++/fplbm.hpp @@ -93,8 +93,8 @@ public: saw::data<sch::Vector<T,Descriptor::D>> vel = vel_f.at(index); // compute_rho_u<T,Descriptor>(dfs,rho,vel); - vel = vel + force * half / rho; auto eq = equilibrium<T,Descriptor>(rho,vel); + vel = vel + force * half / rho; using dfi = df_info<T,Descriptor>; @@ -206,8 +206,7 @@ public: // vel_s is technically time the density of the particle? - force = ( vel_s * rho - vel * rho ) * rho * flip_por; - // force = ( vel_s * rho - vel * 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; |
