portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

TensorRemove not needed now

Browse Source
This commit is contained in:
Azusa Yamaguchi 2015-06-14 00:49:26 +01:00
parent dcfc189f11
commit 2c2112e152
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
tests/Test_gamma.cc
View File

@ -95,7 +95,7 @@ int main (int argc, char ** argv)
for(int mu=0;mu<6;mu++){ for(int mu=0;mu<6;mu++){
result = Gamma(g[mu])* ident * Gamma(g[mu]); result = Gamma(g[mu])* ident * Gamma(g[mu]);
result = result - ident; result = result - ident;
double mag = TensorRemove(norm2(result)); RealD mag = norm2(result);
std::cout << list[mu]<<" " << mag<<std::endl; std::cout << list[mu]<<" " << mag<<std::endl;
} }
@ -103,7 +103,7 @@ int main (int argc, char ** argv)
for(int mu=0;mu<6;mu++){ for(int mu=0;mu<6;mu++){
result = rnd * Gamma(g[mu]); result = rnd * Gamma(g[mu]);
result = result + rnd * Gamma(g[mu+6]); result = result + rnd * Gamma(g[mu+6]);
double mag = TensorRemove(norm2(result)); RealD mag = norm2(result);
std::cout << list[mu]<<" " << mag<<std::endl; std::cout << list[mu]<<" " << mag<<std::endl;
} }
@ -111,7 +111,7 @@ int main (int argc, char ** argv)
for(int mu=0;mu<6;mu++){ for(int mu=0;mu<6;mu++){
result = Gamma(g[mu]) *rnd; result = Gamma(g[mu]) *rnd;
result = result + Gamma(g[mu+6])*rnd; result = result + Gamma(g[mu+6])*rnd;
double mag = TensorRemove(norm2(result)); RealD mag = norm2(result);
std::cout << list[mu]<<" " << mag<<std::endl; std::cout << list[mu]<<" " << mag<<std::endl;
} }