def copula_mv_archimedean(u, transform, args=(), axis=-1):
    """generic multivariate Archimedean copula
    """
    phi = transform.evaluate
    phi_inv = transform.inverse
    cdfv = phi_inv(phi(u, *args).sum(axis), *args)
    return cdfv
