[feat] add det manip control flags for cthyb

these flags where introduced in triqs detmanip exposed in cthyb

Author: the-hampel

python/solid_dmft/dmft_tools/solvers/cthyb_interface.py

@@ -68,6 +68,10 @@ def __init__(
         # Separately stores all params that go into solve() call of solver
         self.triqs_solver_params = {}
         keys_to_pass = (
+            'det_init_size',
+            'det_n_operations_before_check',
+            'det_precision_warning',
+            'det_precision_error',
             'imag_threshold',
             'length_cycle',
             'max_time',

python/solid_dmft/io_tools/default.toml

@@ -68,6 +68,10 @@ crm_dyson_solver = false
 crm_dlr_wmax = "<none>"
 crm_dlr_eps = "<none>"
 delta_interface = false
+det_init_size = 100
+det_n_operations_before_check = 100
+det_precision_warning = 1e-8
+det_precision_error = 1e-5
 diag_delta = false
 fit_max_moment = "<none>"
 fit_max_n = "<none>"

python/solid_dmft/io_tools/documentation.txt

@@ -231,6 +231,14 @@ crm_dlr_eps : float, default = None
         Set to dlr_eps if crm_dlr_eps = None.
 delta_interface : bool, default = False
         use delta interface in cthyb instead of input G0
+det_init_size : int, default = 100
+        The maximum size of the determinant matrix before a resize
+det_n_operations_before_check : int, default = 100
+        Max number of ops before the test of deviation of the det, M^-1 is performed.
+det_precision_warning : float, default = 1e-8
+        Threshold for determinant precision warnings
+det_precision_error : float, default = 1e-5
+        Threshold for determinant precision error
 diag_delta : bool, default = False
         approximate the hybridization function as diagonal when using the delta interface
 fit_max_moment : int, default = None

test/python/svo_cthyb_basic_tf/dmft_config.toml

@@ -45,5 +45,10 @@ fit_max_moment = 4
 fit_min_w = 4
 fit_max_w = 10
 
+det_init_size = 80
+det_n_operations_before_check = 40
+det_precision_warning = 1e-9
+det_precision_error = 1e-6
+
 [advanced]
 dc_orb_shift = [-0.1, 0.2, 0.1]