diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 62e61ac0f8b4fd02be4ceea11c5c8b2590bab5e5..2630adadfc7e29a1d11f80e186a3777d956f9303 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -399,7 +399,6 @@ int main (int argc, char *argv[]) try
         std::cout << "resultant force: " << resultantForce << std::endl;
         std::cout << "resultant torque: " << resultantTorque << std::endl;
 
-#if 1
         VectorType neumannValues(rhs3d.size());
 
         // Using that index 0 is always the left boundary for a uniformly refined OneDGrid
@@ -413,29 +412,6 @@ int main (int argc, char *argv[]) try
                                                         neumannValues,
                                                         rhs3d);
 
-#else
-#ifndef HAVE_LAPACKPP
-#error You need LaPack++ for this!
-#endif
-        // For the time being the Neumann data coming from the rod is a dg function (== not continuous)
-        // Maybe that is not necessary
-        DGIndexSet<GridType> dgIndexSet(grid,grid.maxLevel());
-        dgIndexSet.setup(grid,grid.maxLevel());
-
-        VectorType neumannValues(dgIndexSet.size());
-
-        // Using that index 0 is always the left boundary for a uniformly refined OneDGrid
-        computeAveragePressure<GridType>(resultantForce, resultantTorque, 
-                                              interfaceBoundary[grid.maxLevel()], 
-                                              rodX[0],
-                                              neumannValues);
-
-        rhs3d = 0;
-        assembleAndAddNeumannTerm<GridType, VectorType>(interfaceBoundary[grid.maxLevel()],
-                                                        dgIndexSet,
-                                                        neumannValues,
-                                                        rhs3d);
-#endif
         // ///////////////////////////////////////////////////////////
         //   Solve the Neumann problem for the 3d body
         // ///////////////////////////////////////////////////////////