From 07a668a0b4bb2015592f980b95055d3485ba6a33 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 18 Jul 2014 12:18:37 +0000
Subject: [PATCH] Use Gauss-Seidel instead of IPOpt as the base solver
Gauss-Seidel converges equally good/bad as IPOpt, but it is much faster.
[[Imported from SVN: r9839]]
---
dune/gfe/riemanniantrsolver.cc | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 65347730..130704e6 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -71,13 +71,12 @@ setup(const GridType& grid,
// Create a multigrid solver
// ////////////////////////////////
-#ifdef HAVE_IPOPT
+#if 0//def HAVE_IPOPT
// First create an IPOpt base solver
QuadraticIPOptSolver<MatrixType, CorrectionType>* baseSolver = new QuadraticIPOptSolver<MatrixType,CorrectionType>;
baseSolver->verbosity_ = NumProc::QUIET;
baseSolver->tolerance_ = baseTolerance;
#else
-#warning IPOpt not installed -- falling back onto a Gauss-Seidel base solver
// First create a Gauss-seidel base solver
TrustRegionGSStep<MatrixType, CorrectionType>* baseSolverStep = new TrustRegionGSStep<MatrixType, CorrectionType>;
--
GitLab