diff --git a/examples/diffusion/doc/main.md b/examples/diffusion/doc/main.md
index 28791805adf788e9b36959401d4423087b806ac6..68f2b789ab63d29259a45656378941071d23c301 100644
--- a/examples/diffusion/doc/main.md
+++ b/examples/diffusion/doc/main.md
@@ -139,7 +139,7 @@ We prescribe zero flux over all of the boundary
 
 ```cpp
     NumEqVector neumannAtPos(const GlobalPosition& globalPos) const
-    { return { 0.0 }; }
+    { return NumEqVector(0.0); }
 ```
 
 The diffusion coefficient interface is used in the local residual (see Part 1).
diff --git a/examples/diffusion/doc/model.md b/examples/diffusion/doc/model.md
index b72b22e5c4f08d2d5c2b82311257bde8e5afabd8..29aef2dd638190b8650598d532c9e476249366fd 100644
--- a/examples/diffusion/doc/model.md
+++ b/examples/diffusion/doc/model.md
@@ -123,7 +123,7 @@ F_{K,\sigma} = -D \sum_{B \in \mathcal{B}_K} c_{h,B} \nabla N_B \cdot\boldsymbol
         for (const auto& scv : scvs(fvGeometry))
         {
             const auto& volVars = elemVolVars[scv];
-            // v.axpy(a, w) means v <- v + a*w
+            // v.axpy(a, w) means v += a*w
             gradConcentration.axpy(
                 volVars.priVar(Indices::concentrationIdx),
                 fluxVarCache.gradN(scv.indexInElement())
diff --git a/examples/diffusion/main.cc b/examples/diffusion/main.cc
index 28e4179824c330037dca0555180cd3da70589849..5c188712d277089f2c7b03fb7e1220790ae29ad9 100644
--- a/examples/diffusion/main.cc
+++ b/examples/diffusion/main.cc
@@ -99,7 +99,7 @@ public:
 
     // We prescribe zero flux over all of the boundary
     NumEqVector neumannAtPos(const GlobalPosition& globalPos) const
-    { return { 0.0 }; }
+    { return NumEqVector(0.0); }
 
     // The diffusion coefficient interface is used in the local residual (see Part 1).
     // We can name this interface however we want as long as we adapt the calling site
diff --git a/examples/diffusion/model.hh b/examples/diffusion/model.hh
index 9523db1bfc04904374b7fbc163deddc810925f3d..fa9c0ff4b8d474297a0bc38e06ad56260889afd4 100644
--- a/examples/diffusion/model.hh
+++ b/examples/diffusion/model.hh
@@ -130,7 +130,7 @@ public:
         for (const auto& scv : scvs(fvGeometry))
         {
             const auto& volVars = elemVolVars[scv];
-            // v.axpy(a, w) means v <- v + a*w
+            // v.axpy(a, w) means v += a*w
             gradConcentration.axpy(
                 volVars.priVar(Indices::concentrationIdx),
                 fluxVarCache.gradN(scv.indexInElement())