Main Page | Class Hierarchy | Class List | File List | Class Members | File Members

fermiqcd_staggered_algorithms.h File Reference

Go to the source code of this file.


mdp_matrix Omega4x4 (mdp_site x)
void mul_Q (staggered_field &psi_out, staggered_field &psi_in, gauge_field &U, coefficients &coeff, int parity=EVENODD)
 Executes current Staggered/Asqtad action.
inversion_stats mul_invQ (staggered_field &psi_out, staggered_field &psi_in, gauge_field &U, coefficients &coeff, mdp_real absolute_precision=staggered_inversion_precision, mdp_real relative_precision=0, int max_steps=2000)
 Executes current Staggered/Asqtad inverter.
mdp_array< mdp_real, 1 > lepage_coefficients (mdp_real plaquette, char type[])
void lepage_improved_links (gauge_field &V, gauge_field &U, mdp_array< mdp_real, 1 > c, int project=FALSE)
void staggered_rephase (gauge_field &U, staggered_field &chi)


void(* default_staggered_action )(staggered_field &, staggered_field &, gauge_field &, coefficients &, int) = StaggeredAsqtadActionFast::mul_Q
 Pointer to current Staggered/Asqtad action.
inversion_stats(* default_staggered_inverter )(staggered_field &, staggered_field &, gauge_field &, coefficients &, mdp_real, mdp_real, int) = BiCGStab::inverter<staggered_field,gauge_field>
 Pointer to current Staggered/Asqtad inverter.

Detailed Description

Massimo Di Pierro <>
Various stuff for staggered fermions

This file is copyrighted by Massimo Di Pierro Read attached license in file fermiqcd_license.pdf This file cannot be distributed without file fermiqcd_license.pdf

Function Documentation

mdp_array<mdp_real,1> lepage_coefficients mdp_real  plaquette,
char  type[]

Takes a plaquette and a type of action and returns a 1D array with weights of paths required to build fat links for the action

See also:

void lepage_improved_links gauge_field V,
gauge_field U,
mdp_array< mdp_real, 1 >  c,
int  project = FALSE

Takes a gauge field U and a set of coefficients as computed by lepage_coefficients() and fills the gauge field V with fat links and Long links


    gauge_field U(lattice,nc);
    gauge_field V(lattice,nc);
    float p=1.0; // the average plaquette
    /// now use V instead of U for staggered actions and inverters
Note that the type of action can be
  • "Full" for full as asqtad
  • "Staple+Naik"
  • "Fat3"
  • "Fat5"
  • "Fat7" Also note that if project==TRUE the fat links are projected back to SU(nc)

Generated on Sun Feb 27 15:12:22 2005 by  doxygen 1.4.1