blob: ae6cd808e9009f0175a59a637d96e02484b0ad77 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
|
/* Copyright 2019 Andrew Myers, Luca Fedeli, Maxence Thevenet
* Weiqun Zhang
*
* This file is part of WarpX.
*
* License: BSD-3-Clause-LBNL
*/
#ifndef WARPX_CONST_H_
#define WARPX_CONST_H_
#include <AMReX_REAL.H>
#include <cmath>
// Math constants
namespace MathConst
{
static constexpr amrex::Real pi = 3.14159265358979323846;
}
// Physical constants. Values are the 2018 CODATA recommended values
// https://physics.nist.gov/cuu/Constants/index.html
namespace PhysConst
{
static constexpr amrex::Real c = 299'792'458.;
static constexpr amrex::Real ep0 = 8.8541878128e-12;
static constexpr amrex::Real mu0 = 1.25663706212e-06;
static constexpr amrex::Real q_e = 1.602176634e-19;
static constexpr amrex::Real m_e = 9.1093837015e-31;
static constexpr amrex::Real m_p = 1.67262192369e-27;
static constexpr amrex::Real hbar = 1.054571817e-34;
static constexpr amrex::Real alpha = mu0/(4*MathConst::pi)*q_e*q_e*c/hbar;
static constexpr amrex::Real r_e = 1./(4*MathConst::pi*ep0) * q_e*q_e/(m_e*c*c);
static constexpr double xi = (2.*alpha*alpha*ep0*ep0*hbar*hbar*hbar)/
(45.*m_e*m_e*m_e*m_e*c*c*c*c*c); // SI value is 1.3050122.e-52
static constexpr amrex::Real xi_c2 = xi * c * c; // This should be usable for single precision, though
// very close to smallest number possible: smallest number = 1.2e-38, xi_c2 = 1.1e-35
static constexpr amrex::Real mevpc2_kg = 1.7826619216279e-30; // to convert MeV in kg
}
#endif
|
