portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-10 14:10:46 +01:00
Code Releases Activity

continued with baryons

Browse Source
This commit is contained in:
Felix Erben 2019-02-26 16:41:52 +00:00
parent 50b6db75da
commit f168a9e7ee
1 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
Grid/qcd/utils/A2Autils.h
View File

@ -100,7 +100,6 @@ public:
int orthogdim);
#endif
};
/*
template <class FImpl>
template <typename TensorType>
@ -122,6 +121,8 @@ void A2Autils<FImpl>::BaryonField(TensorType &mat,
typedef iSpinMatrix<vector_type> SpinMatrix_v;
typedef iSpinMatrix<scalar_type> SpinMatrix_s;
typedef iSpinColourMatrix<vector_type> SpinColourMatrix_v;
int oneBlock = mat.dimension(3);
int twoBlock = mat.dimension(4);
int threeBlock = mat.dimension(5);
@ -180,13 +181,16 @@ void A2Autils<FImpl>::BaryonField(TensorType &mat,
auto three_k = three[j]._odata[ss];
SpinMatrix_v vv;
SpinColourMatrix_v vv;
for(int s1=0;s1<Ns;s1++){
for(int s2=0;s2<Ns;s2++){
vv()(s1,s2)() = two_j()(s2)(0) * three_k()(s1)(0) //make this a colorMatrix for the diquark???
+ two_j()(s2)(1) * three_k()(s1)(1)
+ two_j()(s2)(2) * three_k()(s1)(2);
vv()(s1,s2)(0,0) = two_j()(s2)(1) * three_k()(s1)(2) //ideal would be SpinMatrix but ColourVector...
- two_j()(s2)(2) * three_k()(s1)(1); //this is the cross product (two x three)^i
vv()(s1,s2)(1,1) = two_j()(s2)(2) * three_k()(s1)(0)
- two_j()(s2)(0) * three_k()(s1)(2);
vv()(s1,s2)(2,2) = two_j()(s2)(0) * three_k()(s1)(1)
- two_j()(s2)(1) * three_k()(s1)(0);
}}
// After getting the sitewise product do the mom phase loop