force_matsubara_hermitian() [2/2]

template<class GG >

void cntr::force_matsubara_hermitian

(

GG &

G

)

Force the Matsubara component of a two-time contour object to be a hermitian matrix at each \( \tau \).

Purpose

Force the Matsubara component of herm_matrix to be a two-time contour object at each \( \tau \) by \( \hat{G}^M(\tau)=\frac{1}{2} \bigl\{\hat{G}^M(\tau)+[\hat{G}^{M}(\tau) ]^{\dagger} \bigl\}.\)

Parameters

G	[GG] a two-time contour object to be modified

Definition at line 2100 of file cntr_utilities_impl.hpp.

                                      {
    int ntau = G.ntau(), m, s1 = G.size1(), p1, p2;
    cdmatrix Gmat,Gmat_herm;
    for  (m = 0; m <= ntau; m++) {
      G.get_mat(m,Gmat);
      Gmat_herm = 0.5*(Gmat + Gmat.adjoint());
      G.set_mat(m,Gmat_herm);
    }
}