util_compute_8_vectors_s.cc

Go to the documentation of this file.

/*! \file
 * \brief Wrapper code to compute staggered meson correlators.
 *
 * Spectrum calculations
 */
 
#include "handle.h"
//#include "actions/ferm/fermbcs/fermbcs.h"
#include "actions/ferm/fermacts/fermacts_s.h"
#include "meas/hadron/hadron_s.h"
#include "meas/smear/fuzz_smear.h"
#include "meas/sources/srcfil.h"
#include "util/ferm/transf.h"
#include "meas/hadron/vector_meson_s.h"
#include "meas/hadron/stag_propShift_s.h"
 
 
#include "util_compute_quark_prop_s.h"
 
namespace Chroma { 
 
 
  void compute_8_vectors(multi1d<LatticeStaggeredPropagator> & stag_prop,
                         const multi1d<LatticeColorMatrix> & u ,
                         bool gauge_shift, bool sym_shift,
                         XMLWriter & xml_out,
                         int j_decay, int t_length, int t_source,
                         bool binary_meson_dump, std::string binary_name){
 
    Stag_shift_option type_of_shift;
 
    if((gauge_shift) && (sym_shift)){
      type_of_shift=SYM_GAUGE_INVAR;
    }else{
      if((gauge_shift) && (!sym_shift)){
        type_of_shift=GAUGE_INVAR;
      }else{
        if((!gauge_shift) && (sym_shift)){
          type_of_shift=SYM_NON_GAUGE_INVAR;
        }else{
          if((!gauge_shift) && (!sym_shift)){
            type_of_shift=NON_GAUGE_INVAR ;
          }
        }
      }
    }
 
    vector_meson Vector(t_length,  u, type_of_shift) ;
 
 
 
  // ---------- LL ----------
  Vector.compute(stag_prop,j_decay);
 
  if(binary_meson_dump){
 
    std::string tagged_filename_base;
    tagged_filename_base=binary_name+"VT.LL.";
    Vector.binary_dump(t_source,tagged_filename_base);
  }else{
 
    push(xml_out,"Vectors");
    push(xml_out, "Lsink_Lsrc");
    Vector.dump(t_source,xml_out);
    pop(xml_out);
    pop(xml_out); //Vectors
  }
 
  // ------------------------
 
 
  QDPIO::cout << "Computed 8 basic Vector mesons"  << std::endl;
 
}
 
 
 
} // end of namespace