'Adjust quantiles array so that last two axes labels the components of
each data point.'

def _process_quantiles(self, X, dims):

    X = np.asarray(X, dtype=float)
    if (X.ndim == 2):
        X = X[np.newaxis, :]
    if (X.shape[(-2):] != dims):
        raise ValueError('The    shape    of    array    `X`    is    not    compatible    with    the    distribution    parameters.')
    return X
