{4,5,1}};
Matrix L = Matrix.create(dataL);
double D[] = new double[]{1,3,7};
ICholeskyLDUResult ans = CholeskyLDU.decompose(A);
assertNotNull(ans);
Matrix foundL = ans.getL().toMatrix();
assertArrayEquals(L.getElements(),foundL.getElements(), 1e-5);
assertArrayEquals(D, ans.getD().getLeadingDiagonal().toDoubleArray(), 1e-5);
}