Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Interface QED and N3LO #311

Closed
wants to merge 49 commits into from
Closed
Changes from 1 commit
Commits
Show all changes
49 commits
Select commit Hold shift + click to select a range
5bb0389
Init LH23 benchmark
felixhekhorn Aug 9, 2023
36049b3
Pass labels in apply_pdf
felixhekhorn Aug 9, 2023
031527e
Split and enhance lh bench script
felixhekhorn Aug 17, 2023
403afea
Fix and improve LH bench script
felixhekhorn Aug 17, 2023
4680124
Merge branch 'master' into lh-bench-23
giacomomagni Sep 13, 2023
e8a0079
intorduce PTO_match
giacomomagni Sep 20, 2023
9eb1dc4
Interface QED and N3LO
niclaurenti Oct 6, 2023
089391f
Add n3lo_ad_variation key to missing functions
niclaurenti Oct 6, 2023
9db7ac8
Add n3lo_ad_variation key to missing functions. again
niclaurenti Oct 6, 2023
24bbfd3
Fix test_kernels_QEDsinglet.py
niclaurenti Oct 6, 2023
3508713
Add order=4 in fixed_alphaem_exact
niclaurenti Oct 6, 2023
0118d11
Fix test_kernels_QEDns.py
niclaurenti Oct 6, 2023
361afa8
Fix again test
niclaurenti Oct 6, 2023
17404d3
Fix another test
niclaurenti Oct 6, 2023
263f606
Change nf->beta into ns kernels at n3lo
niclaurenti Oct 7, 2023
9393ea2
Add pto = 3 in some qed tests
niclaurenti Oct 8, 2023
da043cc
Merge branch 'master' into lh-bench-23
giacomomagni Oct 9, 2023
66f02cc
add a failing bench
giacomomagni Oct 9, 2023
c712cad
fix bench name
giacomomagni Oct 9, 2023
9d74a98
separate nf black in evolve_pdfs
giacomomagni Oct 9, 2023
df26b33
Update src/eko/io/runcards.py
giacomomagni Oct 9, 2023
6d14529
Update src/eko/runner/parts.py
giacomomagni Oct 9, 2023
8e9cf99
add default matching_order to legacy runner
giacomomagni Oct 9, 2023
e0d74a7
add fix and correct file name
giacomomagni Oct 9, 2023
ce71c5d
run pre-commit
giacomomagni Oct 9, 2023
a05c58c
Tear evolve_pdfs apart
felixhekhorn Oct 9, 2023
a7f5d7a
Add more unit tests for evol_gri
felixhekhorn Oct 9, 2023
e0e4ea6
[pre-commit.ci] pre-commit autoupdate
pre-commit-ci[bot] Oct 9, 2023
3f0ed5f
Merge pull request #313 from NNPDF/evolve_pdfs_by_nfbloc
giacomomagni Oct 10, 2023
e827328
Merge pull request #314 from NNPDF/pre-commit-ci-update-config
felixhekhorn Oct 10, 2023
49b8b12
Merge pull request #299 from NNPDF/lh-bench-23
giacomomagni Oct 10, 2023
fd701a6
Merge branch 'master' into fix_matching_inversion_2
giacomomagni Oct 10, 2023
91e59c7
minor fix on docstring
giacomomagni Oct 10, 2023
569a0fe
restructure self.compute_aem_list
niclaurenti Oct 10, 2023
2f65ee9
Remove unused default parameters and remove qed argument from apply_pdf
niclaurenti Oct 10, 2023
f39b012
Merge pull request #312 from NNPDF/fix_matching_inversion_2
giacomomagni Oct 11, 2023
0d52dac
Interface QED and N3LO
niclaurenti Oct 6, 2023
c6ad7de
Add n3lo_ad_variation key to missing functions
niclaurenti Oct 6, 2023
b7eedfc
Add n3lo_ad_variation key to missing functions. again
niclaurenti Oct 6, 2023
5646753
Fix test_kernels_QEDsinglet.py
niclaurenti Oct 6, 2023
b67c9d0
Add order=4 in fixed_alphaem_exact
niclaurenti Oct 6, 2023
1f95326
Fix test_kernels_QEDns.py
niclaurenti Oct 6, 2023
607d0eb
Fix again test
niclaurenti Oct 6, 2023
0009b3f
Fix another test
niclaurenti Oct 6, 2023
233e1de
Change nf->beta into ns kernels at n3lo
niclaurenti Oct 7, 2023
69ddbde
Add pto = 3 in some qed tests
niclaurenti Oct 8, 2023
ed4dedf
restructure self.compute_aem_list
niclaurenti Oct 10, 2023
b13a260
Remove unused default parameters and remove qed argument from apply_pdf
niclaurenti Oct 10, 2023
9803d22
Merge branch 'n3lo_plus_qed' of github.com:NNPDF/eko into n3lo_plus_qed
niclaurenti Oct 12, 2023
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
Change nf->beta into ns kernels at n3lo
niclaurenti committed Oct 12, 2023
commit 233e1de3ee16562d429c5bcf98f8c9e463d6f477
24 changes: 8 additions & 16 deletions src/eko/kernels/non_singlet.py
Original file line number Diff line number Diff line change
@@ -186,7 +186,7 @@ def nnlo_expanded(gamma_ns, a1, a0, beta):


@nb.njit(cache=True)
def n3lo_expanded(gamma_ns, a1, a0, nf):
def n3lo_expanded(gamma_ns, a1, a0, beta):
"""|N3LO| non-singlet expanded EKO.

Parameters
@@ -206,12 +206,8 @@ def n3lo_expanded(gamma_ns, a1, a0, nf):
|N3LO| non-singlet expanded EKO

"""
beta0 = beta.beta_qcd((2, 0), nf)
b_list = [
beta.b_qcd((3, 0), nf),
beta.b_qcd((4, 0), nf),
beta.b_qcd((5, 0), nf),
]
beta0 = beta[0]
b_list = [betas / beta0 for betas in beta[1:]]
j12 = ei.j12(a1, a0, beta0)
j13 = as4_ei.j13_expanded(a1, a0, beta0, b_list)
j23 = as4_ei.j23_expanded(a1, a0, beta0, b_list)
@@ -225,7 +221,7 @@ def n3lo_expanded(gamma_ns, a1, a0, nf):


@nb.njit(cache=True)
def n3lo_exact(gamma_ns, a1, a0, nf):
def n3lo_exact(gamma_ns, a1, a0, beta):
"""|N3LO| non-singlet exact EKO.

Parameters
@@ -245,12 +241,8 @@ def n3lo_exact(gamma_ns, a1, a0, nf):
|N3LO| non-singlet exact EKO

"""
beta0 = beta.beta_qcd((2, 0), nf)
b_list = [
beta.b_qcd((3, 0), nf),
beta.b_qcd((4, 0), nf),
beta.b_qcd((5, 0), nf),
]
beta0 = beta[0]
b_list = [betas / beta0 for betas in beta[1:]]
roots = as4_ei.roots(b_list)
j12 = ei.j12(a1, a0, beta0)
j13 = as4_ei.j13_exact(a1, a0, beta0, b_list, roots)
@@ -420,6 +412,6 @@ def dispatcher(
"decompose-expanded",
"perturbative-expanded",
]:
return n3lo_expanded(gamma_ns, a1, a0, nf)
return n3lo_exact(gamma_ns, a1, a0, nf)
return n3lo_expanded(gamma_ns, a1, a0, betalist)
return n3lo_exact(gamma_ns, a1, a0, betalist)
raise NotImplementedError("Selected order is not implemented")