def get_D(u, v):
    """Compute the rank-1 dictionary associated with u and v

    Parameters
    ----------
    u: array (n_atoms, n_channels)
    v: array (n_atoms, *atom_support)

    Return
    ------
    D: array (n_atoms, n_channels, *atom_support)
    """
    n_atoms, *atom_support = v.shape
    u = u.reshape(*u.shape, *[1 for _ in atom_support])
    v = v.reshape(n_atoms, 1, *atom_support)
    return u*v