514{
517 auto& points1 = group1.collection->get_points();
518
519 for (size_t i = 0; i < points1.size(); ++i) {
520 const auto& pos_i = points1[i].position;
521 auto& state_i = group1.physics_state[i];
522
525 auto& points2 = group2.collection->get_points();
526
527 size_t start_j = (g1 == g2) ? i + 1 : 0;
528
529 for (size_t j = start_j; j < points2.size(); ++j) {
530 const auto& pos_j = points2[j].position;
531 auto& state_j = group2.physics_state[j];
532
533 glm::vec3 delta = pos_j - pos_i;
534 float distance = glm::length(delta);
535
536 if (distance < m_interaction_radius && distance > 0.001F) {
537 glm::vec3 direction = delta /
distance;
538
540
541 float repulsion_force = 0.0F;
544 }
545
546 glm::vec3 force = direction * (spring_force - repulsion_force);
547
548 state_i.force += force;
549 state_j.force -= force;
550 }
551 }
552 }
553 }
554 }
555}
std::vector< CollectionGroup > m_collections
float m_repulsion_strength
float m_interaction_radius
SpatialField distance(const glm::vec3 &anchor, float radius, DistanceMetric metric=DistanceMetric::EUCLIDEAN)
Normalized distance from an anchor point using the specified metric.