Skip to content

Commit

Permalink
Fix OME patch
Browse files Browse the repository at this point in the history
  • Loading branch information
felixhekhorn committed Sep 16, 2024
1 parent 7d5dac7 commit ddc57bf
Showing 1 changed file with 10 additions and 11 deletions.
21 changes: 10 additions & 11 deletions src/eko/evolution_operator/operator_matrix_element.py.patch
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
diff --git a/src/eko/evolution_operator/operator_matrix_element.py b/src/eko/evolution_operator/operator_matrix_element.py
index 022f8073..a636ea8b 100644
index 098b52db..b56e27d1 100644
--- a/src/eko/evolution_operator/operator_matrix_element.py
+++ b/src/eko/evolution_operator/operator_matrix_element.py
@@ -2,21 +2,18 @@
@@ -3,22 +3,19 @@ evolution."""

import copy
import enum
Expand All @@ -21,13 +21,14 @@ index 022f8073..a636ea8b 100644
from .. import scale_variations as sv
from ..io.types import InversionMethod
from ..matchings import Segment
from ..scale_variations.exponentiated import gamma_variation
-from . import Operator, QuadKerBase
+from . import Operator
+from .quad_ker import cb_quad_ker_ome

logger = logging.getLogger(__name__)

@@ -77,8 +74,6 @@ def build_ome(A, matching_order, a_s, backward_method):
@@ -78,8 +75,6 @@ def build_ome(A, matching_order, a_s, backward_method):
# Print;
# .end
ome = np.eye(len(A[0]), dtype=np.complex_)
Expand All @@ -36,11 +37,10 @@ index 022f8073..a636ea8b 100644
if backward_method is MatchingMethods.BACKWARD_EXPANDED:
# expended inverse
if matching_order[0] >= 1:
@@ -100,107 +95,6 @@ def build_ome(A, matching_order, a_s, backward_method):
ome = np.linalg.inv(ome)
@@ -102,105 +97,6 @@ def build_ome(A, matching_order, a_s, backward_method):
return ome

-

-@nb.njit(cache=True)
-def quad_ker(
- u,
Expand Down Expand Up @@ -97,7 +97,6 @@ index 022f8073..a636ea8b 100644
- -------
- ker : float
- evaluated integration kernel
-
- """
- ker_base = QuadKerBase(u, is_log, logx, mode0)
- integrand = ker_base.integrand(areas)
Expand Down Expand Up @@ -129,7 +128,7 @@ index 022f8073..a636ea8b 100644
-
- # correct for scale variations
- if sv_mode == sv.Modes.exponentiated:
- A = sv.exponentiated.gamma_variation(A, order, nf, Lsv)
- A = gamma_variation(A, order, nf, Lsv)
-
- # build the expansion in alpha_s depending on the strategy
- ker = build_ome(A, order, a_s, backward_method)
Expand All @@ -142,9 +141,9 @@ index 022f8073..a636ea8b 100644
-
-
class OperatorMatrixElement(Operator):
r"""
Internal representation of a single |OME|.
@@ -291,41 +185,14 @@ class OperatorMatrixElement(Operator):
r"""Internal representation of a single |OME|.

@@ -290,41 +186,14 @@ class OperatorMatrixElement(Operator):
)
return labels

Expand Down

0 comments on commit ddc57bf

Please sign in to comment.