summaryrefslogtreecommitdiff
path: root/modules
diff options
context:
space:
mode:
authorClaudius "keldu" Holeksa <mail@keldu.de>2026-08-25 16:49:40 +0200
committerClaudius "keldu" Holeksa <mail@keldu.de>2026-08-25 16:49:40 +0200
commit9a1ee9e0b348f75f5f4f8fb41ea901a0dea02f13 (patch)
treeeebcd9403ac04a0e99b2f57fbfc895b4344a356b /modules
parent29ab0dd49362add54fae1791e1d0776157ee79b0 (diff)
downloadlibs-lbm-9a1ee9e0b348f75f5f4f8fb41ea901a0dea02f13.tar.gz
Weird things are happening?
Diffstat (limited to 'modules')
-rw-r--r--modules/core/c++/fplbm.hpp5
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;