eos_config.hpp

Go to the documentation of this file.

// -------------------------------------------------------//
//
// SHAMROCK code for hydrodynamics
// Copyright (c) 2021-2026 Timothée David--Cléris <tim.shamrock@proton.me>
// SPDX-License-Identifier: CeCILL Free Software License Agreement v2.1
// Shamrock is licensed under the CeCILL 2.1 License, see LICENSE for more information
//
// -------------------------------------------------------//
 
#pragma once
 
#include "shambackends/sycl.hpp"
#include "shamunits/Constants.hpp"
#include "shamunits/UnitSystem.hpp"
 
namespace shamphys {
 
    template<class Tscal>
    struct EOS_Config_Isothermal {
        Tscal cs;
    };
 
    template<class Tscal>
    struct EOS_Config_Adiabatic {
        Tscal gamma;
    };
 
    template<class Tscal>
    inline bool operator==(
        const EOS_Config_Adiabatic<Tscal> &lhs, const EOS_Config_Adiabatic<Tscal> &rhs) {
        return lhs.gamma == rhs.gamma;
    }
 
    template<class Tscal>
    struct EOS_Config_Polytropic {
        Tscal K;
        Tscal gamma;
    };
 
    template<class Tscal>
    inline bool operator==(
        const EOS_Config_Polytropic<Tscal> &lhs, const EOS_Config_Polytropic<Tscal> &rhs) {
        return (lhs.K == rhs.K) && (lhs.gamma == rhs.gamma);
    }
 
    template<class Tscal>
    struct EOS_Config_Fermi {
        Tscal mu_e; 
    };
 
    template<class Tscal>
    inline bool operator==(const EOS_Config_Fermi<Tscal> &lhs, const EOS_Config_Fermi<Tscal> &rhs) {
        return lhs.mu_e == rhs.mu_e;
    }
 
    template<class Tscal>
    struct EOS_Config_LocallyIsothermal_LP07 {
        Tscal cs0;
 
        Tscal q;
 
        Tscal r0;
    };
 
    template<class Tscal>
    inline bool operator==(
        const EOS_Config_LocallyIsothermal_LP07<Tscal> &lhs,
        const EOS_Config_LocallyIsothermal_LP07<Tscal> &rhs) {
        return (lhs.cs0 == rhs.cs0) && (lhs.q == rhs.q) && (lhs.r0 == rhs.r0);
    }
 
    template<class Tscal>
    struct EOS_Config_LocallyIsothermalDisc_Farris2014 {
        Tscal h_over_r = 0.05;
    };
 
    template<class Tscal>
    inline bool operator==(
        const EOS_Config_LocallyIsothermalDisc_Farris2014<Tscal> &lhs,
        const EOS_Config_LocallyIsothermalDisc_Farris2014<Tscal> &rhs) {
        return (lhs.h_over_r == rhs.h_over_r);
    }
 
    template<class Tscal>
    struct EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014 {
        Tscal cs0;
 
        Tscal q;
 
        Tscal r0;
 
        u32 n_sinks;
    };
 
    template<class Tscal>
    inline bool operator==(
        const EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014<Tscal> &lhs,
        const EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014<Tscal> &rhs) {
        return (lhs.cs0 == rhs.cs0) && (lhs.q == rhs.q) && (lhs.r0 == rhs.r0)
               && (lhs.n_sinks == rhs.n_sinks);
    }
 
} // namespace shamphys