Restored identity(n) function, which just calls diagonal(1.0, n) internally.
Restored identity(n) function, which just calls diagonal(1.0, n) internally.

lawmurray/Birch,lawmurray/Birch,lawmurray/Birch