constant.H

Go to the documentation of this file.

/* Copyright 2019-2022 Andrew Myers, Luca Fedeli, Maxence Thevenet,
 *                     Weiqun Zhang, Axel Huebl
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef ABLASTR_CONSTANT_H_
#define ABLASTR_CONSTANT_H_
 
#include <AMReX_REAL.H>
 

namespace ablastr::constant
{
    namespace math
    {
        template<typename T>
        constexpr auto pi_v = static_cast<T>(3.14159265358979323846);

        template<typename T>
        constexpr auto tau_v = static_cast<T>(2.0 * pi_v<double>);

        constexpr auto pi = pi_v<amrex::Real>;

        constexpr auto tau = tau_v<amrex::Real>;

        template<typename T>
        constexpr auto degree2rad_v = pi_v<T> / static_cast<T>(180);

        constexpr auto degree2rad = degree2rad_v<amrex::Real>;
 
    } // namespace math

    namespace SI
    {

        template<typename T>
        constexpr auto c_v = static_cast<T>(299'792'458);

        template<typename T>
        constexpr auto epsilon_0_v = static_cast<T>(8.8541878188e-12);

        template<typename T>
        constexpr auto mu0_v = static_cast<T>(1.2566370612685e-06);

        template<typename T>
        constexpr auto q_e_v = static_cast<T>(1.602176634e-19);

        template<typename T>
        constexpr auto m_e_v = static_cast<T>(9.1093837139e-31);

        template<typename T>
        constexpr auto m_p_v = static_cast<T>(1.67262192595e-27);

        template<typename T>
        constexpr auto m_u_v = static_cast<T>(1.66053906892e-27);

        template<typename T>
        constexpr auto hbar_v = static_cast<T>(1.0545718176461565e-34);

        template<typename T>
        constexpr auto alpha_v = static_cast<T>(0.0072973525643330135);

        template<typename T>
        constexpr auto r_e_v = static_cast<T>(2.8179403205e-15);

        template<typename T>
        constexpr auto xi_c2_v = static_cast<T>(1.1728865075613984e-35);

        template<typename T>
        constexpr auto schwinger_field_v = static_cast<T>(
            (m_e_v<double>*m_e_v<double>)*(c_v<double>*c_v<double>*c_v<double>)/(hbar_v<double>*q_e_v<double>));

        template<typename T>
        constexpr auto kb_v = static_cast<T>(1.380649e-23);

        template<typename T>
        constexpr auto eV_v = q_e_v<T>;

        template<typename T>
        constexpr auto MeV_v = static_cast<T>(q_e_v<double> * 1e6);

        template<typename T>
        constexpr auto eV_inv_c_v = static_cast<T>(eV_v<double> / c_v<double>);

        template<typename T>
        constexpr auto MeV_inv_c_v = static_cast<T>(MeV_v<double>  / c_v<double>);

        template<typename T>
        constexpr auto eV_inv_c2_v = static_cast<T>(eV_v<double> / (c_v<double> * c_v<double>));

        template<typename T>
        constexpr auto MeV_inv_c2_v = static_cast<T>(MeV_v<double> / (c_v<double> * c_v<double>));

        template<typename T>
        constexpr auto sqrt_4_pi_fine_structure_v = static_cast<T>(0.3028221207690299);

        template<typename T>
        constexpr auto c2_v = static_cast<T>(c_v<double>*c_v<double>);

        template<typename T>
        constexpr auto inv_c_v = static_cast<T>(1.0/c_v<double>);

        template<typename T>
        constexpr auto inv_c2_v = static_cast<T>(1.0/c2_v<double>);

        template<typename T>
        constexpr auto reduced_compton_wavelength_v = static_cast<T>(hbar_v<double> / (m_e_v<double> * c_v<double>));

        constexpr auto c = c_v<amrex::Real>;

        constexpr auto epsilon_0 = epsilon_0_v<amrex::Real>;

        constexpr auto mu0 = mu0_v<amrex::Real>;

        constexpr auto q_e = q_e_v<amrex::Real>;

        constexpr auto m_e = m_e_v<amrex::Real>;

        constexpr auto m_p = m_p_v<amrex::Real>;

        constexpr auto m_u = m_u_v<amrex::Real>;

        constexpr auto hbar = hbar_v<amrex::Real>;

        constexpr auto alpha = alpha_v<amrex::Real>;

        constexpr auto r_e = r_e_v<amrex::Real>;

        constexpr double xi = 1.3050122383827227e-52;

        constexpr auto xi_c2 = xi_c2_v<amrex::Real>;

        template<typename T>
        constexpr auto schwinger_field = schwinger_field_v<amrex::Real>;

        constexpr auto kb = kb_v<amrex::Real>;

        constexpr auto eV = eV_v<amrex::Real>;

        constexpr auto MeV = MeV_v<amrex::Real>;

        constexpr auto eV_inv_c = eV_inv_c_v<amrex::Real>;

        constexpr auto MeV_inv_c = MeV_inv_c_v<amrex::Real>;

        constexpr auto eV_inv_c2 = eV_inv_c2_v<amrex::Real>;

        constexpr auto MeV_inv_c2 = MeV_inv_c2_v<amrex::Real>;

        constexpr auto c2 = c2_v<amrex::Real>;

        constexpr auto inv_c = inv_c_v<amrex::Real>;

        constexpr auto inv_c2 = inv_c2_v<amrex::Real>;

        constexpr auto reduced_compton_wavelength = reduced_compton_wavelength_v<amrex::Real>;
 
    } // namespace SI
 
} // namespace ablastr::constant
 
#endif // ABLASTR_CONSTANT_H_