[Lapack] bug in DGEHRD

[Julien Langou]

We have a bug report at:
http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=2&t=1818
It comes from Numpy guys.

The matrix is very pretty to look at with a lots of zero and a lots of 
structure but DGEHRD does not seem to like, not one little bit. See post.

Note that while
    [V,D]=eig(A);
does not work,
    [Q,H]=hess(A);
    [V,D]=eig(H);
    V = Q*V;
works ...

Ideas?

Best wishes,
Julien.