diff --git a/dumux/multidomain/facet/cellcentered/mpfa/couplingmanager.hh b/dumux/multidomain/facet/cellcentered/mpfa/couplingmanager.hh
index 2987fd04df3292d886789d38307dd0bfcb4a31f4..03276ccc5dbcfb38d7af985e7fd40885d5f1ec71 100644
--- a/dumux/multidomain/facet/cellcentered/mpfa/couplingmanager.hh
+++ b/dumux/multidomain/facet/cellcentered/mpfa/couplingmanager.hh
@@ -426,7 +426,7 @@ private:
                     static const int numDiffMethod = getParamFromGroup<int>(this->problem(lowDimId).paramGroup(), "Assembly.NumericDifferenceMethod");
                     static const NumericEpsilon< Scalar<lowDimId>, numEq > eps{this->problem(lowDimId).paramGroup()};
                     NumericDifferentiation::partialDerivative(evalResiduals, origPriVars[pvIdx], partialDerivs,
-                                                              origResidual, eps(origPriVars, pvIdx), numDiffMethod);
+                                                              origResidual, eps(origPriVars[pvIdx], pvIdx), numDiffMethod);
 
                     // update the global stiffness matrix with the current partial derivatives
                     // A[i][col][eqIdx][pvIdx] is the rate of change of the residual of equation