39 if (
mag(
y) >= scalar(1))
42 <<
"The domain of inverse error function argument " 43 <<
"(i.e. y) should be limited to (-1, 1):" <<
nl 47 return std::numeric_limits<scalar>::infinity();
52 constexpr scalar a = 0.147;
57 const scalar
h =
log(scalar(1) -
sqr(
y))/a;
dimensionedScalar log(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
constexpr char nl
The newline '\n' character (0x0a)
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
label k
Boltzmann constant.
scalar erfInv(const scalar y)
Inverse error function of a real-number argument.
constexpr scalar pi(M_PI)
const dimensionedScalar h
Planck constant.
#define WarningInFunction
Report a warning using Foam::Warning.