def test_eig_dyn():
v = 0
for i in xrange(5):
n = 1 + int(mp.rand() * 5)
if mp.rand() > 0.5:
# real
A = 2 * mp.randmatrix(n, n) - 1
if mp.rand() > 0.5:
A *= 10
for x in xrange(n):
for y in xrange(n):
A[x, y] = int(A[x, y])
else:
A = (2 * mp.randmatrix(n, n) -
1) + 1j * (2 * mp.randmatrix(n, n) - 1)
if mp.rand() > 0.5:
A *= 10
for x in xrange(n):
for y in xrange(n):
A[x,
y] = int(mp.re(A[x, y])) + 1j * int(mp.im(A[x, y]))
run_hessenberg(A, verbose=v)
run_schur(A, verbose=v)
run_eig(A, verbose=v)
def test_svd_c_rand():
for i in xrange(5):
full = mp.rand() > 0.5
m = 1 + int(mp.rand() * 10)
n = 1 + int(mp.rand() * 10)
A = (2 * mp.randmatrix(m, n) - 1) + 1j * (2 * mp.randmatrix(m, n) - 1)
if mp.rand() > 0.5:
A *= 10
for x in xrange(m):
for y in xrange(n):
A[x, y] = int(mp.re(A[x, y])) + 1j * int(mp.im(A[x, y]))
run_svd_c(A, full_matrices=full, verbose=False)
def test_svd_r_rand():
for i in xrange(5):
full = mp.rand() > 0.5
m = 1 + int(mp.rand() * 10)
n = 1 + int(mp.rand() * 10)
A = 2 * mp.randmatrix(m, n) - 1
if mp.rand() > 0.5:
A *= 10
for x in xrange(m):
for y in xrange(n):
A[x, y] = int(A[x, y])
run_svd_r(A, full_matrices=full, verbose=False)
def test_svd_c_rand():
for i in xrange(5):
full = mp.rand() > 0.5
m = 1 + int(mp.rand() * 10)
n = 1 + int(mp.rand() * 10)
A = (2 * mp.randmatrix(m, n) - 1) + 1j * (2 * mp.randmatrix(m, n) - 1)
if mp.rand() > 0.5:
A *= 10
for x in xrange(m):
for y in xrange(n):
A[x,y]=int(mp.re(A[x,y])) + 1j * int(mp.im(A[x,y]))
run_svd_c(A, full_matrices=full, verbose=False)
def test_svd_r_rand():
for i in xrange(5):
full = mp.rand() > 0.5
m = 1 + int(mp.rand() * 10)
n = 1 + int(mp.rand() * 10)
A = 2 * mp.randmatrix(m, n) - 1
if mp.rand() > 0.5:
A *= 10
for x in xrange(m):
for y in xrange(n):
A[x,y]=int(A[x,y])
run_svd_r(A, full_matrices = full, verbose = False)
def irandmatrix(n, range=10):
"""
random matrix with integer entries
"""
A = mp.matrix(n, n)
for i in xrange(n):
for j in xrange(n):
A[i, j] = int((2 * mp.rand() - 1) * range)
return A
def irandmatrix(n, range = 10):
"""
random matrix with integer entries
"""
A = mp.matrix(n, n)
for i in xrange(n):
for j in xrange(n):
A[i,j]=int( (2 * mp.rand() - 1) * range)
return A
