more bug fixes in logging

gen_data.py

@@ -76,7 +76,7 @@ def main():
     fourier_coef_group = parser.add_mutually_exclusive_group()
     fourier_coef_group.add_argument("--vac", metavar='float', type=float, nargs='+', default=0., help="Vac, units of Tk")
     fourier_coef_group.add_argument("--fourier_coef", metavar='tuple', type=float, nargs='*', help="Voltage Fourier arguments, units of omega")
-    parser.add_argument("--d", metavar='float', type=float, default=1e9, help="D (UV cutoff), units of Tk")
+    parser.add_argument("--d", metavar='float', type=float, nargs='+', default=1e9, help="D (UV cutoff), units of Tk")
     parser.add_argument("--xL", metavar='float', type=float, nargs='+', default=0.5, help="Asymmetry, 0 < xL < 1")
     parser.add_argument("--compact", metavar='int', type=int, default=0, help="compact FRTRG implementation (0,1, or 2)")
     parser.add_argument("--filename", metavar='file', type=str, default=os.path.join(settings.BASEPATH, settings.FILENAME), help="h5 file to open")

main.py

@@ -266,7 +266,7 @@ def check_convergence(dm, **parameters):
     d_j = table.d[j]
     d_mu = table.d[mu]
     drtol = table.d*table.solver_tol_rel
-    norm = mplcolors.LogNorm(drtol.min(), drtol.max())
+    norm = mplcolors.LogNorm(max(drtol.min(), 0.05), min(drtol.max(), 500))
     c_j = drtol[j]
     c_mu = drtol[mu]
     s_j = 0.05 * table.nmax[j]**2

rtrg.py

@@ -5,9 +5,7 @@ Copyright 2021 Valentin Bruch
 License: MIT
 '''
 
-import sys
 import numpy as np
-from traceback import print_exc
 from numbers import Number
 from scipy.interpolate import interp1d
 import settings
@@ -17,7 +15,7 @@ if settings.USE_CUBLAS:
     try:
         import rtrg_cublas as rtrg_c
     except:
-        settings.logger.warning("failed to load rtrg_cublas, falling back to rtrg_c", file=sys.stderr)
+        settings.logger.warning("failed to load rtrg_cublas, falling back to rtrg_c")
         import rtrg_c
 else:
     import rtrg_c
@@ -53,7 +51,7 @@ class GlobalRGproperties:
         voltage_branches = 0
         vdc = 0 or mu = 0
         '''
-        #settings.logger.info('Created new GlobalRGproperties object', flush=True)
+        settings.logger.debug('Created new GlobalRGproperties object')
         self.__dict__.update(kwargs)
 
     def __getattr__(self, name):
@@ -516,8 +514,7 @@ class RGfunction(RGobj):
                 if settings.VERBOSE:
                     RGfunction.INV_COUNTER[self._values.T.flags.c_contiguous | (self._values.T.flags.f_contiguous << 1)] += 1
             except:
-                settings.logger.warning("padded inversion failed", file=sys.stderr)
-                print_exc()
+                settings.logger.exception("padded inversion failed")
                 res = np.linalg.inv(self._values)
         return RGfunction(self.global_properties, res, self.voltage_shifts, self.symmetry)