I haven’t had my coffee yet, but it seems like following should return true
. Where did I get turned around?
using LinearAlgebra
A = rand(5, 5)
Λ, Y = eigen(A)
Λ = Diagonal(Λ)
A * Y ≈ Λ * Y
I haven’t had my coffee yet, but it seems like following should return true
. Where did I get turned around?
using LinearAlgebra
A = rand(5, 5)
Λ, Y = eigen(A)
Λ = Diagonal(Λ)
A * Y ≈ Λ * Y
Figured it out: You need to right-multiply by Λ
:
A = rand(5, 5)
Λ, Y = eigen(A)
Λ = Diagonal(Λ)
A * Y ≈ Y * Λ # true
If anyone else makes this mistake and needs a good excuse, you can point out that the following works in Python
import numpy as np
A = np.random.rand(5, 5)
Λ, Y = np.linalg.eig(A)
np.allclose(A @ Y, Λ * Y)