kinetic_source.hpp

// SPDX-License-Identifier: MIT
#pragma once
 
#include <cmath>
 
#include "geometry.hpp"
#include "irighthandside.hpp"
 
class KineticSource : public IRightHandSide
{
private:
    double m_amplitude;
    double m_density;
    double m_energy;
    double m_temperature;
    DFieldMemX m_spatial_extent;
    DFieldMemVx m_velocity_shape;
 
public:
    KineticSource(
            IdxRangeX const& gridx,
            IdxRangeVx const& gridv,
            double extent,
            double stiffness,
            double amplitude,
            double density,
            double energy,
            double temperature);
 
    ~KineticSource() override = default;
 
    DFieldSpXVx operator()(DFieldSpXVx allfdistribu, double dt) const override;
};