Function is_identity 

pub fn is_identity<F>(a: &ArrayView2<'_, F>) -> bool
where
    F: Float + NumAssign + Sum + One + Send + Sync + ScalarOperand + 'static,

Check if a matrix is the identity matrix