pub fn L_decomp(l: &mut MatrixF64, tau: &mut VectorF64) -> Value