48 template<typename Matrix, typename Roots>
49 inline void computeRoots(const Matrix& m, Roots& roots)  in computeRoots()  argument
56   // eigenvalues are the roots to this equation, all guaranteed to be  in computeRoots()
62   // Construct the parameters used in classifying the roots of the equation  in computeRoots()
63   // and in solving the equation for the roots in closed form.  in computeRoots()
75   // Compute the eigenvalues by solving for the roots of the polynomial.  in computeRoots()
80   roots(0) = c2_over_3 + Scalar(2)*rho*cos_theta;  in computeRoots()
81   roots(1) = c2_over_3 - rho*(cos_theta + s_sqrt3*sin_theta);  in computeRoots()
82   roots(2) = c2_over_3 - rho*(cos_theta - s_sqrt3*sin_theta);  in computeRoots()
85   if (roots(0) >= roots(1))  in computeRoots()
86     std::swap(roots(0),roots(1));  in computeRoots()
87   if (roots(1) >= roots(2))  in computeRoots()
89     std::swap(roots(1),roots(2));  in computeRoots()
90     if (roots(0) >= roots(1))  in computeRoots()
91       std::swap(roots(0),roots(1));  in computeRoots()