Use Abs() wrapper instead of std::fabs to improve readability
Store the self-energies in real matrices while computing the pole
masses.  This saves memory and avoids temporary objects when calculating
things like "M + M * Re(self_energy) + Re(self_energy) * M".
For now use real loop functions only, because it is faster and
we are using the real ones which Ben has implemented.