def METHOD_NAME(self):
    n = 20
    for k in range(2):
        m = random([n, n])
        for i in range(n):
            m[i, i] = 20*(.1+m[i, i])
        a = dot(transpose(m), m)
        c = cholesky(a)
        a1 = dot(transpose(c), c)
        assert_array_almost_equal(a, a1)
        c = transpose(c)
        a = dot(c, transpose(c))
        assert_array_almost_equal(cholesky(a, lower=1), c)