search

etmc / tmLQCD

tmLQCD is a freely available software suite providing a set of tools to be used in lattice QCD simulations. This is mainly a HMC implementation (including PHMC and RHMC) for Wilson, Wilson Clover and Wilson twisted mass fermions and inverter for different versions of the Dirac operator. The code is fully parallelised and ships with optimisations for various modern architectures, such as commodity PC clusters and the Blue Gene family.
http://www.itkp.uni-bonn.de/~urbach/software.html
GNU General Public License v3.0
32 stars 47 forks source link

Check monomials convergence #531

Closed simone-romiti closed 2 years ago

simone-romiti commented 2 years ago

Reply to the issue https://github.com/etmc/tmLQCD/issues/526 .

I added a global parameter BarrierMonomialsConvergence to be specified in the input file as "yes" or "no". The default is "no". The corresponding variable is g_barrier_monomials_convergence (==1 for "yes" and ==0 for "no").

If set to "yes" the convergence is checked for each monomial. When the returned number of iterations equals -1 , the program is aborted.

For reference: this is implemented with a conditional statement for each solve_*() function in tmLQCD/solver/monomial_solve.c which calls fatal_error() if the solver didn't converge and if g_barrier_monomials_convergence=1

kostrzewa commented 2 years ago

Thanks! (note that I've changed the target branch). Could you also add a description of this option to the TeX-documentation?

simone-romiti commented 2 years ago

Done! Please check if the explanation is clear enough.