Checking the symmetry of a matrix

Are you aware of LinearAlgebra.issymmetric()?

afaik issymmetric() doesn’t have a way of specifying a tolerance though.

There is the build-in function isapprox, which allows you to specify absolute and relative tolerances. You could just broadcast this over the full array (A - A') and reduce with all, e.g you could define:.

isapproxsymmetric(A) = all(isapprox.(A - A', 0; rtol=1e-10))