diff --git a/examples/co2injection_flash_ecfv.cpp b/examples/co2injection_flash_ecfv.cpp
index 2f0564657..3fd5e8cb5 100644
--- a/examples/co2injection_flash_ecfv.cpp
+++ b/examples/co2injection_flash_ecfv.cpp
@@ -50,6 +50,10 @@ SET_TYPE_PROP(
 // else we increase the tolerance of the Newton solver
 #if HAVE_QUAD
 SET_TYPE_PROP(Co2InjectionFlashEcfvProblem, Scalar, quad);
+
+// the default linear solver used for this problem (-> AMG) cannot be used with quadruple
+// precision scalars... (this seems to only apply to Dune >= 2.4)
+SET_TAG_PROP(Co2InjectionFlashEcfvProblem, LinearSolverSplice, ParallelIterativeLinearSolver);
 #else
 SET_SCALAR_PROP(Co2InjectionFlashEcfvProblem, NewtonRawTolerance, 1e-5);
 #endif
diff --git a/examples/co2injection_flash_ni_ecfv.cpp b/examples/co2injection_flash_ni_ecfv.cpp
index 972a2bdf2..b30946b71 100644
--- a/examples/co2injection_flash_ni_ecfv.cpp
+++ b/examples/co2injection_flash_ni_ecfv.cpp
@@ -49,6 +49,10 @@ SET_TYPE_PROP(
 // else we increase the tolerance of the Newton solver
 #if HAVE_QUAD
 SET_TYPE_PROP(Co2InjectionFlashNiEcfvProblem, Scalar, quad);
+
+// the default linear solver used for this problem (-> AMG) cannot be used with quadruple
+// precision scalars... (this seems to only apply to Dune >= 2.4)
+SET_TAG_PROP(Co2InjectionFlashNiEcfvProblem, LinearSolverSplice, ParallelIterativeLinearSolver);
 #else
 SET_SCALAR_PROP(Co2InjectionFlashNiEcfvProblem, NewtonRawTolerance, 1e-5);
 #endif
diff --git a/examples/co2injection_flash_ni_vcfv.cpp b/examples/co2injection_flash_ni_vcfv.cpp
index ed9b720a5..6e62c282e 100644
--- a/examples/co2injection_flash_ni_vcfv.cpp
+++ b/examples/co2injection_flash_ni_vcfv.cpp
@@ -49,6 +49,10 @@ SET_TYPE_PROP(
 // else we increase the tolerance of the Newton solver
 #if HAVE_QUAD
 SET_TYPE_PROP(Co2InjectionFlashNiVcfvProblem, Scalar, quad);
+
+// the default linear solver used for this problem (-> AMG) cannot be used with quadruple
+// precision scalars... (this seems to only apply to Dune >= 2.4)
+SET_TAG_PROP(Co2InjectionFlashNiVcfvProblem, LinearSolverSplice, ParallelIterativeLinearSolver);
 #else
 SET_SCALAR_PROP(Co2InjectionFlashNiVcfvProblem, NewtonRawTolerance, 1e-5);
 #endif
diff --git a/examples/co2injection_flash_vcfv.cpp b/examples/co2injection_flash_vcfv.cpp
index a39f4996d..70af65c32 100644
--- a/examples/co2injection_flash_vcfv.cpp
+++ b/examples/co2injection_flash_vcfv.cpp
@@ -50,6 +50,10 @@ SET_TYPE_PROP(
 // else we increase the tolerance of the Newton solver
 #if HAVE_QUAD
 SET_TYPE_PROP(Co2InjectionFlashVcfvProblem, Scalar, quad);
+
+// the default linear solver used for this problem (-> AMG) cannot be used with quadruple
+// precision scalars... (this seems to only apply to Dune >= 2.4)
+SET_TAG_PROP(Co2InjectionFlashVcfvProblem, LinearSolverSplice, ParallelIterativeLinearSolver);
 #else
 SET_SCALAR_PROP(Co2InjectionFlashVcfvProblem, NewtonRawTolerance, 1e-5);
 #endif