diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 2630adadfc7e29a1d11f80e186a3777d956f9303..0b72c666cd5002bb168ac44ec4319664db688966 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -402,7 +402,7 @@ int main (int argc, char *argv[]) try
         VectorType neumannValues(rhs3d.size());
 
         // Using that index 0 is always the left boundary for a uniformly refined OneDGrid
-        computeAveragePressureIPOpt<GridType>(resultantForce, resultantTorque, 
+        computeAveragePressure<GridType>(resultantForce, resultantTorque, 
                                               interfaceBoundary[grid.maxLevel()], 
                                               rodX[0],
                                               neumannValues);