From 4c7ccc40b550a9428df43eedabdd6a9c02ae47c1 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 23 Mar 2015 09:03:20 +0000
Subject: [PATCH] Make IPOpt use the "mumps" linear solver

Because I am using the Debian version of IPOpt now.  Finally it works!
But the Debian version can only use "mumps", not the default "ma27".

[[Imported from SVN: r10105]]
---
 dune/gfe/trustregionsolver.cc | 1 +
 1 file changed, 1 insertion(+)

diff --git a/dune/gfe/trustregionsolver.cc b/dune/gfe/trustregionsolver.cc
index 50adc5c5..6bb1fe48 100644
--- a/dune/gfe/trustregionsolver.cc
+++ b/dune/gfe/trustregionsolver.cc
@@ -60,6 +60,7 @@ setup(const typename BasisType::GridView::Grid& grid,
     QuadraticIPOptSolver<MatrixType, CorrectionType>* baseSolver = new QuadraticIPOptSolver<MatrixType,CorrectionType>;
     baseSolver->verbosity_ = NumProc::QUIET;
     baseSolver->tolerance_ = baseTolerance;
+    baseSolver->linearSolverType_ = "mumps";
 #else
     // First create a Gauss-seidel base solver
     TrustRegionGSStep<MatrixType, CorrectionType>* baseSolverStep = new TrustRegionGSStep<MatrixType, CorrectionType>;
-- 
GitLab