The LinearAlgebra.jl package’s
inv() function does not preserve matrix symmetry.
while true R = randn((3,3)) R[2,1] = R[1,2]; R[3,1] = R[1,3]; R[2,3]=R[3,2] @assert R == R' @assert inv(R) == inv(R)' end
I’m about to open an issue for this, but in the meantime: is there a way to deal with this?
Note that the original poster on Slack cannot see your response here on Discourse. Consider transcribing the appropriate answer back to Slack, or pinging the poster here on Discourse so they can follow this thread.
(Original message ) (More Info)