UPDATE: bug "sym_conjg" noticed in the last committment has been removed. Testing
the conformity with bugged version and quality of the results.
recall: the bug was in ahm_real_trap having block G_21 not conjg of block G_12
.AND. calling improperly the symmetric matrix inversion routine, where the
matrix G is not symmetric.
small update with respect to previous one, changed a small bug in print_sc_out allocating gap_ij.
fix a misprint in the error file name, missing a /
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.
Merge branch 'devel'