diff --git a/src/scheme/CellbyCellLimitation.hpp b/src/scheme/CellbyCellLimitation.hpp
index 97d11a07e035595735861a173df6210179b4f2d8..a443e71ef6781058e48c093516c2f2b5aaf89a68 100644
--- a/src/scheme/CellbyCellLimitation.hpp
+++ b/src/scheme/CellbyCellLimitation.hpp
@@ -149,7 +149,7 @@ class CellByCellLimitation
     const DiscreteFunctionP0<const double>& epsilon,
     // const DiscreteFunctionDPk<Dimension, double>&
     auto epsilon_R(const CellId cell_id, const Rd& x),
-    CellValue<const double>& lambda_epsilon)   // const
+    CellValue<double>& lambda_epsilon)   // const
   //                                DiscreteFunctionDPk<Dimension, double>& epsilon_R) const
   {
     const auto& cell_to_face_matrix = mesh.connectivity().cellToFaceMatrix();
@@ -171,7 +171,7 @@ class CellByCellLimitation
     // QuadratureManager::instance().getLineFormula(GaussLegendreQuadratureDescriptor(m_quadrature_degree));
 
     parallel_for(
-      mesh.numberOfCells(), PUGS_LAMBDA(const CellId cell_id) {
+      mesh.numberOfCells(), PUGS_LAMBDA(CellId cell_id) {
         const double epsilonj = epsilon[cell_id];
 
         double epsilon_min = epsilonj;