Ejemplo n.º 1
0
        /// Returns a finder to find the next largest eigen value of the receiver's matrix.
        /// @return DhbMatrixAlgebra.LargestEigenvalueFinder
        public LargestEigenvalueFinder NextLargestEigenvalueFinder()
        {
            double    norm = 1.0 / _eigenvector.SecureProduct(_transposedEigenvector);
            DhbVector v1   = _eigenvector * norm;

            return(new LargestEigenvalueFinder(this.DesiredPrecision,
                                               _matrix.SecureProduct(SymmetricMatrix.IdentityMatrix(v1.Dimension)
                                                                     .SecureSubtract(v1.TensorProduct(_transposedEigenvector)))));
        }