#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <type_traits>
#include "quda.h"
#include "test_util.h"
#include "misc.h"
#include "gauge_force_reference.h"

Include dependency graph for gauge_force_reference.cpp:

Classes
struct	fcomplex

struct	dcomplex

struct	fsu3_matrix

struct	fsu3_vector

struct	dsu3_matrix

struct	dsu3_vector

struct	fanti_hermitmat

struct	danti_hermitmat

Macros
#define	CADD(a, b, c)

#define	CMUL(a, b, c)

#define	CSUM(a, b) { (a).real += (b).real; (a).imag += (b).imag; }

#define	CMULJ_(a, b, c)

#define	CMUL_J(a, b, c)

#define	CONJG(a, b) { (b).real = (a).real; (b).imag = -(a).imag; }

Functions
int	neighborIndexFullLattice (int i, int dx4, int dx3, int dx2, int dx1)

template<typename su3_matrix >
void	su3_adjoint (su3_matrix a, su3_matrix b)

template<typename su3_matrix , typename anti_hermitmat >
void	make_anti_hermitian (su3_matrix m3, anti_hermitmat ah3)

template<typename anti_hermitmat , typename su3_matrix >
static void	uncompress_anti_hermitian (anti_hermitmat mat_antihermit, su3_matrix mat_su3)

template<typename su3_matrix , typename Float >
static void	scalar_mult_sub_su3_matrix (su3_matrix a, su3_matrix b, Float s, su3_matrix *c)

template<typename su3_matrix , typename Float >
static void	scalar_mult_add_su3_matrix (su3_matrix a, su3_matrix b, Float s, su3_matrix *c)

template<typename su3_matrix >
static void	mult_su3_nn (su3_matrix a, su3_matrix b, su3_matrix *c)

template<typename su3_matrix >
static void	mult_su3_an (su3_matrix a, su3_matrix b, su3_matrix *c)

template<typename su3_matrix >
static void	mult_su3_na (su3_matrix a, su3_matrix b, su3_matrix *c)

template<typename su3_matrix >
void	print_su3_matrix (su3_matrix *a)

int	gf_neighborIndexFullLattice (int i, int dx4, int dx3, int dx2, int dx1)

template<typename su3_matrix , typename Float >
static void	compute_path_product (su3_matrix staple, su3_matrix sitelink, su3_matrix sitelink_ex_2d, int path, int len, Float loop_coeff, int dir)

template<typename su3_matrix , typename anti_hermitmat , typename Float >
static void	update_mom (anti_hermitmat momentum, int dir, su3_matrix sitelink, su3_matrix staple, Float eb3)

void	gauge_force_reference_dir (void refMom, int dir, double eb3, void sitelink, void sitelink_ex_2d, QudaPrecision prec, int path_dir, int length, void *loop_coeff, int num_paths)

void	gauge_force_reference (void refMom, double eb3, void sitelink, void sitelink_ex_2d, QudaPrecision prec, int *path_dir, int length, void *loop_coeff, int num_paths)

Variables
int	Z [4]

int	V

int	Vh

int	Vh_ex

int	E [4]

Macro Definition Documentation

◆ CADD

#define CADD	(	a,
		b,
		c
	)

Value:

{ (c).real = (a).real + (b).real; \

(c).imag = (a).imag + (b).imag; }

Definition at line 19 of file gauge_force_reference.cpp.

◆ CMUL

#define CMUL	(	a,
		b,
		c
	)

Value:

{ (c).real = (a).real*(b).real - (a).imag*(b).imag; \

(c).imag = (a).real*(b).imag + (a).imag*(b).real; }

Definition at line 21 of file gauge_force_reference.cpp.

Referenced by mult_su3_nn().

◆ CMUL_J

#define CMUL_J	(	a,
		b,
		c
	)

Value:

{ (c).real = (a).real*(b).real + (a).imag*(b).imag; \

(c).imag = (a).imag*(b).real - (a).real*(b).imag; }

Definition at line 30 of file gauge_force_reference.cpp.

Referenced by mult_su3_na().

◆ CMULJ_

#define CMULJ_	(	a,
		b,
		c
	)

Value:

{ (c).real = (a).real*(b).real + (a).imag*(b).imag; \

(c).imag = (a).real*(b).imag - (a).imag*(b).real; }

Definition at line 26 of file gauge_force_reference.cpp.

Referenced by mult_su3_an().

◆ CONJG

#define CONJG	(	a,
		b
	)	{ (b).real = (a).real; (b).imag = -(a).imag; }

Definition at line 33 of file gauge_force_reference.cpp.

Referenced by su3_adjoint().

◆ CSUM

#define CSUM	(	a,
		b
	)	{ (a).real += (b).real; (a).imag += (b).imag; }

Definition at line 23 of file gauge_force_reference.cpp.

Referenced by mult_su3_an(), mult_su3_na(), and mult_su3_nn().

Function Documentation

◆ compute_path_product()

template<typename su3_matrix , typename Float >

static void compute_path_product	(	su3_matrix *	staple,
		su3_matrix **	sitelink,
		su3_matrix **	sitelink_ex_2d,
		int *	path,
		int	len,
		Float	loop_coeff,
		int	dir
	)

static

Definition at line 277 of file gauge_force_reference.cpp.

References su3_matrix< real >::e, gf_neighborIndexFullLattice(), GOES_FORWARDS, memset(), mult_su3_na(), mult_su3_nn(), OPP_DIR, scalar_mult_add_su3_matrix(), su3_adjoint(), and V.

Referenced by gauge_force_reference_dir().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ gauge_force_reference()

void gauge_force_reference	(	void *	refMom,
		double	eb3,
		void **	sitelink,
		void **	sitelink_ex_2d,
		QudaPrecision	prec,
		int ***	path_dir,
		int *	length,
		void *	loop_coeff,
		int	num_paths
	)

Definition at line 408 of file gauge_force_reference.cpp.

References gauge_force_reference_dir().

Referenced by gauge_force_test().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ gauge_force_reference_dir()

void gauge_force_reference_dir	(	void *	refMom,
		int	dir,
		double	eb3,
		void **	sitelink,
		void **	sitelink_ex_2d,
		QudaPrecision	prec,
		int **	path_dir,
		int *	length,
		void *	loop_coeff,
		int	num_paths
	)