def LDLdecomposition(self):
        if not self.is_square:
            raise NonSquareMatrixError("Matrix must be square.")
        if not self.is_symmetric():
            raise ValueError("Matrix must be symmetric.")
        return self._LDLdecomposition()
Returns the LDL Decomposition (L, D) of matrix A,
dataset/ETHPy150Open sympy/sympy/sympy/matrices/matrices.py/MatrixBase.LDLdecomposition