blob: fc40be2908fdc3b7b99d8580062eacb8a1e7d920 (
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
41
42
43
44
45
46
47
48
49
50
|
/* 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>
#include <map>
#include <string>
// Math constants
namespace MathConst
{
static constexpr amrex::Real pi = static_cast<amrex::Real>(3.14159265358979323846);
}
// Physical constants. Values are the 2018 CODATA recommended values
// https://physics.nist.gov/cuu/Constants/index.html
//
// New additions here should also be considered for addition to
// `warpx_constants` in WarpXUtil.cpp's `makeParser`, so that they're
// available in parsing and evaluation of PICMI expressions, as well
// as the corresponding Python definitions
namespace PhysConst
{
static constexpr auto c = static_cast<amrex::Real>( 299'792'458. );
static constexpr auto ep0 = static_cast<amrex::Real>( 8.8541878128e-12 );
static constexpr auto mu0 = static_cast<amrex::Real>( 1.25663706212e-06 );
static constexpr auto q_e = static_cast<amrex::Real>( 1.602176634e-19 );
static constexpr auto m_e = static_cast<amrex::Real>( 9.1093837015e-31 );
static constexpr auto m_p = static_cast<amrex::Real>( 1.67262192369e-27 );
static constexpr auto m_u = static_cast<amrex::Real>( 1.66053906660e-27 );
static constexpr auto hbar = static_cast<amrex::Real>( 1.054571817e-34 );
static constexpr auto alpha = static_cast<amrex::Real>( 0.007297352573748943 );//mu0/(4*MathConst::pi)*q_e*q_e*c/hbar;
static constexpr auto r_e = static_cast<amrex::Real>( 2.817940326204929e-15 );//1./(4*MathConst::pi*ep0) * q_e*q_e/(m_e*c*c);
static constexpr double xi = 1.3050122447005176e-52; //(2.*alpha*alpha*ep0*ep0*hbar*hbar*hbar)/(45.*m_e*m_e*m_e*m_e*c*c*c*c*c);
static constexpr auto xi_c2 = static_cast<amrex::Real>( 1.1728865132395492e-35 ); // This should be usable for single precision, though
// very close to smallest number possible (1.2e-38)
static constexpr auto kb = static_cast<amrex::Real>( 1.380649e-23 ); // Boltzmann's constant, J/K (exact)
}
#endif
|
