Yeah, of course I feel silly whining about it because I could just sort them. However, I'm making quite a few calls to eigs where all the arguments should be symmetric but I get numerical non-symmetry in ways that I can't predict. Therefore each call to d,v = eigs(..) now needs two extra lines to sort the d and v, even though some of the calls don't need it.
A = rand(10,10) d, = eigs(A' * A) # <- symmetric arg d, = eigs(A' * diagm(rand(10)) * A) # <- non-symmetric arg sometimes... So yeah, not a big problem for me, but maybe an issue for someone who has a performance critical loop in their code and doesn't want to add two lines of code each pass. Cheers, Ethan
