In this work, we present a numerical investigation of a novel Lorentz-covariant parametrization to extract $x$-dependent GPDs using off-forward matrix elements of momentum-boosted hadrons coupled to non-local operators. The novelty of the method is the implementation of an asymmetric frame for the momentum transfer between the initial and final hadron state and the parametrization of the matrix elements into Lorentz-invariant amplitudes. The amplitudes can then be related to the standard light-cone GPDs.

GPDs are defined in the symmetric frame, which requires a separate calculation for each value of the momentum transfer, increasing the computational cost significantly. The proposed method is powerful, as one can extract the GPDs at multiple values of the momentum transfer at the computational cost of a single value. For this proof-of-concept calculation, we use one ensemble of $N_f=2+1+1$ twisted mass fermions and a clover improvement with a pion mass of 260 MeV to calculate the proton unpolarized GPDs.