While working on something else, I realized that linspace is not handling requests for returning the sampling spacing consistently:

--

>>> np.linspace(0, 1, 3, retstep=True)

(array([ 0. , 0.5, 1. ]), 0.5)

>>> np.linspace(0, 1, 1, retstep=True)

array([ 0.])

>>> np.linspace(0, 1, 0, retstep=True)

array([], dtype=float64)

Basically, retstep is ignored if the number of samples is 0 or 1. One could argue that it makes sense, because those sequences do not have a spacing defined. But at the very least it should be documented as doing so, and the following inconsistency removed:

>>> np.linspace(0, 1, 1, endpoint=True, retstep=True)

array([ 0.])

>>> np.linspace(0, 1, 1, endpoint=False, retstep=True)

(array([ 0.]), 1.0)

I am personally inclined to think that if a step is requested, then a step should be returned, and if it cannot be calculated in a reasonable manner, then a placeholder such as None, nan, 0 or stop - start should be returned.

What does the collective wisdom think is the best approach for this?

Jaime

(\__/)

( O.o)

( > <) Este es Conejo. Copia a Conejo en tu firma y ayúdale en sus planes de dominación mundial.

( O.o)

( > <) Este es Conejo. Copia a Conejo en tu firma y ayúdale en sus planes de dominación mundial.