Function solve_unit_lower_triangular_in_place_with_conj 

pub fn solve_unit_lower_triangular_in_place_with_conj<T: ComplexField, N: Shape, K: Shape>(
    triangular_unit_lower: MatRef<'_, T, N, N, impl Stride, impl Stride>,
    conj_lhs: Conj,
    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