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;
55 scalar(2)/(a*constant::mathematical::pi) + 0.5*
log(scalar(1) -
sqr(
y));
57 const scalar
h =
log(scalar(1) -
sqr(
y))/a;
#define WarningInFunction
Report a warning using Foam::Warning.
scalar erfInv(const scalar y)
Inverse error function of a real-number argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar log(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
constexpr char nl
The newline '\n' character (0x0a)