tAdd dimensionality check for non-rotating grains during momentum integration - Granular.jl - Julia package for granular dynamics simulation
HTML git clone git://src.adamsgaard.dk/Granular.jl
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit 985ca340ec43f3bc4c581503513756d2bd60083e
DIR parent 0045cc69032ba0ebbbff2f1163130567155291e0
HTML Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Mon, 11 Mar 2019 12:50:20 +0100
Add dimensionality check for non-rotating grains during momentum integration
Diffstat:
M src/interaction.jl | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
---
DIR diff --git a/src/interaction.jl b/src/interaction.jl
t@@ -434,8 +434,8 @@ function interactGrains!(simulation::Simulation, i::Int, j::Int, ic::Int)
simulation.grains[i].torque[3] += -force_t*R_ij + M_t
simulation.grains[j].torque[3] += -force_t*R_ij - M_t
else
- simulation.grains[i].force += force_n*n;
- simulation.grains[j].force -= force_n*n;
+ simulation.grains[i].force += vecTo3d(force_n.*n);
+ simulation.grains[j].force -= vecTo3d(force_n.*n);
end
simulation.grains[i].pressure +=