public static double SpectralRadius_Jacobi(MMatrix A) { Vector Eigenvalues = new Vector(); MMatrix Eigenvectors = new MMatrix(); MMatrix temp = A.Clone(); temp.Jacobi_Cyclic_Method(ref Eigenvalues, ref Eigenvectors); Eigenvalues.Sort(true); return Math.Abs(Eigenvalues[0]); }
public static double SpectralRadius(MMatrix A) { Vector Eigenvalues = new Vector(); if (!A.IsSquared()) throw new MMatrixException("Matrix must be squared."); Eigenvalues = A.QRIterationBasic(GlobalMath.MAXITER).DiagVector(); Eigenvalues.Sort(true); return Math.Abs(Eigenvalues[0]); }