From 2befba239de3c911e97594e864a9d8ef624e107a Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 14 Oct 2014 15:20:54 +0000
Subject: [PATCH] Abort when instantiating a RiemannianTrustRegionSolver class
with a P3 basis
The global indexing infrastructure does not currently support P3 spaces.
[[Imported from SVN: r9924]]
---
dune/gfe/riemanniantrsolver.hh | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/dune/gfe/riemanniantrsolver.hh b/dune/gfe/riemanniantrsolver.hh
index e5f682f7..98e3869f 100644
--- a/dune/gfe/riemanniantrsolver.hh
+++ b/dune/gfe/riemanniantrsolver.hh
@@ -39,7 +39,9 @@ class RiemannianTrustRegionSolver
typedef Dune::BlockVector<Dune::FieldVector<field_type, blocksize> > CorrectionType;
typedef std::vector<TargetSpace> SolutionType;
-#ifdef SECOND_ORDER
+#ifdef THIRD_ORDER
+#error RiemannianTrustRegionSolver cannot currently be used for third-order spaces
+#elif defined SECOND_ORDER
typedef Dune::GlobalP2Mapper<typename GridType::LeafGridView> GUIndex;
#else
typedef GlobalUniqueIndex<typename GridType::LeafGridView, gridDim> GUIndex;
--
GitLab