================================
  Test of eigenvalue routines
================================
Real symmetric matrix
A = [[5.65396 2.97436 0.781578 0.86311 2.85983]
 [2.97436 4.04254 1.71215 0.192244 1.26487]
 [0.781578 1.71215 1.71799 -0.0816675 0.307146]
 [0.86311 0.192244 -0.0816675 4.82258 -1.03047]
 [2.85983 1.26487 0.307146 -1.03047 5.0802]]
V = [[-0.281161 0.636489 -0.137997 -0.174451 0.682901]
 [0.593213 -0.172618 0.598425 -0.145853 0.48879]
 [-0.750183 -0.438909 0.457935 -0.0475343 0.180611]
 [0.0336957 -0.330722 -0.362581 -0.87024 0.0265416]
 [0.071676 -0.5129 -0.530729 0.434415 0.511279]]
d = [0.631351 1.55534 3.42858 5.53776 10.1642]
only d = [0.631351 1.55534 3.42858 5.53776 10.1642]
norm(A*V-V*diag(d)) = 1.82327e-15

Real non-symmetric matrix
A = [[-0.516334 -1.63085 0.352447 -0.899812 -0.837553]
 [-1.96168 -1.36627 -1.02175 0.173576 -0.518217]
 [-2.85784 -0.45695 1.11508 -0.970581 0.710777]
 [1.13228 0.35137 0.67581 -0.346552 1.05022]
 [0.887216 0.552351 -0.804102 0.860294 -1.19285]]
V = [[0.0267734-0.393232i 0.0267734+0.393232i 0.361951+0i 0.0136296+0i 0.00582372+0i]
 [-0.0645099+0.366949i -0.0645099-0.366949i 0.798302+0i 0.283819+0i -0.324954+0i]
 [0.769948+0i 0.769948+-0i 0.319263+0i 0.432832+0i 0.550581+0i]
 [-0.028264-0.177314i -0.028264+0.177314i -0.32406+0i 0.38122+0i 0.768105+0i]
 [-0.277587+0.0610262i -0.277587-0.0610262i -0.15736+0i -0.765891+0i 0.0353134+0i]]
d = [0.833361+1.52165i 0.833361-1.52165i -2.63264+0i -1.38711+0i 0.0460897+0i]
only d = [0.833361+1.52165i 0.833361-1.52165i -2.63264+0i -1.38711+0i 0.0460897+0i]
norm(A*V-V*diag(d)) = 2.36814e-15

Complex hermitian matrix
A = [[2.17458+6.1664e-19i 0.790511+2.46594i -0.770475+0.642747i -0.193795+0.79802i -0.155882+0.932205i]
 [0.790511-2.46594i 8.1106-1.9868e-17i 2.82259-3.14959i -2.14686+2.08479i 0.718449-2.7678i]
 [-0.770475-0.642747i 2.82259+3.14959i 7.2953-3.78387e-17i -2.16851-0.472972i 2.62224-2.85262i]
 [-0.193795-0.79802i -2.14686-2.08479i -2.16851+0.472972i 5.83339-3.03035e-17i -1.01055+2.83016i]
 [-0.155882-0.932205i 0.718449+2.7678i 2.62224+2.85262i -1.01055-2.83016i 4.35633-3.63208e-17i]]
V = [[-0.609145-0.499627i 0.016469-0.268674i 0.0730936+0.0894384i -0.343116+0.410956i 0.0807966+0.0193406i]
 [0.425651-0.223503i -0.114141+0.0272727i 0.296254+0.161969i 0.20995+0.493377i 0.0180506-0.594365i]
 [-0.19231-0.254923i -0.180642+0.253951i -0.640168-0.0532483i 0.217686-0.175364i 0.401801-0.385584i]
 [0.180073-0.0627541i 0.0948081-0.383457i -0.529543+0.383122i 0.374558+0.289377i -0.228432+0.322792i]
 [-0.0992416+0i 0.812786+0i 0.181131+0i 0.351731+0i 0.415942+0i]]
d = [0.154865 0.924081 3.96533 6.44143 16.2845]
only d = [0.154865 0.924081 3.96533 6.44143 16.2845]
norm(A*V-V*diag(d)) = 5.4507e-15

Complex non-hermitian matrix
A = [[-0.10322+1.13462i -0.171829+0.696077i -0.446256-0.266404i -0.245456+0.387019i 0.772273-0.484054i]
 [-0.535805+0.546278i -0.266859-0.438964i -0.558091-0.0390739i -0.537047-0.479537i 0.361024+0.924574i]
 [-0.11788-1.15986i 0.662834-0.648018i -0.276296-0.0400101i -0.238933-0.742285i -0.53522-0.0752561i]
 [-0.112724-1.15655i -0.953276+0.552194i 0.310217+1.42559i 0.518141-0.586608i -0.149551-0.433852i]
 [0.272566-0.548203i 0.260513+1.41633i -0.187267-0.253276i -1.18154+0.990668i -0.547819-0.238739i]]
V = [[-0.162457+0.142576i 0.199879+0.216899i 0.11755+0.34397i -0.561358+0.171473i 0.349666+0.189041i]
 [-0.312734-0.272181i -0.360804-0.0284117i -0.0173327+0.0902266i 0.064481-0.00401606i 0.533983+0i]
 [-0.0723421-0.228206i 0.115117-0.0189655i -0.0248118+0.513119i 0.325984+0.0851569i -0.138624-0.412258i]
 [0.786521+0i -0.0289235+0.433014i 0.270784+0.288182i 0.293296-0.257866i -0.364173+0.00406151i]
 [-0.273494+0.175124i 0.761621+0i 0.662681+0i 0.620698+0i 0.428918-0.226004i]]
d = [1.8549-0.361396i -0.943723-1.72866i -1.12506-0.440609i -0.816995+1.27807i 0.354823+1.08289i]
only d = [1.8549-0.361396i -0.943723-1.72866i -1.12506-0.440609i -0.816995+1.27807i 0.354823+1.08289i]
norm(A*V-V*diag(d)) = 1.66924e-15