Go to the documentation of this file.
37 if (
mag(
y) >= scalar(1))
40 <<
"The domain of inverse error function argument "
41 <<
"(i.e. y) should be limited to (-1, 1):" <<
nl
45 return std::numeric_limits<scalar>::infinity();
50 constexpr scalar a = 0.147;
55 const scalar
h =
log(scalar(1) -
sqr(
y))/a;
Ostream & endl(Ostream &os)
Add newline and flush stream.
const dimensionedScalar h
Planck constant.
scalar erfInv(const scalar y)
Inverse error function of a real-number argument.
dimensionedScalar log(const dimensionedScalar &ds)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
constexpr scalar pi(M_PI)
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label k
Boltzmann constant.
#define WarningInFunction
Report a warning using Foam::Warning.