Function identity 

pub fn identity<T: Zero + One>(n: usize) -> DTensor<T, 2>

Creates an identity matrix of size n x n.

Examples

use mdarray::tensor;
use mdarray_linalg::identity;

let i3 = identity::<f64>(3);
assert_eq!(i3, tensor![[1.,0.,0.],[0.,1.,0.],[0.,0.,1.]]);