Skip to content

Commit

Permalink
Force argument to self.nf
Browse files Browse the repository at this point in the history
  • Loading branch information
niclaurenti committed Jan 10, 2024
1 parent 75bfca0 commit b623f97
Showing 1 changed file with 6 additions and 2 deletions.
8 changes: 6 additions & 2 deletions src/eko/evolution_operator/__init__.py
Original file line number Diff line number Diff line change
Expand Up @@ -707,13 +707,17 @@ def compute_aem_list(self):
couplings = self.managers["couplings"]
mu2_steps = utils.geomspace(self.q2_from, self.q2_to, 1 + ev_op_iterations)
mu2_l = mu2_steps[0]
as_list = np.array([couplings.a_s(scale_to=mu2) for mu2 in mu2_steps])
as_list = np.array(
[couplings.a_s(scale_to=mu2, nf_to=self.nf) for mu2 in mu2_steps]
)
a_half = np.zeros((ev_op_iterations, 2))
for step, mu2_h in enumerate(mu2_steps[1:]):
mu2_half = (mu2_h + mu2_l) / 2.0
a_s, aem = couplings.a(scale_to=mu2_half)
a_s, aem = couplings.a(scale_to=mu2_half, nf_to=self.nf)
a_half[step] = [a_s, aem]
mu2_l = mu2_h
print(as_list.tolist())
print(a_half.tolist())
return as_list, a_half

@property
Expand Down

0 comments on commit b623f97

Please sign in to comment.