Function is_diagonal 

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

Check if a matrix is diagonal