search

andrenarchy / pseudopy

PseudoPy computes and visualizes the pseudospectrum of a matrix
MIT License
24 stars 5 forks source link

PseudoPsectrum of Random Matrix. #9

Open nikosmatsa opened 3 years ago

nikosmatsa commented 3 years ago

why my kernel dies when I put : 

n  = 2000
mu = 0  
sigma = 1
A  = np.random.normal(mu, sigma, (n,n))

pseudo = NonnormalAuto(A, 1e-5, 1)
pseudo.plot([10**k for k in range(-4, 3)], spectrum=eigvals(A))
pyplot.show()