Handles global updates and parallel tempering. More...

Functions/Subroutines
subroutine	exchange_step (phase, gr, udvr, udvl, stab_nt, udvst, n_exchange_steps, tempering_calc_det)
	Handles parrallel tempering This subroutine is called only if the tempering flag is switched on Requires MPI.
subroutine	global_updates (phase, gr, udvr, udvl, stab_nt, udvst, n_global)
	Carries out N_global global updates as defeined in the Global_move routine of the Hamiltonian module.
complex(kind=kind(0.d0)) function	compute_ratio_global (phase_det_old, phase_det_new, det_vec_old, det_vec_new, nsigma_old, t0_proposal_ratio, ratio)
	This fucntion computes ratio of weights for global moves.
subroutine	compute_fermion_det (phase_det, det_vec, udvl, udvst, stab_nt, storage)
	Computes the fermion determinant from scratch or from the storage.
integer function	npbc_tempering (n, isize)
	Periodic boundary conditions required to defined master and slave for the tempering.
subroutine	global_tempering_setup
	Allocates the Tempering_acceptance (Type Obser_vec) variable that monitors the exchange acceptance.
subroutine	global_tempering_init_obs
	Initializes Tempering_acceptance.
subroutine	global_tempering_obser (toggle)
	Measures Tempering_acceptance.
subroutine	global_tempering_pr
	Prints Tempering_acceptance.

Variables
type(obser_vec), private	tempering_acceptance

Detailed Description

Handles global updates and parallel tempering.

Author: ALF-project

Function/Subroutine Documentation

◆ compute_fermion_det()

subroutine global_mod::compute_fermion_det	(	complex (kind=kind(0.d0)), dimension(:), intent(out)	phase_det,
		real (kind=kind(0.d0)), dimension(:,:), intent(out)	det_vec,
		class(udv_state), dimension(:), intent(inout), allocatable	udvl,
		class(udv_state), dimension(:,:), intent(inout), allocatable	udvst,
		integer, dimension(:), intent(in), allocatable	stab_nt,
		character (len=64), intent(in)	storage )

Computes the fermion determinant from scratch or from the storage.

Author

!> Finite temperature: If the storage is full then computes 
!>                     det( 1 +  VL*DL*UL^{dag} )= Phase_det*e^{ \sum_{n=1}^{ndim} Det_vec(n)}
!>                     Note that  Udvl%U = UL, Udvl%D = DL and  Udvl%V = VL^{dag}
!>                     If storage is empty one first computes  Udvl and in doing so fills up the storage.
!> Zero temperature:   If storage is full then computes  det \Psi_L U(\theta_tot,0) \Psi_R  =  Phase_det*e^{ \sum_{n=1}^{N_part} Det_vec(n) }
!>                     For the projective code the required scales which one can throw away for the Green function are in udvst.
!>                     If the storage is not full, then things will be computed from scratch and intermediate results will be stored in udvst.
!>

Parameters

[in,out]	udvl,udvr	Class(UDV_State) !> If storage=Full : On entry and exit udvl contains the last udv decomposition such that G=(1 + VLDLUL^{dag})^{-1} !> If storage=Empty : On exit udv1 will contain the last udv decomposition such that G=(1 + VLDLUL^{dag})^{-1} !>
[out]	Phase_det	Complex, Dimension(N_FL) !> The phase, per flavor !>
[out]	Det_vec	Real, Dimension(:,N_FL)
[in]	Stab_nt	Integer, Dimension(:) !> List of time slices for stabilization. !>
[in]	Storage	Character !> storage = . Compute the fermion det with the use of the storage. !> storage = . Compute the fermion det from scratch and in doing so, fill the storage. !>

Local variables

Definition at line 780 of file Global_mod.F90.

◆ compute_ratio_global()

complex (kind=kind(0.d0)) function global_mod::compute_ratio_global	(	complex (kind=kind(0.d0)), dimension(:), intent(in), allocatable	phase_det_old,
		complex (kind=kind(0.d0)), dimension(:), intent(in), allocatable	phase_det_new,
		real (kind=kind(0.d0)), dimension(:,:), intent(in), allocatable	det_vec_old,
		real (kind=kind(0.d0)), dimension(:,:), intent(in), allocatable	det_vec_new,
		type (fields), intent(in)	nsigma_old,
		real (kind=kind(0.d0)), intent(in)	t0_proposal_ratio,
		complex (kind=kind(0.d0)), dimension(2), intent(out)	ratio )

This fucntion computes ratio of weights for global moves.

Author: The ALF Project contributors

!>  Ratio_Global = T0(nsigma--> nsigma_old) W(nsigma)/T0(nsigma_old--> nsigma) W(nsigma_old) =
!>                 T0_Proposal_ratio   W(nsigma)/  W(nsigma_old)
!>  On entry the old and new fermion determinants read: Phase_det*e^{ \sum_{n=1}^{ndim} Det_vec(n) }
!>  as obtained from the Compute_Fermion_det routine.
!>  On exit Ratio_Global = Ratio(1)*exp(Ratio(2))
!>

Parameters

[in]	Phase_det_new	Complex, Dimension(N_FL)
[in]	Phase_det_old	Complex, Dimension(N_FL)
[in]	Det_vec_new	Real, Dimension(:,N_FL)
[in]	Det_vec_old	Real, Dimension(:,N_FL)
[in]	T0_proposal_ratio	Real
[in]	nsigma_old	Type(Fields) !> Old configuration. The new configuration is stored in nsigma. nsigma is a globale variable !> contained in the Hamiltonian module. !>

Definition at line 669 of file Global_mod.F90.

◆ exchange_step()

subroutine global_mod::exchange_step	(	complex (kind=kind(0.d0)), intent(inout)	phase,
		complex (kind=kind(0.d0)), dimension(:,:,:), intent(inout), allocatable	gr,
		class(udv_state), dimension(:), intent(inout), allocatable	udvr,
		class(udv_state), dimension(:), intent(inout), allocatable	udvl,
		integer, dimension(:), intent(in), allocatable	stab_nt,
		class(udv_state), dimension(:, :), intent(inout), allocatable	udvst,
		integer, intent(in)	n_exchange_steps,
		logical, intent(in)	tempering_calc_det )

Handles parrallel tempering
This subroutine is called only if the tempering flag is switched on
Requires MPI.

Author: ALF-project

On entry and on exit the left storage is full, the Green function is on time slice 0 and the phase is set. The same holds on exit but with updated quantities.

Parameters

[in,out]	Phase	Complex !> Is updated upon acceptance !>
[in,out]	udvl,udvr	Class(UDV_State) !> On entry udvl contains the last udv decomposition such that G=(1 + VLDLUL^{dag})^{-1} !> udvr is used as working space !>
[in,out]	GR	Complex !> Green function. Is updated upon acceptance. !>
[in,out]	udvst	Class(UDV_State) !> Storage. Is updated with left propagation upon acceptance !>
[in]	Stab_nt	Integer !> List of time slices for stabilization. !>
[in]	N_exchange_steps	Integer !> Number of exchange steps !>
[in]	Tempering_calc_det	Logical !> If true then the fermion determinant is computed in the calulcation of the exchange weight. !>

Local variables.

Store old configuration

Send result of acceptance/rejection decision form master to slave

Definition at line 108 of file Global_mod.F90.

◆ global_tempering_init_obs()

subroutine global_mod::global_tempering_init_obs

Initializes Tempering_acceptance.

Author: The ALF Project contributors

Definition at line 995 of file Global_mod.F90.

◆ global_tempering_obser()

subroutine global_mod::global_tempering_obser ( logical, intent(in) toggle )

Measures Tempering_acceptance.

Author: The ALF Project contributors

Definition at line 1007 of file Global_mod.F90.

◆ global_tempering_pr()

subroutine global_mod::global_tempering_pr

Prints Tempering_acceptance.

Author: The ALF Project contributors

Definition at line 1023 of file Global_mod.F90.

◆ global_tempering_setup()

subroutine global_mod::global_tempering_setup

Allocates the Tempering_acceptance (Type Obser_vec) variable that monitors the exchange acceptance.

Author: The ALF Project contributors

Definition at line 979 of file Global_mod.F90.

◆ global_updates()

subroutine global_mod::global_updates	(	complex (kind=kind(0.d0)), intent(inout)	phase,
		complex (kind=kind(0.d0)), dimension(:,:,:), intent(inout), allocatable	gr,
		class (udv_state), dimension(:), intent(inout), allocatable	udvr,
		class (udv_state), dimension(:), intent(inout), allocatable	udvl,
		integer, dimension(:), intent(in), allocatable	stab_nt,
		class(udv_state), dimension(:,:), intent(inout), allocatable	udvst,
		integer, intent(in)	n_global )

Carries out N_global global updates as defeined in the Global_move routine of the Hamiltonian module.

Author: The ALF Project contributors

On entry and on exit the left storage is full, the Green function is on time slice 0 and the phase is set. The same holds on exit but with updated quantities.

Parameters

[in,out]	Phase	Complex !> Is updated upon acceptance !>
[in,out]	udvl,udvr	Class(UDV_State) !> On entry udvl contains the last udv decomposition such that G=(1 + VLDLUL^{dag})^{-1} !> udvr is used as working space !>
[in,out]	GR	Complex !> Green function. Is updated upon acceptance. !>
[in,out]	udvst	Class(UDV_State) !> Storage. Is updated with left propagation upon acceptance !>
[in]	Stab_nt	Integer !> List of time slices for stabilization. !>
[in]	N_Global	Integer !> Number of global moves that will be caried out !>

Definition at line 450 of file Global_mod.F90.

◆ npbc_tempering()

integer function global_mod::npbc_tempering	(	integer, intent(in)	n,
		integer, intent(in)	isize )

Periodic boundary conditions required to defined master and slave for the tempering.

Author: The ALF Project contributors

Definition at line 961 of file Global_mod.F90.

Variable Documentation

◆ tempering_acceptance

type (obser_vec ), private global_mod::tempering_acceptance

private

Definition at line 59 of file Global_mod.F90.

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ compute_fermion_det()

◆ compute_ratio_global()

◆ exchange_step()

◆ global_tempering_init_obs()

◆ global_tempering_obser()

◆ global_tempering_pr()

◆ global_tempering_setup()

◆ global_updates()

◆ npbc_tempering()

Variable Documentation

◆ tempering_acceptance