Bug Fixed in const template in limitor expression

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;