diff --git a/src/eko/evolution_operator/operator_matrix_element.py.patch b/src/eko/evolution_operator/operator_matrix_element.py.patch
index fbd1665d8..d1aa3861e 100644
--- a/src/eko/evolution_operator/operator_matrix_element.py.patch
+++ b/src/eko/evolution_operator/operator_matrix_element.py.patch
@@ -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
@@ -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_)
@@ -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,
@@ -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)
@@ -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)
@@ -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