velocity_verlet.hpp

Go to the documentation of this file.

#pragma once

 
#include <string_view>
 
#include "orrery/core/particle_data.hpp"
#include "orrery/core/types.hpp"
#include "orrery/integrators/acceleration_field.hpp"
#include "orrery/integrators/integrator.hpp"
 
namespace orrery::integrators {

void velocity_verlet_step(core::ParticleData& data, core::Real timestep, AccelerationField& field);

class VelocityVerlet final : public Integrator {
public:
    void step(core::ParticleData& data, core::Real timestep, AccelerationField& field) override;
 
    [[nodiscard]] std::string_view name() const noexcept override { return "velocity Verlet"; }
 
    [[nodiscard]] int order() const noexcept override { return 2; }
 
    [[nodiscard]] bool is_symplectic() const noexcept override { return true; }
 
    [[nodiscard]] core::Index force_evaluations_per_step() const noexcept override { return 1; }
};
 
} // namespace orrery::integrators