Skip to content

Added code for u- and d- flavour sequential sources, supporting multiple t-sink and projector input. #2

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
ckallidonis wants to merge 19 commits into
base: master
Choose a base branch
from
Open

Added code for u- and d- flavour sequential sources, supporting multiple t-sink and projector input. #2

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
19 commits
Select commit Hold shift + click to select a range
71255ad
Added code for u- and d- flavour sequential source
Nov 27, 2014
d0cf091
Added a Makefile for FERMI
Nov 27, 2014
928aa25
Merge branch 'master' of https://github.com/g-koutsou/qcd
ckallidonis Jun 4, 2015
4e71d8a
Added a Makefile for Juropa when compiling with gnu compilers
ckallidonis Jun 4, 2015
cd50bea
Moved makefile for Juropa to correct directory
ckallidonis Jun 4, 2015
f3101b4
Fixed a bug for 2pt and 3pt contractions
ckallidonis Jun 4, 2015
fde1b79
Fixed 2pt contractions for pseudoscalar mesons
ckallidonis Jun 4, 2015
507b691
Added a function for computing Propagator*gamma5 in ./lib/qcd_gamma.c
ckallidonis Jun 4, 2015
bdeec67
Renamed: twop_all_opt.c -> twop_all_hadr.c. Added code for scalar, ve…
ckallidonis Jun 8, 2015
8780e5b
Various tweaks
ckallidonis Jun 8, 2015
5c37e68
Added contractions for heavy vector and ax.vector mesons
ckallidonis Jun 11, 2015
34aa8f2
Added a Makefile.in for JURECA
ckallidonis Jul 13, 2015
1cae79f
Corrected a bug for the 3pt-function contractions
ckallidonis Sep 2, 2015
b409206
Added a program to perform 4D APE smearing on configurations
ckallidonis Dec 17, 2015
1075fd5
Fixed a bug concerning lex. definition of position coordinates
ckallidonis Jan 19, 2016
ea28bd2
Added a Makefile for the PHI cluster
ckallidonis Jul 21, 2016
2c45a76
Added a program to create Unity Gauge Fields
ckallidonis Jul 21, 2016
d0b510f
Added a program to read LIME propagators and write them in ASCII format
ckallidonis Jul 21, 2016
57498dc
Added a program to read LIME vectors and write them in ASCII format
ckallidonis Jul 21, 2016
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
7 changes: 7 additions & 0 deletions Makefile.in.FERMI
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
AR=ar
RANLIB=ranlib
CC=mpixlc_r -DINLINE=""
CFLAGS=-O3 -qstrict -I./ -I$(HOME)/install/include/ -I$(GSL_LIB)/../include -qsmp=omp -qthreaded
LDFLAGS=-L$(HOME)/install/lib -L$(GSL_LIB)/
LIBS=-llime -lgslcblas -lgsl -lm -qsmp=omp -qthreaded

6 changes: 6 additions & 0 deletions Makefile.in.JURECA
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
AR=ar
RANLIB=ranlib
CC=mpicc -std=c99
CFLAGS=-O3 -I$(HOME)/install_gnu/include -fopenmp -I./ -I$(HOME)/install_gnu/include -DINLINE="static inline"
LDFLAGS=-L$(HOME)/install_gnu/lib
LIBS=-llime -lgslcblas -lgsl -lxml2 -lm -openmp
6 changes: 6 additions & 0 deletions Makefile.in.JUROPA_gnu
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
AR=ar
RANLIB=ranlib
CC=mpicc -std=c99
CFLAGS=-O3 -fopenmp -m64 -override-limits -I./ -I$(HOME)/install_gnu/include -DINLINE="static inline"
LDFLAGS=-L$(HOME)/install_gnu/lib
LIBS=-llime -lgslcblas -lgsl -lm
7 changes: 7 additions & 0 deletions Makefile.in.PHI
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
AR=ar
RANLIB=ranlib
CC=mpiicc
CFLAGS=-O3 -std=c99 -openmp -I./ -I$(HOME)/install/include/
LDFLAGS=-L$(HOME)/install/lib
LIBS=-llime -lgslcblas -lgsl -lxml2 -lm -liomp5

6 changes: 6 additions & 0 deletions Makefile.in.PIZ_DAINT
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
AR=ar
RANLIB=ranlib
CC=cc -std=c99
CFLAGS=-O3 -openmp -unroll-aggressive -opt-prefetch -finline-functions -I./ -I$(HOME)/install/include/
LDFLAGS=-L$(HOME)/install/lib
LIBS=-llime -lgslcblas -lgsl -lxml2 -lm -openmp
2 changes: 2 additions & 0 deletions include/qcd.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -185,6 +185,8 @@ void qcd_destroyVector3d(qcd_vector3d *vec);
int qcd_initPropagator(qcd_propagator *prp, qcd_geometry *geo);
void qcd_destroyPropagator(qcd_propagator *prp);
int qcd_initGaugeField(qcd_gaugeField *u, qcd_geometry *geo);
int qcd_setZeroGaugeField(qcd_gaugeField *u, qcd_geometry *geo);
int qcd_setUnityGaugeField(qcd_gaugeField *u, qcd_geometry *geo);
void qcd_destroyGaugeField(qcd_gaugeField *u);
int qcd_initGaugeTransformation(qcd_gaugeTransformation *u, qcd_geometry *geo);
void qcd_destroyGaugeTransformation(qcd_gaugeTransformation *u);
Expand Down
19 changes: 19 additions & 0 deletions include/qcd_contract_2pt_pr.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,6 +6,10 @@
*
*/

#ifndef H_QCD_CONTRACT_2PT_PR
#define H_QCD_CONTRACT_2PT_PR 1


int qcd_projector12_2pt(qcd_complex_16 gamma12[4][4], qcd_complex_16 gamma13[4][4], qcd_complex_16 gamma23[4][4],
qcd_complex_16 *block[4][4], qcd_complex_16 *block_pr[9][4][4], qcd_geometry *geo);

Expand All @@ -20,6 +24,18 @@ int qcd_f1f1f1_2pt_pr(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16

int qcd_deltas_xistar_omegastar_2pt_pr(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra11_2pt_pr(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra12_2pt_pr(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra21_2pt_pr(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra22_2pt_pr(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_f1f2f1_2pt_pr(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_geometry *geo, qcd_uint_4 lt,qcd_uint_4 elem);
Expand All @@ -37,3 +53,6 @@ int qcd_sigmas2_xistar_2pt_pr(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_co
void qcd_contractions2pt_pr(qcd_uint_4 particle_id, qcd_complex_16 *block12[4][4], qcd_complex_16 *block32[4][4], qcd_complex_16 *blocknp[4][4],
qcd_propagator *uprop,qcd_propagator *dprop,qcd_propagator *sprop,qcd_propagator *cprop,
qcd_geometry *geo, qcd_uint_4 lt);


#endif
114 changes: 68 additions & 46 deletions include/qcd_contract_3pt_all.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,71 +5,93 @@
* Header file for qcd_contract_3pt_all.c
*
*/

int qcd_projector_3pt(qcd_complex_16 *block[5], qcd_complex_16 *block_pr[10][4][4], qcd_geometry *geo, qcd_uint_4 prid);

int qcd_projector32_3pt(qcd_complex_16 gamma12[4][4], qcd_complex_16 gamma13[4][4], qcd_complex_16 gamma23[4][4],
qcd_complex_16 *block_pr[10][4][4], qcd_geometry *geo);

int qcd_projector_3pt_spin12(qcd_complex_16 *block[5], qcd_complex_16 *block_pr[9][4][4], qcd_geometry *geo);

int qcd_projector_pr32_3pt_spin32(qcd_complex_16 *block[5], qcd_complex_16 *block_pr[9][4][4],
qcd_complex_16 gamma12[4][4], qcd_complex_16 gamma13[4][4], qcd_complex_16 gamma23[4][4], qcd_geometry *geo);
//-----------------------------------

int qcd_f1f1f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *prop, qcd_propagator *seqprop, qcd_geometry *geo, qcd_uint_4 lt,qcd_uint_4 elem);
int qcd_f1f1f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *prop, qcd_propagator *seqprop, qcd_geometry *geo, qcd_uint_4 lt,qcd_uint_4 elem);

int qcd_f1f2f1_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_f1f2f1_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_pnextra_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_pnextra_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_f1f2f1_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_f1f2f1_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_pnextra_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_pnextra_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_deltas_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_deltas_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_deltas_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_f123f321_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf2, qcd_propagator *propf3, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_deltas_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra11_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra21_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra12_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra22_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra11_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra12_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra21_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_xistar_extra22_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_f123f321_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf2, qcd_propagator *propf3, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_f123f321_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf3, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_f123f321_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf3, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_f123f321_f3_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf3, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_f123f321_f3_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf3, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_f1f2f3_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *propf3, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_f1f2f3_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *propf3, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_sigmas4_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_sigmas4_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_sigmas4_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_sigmas4_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_lambdas_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf3, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_lambdas_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf3, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_lambdas_f3_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf3, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);
int qcd_lambdas_f3_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf3, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem);

int qcd_sigmas2_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf2, qcd_propagator *propf3, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem, qcd_uint_4 xis);
int qcd_sigmas2_f1_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf2, qcd_propagator *propf3, qcd_propagator *seqpropf1, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem, qcd_uint_4 xis);

int qcd_sigmas2_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf3, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem, qcd_uint_4 xis);

int qcd_sigmas2_f3_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[10][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf3, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem, qcd_uint_4 xis);
int qcd_sigmas2_f2_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf3, qcd_propagator *seqpropf2, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem, qcd_uint_4 xis);
int qcd_sigmas2_f3_3pt(qcd_int_4 ctr, qcd_int_2 cgcg_ind[16*16][4],qcd_complex_16 cgcg_val[16*16], qcd_complex_16 *block[9][4][4],
qcd_propagator *propf1, qcd_propagator *propf2, qcd_propagator *seqpropf3, qcd_geometry *geo, qcd_uint_4 lt, qcd_uint_4 elem, qcd_uint_4 xis);

//=====================================================================

void qcd_contractions3pt_new(qcd_uint_4 p_id, qcd_uint_4 np, qcd_complex_16 *block[5],
qcd_propagator *uprop, qcd_propagator *dprop, qcd_propagator *sprop, qcd_propagator *cprop,
qcd_propagator *sequprop, qcd_propagator *seqdprop, qcd_propagator *seqsprop, qcd_propagator *seqcprop,
qcd_geometry *geo, qcd_uint_4 lt);
qcd_propagator *uprop, qcd_propagator *dprop, qcd_propagator *sprop, qcd_propagator *cprop,
qcd_propagator *sequprop, qcd_propagator *seqdprop, qcd_propagator *seqsprop, qcd_propagator *seqcprop,
qcd_geometry *geo, qcd_uint_4 lt);
1 change: 1 addition & 0 deletions include/qcd_gamma.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -25,5 +25,6 @@ extern const qcd_complex_16 qcd_BAR_C[4][4];
/* prototypes */
void qcd_gamma5Vector(qcd_vector *v);
void qcd_gamma5Propagator(qcd_propagator *p);
void qcd_Propagatorgamma5(qcd_propagator *p);

#endif
2 changes: 2 additions & 0 deletions include/qcd_io.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,6 +39,8 @@ int qcd_getVectorHMCV(char *fname, qcd_uint_2 nu, qcd_uint_2 c2, qcd_vector *v);
int qcd_getVector(char *fname, int type, qcd_uint_2 nu, qcd_uint_2 c2, qcd_vector *v);
int qcd_writePropagatorCMI(char *fname, qcd_propagator *p);
int qcd_writePropagator(char *fname, int type, qcd_propagator *p);
int qcd_writePropagatorASCII(char *fname_basis, qcd_propagator *prop);
int qcd_writeVectorASCII(char *fname_basis, int mu, qcd_vector *vec);
int qcd_writeVectorCMI(char *fname, qcd_uint_2 nu, qcd_uint_2 c2, qcd_vector *v);
int qcd_writeVectorLime(char *fname, int type, qcd_vector *v);
int qcd_writeVector(char *fname, int type, qcd_uint_2 nu, qcd_uint_2 c2, qcd_vector *v);
Expand Down
10 changes: 5 additions & 5 deletions include/qcd_particle_data.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -9,10 +9,10 @@
#ifndef H_QCD_PARTICLE_DATA
#define H_QCD_PARTICLE_DATA 1

extern char *particles_parts[53][4];
extern int particles_pnum[53][4];
extern char *particle_names[53];
extern char *meson_names[4];
extern int particles32[53];
extern char *particles_parts[61][4];
extern int particles_pnum[61][4];
extern char *particle_names[61];
extern char *meson_names[20];
extern int particles32[61];

#endif
1 change: 1 addition & 0 deletions include/qcd_smearing.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,5 +39,6 @@
int qcd_gaussIteration3d(qcd_vector *v, qcd_gaugeField *u, qcd_real_8 alpha, qcd_uint_4 t);
int qcd_gaussIteration3dAll(qcd_vector *v, qcd_gaugeField *u, qcd_real_8 alpha, qcd_uint_2 gaugeCom);
int qcd_apeSmear3d(qcd_gaugeField *apeu, qcd_gaugeField *u, qcd_real_8 alpha);
int qcd_apeSmear4d(qcd_gaugeField *apeu, qcd_gaugeField *u, qcd_real_8 alpha);

#endif
Loading