def testIncompleteCholesky(self): numpy.random.seed(21) A = numpy.random.rand(5, 5) B = A.T.dot(A) k = 4 R = Util.incompleteCholesky2(B, k) R2 = numpy.linalg.cholesky(B)