diff --git a/dune/gdt/test/stationary-eocstudies/diffusion-ipdg.hh b/dune/gdt/test/stationary-eocstudies/diffusion-ipdg.hh
index e3186781447fecb556bb1d6bc22ef5e63b01bf6c..90ea966d884aedbf14731bc56b245c50e9c734ea 100644
--- a/dune/gdt/test/stationary-eocstudies/diffusion-ipdg.hh
+++ b/dune/gdt/test/stationary-eocstudies/diffusion-ipdg.hh
@@ -27,9 +27,9 @@
 #include <dune/gdt/functionals/vector-based.hh>
 #include <dune/gdt/local/functionals/integrals.hh>
 #include <dune/gdt/local/bilinear-forms/integrals.hh>
-#include <dune/gdt/local/integrands/elliptic-ipdg.hh>
+#include <dune/gdt/local/integrands/laplace.hh>
+#include <dune/gdt/local/integrands/laplace-ipdg.hh>
 #include <dune/gdt/local/integrands/ipdg.hh>
-#include <dune/gdt/local/integrands/elliptic.hh>
 #include <dune/gdt/local/integrands/product.hh>
 #include <dune/gdt/local/integrands/conversion.hh>
 #include <dune/gdt/norms.hh>
@@ -265,10 +265,10 @@ protected:
         (space_type_.size() >= 2 && space_type_.substr(0, 2) == "cg") ? Stencil::element
                                                                       : Stencil::element_and_intersection));
     // - volume term
-    lhs_op->append(LocalElementIntegralBilinearForm<E>(LocalEllipticIntegrand<E>(one, this->diffusion())));
+    lhs_op->append(LocalElementIntegralBilinearForm<E>(LocalLaplaceIntegrand<E>(this->diffusion())));
     // - inner faces
     lhs_op->append(LocalIntersectionIntegralBilinearForm<I>(
-                       LocalEllipticIPDGIntegrands::InnerCoupling<I>(1., this->diffusion(), this->diffusion())),
+                       LocalLaplaceIPDGIntegrands::InnerCoupling<I>(1., this->diffusion(), this->diffusion())),
                    {},
                    XT::Grid::ApplyOn::InnerIntersectionsOnce<GV>());
     lhs_op->append(
@@ -279,7 +279,7 @@ protected:
     // - Dirichlet faces
     lhs_op->append(
         LocalIntersectionIntegralBilinearForm<I>(
-            LocalEllipticIPDGIntegrands::DirichletCoupling<I>(1., this->diffusion())),
+            LocalLaplaceIPDGIntegrands::DirichletCoupling<I>(1., this->diffusion())),
         {},
         XT::Grid::ApplyOn::CustomBoundaryIntersections<GV>(this->boundary_info(), new XT::Grid::DirichletBoundary()));
     lhs_op->append(
diff --git a/dune/gdt/tools/grid-quality-estimates.hh b/dune/gdt/tools/grid-quality-estimates.hh
index 48691f1ac5beb570a8b69f570196930538d10dbc..0c1b476389c352482d8022a9ada692ff527152a3 100644
--- a/dune/gdt/tools/grid-quality-estimates.hh
+++ b/dune/gdt/tools/grid-quality-estimates.hh
@@ -24,7 +24,7 @@
 #include <dune/xt/grid/type_traits.hh>
 
 #include <dune/gdt/local/bilinear-forms/integrals.hh>
-#include <dune/gdt/local/integrands/elliptic.hh>
+#include <dune/gdt/local/integrands/laplace.hh>
 #include <dune/gdt/local/integrands/product.hh>
 #include <dune/gdt/spaces/interface.hh>
 
@@ -43,9 +43,9 @@ double estimate_inverse_inequality_constant(const SpaceInterface<GV, r>& space)
     basis->bind(element);
     const double h = XT::Grid::diameter(element);
     auto H1_product_matrix = XT::LA::convert_to<XT::LA::CommonDenseMatrix<double>>(
-        LocalElementIntegralBilinearForm<E, r>(LocalEllipticIntegrand<E, r>(1.)).apply2(*basis, *basis));
+        LocalElementIntegralBilinearForm<E, r>(LocalLaplaceIntegrand<E, r>()).apply2(*basis, *basis));
     auto L2_product_matrix = XT::LA::convert_to<XT::LA::CommonDenseMatrix<double>>(
-        LocalElementIntegralBilinearForm<E, r>(LocalElementProductIntegrand<E, r>(1.)).apply2(*basis, *basis));
+        LocalElementIntegralBilinearForm<E, r>(LocalElementProductIntegrand<E, r>()).apply2(*basis, *basis));
     auto evs =
         XT::LA::make_generalized_eigen_solver(H1_product_matrix,
                                               L2_product_matrix,