This test was wrong, it was computing A * A^T instead of A^T * A.
The latter is correct since our implementation of Cholesky decomposition
stores the triangular matrix in the upper triangle
The test was passing only because the Cholesky decomposition was tested
with diagonal matrix only, than this does not matter.