Skip to content

Commit

Permalink
small fix
Browse files Browse the repository at this point in the history
  • Loading branch information
tgiani committed Jan 25, 2024
1 parent e082fb1 commit c201f85
Showing 1 changed file with 7 additions and 7 deletions.
14 changes: 7 additions & 7 deletions src/ekore/scet_I/as1.py
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
import numba as nb
import numpy as np

from eko.constants import CF
from eko.constants import CF, zeta2
from ..harmonics import cache as c


Expand Down Expand Up @@ -179,11 +179,11 @@ def Agg10(n, cache):
+ 5 * np.power(n, 2)
+ 2 * np.power(n, 3)
+ 2 * np.power(n, 4)
+ 2 * np.power(n, 2) * z2
- np.power(n, 3) * z2
- 3 * np.power(n, 4) * z2
+ np.power(n, 5) * z2
+ np.power(n, 6) * z2
+ 2 * np.power(n, 2) * zeta2
- np.power(n, 3) * zeta2
- 3 * np.power(n, 4) * zeta2
+ np.power(n, 5) * zeta2
+ np.power(n, 6) * zeta2
)
)
/ (np.power(-1 + n, 2) * np.power(n, 2) * (1 + n) * (2 + n))
Expand Down Expand Up @@ -254,7 +254,7 @@ def Aqq10(n, cache):
S1 = c.get(c.S1, cache, n)
S2 = c.get(c.S2, cache, n)
res = (
(2 * (1 + 2 * n + 2 * np.power(n, 2) * z2 + 2 * np.power(n, 3) * z2))
(2 * (1 + 2 * n + 2 * np.power(n, 2) * zeta2 + 2 * np.power(n, 3) * zeta2))
/ (3.0 * np.power(n, 2) * (1 + n))
- (2 * S1) / (3.0 * n * (1 + n))
+ (2 * np.power(S1, 2)) / 3.0
Expand Down

0 comments on commit c201f85

Please sign in to comment.