[feat] add print of avg order for ctseg

TRIQS · Jun 10, 2024 · e3a85d1 · e3a85d1
1 parent 9637825
commit e3a85d1
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 0 deletions.
diff --git a/python/solid_dmft/dmft_tools/results_to_archive.py b/python/solid_dmft/dmft_tools/results_to_archive.py
@@ -110,6 +110,7 @@ def _compile_information(sum_k, general_params, solver_params, solvers, map_imp_
                 write_to_h5['F_time_{}'.format(icrsh)] = solvers[icrsh].F_time
             if solver_params[isolvsec]['measure_pert_order']:
                 write_to_h5['pert_order_histo_imp_{}'.format(icrsh)] = solvers[icrsh].perturbation_order_histo
+                write_to_h5['avg_order_imp_{}.format(icrsh)'] = solvers[icrsh].avg_pert_order
             if solver_params[isolvsec]['measure_nnt']:
                 write_to_h5['O_NN_{}'.format(icrsh)] = solvers[icrsh].triqs_solver.results.nn_tau
             if solver_params[isolvsec]['measure_statehist']:

diff --git a/python/solid_dmft/dmft_tools/solver.py b/python/solid_dmft/dmft_tools/solver.py
@@ -1382,6 +1382,7 @@ def set_Gs_from_G_l():
         # first print average sign
         if mpi.is_master_node():
             print('\nAverage sign: {}'.format(self.triqs_solver.results.sign))
+
         # get Delta_time from solver
         self.Delta_time << self.triqs_solver.Delta_tau
 
@@ -1511,5 +1512,9 @@ def set_Gs_from_G_l():
 
         if self.solver_params['measure_pert_order']:
             self.perturbation_order_histo  = self.triqs_solver.results.perturbation_order_histo_Delta
+            bin_vec = np.arange(0, self.perturbation_order_histo.data.shape[0])
+            self.avg_pert_order = np.sum(bin_vec * self.perturbation_order_histo.data[:])
+            if mpi.is_master_node():
+                print(f'Average perturbation order: {self.avg_pert_order}')
 
         return