On Wed, Aug 12, 2009 at 9:29 AM, Scott Sinclair <scott.sinclair.za@gmail.com
wrote:
2009/8/12 Keith Goodman <kwgoodman@gmail.com>: On Wed, Aug 12, 2009 at 7:24 AM, Keith Goodman<kwgoodman@gmail.com> wrote:
On Wed, Aug 12, 2009 at 1:31 AM, Lars Bittrich<lars.bittrich@googlemail.com> wrote:
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). The presented solution was much faster. Someone wondered if the change was already made in the svn.
Things tend to get lost on the mailing list. The next step would be to file a ticket on the numpy trac. (I've never done that) That would increase the chance of someone important taking a look at it.
Here's the ticket:
A patch against recent SVN trunk is attached to the ticket. Please review...
Already done. Thanks. Chuck