portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-09-20 09:15:38 +01:00
Code Releases Activity

changed debug output to variable direct

Browse Source
This commit is contained in:
nmeyer-ur 2020-04-10 12:18:30 +02:00
parent a28bc0de90
commit e699b7e9f9
1 changed files with 204 additions and 408 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

612
Grid/qcd/action/fermion/implementation/WilsonKernelsHandImplementation.h
View File

@ -91,30 +91,18 @@ Author: paboyle <paboyle@ph.ed.ac.uk>
Chimu_31=ref()(3)(1);\
Chimu_32=ref()(3)(2);\
std::cout << std::endl << "DEBUG -- LOAD_CHIMU" << std::endl; \
svst1(pg1, &debugreg.v, Chimu_00); \
std::cout << "Chimu_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_01); \
std::cout << "Chimu_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_02); \
std::cout << "Chimu_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_10); \
std::cout << "Chimu_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_11); \
std::cout << "Chimu_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_12); \
std::cout << "Chimu_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_20); \
std::cout << "Chimu_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_21); \
std::cout << "Chimu_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_22); \
std::cout << "Chimu_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_30); \
std::cout << "Chimu_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_31); \
std::cout << "Chimu_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chimu_32); \
std::cout << "Chimu_32 -- " << debugreg << std::endl; \
std::cout << "Chimu_00 -- " << Chimu_00 << std::endl; \
std::cout << "Chimu_01 -- " << Chimu_01 << std::endl; \
std::cout << "Chimu_02 -- " << Chimu_02 << std::endl; \
std::cout << "Chimu_10 -- " << Chimu_10 << std::endl; \
std::cout << "Chimu_11 -- " << Chimu_11 << std::endl; \
std::cout << "Chimu_12 -- " << Chimu_12 << std::endl; \
std::cout << "Chimu_20 -- " << Chimu_20 << std::endl; \
std::cout << "Chimu_21 -- " << Chimu_21 << std::endl; \
std::cout << "Chimu_22 -- " << Chimu_22 << std::endl; \
std::cout << "Chimu_30 -- " << Chimu_30 << std::endl; \
std::cout << "Chimu_31 -- " << Chimu_31 << std::endl; \
std::cout << "Chimu_32 -- " << Chimu_32 << std::endl; \
}
#define LOAD_CHI\
@ -126,18 +114,12 @@ Author: paboyle <paboyle@ph.ed.ac.uk>
Chi_11 = ref()(1)(1);\
Chi_12 = ref()(1)(2);\
std::cout << std::endl << "DEBUG -- LOAD_CHI" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl; \
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl; \
}
// To splat or not to splat depends on the implementation
@ -171,35 +153,23 @@ Author: paboyle <paboyle@ph.ed.ac.uk>
UChi_02+= U_20*Chi_02;\
UChi_12+= U_20*Chi_12;\
std::cout << std::endl << "DEBUG -- MULT_2SPIN" << std::endl; \
svst1(pg1, &debugreg.v, UChi_00); \
std::cout << "UChi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, UChi_01); \
std::cout << "UChi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, UChi_02); \
std::cout << "UChi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, UChi_10); \
std::cout << "UChi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, UChi_11); \
std::cout << "UChi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, UChi_12); \
std::cout << "UChi_12 -- " << debugreg << std::endl; \
std::cout << "UChi_00 -- " << UChi_00 << std::endl; \
std::cout << "UChi_01 -- " << UChi_01 << std::endl; \
std::cout << "UChi_02 -- " << UChi_02 << std::endl; \
std::cout << "UChi_10 -- " << UChi_10 << std::endl; \
std::cout << "UChi_11 -- " << UChi_11 << std::endl; \
std::cout << "UChi_12 -- " << UChi_12 << std::endl; \
}
#define PERMUTE_DIR(dir) \
std::cout << std::endl << "DEBUG -- PERM PRE" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl; \
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl; \
permute##dir(Chi_00,Chi_00);\
permute##dir(Chi_01,Chi_01);\
permute##dir(Chi_02,Chi_02);\
@ -207,18 +177,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
permute##dir(Chi_11,Chi_11);\
permute##dir(Chi_12,Chi_12);\
std::cout << std::endl << "DEBUG -- PERM POST" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
// hspin(0)=fspin(0)+timesI(fspin(3));
// hspin(1)=fspin(1)+timesI(fspin(2));
@ -230,18 +194,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
Chi_11 = Chimu_11+timesI(Chimu_21);\
Chi_12 = Chimu_12+timesI(Chimu_22);\
std::cout << std::endl << "DEBUG -- XP_PROJ" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
#define YP_PROJ \
Chi_00 = Chimu_00-Chimu_30;\
@ -251,18 +209,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
Chi_11 = Chimu_11+Chimu_21;\
Chi_12 = Chimu_12+Chimu_22;\
std::cout << std::endl << "DEBUG -- YP_PROJ" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
#define ZP_PROJ \
Chi_00 = Chimu_00+timesI(Chimu_20); \
@ -272,18 +224,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
Chi_11 = Chimu_11-timesI(Chimu_31); \
Chi_12 = Chimu_12-timesI(Chimu_32);\
std::cout << std::endl << "DEBUG -- ZP_PROJ" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
#define TP_PROJ \
Chi_00 = Chimu_00+Chimu_20; \
@ -293,18 +239,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
Chi_11 = Chimu_11+Chimu_31; \
Chi_12 = Chimu_12+Chimu_32;\
std::cout << std::endl << "DEBUG -- TP_PROJ" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
// hspin(0)=fspin(0)-timesI(fspin(3));
@ -317,18 +257,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
Chi_11 = Chimu_11-timesI(Chimu_21);\
Chi_12 = Chimu_12-timesI(Chimu_22);\
std::cout << std::endl << "DEBUG -- XM_PROJ" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
#define YM_PROJ \
Chi_00 = Chimu_00+Chimu_30;\
@ -338,18 +272,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
Chi_11 = Chimu_11-Chimu_21;\
Chi_12 = Chimu_12-Chimu_22;\
std::cout << std::endl << "DEBUG -- YM_PROJ" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
#define ZM_PROJ \
Chi_00 = Chimu_00-timesI(Chimu_20); \
@ -359,18 +287,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
Chi_11 = Chimu_11+timesI(Chimu_31); \
Chi_12 = Chimu_12+timesI(Chimu_32);\
std::cout << std::endl << "DEBUG -- ZM_PROJ" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
#define TM_PROJ \
Chi_00 = Chimu_00-Chimu_20; \
@ -380,18 +302,12 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
Chi_11 = Chimu_11-Chimu_31; \
Chi_12 = Chimu_12-Chimu_32;\
std::cout << std::endl << "DEBUG -- TM_PROJ" << std::endl; \
svst1(pg1, &debugreg.v, Chi_00); \
std::cout << "Chi_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_01); \
std::cout << "Chi_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_02); \
std::cout << "Chi_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_10); \
std::cout << "Chi_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_11); \
std::cout << "Chi_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, Chi_12); \
std::cout << "Chi_12 -- " << debugreg << std::endl;
std::cout << "Chi_00 -- " << Chi_00 << std::endl; \
std::cout << "Chi_01 -- " << Chi_01 << std::endl; \
std::cout << "Chi_02 -- " << Chi_02 << std::endl; \
std::cout << "Chi_10 -- " << Chi_10 << std::endl; \
std::cout << "Chi_11 -- " << Chi_11 << std::endl; \
std::cout << "Chi_12 -- " << Chi_12 << std::endl;
// fspin(0)=hspin(0);
// fspin(1)=hspin(1);
@ -411,30 +327,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31 = timesMinusI(UChi_01);\
result_32 = timesMinusI(UChi_02);\
std::cout << std::endl << "DEBUG -- XP_RECON" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define XP_RECON_ACCUM\
result_00+=UChi_00;\
@ -450,30 +354,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31-=timesI(UChi_01);\
result_32-=timesI(UChi_02);\
std::cout << std::endl << "DEBUG -- XP_RECON_ACCUM" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define XM_RECON\
result_00 = UChi_00;\
@ -489,30 +381,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31 = timesI(UChi_01);\
result_32 = timesI(UChi_02);\
std::cout << std::endl << "DEBUG -- XM_RECON" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define XM_RECON_ACCUM\
result_00+= UChi_00;\
@ -528,30 +408,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31+= timesI(UChi_01);\
result_32+= timesI(UChi_02);\
std::cout << std::endl << "DEBUG -- XM_RECON_ACCUM" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define YP_RECON_ACCUM\
result_00+= UChi_00;\
@ -567,30 +435,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31-= UChi_01;\
result_32-= UChi_02;\
std::cout << std::endl << "DEBUG -- YP_RECON_ACCUM" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define YM_RECON_ACCUM\
result_00+= UChi_00;\
@ -606,30 +462,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31+= UChi_01;\
result_32+= UChi_02;\
std::cout << std::endl << "DEBUG -- YM_RECON_ACCUM" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define ZP_RECON_ACCUM\
result_00+= UChi_00;\
@ -645,30 +489,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31+= timesI(UChi_11); \
result_32+= timesI(UChi_12);\
std::cout << std::endl << "DEBUG -- ZP_RECON_ACCUM" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define ZM_RECON_ACCUM\
result_00+= UChi_00;\
@ -684,30 +516,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31-= timesI(UChi_11); \
result_32-= timesI(UChi_12);\
std::cout << std::endl << "DEBUG -- ZM_RECON_ACCUM" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define TP_RECON_ACCUM\
result_00+= UChi_00;\
@ -723,30 +543,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31+= UChi_11; \
result_32+= UChi_12;\
std::cout << std::endl << "DEBUG -- TP_RECON_ACCUM" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define TM_RECON_ACCUM\
result_00+= UChi_00;\
@ -762,30 +570,18 @@ std::cout << "Chi_12 -- " << debugreg << std::endl; \
result_31-= UChi_11; \
result_32-= UChi_12;\
std::cout << std::endl << "DEBUG -- TM_RECON_ACCUM" << std::endl; \
svst1(pg1, &debugreg.v, result_00); \
std::cout << "result_00 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_01); \
std::cout << "result_01 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_02); \
std::cout << "result_02 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_10); \
std::cout << "result_10 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_11); \
std::cout << "result_11 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_12); \
std::cout << "result_12 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_20); \
std::cout << "result_20 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_21); \
std::cout << "result_21 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_22); \
std::cout << "result_22 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_30); \
std::cout << "result_30 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_31); \
std::cout << "result_31 -- " << debugreg << std::endl; \
svst1(pg1, &debugreg.v, result_32); \
std::cout << "result_32 -- " << debugreg << std::endl;
std::cout << "result_00 -- " << result_00 << std::endl; \
std::cout << "result_01 -- " << result_01 << std::endl; \
std::cout << "result_02 -- " << result_02 << std::endl; \
std::cout << "result_10 -- " << result_10 << std::endl; \
std::cout << "result_11 -- " << result_11 << std::endl; \
std::cout << "result_12 -- " << result_12 << std::endl; \
std::cout << "result_20 -- " << result_20 << std::endl; \
std::cout << "result_21 -- " << result_21 << std::endl; \
std::cout << "result_22 -- " << result_22 << std::endl; \
std::cout << "result_30 -- " << result_30 << std::endl; \
std::cout << "result_31 -- " << result_31 << std::endl; \
std::cout << "result_32 -- " << result_32 << std::endl;
#define HAND_STENCIL_LEG(PROJ,PERM,DIR,RECON) \
SE=st.GetEntry(ptype,DIR,ss); \