def METHOD_NAME(dist, x):
    x = dist.check_x_dimension(x)
    # Compute pdf of independent marginals
    pdf_val = np.ones((x.shape[0],))
    for ind_m in range(len(self.marginals)):
        if hasattr(self.marginals[ind_m], "pdf"):
            pdf_val *= marginals[ind_m].pdf(x[:, ind_m])
        else:
            pdf_val *= marginals[ind_m].pmf(x[:, ind_m])
    return pdf_val