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