collisioninfo.hpp

// SPDX-License-Identifier: MIT
 
#pragma once
 
#include <paraconf.h>
 
#include "paraconfpp.hpp"
#include "species_info.hpp"
 
class CollisionInfo
{
    // value of nustar0 = to nustar0_rpeak read in the YAML input file
    double const m_nustar0;
 
    // boolean that is equal to true if inter-species collisions are taken into account
    // read in the YAML input file
    std::int8_t m_collisions_interspecies;
 
    // AD coefficient almost equivalent to the collision frequency
    double m_coeff_AD;
 
public:
    using radial_chunk_type = double;
 
    explicit CollisionInfo(PC_tree_t const& yaml_input_file)
        : m_nustar0 {PCpp_double(yaml_input_file, ".CollisionsInfo.nustar0_rpeak")}
        , m_collisions_interspecies {
                  PCpp_bool(yaml_input_file, ".CollisionsInfo.collisions_interspecies")}
    {
        // nustar fixed to 1. => normalisation of time is equal to 1./nustar
        if (m_nustar0 != 1.0) {
            throw std::invalid_argument("nustar0 must be equal to 1");
        }
 
        m_coeff_AD = Kokkos::sqrt(2.);
    };
 
    ~CollisionInfo() = default;
 
    std::int8_t collisions_interspecies() const
    {
        return m_collisions_interspecies;
    }
 
    double coeff_AD() const
    {
        return m_coeff_AD;
    }
};