Function schur_ndarray 

pub fn schur_ndarray<T>(a: &Array2<T>) -> LapackResult<SchurResult<T>>
where
    T: Field + Clone + Zeroable + Real,

Computes the real Schur decomposition of a square matrix.

A = Q T Q^T where:

• Q is orthogonal (Q^T Q = I)
• T is quasi-upper triangular (upper triangular with possible 2×2 blocks on the diagonal for complex eigenvalue pairs)

Arguments

• a - The input square matrix (n×n)

Returns

Schur decomposition with Q, T, and eigenvalues