maxwellianequilibrium.hpp

// SPDX-License-Identifier: MIT
 
#pragma once
 
#include <paraconf.h>
 
#include "geometry.hpp"
#include "iequilibrium.hpp"
#include "paraconfpp.hpp"
#include "species_info.hpp"
 
class MaxwellianEquilibrium : public IEquilibrium
{
    // equilibrium density of all kinetic species
    host_t<DFieldMemSp> m_density_eq;
 
    // equilibrium temperature of all kinetic species
    host_t<DFieldMemSp> m_temperature_eq;
 
    // equilibrium mean velocity of all kinetic species
    host_t<DFieldMemSp> m_mean_velocity_eq;
 
public:
    MaxwellianEquilibrium(
            host_t<DFieldMemSp> density_eq,
            host_t<DFieldMemSp> temperature_eq,
            host_t<DFieldMemSp> mean_velocity_eq);
 
    ~MaxwellianEquilibrium() override = default;
 
    static MaxwellianEquilibrium init_from_input(
            IdxRangeSp idx_range_kinsp,
            PC_tree_t const& yaml_input_file);
 
    DFieldSpVxVy operator()(DFieldSpVxVy allfequilibrium) const override;
 
 
    static void compute_maxwellian(
            DFieldVxVy const fMaxwellian,
            double const density,
            double const temperature,
            double const mean_velocity);
 
    host_t<ConstFieldSp<double>> density_eq() const
    {
        return get_const_field(m_density_eq);
    }
 
    host_t<ConstFieldSp<double>> temperature_eq() const
    {
        return get_const_field(m_temperature_eq);
    }
 
    host_t<ConstFieldSp<double>> mean_velocity_eq() const
    {
        return get_const_field(m_mean_velocity_eq);
    }
};