On Oct 6, 2011, at 9:26 AM, Pauli Virtanen wrote: > I think adding a new function called `diags` could be justified here. > Numpy already has a function called `diag`, and with the new semantics > one would have > > diags([a, b, c], [i, j, k]) == diag(a, i) + diag(b, j) + diag(c, k) > > which hangs together nicely. +1