diff --git a/src/geodesicfeassembler.hh b/src/geodesicfeassembler.hh
index ce2b4519c384518c9aacba6a44611a989a29fe8e..f657304672684eb7056261c5299c637e59b2013b 100644
--- a/src/geodesicfeassembler.hh
+++ b/src/geodesicfeassembler.hh
@@ -34,8 +34,10 @@ class GeodesicFEAssembler {
 public:
     
     /** \brief Constructor for a given grid */
-    GeodesicFEAssembler(const GridView& gridView) : 
-        gridView_(gridView)
+    GeodesicFEAssembler(const GridView& gridView,
+                        LocalGeodesicFEStiffness<GridView,TargetSpace>* localStiffness)
+        : gridView_(gridView),
+          localStiffness_(localStiffness)
     {}
     
     /** \brief Assemble the tangent stiffness matrix
diff --git a/src/rodassembler.hh b/src/rodassembler.hh
index 8928383cf5747e8a27ef578dfd184b9f07641e4c..407a52a78e40e5fc2290b449483cdefa1c1fd7a4 100644
--- a/src/rodassembler.hh
+++ b/src/rodassembler.hh
@@ -9,6 +9,7 @@
 #include <dune/ag-common/boundarypatch.hh>
 
 #include "rigidbodymotion.hh"
+#include "rodlocalstiffness.hh"
 #include "geodesicfeassembler.hh"
 
 /** \brief The FEM operator for an extensible, shearable rod
@@ -47,9 +48,11 @@ class RodAssembler : public GeodesicFEAssembler<typename GridType::LeafGridView,
     public:
         
         //! ???
-        RodAssembler(const GridType &grid) : 
-            GeodesicFEAssembler<typename GridType::LeafGridView, RigidBodyMotion<3> >(grid.leafView()),
-            grid_(&grid)
+    RodAssembler(const GridType &grid,
+                 RodLocalStiffness<typename GridType::LeafGridView,double>* localStiffness) : 
+        GeodesicFEAssembler<typename GridType::LeafGridView, RigidBodyMotion<3> >(grid.leafView(),
+                                                                                  localStiffness),
+        grid_(&grid)
         { 
             // Set dummy material parameters
             K_[0] = K_[1] = K_[2] = 1;