diff --git a/dumux/material/constraintsolvers/compositionfromfugacities.hh b/dumux/material/constraintsolvers/compositionfromfugacities.hh
index 419ca4ab6592eeb65148c9f1573c9d7b32acd53d..cd5aad58ff79ebe6a9bcbe027a817dae60d23a3b 100644
--- a/dumux/material/constraintsolvers/compositionfromfugacities.hh
+++ b/dumux/material/constraintsolvers/compositionfromfugacities.hh
@@ -106,9 +106,9 @@ public:
             solveIdealMix_(fluidState, paramCache, phaseIdx, targetFug);
             return;
         }
-
+#if !DUNE_VERSION_NEWER(DUNE_COMMON,2,6)
         Dune::FMatrixPrecision<Scalar>::set_singular_limit(1e-25);
-
+#endif
         // save initial composition in case something goes wrong
         Dune::FieldVector<Scalar, numComponents> xInit;
         for (int i = 0; i < numComponents; ++i) {
diff --git a/dumux/material/constraintsolvers/immiscibleflash.hh b/dumux/material/constraintsolvers/immiscibleflash.hh
index 3f2aa31884e14ebf44c228bad56a329c8b333072..93183617fc6d9eb3dc03dc3b01a139a5d30c5efd 100644
--- a/dumux/material/constraintsolvers/immiscibleflash.hh
+++ b/dumux/material/constraintsolvers/immiscibleflash.hh
@@ -120,7 +120,9 @@ public:
                       const typename MaterialLaw::Params &matParams,
                       const ComponentVector &globalMolarities)
     {
+#if !DUNE_VERSION_NEWER(DUNE_COMMON,2,6)
         Dune::FMatrixPrecision<Scalar>::set_singular_limit(1e-25);
+#endif
 
         /////////////////////////
         // Check if all fluid phases are incompressible
diff --git a/dumux/material/constraintsolvers/ncpflash.hh b/dumux/material/constraintsolvers/ncpflash.hh
index 4ffcc557484e360102fcf52c51c5f957fe0fcf4d..4bcb1d7d95d77bb0facee04ce05d4ce8de22184f 100644
--- a/dumux/material/constraintsolvers/ncpflash.hh
+++ b/dumux/material/constraintsolvers/ncpflash.hh
@@ -146,7 +146,9 @@ public:
                       const typename MaterialLaw::Params &matParams,
                       const ComponentVector &globalMolarities)
     {
+#if !DUNE_VERSION_NEWER(DUNE_COMMON,2,6)
         Dune::FMatrixPrecision<Scalar>::set_singular_limit(1e-25);
+#endif
 
         /////////////////////////
         // Newton method