diff --git a/CHANGELOG.md b/CHANGELOG.md
index 23d757975e52f53c604e3b8c11f6f0ae4d5ae516..f8933304a325ee89ae5d92c32e39335b9fbe4864 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -4,6 +4,9 @@
 - Introduce class `LocalIntegralEnergy` to work with the densities
 - Local energies and FE assemblers use now `dune-functions` power bases instead of scalar `dune-fufem` bases; the key element is the `LocalView` which contains the information for each element
 - Introduce class `LocalHyperDualStiffness` and `hyperdual` to calculate gradient and hessian using hyper-dual numbers
+- The class `FEAssembler` now stores a `const` reference to the basis, rather than a value.
+  This allows to use the assembler with bases that are not copyable
+  (as, e.g., the `RaviartThomasBasis`, see [dune-functions issue 58](https://gitlab.dune-project.org/staging/dune-functions/-/issues/58)).
 
 ## Deprecations and removals
 
diff --git a/dune/elasticity/assemblers/feassembler.hh b/dune/elasticity/assemblers/feassembler.hh
index 73ee29616d032bb4e2de33cc4840cf20395c6b39..43a50ac9198a6df0d8f6bfb9fea38d8513b11b73 100644
--- a/dune/elasticity/assemblers/feassembler.hh
+++ b/dune/elasticity/assemblers/feassembler.hh
@@ -33,7 +33,7 @@ class FEAssembler {
 
 public:
 
-    const Basis basis_;
+    const Basis& basis_;
 
     /** \brief Partition type on which to assemble
      *