diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 7344f59f39e541016e8067daff3bc1f7606540c4..ae34fe82ce8cb690d71b8a9dec839e5a0fb7b8a3 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -242,7 +242,7 @@ int main (int argc, char *argv[]) try
     RodLocalStiffness<RodGridType::LeafGridView,double> rodLocalStiffness(rodGrid.leafView(),
                                                                        rodA, rodJ1, rodJ2, rodE, rodNu);
 
-    RodAssembler<RodGridType::LeafGridView> rodAssembler(rodGrid.leafView(), &rodLocalStiffness);
+    RodAssembler<RodGridType::LeafGridView,3> rodAssembler(rodGrid.leafView(), &rodLocalStiffness);
 
     RiemannianTrustRegionSolver<RodGridType,RigidBodyMotion<3> > rodSolver;
     rodSolver.setup(rodGrid,