integrator.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"
 
namespace orrery::integrators {

inline void refresh_accelerations(core::ParticleData& data, AccelerationField& field) {
    field.evaluate(data.positions(), data.masses(), data.accelerations());
}

class Integrator {
public:
    virtual ~Integrator() = default;

    virtual void step(core::ParticleData& data, core::Real timestep, AccelerationField& field) = 0;

    [[nodiscard]] virtual std::string_view name() const noexcept = 0;

    [[nodiscard]] virtual int order() const noexcept = 0;

    [[nodiscard]] virtual bool is_symplectic() const noexcept = 0;

    [[nodiscard]] virtual core::Index force_evaluations_per_step() const noexcept = 0;
 
protected:
    // As on `AccelerationField`, and for the same reasons.
    Integrator() = default;
    Integrator(const Integrator&) = default;
    Integrator(Integrator&&) = default;
    Integrator& operator=(const Integrator&) = default;
    Integrator& operator=(Integrator&&) = default;
};
 
} // namespace orrery::integrators