elif ddef == -1:
            cp1 = np.clip(np.ceil(cp1), 1, n - 1)
            cp2 = np.clip(np.floor(cp2), 0, n - 2)
        else:
            raise ValueError('Undefined discretization definition, ddef = {}'.format(ddef))