diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc
index 1d08d09c2506074cf0c0e701feb2572375935102..4a9e486a63b6981f2bc3e42847f1309a73be12c9 100644
--- a/cosserat-continuum.cc
+++ b/cosserat-continuum.cc
@@ -267,17 +267,6 @@ int main (int argc, char *argv[]) try
     std::cout << "Material parameters:" << std::endl;
     materialParameters.report();
 
-#if 0
-    CosseratEnergyLocalStiffness<GridType::LeafGridView,
-                                 P1Basis::LocalFiniteElement,
-                                 3> cosseratEnergyLocalStiffness(materialParameters,
-                                                                 &neumannBoundary,
-                                                                 neumannFunction.get());
-
-    GeodesicFEAssembler<P1Basis,TargetSpace> assembler(grid->leafGridView(),
-                                                                      &cosseratEnergyLocalStiffness);
-
-#else
     // Assembler using ADOL-C
     CosseratEnergyLocalStiffness<GridType::LeafGridView,
                                  FEBasis::LocalFiniteElement,
@@ -290,7 +279,7 @@ int main (int argc, char *argv[]) try
 
     GeodesicFEAssembler<FEBasis,TargetSpace> assembler(grid->leafGridView(),
                                                             &localGFEADOLCStiffness);
-#endif
+
     // /////////////////////////////////////////////////
     //   Create a Riemannian trust-region solver
     // /////////////////////////////////////////////////