contract_4spinor.cpp File Reference

#include "Measurements/Fermion/contract_4spinor.h"
#include "Field/index_lex.h"
#include <cassert>

Go to the source code of this file.

Macros
#define	NDIM   4

Functions
void	contract_at_t (dcomplex &corr, const GammaMatrix &gm_sink, const Field_F &v1, const Field_F &v2, const int time)
	Contraction of hadron for 4-spinor fermion. More...
void	contract_at_t (std::vector< dcomplex > &corr_global, const GammaMatrix &gm_sink, const Field_F &v1, const Field_F &v2)
void	contract_at_t (dcomplex &corr, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2, const int time)
void	contract_at_t (std::vector< dcomplex > &corr_global, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2)
void	contract_at_t_cos (dcomplex &corr, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2, const int time)
void	contract_at_t_cos (std::vector< dcomplex > &corr_global, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2)
void	contract_at_t (dcomplex &corr, const GammaMatrix &gm_sink, const int i_alpha, const Field_F &v1, const Field_F &v2, const Field_F &v3, const int time)
	contraction for baryon (Nc=3 case only) at a given time t. More...
void	contract_at_x (dcomplex &corr, const GammaMatrix &gm_sink, const Field_F &v1, const Field_F &v2, int x)
void	contract_at_x (std::vector< dcomplex > &corr_global, const GammaMatrix &gm_sink, const Field_F &v1, const Field_F &v2)
void	contract_at_x (dcomplex &corr, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2, const int x)
	contraction for meson at a given x with Fourier transformation, where (p_y,p_z,p_t) is given by momentum_sink. More...
void	contract_at_x (std::vector< dcomplex > &corr_global, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2)
void	contract_at_x_cos (dcomplex &corr, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2, const int x)
	contraction for meson at a given x with Fourier transformation, where (p_y,p_z,p_t) is given by momentum_sink. More...
void	contract_at_x_cos (std::vector< dcomplex > &corr_global, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2)
void	contract_at_x (dcomplex &corr, const GammaMatrix &gm_sink, const int i_alpha, const Field_F &v1, const Field_F &v2, const Field_F &v3, const int x)
	contraction for baryon (Nc=3 case only) at a given x. More...
void	contract_at_y (dcomplex &corr, const GammaMatrix &gm_sink, const Field_F &v1, const Field_F &v2, const int y)
	contraction for meson at a given y. More...
void	contract_at_y (std::vector< dcomplex > &corr_global, const GammaMatrix &gm_sink, const Field_F &v1, const Field_F &v2)
void	contract_at_y (dcomplex &corr, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2, const int y)
	contraction for meson at a given y with Fourier transformation, where (p_x,p_z,p_t) is given by momentum_sink. More...
void	contract_at_y (std::vector< dcomplex > &corr_global, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2)
void	contract_at_y_cos (dcomplex &corr, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2, const int y)
	contraction for meson at a given y with Fourier transformation, where (p_x,p_z,p_t) is given by momentum_sink. More...
void	contract_at_y_cos (std::vector< dcomplex > &corr_global, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2)
void	contract_at_z (dcomplex &corr, const GammaMatrix &gm_sink, const Field_F &v1, const Field_F &v2, const int z)
	contraction for meson at a given z. More...
void	contract_at_z (std::vector< dcomplex > &corr_global, const GammaMatrix &gm_sink, const Field_F &v1, const Field_F &v2)
void	contract_at_z (dcomplex &corr, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2, const int z)
	contraction for meson at a given z with Fourier transformation, where (p_x,p_y,p_t) is given by momentum_sink. More...
void	contract_at_z (std::vector< dcomplex > &corr_global, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2)
void	contract_at_z_cos (dcomplex &corr, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2, const int z)
	contraction for meson at a given z with Fourier transformation, where (p_x,p_y,p_t) is given by momentum_sink. More...
void	contract_at_z_cos (std::vector< dcomplex > &corr_global, const std::vector< int > &momentum_sink, const GammaMatrix &gm_sink, const std::vector< int > &source_position, const Field_F &v1, const Field_F &v2)
void	global_corr_x (std::vector< dcomplex > &corr_global, std::vector< dcomplex > &corr_local)
	transform node-local correlator in x to global. More...
void	global_corr_y (std::vector< dcomplex > &corr_global, std::vector< dcomplex > &corr_local)
	transform node-local correlator in y to global. More...
void	global_corr_z (std::vector< dcomplex > &corr_global, std::vector< dcomplex > &corr_local)
	transform node-local correlator in z to global. More...
void	global_corr_t (std::vector< dcomplex > &corr_global, std::vector< dcomplex > &corr_local)
	transform node-local correlator in t to global. More...

Macro Definition Documentation

#define NDIM   4

Definition at line 19 of file contract_4spinor.cpp.

Function Documentation

void contract_at_t	(	dcomplex &	corr,
		const GammaMatrix &	gm_sink,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	time
	)

Contraction of hadron for 4-spinor fermion.

This class calculates contraction of hadron for 4-spinor fermion. 15 Mar 2012 H.Matsufuru Add momentum of sink. [21 Mar 2015 Y.Namekawa] Add functions contract_at_x(), contract_at_y(), contract_at_z() to measure spatial correlation functions. Add functions contract_at_[x-z]() to measure correlation functions directly by corr_global(t)=sum_{xyz}(v2^*)_alpha(t)(gm_sink)_{alpha,beta}(v1)_beta(t). Add functions contract_at_[x-z]_cos() to measure Fourier transformed correlation functions in terms of cos by corr=sum cos(px)(v2^*)_alpha(t)(gm_sink)_{alpha,beta}(v1)_beta(t). Add functions global_corr_x(), global_corr_y(), global_corr_z() which should always be used with contract_at_[x-z](). A function global_corr_t() is moved here from corr2pt_4spinor.h . [25 May 2017 Y.Taniguchi] Fix assert in Org/contact_4spinor.cpp, kindly reported by Hiroki Ohata. [ 7 May 2021 Y.Namekawa]contraction of two Field_F's at a given time t. corr=sum_{xyz}(v2^*)_alpha(t) (gm_sink)_{alpha,beta} (v1)_beta(t)

Definition at line 39 of file contract_4spinor.cpp.

void contract_at_t	(	std::vector< dcomplex > &	corr_global,
		const GammaMatrix &	gm_sink,
		const Field_F &	v1,
		const Field_F &	v2
	)

contraction of two Field_F's and makes correlator in t. corr_global(t)=sum_{xyz}(v2^*)_alpha(t) (gm_sink)_{alpha,beta} (v1)_beta(t)

Definition at line 94 of file contract_4spinor.cpp.

void contract_at_t	(	dcomplex &	corr,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	time
	)

contraction for meson at a given time t with Fourier transformation. corr=sum_{xyz}e^{ip_xx+ip_yy+ip_zz}(v2^*)_alpha(t) (gm_sink)_{alpha,beta} (v1)_beta(t), where (p_x,p_y,p_z) is given by momentum_sink.

Definition at line 114 of file contract_4spinor.cpp.

void contract_at_t	(	std::vector< dcomplex > &	corr_global,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 205 of file contract_4spinor.cpp.

void contract_at_t	(	dcomplex &	corr,
		const GammaMatrix &	gm_sink,
		const int	i_alpha,
		const Field_F &	v1,
		const Field_F &	v2,
		const Field_F &	v3,
		const int	time
	)

contraction for baryon (Nc=3 case only) at a given time t.

Definition at line 352 of file contract_4spinor.cpp.

void contract_at_t_cos	(	dcomplex &	corr,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	time
	)

contraction for meson at a given time t with Fourier transformation. corr=sum_{xyz}cos(p_xx+p_yy+p_zz)(v2^*)_alpha(t) (gm_sink)_{alpha,beta} (v1)_beta(t), where (p_x,p_y,p_z) is given by momentum_sink. This funciton is intended to be used for a Fourier transformation of the real part Re((v2^*)(gm_sink)(v1)) only when imaginary part is unphysical.

Definition at line 227 of file contract_4spinor.cpp.

void contract_at_t_cos	(	std::vector< dcomplex > &	corr_global,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 330 of file contract_4spinor.cpp.

void contract_at_x	(	dcomplex &	corr,
		const GammaMatrix &	gm_sink,
		const Field_F &	v1,
		const Field_F &	v2,
		int	x
	)

contraction for meson at a given x. corr=sum_{yzt}(v2^*)_alpha(x) (gm_sink)_{alpha,beta} (v1)_beta(x)

Definition at line 460 of file contract_4spinor.cpp.

void contract_at_x	(	std::vector< dcomplex > &	corr_global,
		const GammaMatrix &	gm_sink,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 531 of file contract_4spinor.cpp.

void contract_at_x	(	dcomplex &	corr,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	x
	)

contraction for meson at a given x with Fourier transformation, where (p_y,p_z,p_t) is given by momentum_sink.

Definition at line 551 of file contract_4spinor.cpp.

void contract_at_x	(	std::vector< dcomplex > &	corr_global,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 655 of file contract_4spinor.cpp.

void contract_at_x	(	dcomplex &	corr,
		const GammaMatrix &	gm_sink,
		const int	i_alpha,
		const Field_F &	v1,
		const Field_F &	v2,
		const Field_F &	v3,
		const int	x
	)

contraction for baryon (Nc=3 case only) at a given x.

Definition at line 806 of file contract_4spinor.cpp.

void contract_at_x_cos	(	dcomplex &	corr,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	x
	)

contraction for meson at a given x with Fourier transformation, where (p_y,p_z,p_t) is given by momentum_sink.

Definition at line 677 of file contract_4spinor.cpp.

void contract_at_x_cos	(	std::vector< dcomplex > &	corr_global,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 784 of file contract_4spinor.cpp.

void contract_at_y	(	dcomplex &	corr,
		const GammaMatrix &	gm_sink,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	y
	)

contraction for meson at a given y.

Definition at line 920 of file contract_4spinor.cpp.

void contract_at_y	(	std::vector< dcomplex > &	corr_global,
		const GammaMatrix &	gm_sink,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 992 of file contract_4spinor.cpp.

void contract_at_y	(	dcomplex &	corr,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	y
	)

contraction for meson at a given y with Fourier transformation, where (p_x,p_z,p_t) is given by momentum_sink.

Definition at line 1012 of file contract_4spinor.cpp.

void contract_at_y	(	std::vector< dcomplex > &	corr_global,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 1116 of file contract_4spinor.cpp.

void contract_at_y_cos	(	dcomplex &	corr,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	y
	)

contraction for meson at a given y with Fourier transformation, where (p_x,p_z,p_t) is given by momentum_sink.

Definition at line 1138 of file contract_4spinor.cpp.

void contract_at_y_cos	(	std::vector< dcomplex > &	corr_global,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 1245 of file contract_4spinor.cpp.

void contract_at_z	(	dcomplex &	corr,
		const GammaMatrix &	gm_sink,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	z
	)

contraction for meson at a given z.

Definition at line 1268 of file contract_4spinor.cpp.

void contract_at_z	(	std::vector< dcomplex > &	corr_global,
		const GammaMatrix &	gm_sink,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 1340 of file contract_4spinor.cpp.

void contract_at_z	(	dcomplex &	corr,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	z
	)

contraction for meson at a given z with Fourier transformation, where (p_x,p_y,p_t) is given by momentum_sink.

Definition at line 1360 of file contract_4spinor.cpp.

void contract_at_z	(	std::vector< dcomplex > &	corr_global,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 1464 of file contract_4spinor.cpp.

void contract_at_z_cos	(	dcomplex &	corr,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2,
		const int	z
	)

contraction for meson at a given z with Fourier transformation, where (p_x,p_y,p_t) is given by momentum_sink.

Definition at line 1486 of file contract_4spinor.cpp.

void contract_at_z_cos	(	std::vector< dcomplex > &	corr_global,
		const std::vector< int > &	momentum_sink,
		const GammaMatrix &	gm_sink,
		const std::vector< int > &	source_position,
		const Field_F &	v1,
		const Field_F &	v2
	)

Definition at line 1593 of file contract_4spinor.cpp.

void global_corr_t	(	std::vector< dcomplex > &	corr_global,
		std::vector< dcomplex > &	corr_local
	)

transform node-local correlator in t to global.

Definition at line 1696 of file contract_4spinor.cpp.

void global_corr_x	(	std::vector< dcomplex > &	corr_global,
		std::vector< dcomplex > &	corr_local
	)

transform node-local correlator in x to global.

Definition at line 1615 of file contract_4spinor.cpp.

void global_corr_y	(	std::vector< dcomplex > &	corr_global,
		std::vector< dcomplex > &	corr_local
	)

transform node-local correlator in y to global.

Definition at line 1642 of file contract_4spinor.cpp.

void global_corr_z	(	std::vector< dcomplex > &	corr_global,
		std::vector< dcomplex > &	corr_local
	)

transform node-local correlator in z to global.

Definition at line 1669 of file contract_4spinor.cpp.