test_self_adjoint_eigen_solver.py

Go to the documentation of this file.

import eigenpy
 
import numpy as np
 
dim = 100
A = np.random.rand(dim, dim)
A = (A + A.T) * 0.5
 
es = eigenpy.SelfAdjointEigenSolver(A)
 
V = es.eigenvectors()
D = es.eigenvalues()
 
assert eigenpy.is_approx(A.dot(V), V.dot(np.diag(D)), 1e-6)