diff --git a/src/averageinterface.hh b/src/averageinterface.hh
index 816f6279e4aa32dd02029f16e046f849444d71c6..ce94d518c88d55bcb5f0d94d6a1f4e4057a084e9 100644
--- a/src/averageinterface.hh
+++ b/src/averageinterface.hh
@@ -264,9 +264,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
 
             const typename NeighborIterator::Geometry& segmentGeometry = nIt.intersectionGlobal();
 
-            const ReferenceElement<double,dim>& refElement = ReferenceElements<double, dim>::general(eIt->type());
-            int nDofs = refElement.size(nIt.numberInSelf(),1,dim);
-
             // Get quadrature rule
             const QuadratureRule<double, dim-1>& quad = QuadratureRules<double, dim-1>::rule(segmentGeometry.type(), dim-1);
 
@@ -295,9 +292,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
 
                 const FieldMatrix<double,dim,dim>& inv = eIt->geometry().jacobianInverseTransposed(quadPos);
                 
-                /* Compute the weight of the current integration point */
-                double weight = quad[ip].weight() * integrationElement;
-                
                 /**********************************************/
                 /* compute gradients of the shape functions   */
                 /**********************************************/
@@ -381,7 +375,7 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
     deformationGradient = U;
     deformationGradient.rightmultiply(VT);
 #endif
-    std::cout << deformationGradient << std::endl;
+    //std::cout << deformationGradient << std::endl;
 
     // deformationGradient now contains the orthogonal part of the polar decomposition
     assert( std::abs(1-deformationGradient.determinant()) < 1e-3);