diff --git a/src/spatial-solving/fixedpointiterator.cc b/src/spatial-solving/fixedpointiterator.cc
index 7b70ca59b11540d794282b136a344d03be313d8c..1d4756840593d9c537e373eb43f57a9a112b66d2 100644
--- a/src/spatial-solving/fixedpointiterator.cc
+++ b/src/spatial-solving/fixedpointiterator.cc
@@ -125,8 +125,6 @@ void FixedPointIterator<Factory, Updaters, ErrorNorm>::relativeVelocities(std::v
   using field_type = typename Factory::Matrix::field_type;
 
 
-
-
   // adaptation of DualMortarCoupling::setup()
 
   const size_t dim = DeformedGrid::dimension;
@@ -147,10 +145,14 @@ void FixedPointIterator<Factory, Updaters, ErrorNorm>::relativeVelocities(std::v
     gridFunctions[i] = new BasisGridFunction<MyAssembler::VertexBasis, Vector>(assemblers[i]->vertexBasis, v_m[i]);
   }
 
-  const auto& contactCouplings = nBodyAssembler_.getContactCouplings();
-  for (size_t i=0; i<contactCouplings.size(); i++) {
-    auto contactCoupling = contactCouplings[i];
-    auto glue = contactCoupling->getGlue();
+  /*
+  for (size_t i=0; i<nBodyAssembler_.nCouplings(); i++) {
+    const auto& coupling = nBodyAssembler_.getCoupling(i);
+    auto glue = coupling.backend();
+
+    const std::array<int, 2> gridIdx = coupling.gridIdx_;
+    const int nonmortarGridIdx = ;
+    const int mortarGridIdx = ;
 
     // loop over all intersections
     for (const auto& rIs : intersections(glue)) {
@@ -288,7 +290,7 @@ void FixedPointIterator<Factory, Updaters, ErrorNorm>::relativeVelocities(std::v
               }
           }
       }
-
+    */
 
 }