pub fn solve_in_place<I: Index, T: ComplexField, C: Conjugate<Canonical = T>>(
L: MatRef<'_, C>,
diagonal: DiagRef<'_, C>,
subdiagonal: DiagRef<'_, C>,
perm: PermRef<'_, I>,
rhs: MatMut<'_, T>,
par: Par,
stack: &mut MemStack,
)