diff --git a/dune/mmdg/mmdg.hh b/dune/mmdg/mmdg.hh
index c830287cfacdf8542b50acebc5467ac614e7a03e..1076ec954a92b74e5bf25a29daca35ac56736a6d 100644
--- a/dune/mmdg/mmdg.hh
+++ b/dune/mmdg/mmdg.hh
@@ -350,10 +350,10 @@ private:
 
         const Scalar couplingUpdate1 = quadratureFactor * betaEvaluation2;
 
-        (*Base::A)[elemInIdxSLE][iElemIdx] -= couplingUpdate1;
-        (*Base::A)[iElemIdx][elemInIdxSLE] -= couplingUpdate1;
-        (*Base::A)[elemOutIdxSLE][iElemIdx] -= couplingUpdate1;
-        (*Base::A)[iElemIdx][elemOutIdxSLE] -= couplingUpdate1;
+        (*Base::A)[elemInIdxSLE][iElemIdxSLE] -= couplingUpdate1;
+        (*Base::A)[iElemIdxSLE][elemInIdxSLE] -= couplingUpdate1;
+        (*Base::A)[elemOutIdxSLE][iElemIdxSLE] -= couplingUpdate1;
+        (*Base::A)[iElemIdxSLE][elemOutIdxSLE] -= couplingUpdate1;
 
         for (int i = 0; i < dim - 1; i++)
         {
@@ -361,10 +361,10 @@ private:
           const Scalar couplingUpdate3 =
             (iFrame[i] * qpGlobal) * couplingUpdate1;
 
-          (*Base::A)[elemInIdxSLE + i + 1][iElemIdx] -= couplingUpdate2;
-          (*Base::A)[iElemIdx][elemInIdxSLE + i + 1] -= couplingUpdate2;
-          (*Base::A)[elemOutIdxSLE + i + 1][iElemIdx] -= couplingUpdate2;
-          (*Base::A)[iElemIdx][elemOutIdxSLE + i + 1] -= couplingUpdate2;
+          (*Base::A)[elemInIdxSLE + i + 1][iElemIdxSLE] -= couplingUpdate2;
+          (*Base::A)[iElemIdxSLE][elemInIdxSLE + i + 1] -= couplingUpdate2;
+          (*Base::A)[elemOutIdxSLE + i + 1][iElemIdxSLE] -= couplingUpdate2;
+          (*Base::A)[iElemIdxSLE][elemOutIdxSLE + i + 1] -= couplingUpdate2;
 
           (*Base::A)[elemInIdxSLE][iElemIdx + i + 1] -= couplingUpdate3;
           (*Base::A)[iElemIdx + i + 1][elemInIdxSLE] -= couplingUpdate3;
@@ -374,10 +374,10 @@ private:
 
         //i = dim - 1
         const Scalar couplingUpdate4 = qpGlobal[dim] * couplingUpdate1;
-        (*Base::A)[elemInIdxSLE + dim + 1][iElemIdx] -= couplingUpdate4;
-        (*Base::A)[iElemIdx][elemInIdxSLE + dim + 1] -= couplingUpdate4;
-        (*Base::A)[elemOutIdxSLE + dim + 1][iElemIdx] -= couplingUpdate4;
-        (*Base::A)[iElemIdx][elemOutIdxSLE + dim + 1] -= couplingUpdate4;
+        (*Base::A)[elemInIdxSLE + dim + 1][iElemIdxSLE] -= couplingUpdate4;
+        (*Base::A)[iElemIdxSLE][elemInIdxSLE + dim + 1] -= couplingUpdate4;
+        (*Base::A)[elemOutIdxSLE + dim + 1][iElemIdxSLE] -= couplingUpdate4;
+        (*Base::A)[iElemIdxSLE][elemOutIdxSLE + dim + 1] -= couplingUpdate4;
       }
 
       // === interface entries ===