On Mon, Jul 20, 2009 at 1:44 PM, Citi, Luca<lciti@essex.ac.uk> wrote: > Just my 2 cents. > It is duplicated code. > But it is only 3 lines. > "identity" does not need to handle rectangular matrices and non-principal diagonals, > therefore it can be reasonably faster (especially for small matrices, I guess). I think that's the best trade-off too.