diff options
Diffstat (limited to 'modules/core/c++/fplbm.hpp')
| -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 4acdd36..1f38f51 100644 --- a/modules/core/c++/fplbm.hpp +++ b/modules/core/c++/fplbm.hpp @@ -208,7 +208,7 @@ public: 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,sss); + verlet_step_lambda<T,Descriptor::D>(pi,sss); } } }; @@ -286,8 +286,7 @@ public: // vel_s is technically time the density of the particle? - force = ( vel_s * rho - vel * rho ) * flip_por / (one + flip_por / two); - // force = ( vel_s * rho - vel * rho ) * two * flip_por; + force = ( vel_s * rho - vel * rho ) * two * flip_por / (one + flip_por); force_p = force_p - force; }, aabb.template get<"a">(), aabb.template get<"b">()); |