diff --git a/dune/tectonic/pointtractionboundaryassembler.hh b/dune/tectonic/pointtractionboundaryassembler.hh
index 648ec7e2ae9f42e521b92fba4539853edfcf3434..c64b92276c80f158f78af3e68faa168a2129fca1 100644
--- a/dune/tectonic/pointtractionboundaryassembler.hh
+++ b/dune/tectonic/pointtractionboundaryassembler.hh
@@ -39,7 +39,7 @@ class PointTractionBoundaryAssembler
   template <class TrialLocalFE, class BoundaryIterator>
   void assemble(const BoundaryIterator& it, LocalVector& localVector,
                 const TrialLocalFE& tFE) {
-
+    auto const inside = it->inside();
     localVector = 0.0;
 
     // geometry of the boundary face
@@ -47,9 +47,8 @@ class PointTractionBoundaryAssembler
         it->geometry();
 
     // get quadrature rule
-    const int order = (it->inside()->type().isSimplex())
-                          ? 2 * (order_ - 1)
-                          : 2 * (order_ * dim - 1);
+    const int order =
+        (inside.type().isSimplex()) ? 2 * (order_ - 1) : 2 * (order_ * dim - 1);
 
     // get quadrature rule
     const Dune::QuadratureRule<ctype, dim - 1>& quad =
@@ -75,8 +74,8 @@ class PointTractionBoundaryAssembler
 
       // evaluate the displacement gradient at quad point of the element
       typename GridFunction::DerivativeType localDispGradient;
-      if (displacement_->isDefinedOn(*it->inside()))
-        displacement_->evaluateDerivativeLocal(*it->inside(), elementQuadPos,
+      if (displacement_->isDefinedOn(inside))
+        displacement_->evaluateDerivativeLocal(inside, elementQuadPos,
                                                localDispGradient);
       else
         displacement_->evaluateDerivative(segmentGeometry.global(quadPos),