def METHOD_NAME(w, J, sf):
    """
    Inverse 3-D Discrete Wavelet Transform
    Parameters
    ----------
    w : cell array
        wavelet coefficient
    J : int
        number of stages
    sf : 2D ndarray
        synthesis filters
    Returns
    -------
    y : 3D ndarray
        output array
    """
    y = w[J]
    for k in range(J)[::-1]:
        y = sfb3D(y, w[k], sf, sf, sf)
    return y