def METHOD_NAME(self):
    # Hermitian positive definite banded matrix `a`
    a = array([[4.0, 1.0, 0.0, 0.0],
               [1.0, 4.0, 0.5, 0.0],
               [0.0, 0.5, 4.0, -0.2j],
               [0.0, 0.0, 0.2j, 4.0]])
    # Banded storage form of `a`.
    ab = array([[-1.0, 1.0, 0.5, -0.2j],
                [4.0, 4.0, 4.0, 4.0]])
    c = cholesky_banded(ab, lower=False)
    ufac = zeros_like(a)
    ufac[list(range(4)), list(range(4))] = c[-1]
    ufac[(0, 1, 2), (1, 2, 3)] = c[0, 1:]
    assert_array_almost_equal(a, dot(ufac.conj().T, ufac))
    b = array([0.0, 0.5, 4.0-0.2j, 0.2j + 4.0])
    x = cho_solve_banded((c, False), b)
    assert_array_almost_equal(x, [0.0, 0.0, 1.0, 1.0])