Function tridiag_in_place 

pub fn tridiag_in_place<T: ComplexField>(
    A: MatMut<'_, T>,
    householder: MatMut<'_, T>,
    par: Par,
    stack: &mut MemStack,
    params: Spec<TridiagParams, T>,
)

computes a self-adjoint matrix $A$’s tridiagonalization such that $A = Q T Q^H$

$T$ is a self-adjoint tridiagonal matrix stored in $A$’s diagonal and subdiagonal

$Q$ is a sequence of householder reflections stored in the unit lower triangular half of $A$ (excluding the diagonal), with the householder coefficients being stored in householder