anchor.diff:     INTEGER,    ALLOCATABLE                   :: KEY(:)
      INTEGER(IntKi)                            :: ErrStat2
      CHARACTER(ErrMsgLen)                      :: ErrMsg2
  -   REAL(FEKi) :: normA
  +   REAL(FEKi) :: normA, Om2
  +   REAL(FEKi) :: Omega2(N)  !< Squared eigenvalues
         
      ErrStat = ErrID_None
      ErrMsg  = ''
  
positive.code:    REAL(FEKi),              INTENT(INOUT)    :: M(N, N)       !< Mass matrix 
     LOGICAL,                 INTENT(IN   )    :: bCheckSingularity                  ! If True, the solver will fail if rigid modes are present 
     REAL(FEKi),              INTENT(INOUT)    :: EigVect(N, N) !< Returned Eigenvectors
     REAL(FEKi),              INTENT(INOUT)    :: Omega2(N)      !< Returned Eigenvalues
     INTEGER(IntKi),          INTENT(  OUT)    :: ErrStat       !< Error status of the operation
     CHARACTER(*),            INTENT(  OUT)    :: ErrMsg        !< Error message if ErrStat /= ErrID_None
     ! LOCALS         
  
negative.code:    ! LOCALS         
     REAL(FEKi), ALLOCATABLE                   :: K_FEKi(:,:), M_FEKi(:,:) 
     REAL(FEKi), ALLOCATABLE                   :: EigVect_FEKi(:,:), Omega2_FEKi(:) 
     INTEGER(IntKi)                            :: N, i
     INTEGER(IntKi)                            :: ErrStat2
     CHARACTER(ErrMsgLen)                      :: ErrMsg2
  