UPDATE. *symmcong* bug removed. Real-axis codes have been updated to remove the symmconjg bug.
Initially the code was setting the block 21 of the real-space matrix G(w) to the
wrong function:
G_21=-S(-w)
The correct function has been restored, in accord with the math, IPT codes and solver part:
G_21=-S^*(-w).
In any case this was not introducing a real bug as we used symmetric matrices inversion routine
from Lapack. This routine consider only the upper triangle block matrix (G_11,12,22), so that
wrong equality was disregarded.
The symmetric matrix routine enforce the symmetry:
S^*(-w)=S(w).
We keep using this routine. Test with generic matrix routine shows small error propagation
caused by numerical issues in the inversions.
So the situation now is:
correct equality (anyway useless with respect to sym_matrix_inversion)
symmetric matrix inversion.