This bugzilla service is closed. All entries have been migrated to https://gitlab.com/libeigen/eigen
View | Details | Raw Unified | Return to bug 1013
Collapse All | Expand All

(-)eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h~ (+6 lines)
Lines 708-713 Link Here
708
    // compute the eigen vectors
708
    // compute the eigen vectors
709
    if(computeEigenvectors)
709
    if(computeEigenvectors)
710
    {
710
    {
711
      if((eivals(1)-eivals(0))<=Eigen::NumTraits<Scalar>::epsilon())
712
      {
713
        eivecs.setIdentity();
714
      }
715
      else {
711
      scaledMat.diagonal().array () -= eivals(1);
716
      scaledMat.diagonal().array () -= eivals(1);
712
      Scalar a2 = numext::abs2(scaledMat(0,0));
717
      Scalar a2 = numext::abs2(scaledMat(0,0));
713
      Scalar c2 = numext::abs2(scaledMat(1,1));
718
      Scalar c2 = numext::abs2(scaledMat(1,1));
Lines 724-729 Link Here
724
      }
729
      }
725
730
726
      eivecs.col(0) << eivecs.col(1).unitOrthogonal();
731
      eivecs.col(0) << eivecs.col(1).unitOrthogonal();
732
      }
727
    }
733
    }
728
    
734
    
729
    // Rescale back to the original size.
735
    // Rescale back to the original size.

Return to bug 1013