Dear ROOT team,
I just stumbled on a feature of TMatrixDSym::Similarity() that I didn’t expect…
it turns out that this method changes the object itself and returns it, rather than creating a new matrix. There is nothing in the documentation clarifying this. It would be good to be explicit about this.