From e699b7e9f9c2643f156fd1cba1b6a01131f7107d Mon Sep 17 00:00:00 2001 From: nmeyer-ur Date: Fri, 10 Apr 2020 12:18:30 +0200 Subject: [PATCH] changed debug output to variable direct --- .../WilsonKernelsHandImplementation.h | 612 ++++++------------ 1 file changed, 204 insertions(+), 408 deletions(-) diff --git a/Grid/qcd/action/fermion/implementation/WilsonKernelsHandImplementation.h b/Grid/qcd/action/fermion/implementation/WilsonKernelsHandImplementation.h index cdf1ade5..6eec803a 100644 --- a/Grid/qcd/action/fermion/implementation/WilsonKernelsHandImplementation.h +++ b/Grid/qcd/action/fermion/implementation/WilsonKernelsHandImplementation.h @@ -91,30 +91,18 @@ Author: paboyle 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 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 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); \