Hi,
a colleague made me aware of a speed issue with numpy.identity. Since he was
using numpy.diag(numpy.ones(N)) before, he expected identity to be at least as
fast as diag. But that is not the case.
We found that there was a discussion on the list (July, 20th; "My identity" by
Keith Goodman).
On Monday 05 March 2007 08:01, Steffen Loeck wrote:
> Has there been any progress in solving this problem? I get the same error
> message and have no idea how to solve it.
I do not understand those code parts very well but I think the values passed
to the lapack routines must be integer and not l