I think that the name of this function can be a bit misleading: it is not garanteed to return a f-contiguous array: X may be neither C nor F contiguous.

You probably need tests with such an array, by the way:

<pre>
X = np.ones((10, 10))[::2]
</pre>