Computational Thinking, Linear Algebra, Python, Jupyter Notebook, Math.
possible to just use the QR algorithm, but ridiculously slow. Arnoldi Iteration We can use the Arnoldi iteration for phase 1 (and the QR algorithm for phase 2). Initialization import numpy as np n = 5 A0 = np.random.rand(n,n) #.astype(np.float64) A = A0 @ A0.T np.set_printoptions(precision=5,