diff --git a/rod3d.cc b/rod3d.cc
index f7dfc75649d062c209ede50925944584adce705e..9be87d0d49b38649ac4268bb2c369358aef7846d 100644
--- a/rod3d.cc
+++ b/rod3d.cc
@@ -117,7 +117,10 @@ int main (int argc, char *argv[]) try
     // ///////////////////////////////////////////
     //   Create a solver for the rod problem
     // ///////////////////////////////////////////
-    RodAssembler<GridType> rodAssembler(grid);
+
+    RodLocalStiffness<GridType::LeafGridView,double> localStiffness;
+
+    RodAssembler<GridType> rodAssembler(grid, &localStiffness);
     rodAssembler.setShapeAndMaterial(A, J1, J2, E, nu);
 
     RiemannianTrustRegionSolver<GridType,RigidBodyMotion<3> > rodSolver;