I have been using some stuff from Distributions.jl and ConjugatePriors.jl to (amongst other things) sample from a NormalWishart distribution. I am getting a whole bunch of problems that trace back to Base.LinAlg.ishermitian detecting a matrix that is definitely Hermitian as being asymmetric. When I investigate it I see something like:
Are there no checks on relative error in ishermitian? This seems like something that would cause problems elsewhere, and be relatively simple to fix.