diff --git a/examples/CuboidMesh.cc b/examples/CuboidMesh.cc
index 9348926af3891efe2f5be4cee553bb0e2bd84664..2b5336a71ed57b5a5d7407bbe49f4c5b1b6fa97b 100644
--- a/examples/CuboidMesh.cc
+++ b/examples/CuboidMesh.cc
@@ -60,7 +60,8 @@ constexpr auto dof_per_vertex = Mesh::size_type{3};
 
 // First, we select a Hookean (linear elastic) material model in 3D.
 
-using MaterialModel = materialmodels::Hookean<dimension>;
+using MaterialModel =
+    materialmodels::Hookean<dimension, Vector::value_type, Vector::real_type>;
 
 // We'll choose the Hexa8 shape functions (8 shape functions) that are defined
 // in 3D reference space.
@@ -80,11 +81,14 @@ using Quadrature = quadrature::Quadrature<quadrature::QuadratureType::Cube,
 // on the selected quadrature rule and the embedding. In the case of the
 // isoparametric embedding, this embedding is specified by the "Shape".
 
-using Integrator = integrator::IsoparametricIntegrator<Shape, Quadrature>;
+using Integrator =
+    integrator::IsoparametricIntegrator<Shape, Quadrature, Vector::value_type,
+                                        Vector::real_type>;
 
 // Now we are ready to select our element type.
 
-using Element = elements::CoreElement<MaterialModel, Integrator>;
+using Element = elements::CoreElement<MaterialModel, Integrator,
+                                      Vector::value_type, Vector::real_type>;
 
 // We will assemble e.g. energy using a collection of elements. This is done by
 // the assembler. (The list DefaultFeaturePlugins contain the features (e.g.