smaller interaction radius for cylinders

BioDynaMo · Jan 9, 2024 · f673903 · f673903
1 parent 60e9619
commit f673903
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/core/operation/mechanical_forces_op.h b/src/core/operation/mechanical_forces_op.h
@@ -93,6 +93,10 @@ class MechanicalForcesOp : public AgentOperationImpl {
       last_time_run_[tid] = current_time;
     }
 
+    if(agent->GetShape() == Shape::kCylinder) {
+      squared_radius_ = 0.0001;
+    }
+
     const auto& displacement =
         agent->CalculateDisplacement(force_, squared_radius_, delta_time_[tid]);
     agent->ApplyDisplacement(displacement);