grassmann/functions/
lps.rs

1use crate::{Number, core::matrix::Matrix};
2
3
4
5pub fn lps<T: Number>(M: &Matrix<T>, n: usize) -> Matrix<T> {
6    let mut A = Matrix::new(n, n);
7
8    for i in 0..n {
9        for j in 0..n {
10            A[[i,j]] = M[[i,j]];
11        }
12    }
13
14    A
15}