From cb467ee59f5a4de05fd8412106366e1cc30607b7 Mon Sep 17 00:00:00 2001
From: Carl-Martin Pfeiler <carl-martin.pfeiler@asc.tuwien.ac.at>
Date: Fri, 8 Feb 2019 16:16:51 +0100
Subject: [PATCH] TUNE TPS2; f) Use initial guess x0 in gmres
---
integrators/tps2.py | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/integrators/tps2.py b/integrators/tps2.py
index ad55678..5af8bd5 100644
--- a/integrators/tps2.py
+++ b/integrators/tps2.py
@@ -234,9 +234,10 @@ class TPS2(_tpsx.TPSX, \
RHS = self._Q.dot( f )
- w, succ = scipy.sparse.linalg.gmres(self._LHS, RHS , x0=None \
+ self._UpdateGuess()
+ self._solution, succ = scipy.sparse.linalg.gmres(self._LHS, RHS , x0=self._guess \
, tol=self._solvetol, maxiter=4000, M=self._preconditioner)
- self._Mag.v.gf.vec.FV().NumPy()[:] = self._Q.transpose().dot(w)
+ self._Mag.v.gf.vec.FV().NumPy()[:] = self._Q.transpose().dot(self._solution)
diff_gf.vec.FV().NumPy()[:] = self._Mag.v.gf.vec.FV().NumPy()-oldIterate
itError = sqrt( Integrate(diff_cf*diff_cf, self._GetMesh()) )
--
GitLab