Set tail of Sigma_iw at all zero

The high-frequency tail of Sigma will be set properly, if both G and G0 have correct tail. At present, the tail of G is not computed. This should be improved.

python/hubbard_one_solver.py

@@ -17,6 +17,7 @@ def __init__(self, beta, gf_struct, n_iw, spin_orbit=False, verbose=True):
         self.beta = beta
         self.gf_struct = gf_struct
         self.n_iw = n_iw
+        self.spin_orbit = spin_orbit
         self.verbose = verbose
 
         if self.verbose and mpi.is_master_node():
@@ -81,6 +82,14 @@ def solve(self, H_int, E_levels, const_of_motion=None, density_matrix_cutoff=1e-
         self.Sigma_iw << self.G0_iw - inverse(self.G_iw)
         self.G0_iw.invert()
 
+        # *********************************************************************
+        # TODO: tail
+        # set tail of Sigma at all zero in the meantime, because tail of G is
+        #  not computed in pomerol solver. This part should be improved.
+        # *********************************************************************
+        for s, sig in self.Sigma_iw:
+            sig.tail.zero()
+
     def __copy_E_levels(self, E_levels):
         """Copy E_levels after checking the data structure"""