bumpontailequilibrium.hpp

// SPDX-License-Identifier: MIT
 
#pragma once
 
#include <paraconf.h>
 
#include "geometry.hpp"
#include "iequilibrium.hpp"
#include "paraconfpp.hpp"
#include "species_info.hpp"
 
class BumpontailEquilibrium : public IEquilibrium
{
    host_t<DFieldMemSp> m_epsilon_bot;
 
    host_t<DFieldMemSp> m_temperature_bot;
 
    host_t<DFieldMemSp> m_mean_velocity_bot;
 
public:
    void compute_twomaxwellian(
            DFieldVx fMaxwellian,
            double epsilon_bot,
            double temperature_bot,
            double mean_velocity_bot) const;
    BumpontailEquilibrium(
            host_t<DFieldMemSp> epsilon_bot,
            host_t<DFieldMemSp> temperature_bot,
            host_t<DFieldMemSp> mean_velocity_bot);
 
    ~BumpontailEquilibrium() override = default;
 
    static BumpontailEquilibrium init_from_input(
            IdxRangeSp idx_range_kinsp,
            PC_tree_t const& yaml_input_file);
 
    DFieldSpVx operator()(DFieldSpVx allfequilibrium) const override;
 
    host_t<DConstFieldSp> epsilon_bot() const
    {
        return get_const_field(m_epsilon_bot);
    }
 
    host_t<DConstFieldSp> temperature_bot() const
    {
        return get_const_field(m_temperature_bot);
    }
 
    host_t<DConstFieldSp> mean_velocity_bot() const
    {
        return get_const_field(m_mean_velocity_bot);
    }
};