Skip to content
Fetching contributors…
Cannot retrieve contributors at this time
13 lines (9 sloc) 212 Bytes
load("../extras/arpack.jl")
# arpack
n = 10
a = rand(n,n)
asym = a+a'+n*eye(n)
(d,v) = eigs(asym, 3)
@assert sum(asym*v[:,1]-d[1]*v[:,1]) < 1e-8
(d,v) = eigs(a,3)
@assert abs(sum(a*v[:,2]-d[2]*v[:,2])) < 1e-8
Something went wrong with that request. Please try again.