36 const scalar
a = this->
a();
37 const scalar
b = this->
b();
38 const scalar
c = this->
c();
48 const scalar numDiscr = fma(-
a,
c, w) + fma(
b,
b/4, -w);
49 const scalar discr = (
mag(numDiscr) > VSMALL) ? numDiscr : 0;
52 const bool twoReal = discr > 0;
53 const bool twoComplex = discr < 0;
dimensionedScalar sign(const dimensionedScalar &ds)
Roots< 2 > roots() const
Return the roots of the quadratic equation with no particular order.
scalar a() const noexcept
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
Templated storage for the roots of polynomial equations, plus flags to indicate the nature of the roo...
scalar b() const noexcept
scalar c() const noexcept
Container to encapsulate various operations for linear equation of the forms with real coefficients: ...
Roots< 1 > roots() const
Return the real root of the linear equation.