diff --git a/src/OpenFOAM/primitives/polynomialEqns/cubicEqn/cubicEqn.C b/src/OpenFOAM/primitives/polynomialEqns/cubicEqn/cubicEqn.C index 9633937f6..5e3b8899f 100644 --- a/src/OpenFOAM/primitives/polynomialEqns/cubicEqn/cubicEqn.C +++ b/src/OpenFOAM/primitives/polynomialEqns/cubicEqn/cubicEqn.C @@ -96,7 +96,7 @@ Foam::Roots<3> Foam::cubicEqn::roots() const if (oneReal) { - const Roots<1> r = linearEqn(- a, b/3).roots(); + const Roots<1> r = linearEqn(a, b/3).roots(); return Roots<3>(r.type(0), r[0]); } else if (twoReal) diff --git a/src/OpenFOAM/primitives/polynomialEqns/quadraticEqn/quadraticEqn.C b/src/OpenFOAM/primitives/polynomialEqns/quadraticEqn/quadraticEqn.C index 0fea8e6e9..cfcfceafa 100644 --- a/src/OpenFOAM/primitives/polynomialEqns/quadraticEqn/quadraticEqn.C +++ b/src/OpenFOAM/primitives/polynomialEqns/quadraticEqn/quadraticEqn.C @@ -71,7 +71,7 @@ Foam::Roots<2> Foam::quadraticEqn::roots() const if (oneReal) { - const Roots<1> r = linearEqn(- a, b/2).roots(); + const Roots<1> r = linearEqn(a, b/2).roots(); return Roots<2>(r, r); } else if (twoReal)