Skip to content

Commit

Permalink
Change HPL
Browse files Browse the repository at this point in the history
  • Loading branch information
@niclaurenti
niclaurenti committed Mar 25, 2024
1 parent 376f933 commit 37c53cf
Show file tree
Hide file tree
Showing 11 changed files with 281 additions and 313 deletions.
82 changes: 41 additions & 41 deletions inc/adani/SpecialFunctions.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -45,47 +45,47 @@ double S12(double x);
// Harmonic polylogarithms up to weight 3.
//------------------------------------------------------------------------------------------//

double H0(double x);
double Hm1(double x);
double H1(double x);

double Hm1m1(double x);
double Hm10(double x);
double Hm11(double x);
double H0m1(double x);
double H00(double x);
double H01(double x);
double H1m1(double x);
double H10(double x);
double H11(double x);

double Hm1m1m1(double x);
double Hm1m10(double x);
double Hm1m11(double x);
double Hm10m1(double x);
double Hm100(double x);
double Hm101(double x);
double Hm11m1(double x);
double Hm110(double x);
double Hm111(double x);
double H0m1m1(double x);
double H0m10(double x);
double H0m11(double x);
double H00m1(double x);
double H000(double x);
double H001(double x);
double H01m1(double x);
double H010(double x);
double H011(double x);
double H1m1m1(double x);
double H1m10(double x);
double H1m11(double x);
double H10m1(double x);
double H100(double x);
double H101(double x);
double H11m1(double x);
double H110(double x);
double H111(double x);
double H_0(double x);
double H_m1(double x);
double H_1(double x);

double H_m1m1(double x);
double H_m10(double x);
double H_m11(double x);
double H_0m1(double x);
double H_00(double x);
double H_01(double x);
double H_1m1(double x);
double H_10(double x);
double H_11(double x);

double H_m1m1m1(double x);
double H_m1m10(double x);
double H_m1m11(double x);
double H_m10m1(double x);
double H_m100(double x);
double H_m101(double x);
double H_m11m1(double x);
double H_m110(double x);
double H_m111(double x);
double H_0m1m1(double x);
double H_0m10(double x);
double H_0m11(double x);
double H_00m1(double x);
double H_000(double x);
double H_001(double x);
double H_01m1(double x);
double H_010(double x);
double H_011(double x);
double H_1m1m1(double x);
double H_1m10(double x);
double H_1m11(double x);
double H_10m1(double x);
double H_100(double x);
double H_101(double x);
double H_11m1(double x);
double H_110(double x);
double H_111(double x);

//==========================================================================================//
// Harmonic polylogarithms up to weight 5.
Expand Down
36 changes: 18 additions & 18 deletions src/HighEnergyCoefficientFunction.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -237,7 +237,7 @@ double HighEnergyCoefficientFunction::C2_g2_highenergy(

double L = log((1. + z) / (1. - z));

double Hmp = H11(z) + H1m1(z) - Hm11(z) - Hm1m1(z);
double Hmp = H_11(z) + H_1m1(z) - H_m11(z) - H_m1m1(z);

double I = 4. * z * Hmp;
double J = 4. * z * L;
Expand Down Expand Up @@ -277,7 +277,7 @@ double HighEnergyCoefficientFunction::CL_g2_highenergy(

double L = log((1. + z) / (1. - z));

double Hmp = H11(z) + H1m1(z) - Hm11(z) - Hm1m1(z);
double Hmp = H_11(z) + H_1m1(z) - H_m11(z) - H_m1m1(z);

double I = 4. * z * Hmp;
double J = 4. * z * L;
Expand Down Expand Up @@ -440,11 +440,11 @@ double HighEnergyCoefficientFunction::C2_g3_highenergyLL(

double L = log((1. + z) / (1. - z));

double Hmp = H11(z) + H1m1(z) - Hm11(z) - Hm1m1(z);
double Hmp = H_11(z) + H_1m1(z) - H_m11(z) - H_m1m1(z);

double Hmpm = H111(z) - H11m1(z) + H1m11(z)
- H1m1m1(z) - Hm111(z) + Hm11m1(z)
- Hm1m11(z) + Hm1m1m1(z);
double Hmpm = H_111(z) - H_11m1(z) + H_1m11(z)
- H_1m1m1(z) - H_m111(z) + H_m11m1(z)
- H_m1m11(z) + H_m1m1m1(z);

double I = 4 * z * Hmp;
double J = 4 * z * L;
Expand Down Expand Up @@ -483,11 +483,11 @@ Value HighEnergyCoefficientFunction::C2_g3_highenergyNLL(

double L = log((1. + z) / (1. - z));

double Hmp = H11(z) + H1m1(z) - Hm11(z) - Hm1m1(z);
double Hmp = H_11(z) + H_1m1(z) - H_m11(z) - H_m1m1(z);

double Hmpm = H111(z) - H11m1(z) + H1m11(z)
- H1m1m1(z) - Hm111(z) + Hm11m1(z)
- Hm1m11(z) + Hm1m1m1(z);
double Hmpm = H_111(z) - H_11m1(z) + H_1m11(z)
- H_1m1m1(z) - H_m111(z) + H_m11m1(z)
- H_m1m11(z) + H_m1m1m1(z);

double II = 4 * z * Hmp;
double J = 4 * z * L;
Expand Down Expand Up @@ -634,11 +634,11 @@ double HighEnergyCoefficientFunction::CL_g3_highenergyLL(

double L = log((1. + z) / (1. - z));

double Hmp = H11(z) + H1m1(z) - Hm11(z) - Hm1m1(z);
double Hmp = H_11(z) + H_1m1(z) - H_m11(z) - H_m1m1(z);

double Hmpm = H111(z) - H11m1(z) + H1m11(z)
- H1m1m1(z) - Hm111(z) + Hm11m1(z)
- Hm1m11(z) + Hm1m1m1(z);
double Hmpm = H_111(z) - H_11m1(z) + H_1m11(z)
- H_1m1m1(z) - H_m111(z) + H_m11m1(z)
- H_m1m11(z) + H_m1m1m1(z);

double II = 4 * z * Hmp;
double J = 4 * z * L;
Expand Down Expand Up @@ -683,11 +683,11 @@ Value HighEnergyCoefficientFunction::CL_g3_highenergyNLL(

double L = log((1. + z) / (1. - z));

double Hmp = H11(z) + H1m1(z) - Hm11(z) - Hm1m1(z);
double Hmp = H_11(z) + H_1m1(z) - H_m11(z) - H_m1m1(z);

double Hmpm = H111(z) - H11m1(z) + H1m11(z)
- H1m1m1(z) - Hm111(z) + Hm11m1(z)
- Hm1m11(z) + Hm1m1m1(z);
double Hmpm = H_111(z) - H_11m1(z) + H_1m11(z)
- H_1m1m1(z) - H_m111(z) + H_m11m1(z)
- H_m1m11(z) + H_m1m1m1(z);

double II = 4 * z * Hmp;
double J = 4 * z * L;
Expand Down
64 changes: 32 additions & 32 deletions src/HighScaleCoefficientFunction.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -220,13 +220,13 @@ double HighScaleCoefficientFunction::D2_ps2_highscale(
double L_Q = log(1. / m2Q2) + L_M;
double L_Q2 = L_Q * L_Q;

double H0 = H0(z);
double H1 = H1(z);
double Hm1 = H(z, -1);
double H01 = H(z, 0, 1);
double H0m1 = H(z, 0, -1);
double H001 = H(z, 0, 0, 1);
double H011 = H(z, 0, 1, 1);
double H0 = H_0(z);
double H1 = H_1(z);
double Hm1 = H_m1(z);
double H01 = H_01(z);
double H0m1 = H_0m1(z);
double H001 = H_001(z);
double H011 = H_011(z);

return CF * TR
* ((32. / 3 * H0m1 - 32. / 3 * Hm1 * H0) * (z + 1) * (z + 1)
Expand Down Expand Up @@ -297,23 +297,23 @@ double HighScaleCoefficientFunction::D2_g2_highscale(

double x2 = x * x;

double H0 = H(x, 0);
double H1 = H(x, 1);
double Hm1 = H(x, -1);
double H00 = H(x, 0, 0);
double H01 = H(x, 0, 1);
double H11 = H(x, 1, 1);
double H10 = H(x, 1, 0);
double H011 = H(x, 0, 1, 1);
double H111 = H(x, 1, 1, 1);
double H100 = H(x, 1, 0, 0);
double H000 = H(x, 0, 0, 0);
double H010 = H(x, 0, 1, 0);
double H101 = H(x, 1, 0, 1);
double Hm1m10 = H(x, -1, -1, 0);
double Hm10 = H(x, -1, 0);
double Hm100 = H(x, -1, 0, 0);
double H110 = H(x, 1, 1, 0);
double H0 = H_0(x);
double H1 = H_1(x);
double Hm1 = H_m1(x);
double H00 = H_00(x);
double H01 = H_01(x);
double H11 = H_11(x);
double H10 = H_10(x);
double H011 = H_011(x);
double H111 = H_111(x);
double H100 = H_100(x);
double H000 = H_000(x);
double H010 = H_010(x);
double H101 = H_101(x);
double Hm1m10 = H_m1m10(x);
double Hm10 = H_m10(x);
double Hm100 = H_m100(x);
double H110 = H_110(x);

double LQm = log(1. / m2Q2);
double LQm2 = LQm * LQm;
Expand Down Expand Up @@ -416,11 +416,11 @@ double HighScaleCoefficientFunction::DL_g2_highscale(
double L_M = log(m2mu2);
double L_Q = log(1. / m2Q2) + L_M;

double H0 = H0(z);
double H1 = H1(z);
double Hm1 = H(z, -1);
double H01 = H(z, 0, 1);
double H0m1 = H(z, 0, -1);
double H0 = H_0(z);
double H1 = H_1(z);
double Hm1 = H_m1(z);
double H01 = H_01(z);
double H0m1 = H_0m1(z);

return TR * TR * (-64. / 3 * (z - 1.) * z * L_M)
+ CA * TR
Expand Down Expand Up @@ -461,9 +461,9 @@ double HighScaleCoefficientFunction::DL_ps2_highscale(
double L_M = log(m2mu2);
double L_Q = log(1. / m2Q2) + L_M;

double H0 = H0(z);
double H1 = H1(z);
double H01 = H(z, 0, 1);
double H0 = H_0(z);
double H1 = H_1(z);
double H01 = H_01(z);

return CF * TR
* (32. / 9. / z * (z - 1.) * (10. * z2 - 2. * z + 1.)
Expand Down
Loading

0 comments on commit 37c53cf

Please sign in to comment.