Function solve_unit_lower_triangular_in_place 

pub fn solve_unit_lower_triangular_in_place<T: ComplexField, LhsT: Conjugate<Canonical = T>, N: Shape, K: Shape>(
    triangular_unit_lower: MatRef<'_, LhsT, N, N, impl Stride, impl Stride>,
    rhs: MatMut<'_, T, N, K, impl Stride, impl Stride>,
    par: Par,
)

solves $L x = b$, replacing the diagonal of $L$ with ones, and implicitly conjugating $L$ if needed, and stores the result in rhs