apply_interaction_forces()

void MayaFlux::Nodes::Network::ParticleNetwork::apply_interaction_forces ( )

private

Apply interaction forces based on topology.

Definition at line 366 of file ParticleNetwork.cpp.

{
    if (m_topology == Topology::SPATIAL) {
        if (m_neighbor_map_dirty) {
            rebuild_spatial_neighbors();
        }
 
        for (const auto& [idx, neighbors] : m_neighbor_map) {
            auto& particle_a = m_particles[idx];
            glm::vec3 pos_a = particle_a.point->get_position();
 
            for (size_t neighbor_idx : neighbors) {
                auto& particle_b = m_particles[neighbor_idx];
                glm::vec3 pos_b = particle_b.point->get_position();
 
                glm::vec3 delta = pos_b - pos_a;
                float distance = glm::length(delta);
 
                if (distance > 0.001F) {
                    glm::vec3 direction = delta / distance;
 
                    float spring_force = m_spring_stiffness * (distance - m_interaction_radius);
 
                    float repulsion_force = 0.0F;
                    if (distance < m_interaction_radius * 0.5F) {
                        repulsion_force = m_repulsion_strength / (distance * distance);
                    }
 
                    glm::vec3 force = direction * (spring_force - repulsion_force);
                    particle_a.force += force;
                }
            }
        }
    } else if (m_topology == Topology::GRID_2D || m_topology == Topology::GRID_3D || m_topology == Topology::RING || m_topology == Topology::CHAIN) {
        for (const auto& [idx, neighbors] : m_neighbor_map) {
            auto& particle_a = m_particles[idx];
            glm::vec3 pos_a = particle_a.point->get_position();
 
            for (size_t neighbor_idx : neighbors) {
                auto& particle_b = m_particles[neighbor_idx];
                glm::vec3 pos_b = particle_b.point->get_position();
 
                glm::vec3 delta = pos_b - pos_a;
                float distance = glm::length(delta);
 
                if (distance > 0.001F) {
                    glm::vec3 direction = delta / distance;
                    float spring_force = m_spring_stiffness * (distance - 1.0F); // Rest length = 1.0
                    particle_a.force += direction * spring_force;
                }
            }
        }
    }
}

References MayaFlux::Nodes::Network::NodeNetwork::CHAIN, MayaFlux::Nodes::Network::NodeNetwork::GRID_2D, MayaFlux::Nodes::Network::NodeNetwork::GRID_3D, m_interaction_radius, m_neighbor_map, m_neighbor_map_dirty, m_particles, m_repulsion_strength, m_spring_stiffness, MayaFlux::Nodes::Network::NodeNetwork::m_topology, rebuild_spatial_neighbors(), MayaFlux::Nodes::Network::NodeNetwork::RING, and MayaFlux::Nodes::Network::NodeNetwork::SPATIAL.

Referenced by process_batch().

Here is the call graph for this function:

Here is the caller graph for this function: