From 03faf2e067b23191803d05fdd3e74a2664f35179 Mon Sep 17 00:00:00 2001
From: Carl-Martin Pfeiler <carl-martin.pfeiler@asc.tuwien.ac.at>
Date: Wed, 16 Jan 2019 11:04:52 +0100
Subject: [PATCH] integrator tps2ab + tps2ee in first step

---
 integrators/__init__.py              |  3 +++
 integrators/tps2ab_initial_tps2ee.py | 25 +++++++++++++++++++++++++
 2 files changed, 28 insertions(+)
 create mode 100644 integrators/tps2ab_initial_tps2ee.py

diff --git a/integrators/__init__.py b/integrators/__init__.py
index f8d59cd..2a28c5e 100644
--- a/integrators/__init__.py
+++ b/integrators/__init__.py
@@ -27,6 +27,8 @@ from commics.integrators.kw1_mps import KW1_MPS
 from commics.integrators.tps1_kw2 import TPS1_KW2
 from commics.integrators.tps1_mps import TPS1_MPS
 
+from commics.integrators.tps2ab_initial_tps2ee import TPS2AB_INITIAL_TPS2EE
+
 __all__ = \
   [ "Integrator" \
   , "DUMMY" \
@@ -36,4 +38,5 @@ __all__ = \
   , "TPS1HEE" \
   , "TPS1HEE_MPS" \
   , "KW1_MPS", "TPS1_KW2", "TPS1_MPS" \
+  , "TPS2AB_INITIAL_TPS2EE"
   ]
diff --git a/integrators/tps2ab_initial_tps2ee.py b/integrators/tps2ab_initial_tps2ee.py
new file mode 100644
index 0000000..2cb36ba
--- /dev/null
+++ b/integrators/tps2ab_initial_tps2ee.py
@@ -0,0 +1,25 @@
+###  TPS2+AB + TPS2EE for first time step ###
+from . import tps2ab
+
+
+################################################################################
+
+
+class TPS2AB_INITIAL_TPS2EE(tps2ab.TPS2AB):
+  '''
+  implements TPS2AB from second time-step on,
+  BUT: At first time step use TPS2EE = EE approximation of pi, i.e., pi(v) := 0
+  '''
+
+
+#------------------------------------------------------------------------------#
+
+
+  def _DerivLinearTerms_pi_Implicit_RHS(self):
+    return self._ZeroCF()
+
+
+#------------------------------------------------------------------------------#
+
+
+################################################################################
-- 
GitLab