qdp_generic_fused_spin_recon_evaluates_wrapper.h

Go to the documentation of this file.

#ifndef QDP_GENERIC_FUSED_SPIN_RECON_EVALUATES_WRAPPER_H
#define QDP_GENERIC_FUSED_SPIN_RECON_EVALUATES_WRAPPER_H

// Threading evaluates wrappers
//
// by Xu Guo, EPCC, 28 August, 2008
 
// user arg for evaluate having order
struct ordered_fused_spin_recon_user_arg{
  const OLattice< SU3Mat >& u;
  const OLattice< HVec >& a;
  OLattice< FVec >& d;
  int base;
  void (*func)(const REAL*, REAL*, unsigned int);
};
 
 
// user func for evaluate having order
inline
void ordered_fused_spin_recon_evaluate_function (int lo, int hi, int myId, ordered_fused_spin_recon_user_arg* arg){
 
  const OLattice< SU3Mat >& u = arg->u;
  const OLattice< HVec >& a = arg->a ;
  OLattice< FVec >& d = arg->d;
  int base = arg->base;
  void (*func)(const REAL*, REAL*, unsigned int) = arg->func; 
 
  int low = lo + base;
  int high = hi + base;
 
  for (int site = low; site < high; site++){
      HVec tmp;
      _inline_mult_su3_mat_vec(u.elem(site).elem(), a.elem(site).elem(0), tmp.elem(0));
      _inline_mult_su3_mat_vec(u.elem(site).elem(), a.elem(site).elem(1), tmp.elem(1));
      
      
      func( (REAL *)&(tmp.elem(0).elem(0).real()), (REAL *)&(d.elem(site).elem(0).elem(0).real()), 1);
  }
 
}
 
 
// user arg for evaluate NOT having order
struct unordered_fused_spin_recon_user_arg{
  const OLattice< SU3Mat >& u;
  const OLattice< HVec >& a;
  OLattice< FVec >& d;
  const int *tab;
  void (*func)(const REAL*, REAL*, unsigned int);
};
 
 
// user func for evaluate NOT having order
inline
void unordered_fused_spin_recon_evaluate_function (int lo, int hi, int myId, unordered_fused_spin_recon_user_arg* arg){
 
  const OLattice< SU3Mat >& u = arg->u;
  const OLattice< HVec >& a = arg->a ;
  OLattice< FVec >& d = arg->d;
  const int *tab = arg->tab;
  void (*func)(const REAL*, REAL*, unsigned int) = arg->func; 
 
  for (int j = lo; j < hi; j++){
 
    int site = tab[j];
      
    HVec tmp;
    _inline_mult_su3_mat_vec(u.elem(site).elem(), a.elem(site).elem(0), tmp.elem(0));
    _inline_mult_su3_mat_vec(u.elem(site).elem(), a.elem(site).elem(1), tmp.elem(1));
      
      
    func( (REAL *)&(tmp.elem(0).elem(0).real()), (REAL *)&(d.elem(site).elem(0).elem(0).real()),1);
  }
 
}
 
 
 
 
 
#endif