Skip to content

Commit

Permalink
restore tests and previous Li2 mellin
Browse files Browse the repository at this point in the history
  • Loading branch information
giacomomagni committed May 14, 2024
1 parent 7f281e8 commit 1bf6ae5
Show file tree
Hide file tree
Showing 2 changed files with 13 additions and 14 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
import numba as nb
import numpy as np

from eko.constants import CA, CF, TR, zeta2, zeta3, log2
from eko.constants import CA, CF, TR, log2, zeta2, zeta3

from ....anomalous_dimensions.polarized.space_like.as1 import gamma_gq as gamma0_gq
from ....anomalous_dimensions.polarized.space_like.as1 import gamma_ns as gamma0_qq
Expand Down Expand Up @@ -159,12 +159,11 @@ def A_hg(n, cache, L):
S1 = c.get(c.S1, cache, n)
S2 = c.get(c.S2, cache, n)
S3 = c.get(c.S3, cache, n)
Sm21 = c.get(c.Sm21, cache, n, is_singlet=True)
Sm21 = c.get(c.Sm21, cache, n, is_singlet=False)
S2h = c.get(c.S2h, cache, n)
S3h = c.get(c.S3h, cache, n)
S2ph = c.get(c.S2ph, cache, n)
S3ph = c.get(c.S3ph, cache, n)
Sm1 = c.get(c.Sm1, cache, n, is_singlet=True)

a_hg = (
(1 / (6 * n**4 * (1 + n) ** 4 * (2 + n)))
Expand Down Expand Up @@ -239,7 +238,7 @@ def A_hg(n, cache, L):
+ 6 * (1 + n) ** 3 * S3ph
- 6 * (1 + n) ** 3 * (S3h - zeta3)
+ 18 * (1 + n) ** 3 * zeta3
- 12 * (1 + n) ** 3 * (8 * Sm21 + 5 * zeta3 - 8 * zeta2 * log2 - 8 * zeta2 * Sm1)
- 12 * (1 + n) ** 3 * (8 * Sm21 + 5 * zeta3)
)
)
)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -27,18 +27,18 @@ def test_quark_number_conservation():
def test_hg():
refs = {
0: [
-242.3869306886845,
-480.6122566782359,
-814.6150038529145,
-1244.9445223204148,
-1771.5995558739958,
-14.120648885210752,
-15.093091588201201,
-14.014668171450229,
-12.815826533022285,
-11.761139164575077,
],
10: [
-21.21409118251165,
-295.28287396218656,
-664.653836052461,
-1121.4595082763883,
-1667.8336593381266,
207.0521906209621,
170.2362911278482,
135.9464996290033,
110.6691875110041,
92.0047573712939,
],
}
for L, vals in refs.items():
Expand Down

0 comments on commit 1bf6ae5

Please sign in to comment.