433{
436 auto& points1 = group1.collection->get_points();
437
438 for (size_t i = 0; i < points1.size(); ++i) {
439 const auto& pos_i = points1[i].position;
440 auto& state_i = group1.physics_state[i];
441
444 auto& points2 = group2.collection->get_points();
445
446 size_t start_j = (g1 == g2) ? i + 1 : 0;
447
448 for (size_t j = start_j; j < points2.size(); ++j) {
449 const auto& pos_j = points2[j].position;
450 auto& state_j = group2.physics_state[j];
451
452 glm::vec3 delta = pos_j - pos_i;
453 float distance = glm::length(delta);
454
455 if (distance < m_interaction_radius && distance > 0.001F) {
456 glm::vec3 direction = delta /
distance;
457
459
460 float repulsion_force = 0.0F;
463 }
464
465 glm::vec3 force = direction * (spring_force - repulsion_force);
466
467 state_i.force += force;
468 state_j.force -= force;
469 }
470 }
471 }
472 }
473 }
474}
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.