'Evaluate the interpolant.
Parameters
t : float or array_like with shape (n_points,)
Points to evaluate the solution at.
Returns
y : ndarray, shape (n,) or (n, n_points)
Computed values. Shape depends on whether `t` was a scalar or a
1-d array.'

def __call__(self, t):

    t = np.asarray(t)
    if (t.ndim > 1):
        raise ValueError('`t`    must    be    float    or    1-d    array.')
    return self._call_impl(t)
