VerletFunctor Class Reference

Implementation of the Strömer-Verlet integration method. Implements the interface IntegrationFunctor. More...

#include <VerletFunctor.h>

Inheritance diagram for VerletFunctor:

Collaboration diagram for VerletFunctor:

Public Member Functions
void	step (std::unique_ptr< ParticleContainer > &particle_container, const std::vector< std::shared_ptr< SimpleForceSource >> &simple_force_sources, const std::vector< std::shared_ptr< PairwiseForceSource >> &pairwise_force_sources, const std::vector< std::shared_ptr< TargettedForceSource >> &targetted_force_sources, double delta_t, double curr_simulation_time) const override
	Performs one step with the respective integration method. More...
Public Member Functions inherited from IntegrationFunctor
virtual	~IntegrationFunctor ()=default
	Virtual default destructor for the interface. More...

Detailed Description

Implementation of the Strömer-Verlet integration method. Implements the interface IntegrationFunctor.

Implements the IntegrationFunctor interface, and therefore updates all particles in the particle_container according to the Verlet integration method.

Definition at line 11 of file VerletFunctor.h.

Member Function Documentation

step()

void VerletFunctor::step ( std::unique_ptr< ParticleContainer > &  particle_container, const std::vector< std::shared_ptr< SimpleForceSource >> &  simple_force_sources, const std::vector< std::shared_ptr< PairwiseForceSource >> &  pairwise_force_sources, const std::vector< std::shared_ptr< TargettedForceSource >> &  targetted_force_sources, double  delta_t, double  curr_simulation_time  ) const

overridevirtual

Performs one step with the respective integration method.

Parameters

particle_container	Container of particles on which the integration step is applied
simple_force_sources	Vector of simple force sources which are used to calculate the new forces
pairwise_force_sources	Vector of pairwise force sources which are used to calculate the new forces
targetted_force_sources	Vector of targetted force sources which are used to calculate the new forces
delta_t	Time step
curr_simulation_time	Current simulation time

Implements IntegrationFunctor.

Definition at line 5 of file VerletFunctor.cpp.

                                                            {
    for (auto& p : *particle_container) {
        // update position
        if (!p.isLocked()) {
            const std::array<double, 3> new_x = p.getX() + delta_t * p.getV() + (delta_t * delta_t / (2 * p.getM())) * p.getF();
            p.setX(new_x);
        }
 
        // reset forces
        p.setOldF(p.getF());
        p.setF({0, 0, 0});
    }
 
    // calculate new forces
    particle_container->prepareForceCalculation();
    particle_container->applySimpleForces(simple_force_sources);
    particle_container->applyPairwiseForces(pairwise_force_sources);
    particle_container->applyTargettedForces(targetted_force_sources, curr_simulation_time);
 
    // update velocity
    for (auto& p : *particle_container) {
        if (!p.isLocked()) {
            const std::array<double, 3> new_v = p.getV() + (delta_t / (2 * p.getM())) * (p.getF() + p.getOldF());
            p.setV(new_v);
        } else {
            p.setV({0, 0, 0});
        }
    }
}

---

The documentation for this class was generated from the following files:

• VerletFunctor.h
• VerletFunctor.cpp