jakevdp/AR_crash.py

## AR_crash.py
import numpy as np
from scipy.sparse.linalg import eigs

N = 6
k = 2

# with this random seed, I get a memory error on the third iteration below
np.random.seed(2301)
A = np.random.random((N,N))

for i in range(3):
    print i
    evals2, evecs2 = eigs(A, k)
	import numpy as np
	from scipy.sparse.linalg import eigs

	N = 6
	k = 2

	# with this random seed, I get a memory error on the third iteration below
	np.random.seed(2301)
	A = np.random.random((N,N))

	for i in range(3):
	print i
	evals2, evecs2 = eigs(A, k)