tm_sub_Hopping_Matrix.c

/**********************************************************************
 *
 * Copyright (C) 2012 Carsten Urbach
 *
 * This file is based on an implementation of the Dirac operator
 * written by Martin Luescher, modified by Martin Hasenbusch in 2002
 * and modified and extended by Carsten Urbach from 2003-2008
 *
 * This file is part of tmLQCD.
 *
 * tmLQCD is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * tmLQCD is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with tmLQCD.  If not, see <http://www.gnu.org/licenses/>.
 *
 *
 * Hopping_Matrix is the conventional Wilson
 * hopping matrix
 *
 ****************************************************************/

#ifdef HAVE_CONFIG_H
#include <tmlqcd_config.h>
#endif
#include <stdio.h>
#include <stdlib.h>
#ifdef TM_USE_OMP
#include <omp.h>
#endif
#include <complex.h>
#include "global.h"
#include "su3.h"
#ifdef TM_USE_MPI
#include "xchange/xchange.h"
#endif
#include "boundary.h"
#include "init/init_dirac_halfspinor.h"
#include "tm_sub_Hopping_Matrix.h"
#include "update_backward_gauge.h"

// now comes the definition of tm_times_Hopping_Matrix
// which does (a + g5 i b - Hopping_Matrix)
// where cfactor = a + i b
//

#if (defined TM_USE_HALFSPINOR)
#include "operator/halfspinor_hopping.h"

void tm_sub_Hopping_Matrix(const int ieo, spinor* const l, spinor* const p, spinor* const k,
                           complex double const cfactor) {
#ifdef TM_USE_GAUGE_COPY
  if (g_update_gauge_copy) {
    update_backward_gauge(g_gauge_field);
  }
#endif

#ifdef TM_USE_OMP
#pragma omp parallel
  {
    su3* restrict u0 ALIGN;
#endif

#define _TM_SUB_HOP
    spinor* pn;
#include "operator/halfspinor_body.inc"
#undef _TM_SUB_HOP
#ifdef TM_USE_OMP
  } /* OpenMP closing brace */
#endif
  return;
}

#elif (!defined _NO_COMM && !defined TM_USE_HALFSPINOR)
#include "hopping.h"

void tm_sub_Hopping_Matrix(const int ieo, spinor* const l, spinor* p, spinor* const k,
                           complex double const cfactor) {
#ifdef TM_USE_GAUGE_COPY
  if (g_update_gauge_copy) {
    update_backward_gauge(g_gauge_field);
  }
#endif

#if (defined TM_USE_MPI)
  xchange_field(k, ieo);
#endif

#ifdef TM_USE_OMP
#pragma omp parallel
  {
#endif
#define _TM_SUB_HOP
    spinor* pn;

#include "operator/hopping_body_dbl.inc"

#undef _TM_SUB_HOP
#ifdef TM_USE_OMP
  } /* OpenMP closing brace */
#endif
  return;
}
#endif