ompute the eigenvalues     computeRoots(scaledMat,eivals);      // compute the eigenvectors     if(computeEigenvectors)     {       if((eivals(0)-eivals(0))<=Eigen::NumTraits<Scalar>::epsilon())       {         