Function reflection_action 

fn reflection_action(
    lower: MatRef<'_, Complex>,
    upper: MatRef<'_, Complex>,
    x: MatRef<'_, Complex>,
    adjoint: bool,
) -> DMatrix<Complex>

Apply I - L U, or its adjoint, using only thin products.