`linalg.hankel` returns unexpected results

A Hankel matrix is supposed to be a square matrix, but the scipy implementation can sometimes return non-square matrices, which I think is misleading/unexpected:

In [1]: from scipy import linalg

In [2]: linalg.hankel([1, 2, 3], [4, 5, 6, 7, 8, 9])
Out[2]:
array([[1, 2, 3, 5, 6, 7],
       [2, 3, 5, 6, 7, 8],
       [3, 5, 6, 7, 8, 9]])

Note that 4 is not included in the output, which may be unexpected (I know that the documentation says that r[0] is ignored, but the initial description of r is “the last row” of the matrix, so at the very least the documentation could do a better job clarifying). This output seems much more reasonable as it captures all elements from both arrays and returns a square matrix:

array([[1, 2, 3, 4, 5],
       [2, 3, 4, 5, 6],
       [3, 4, 5, 6, 7],
       [4, 5, 6, 7, 8],
       [5, 6, 7, 8, 9]])

[ilayn]

There are a couple of things to unpack here. First is that when given two arguments hence working in the context of nonsymmetric Hankel matrix the result is exactly like as required. It has to be a rectangle matrix otherwise the result is not a Hankel matrix.

The data problem here is that the user supplied arrays are not compatible for forming a Hankel matrix, that is to say, the last element of the c must match the first element of r. Because that is the property of being a Hankel matrix; constant anti-diagonals.

Matlab handles the violation with a warning in scipy historically it is chosen to ignore the first element since it is supposed to be the last element of the column anyways otherwise Hankelness is not preserved.

The point can be emphasized better and I’d appreciate a PR for that but this property of the Hankel matrix can be expected also from the user since the result should be familiar.

[ilayn]

No problem at all. The difference is the way the function is called. Suppose you are trying to form the array you got from the current implementation

array([[1, 2, 3, 5, 6, 7],
       [2, 3, 5, 6, 7, 8],
       [3, 5, 6, 7, 8, 9]])

what should be the input in your proposed syntax?