krook_source_adaptive.hpp

// SPDX-License-Identifier: MIT
#pragma once
 
#include "geometry.hpp"
#include "irighthandside.hpp"
 
class KrookSourceAdaptive : public IRightHandSide
{
private:
    RhsType m_type;
    double m_extent;
    double m_stiffness;
    double m_amplitude;
    double m_density;
    double m_temperature;
    DFieldMemX m_mask;
    DFieldMemVx m_ftarget;
 
public:
    KrookSourceAdaptive(
            IdxRangeX const& gridx,
            IdxRangeVx const& gridvx,
            RhsType const type,
            double extent,
            double stiffness,
            double amplitude,
            double density,
            double temperature);
 
    KrookSourceAdaptive(KrookSourceAdaptive&&) = default;
 
    ~KrookSourceAdaptive() override = default;
 
    DFieldSpXVx operator()(DFieldSpXVx allfdistribu, double dt) const override;
 
public:
    void get_amplitudes(DFieldSpX amplitudes, DConstFieldSpXVx allfdistribu) const;
 
    void get_derivative(DFieldSpXVx df, DConstFieldSpXVx f, DConstFieldSpXVx f0) const;
};