pub fn linear_Lk(p: usize, k: usize, L: &mut MatrixF64) -> Value