TUNE TPS2; d) Save 1 RHS computation at start

+ is done in fixpoint loop anyhow

integrators/tps2.py

@@ -258,6 +258,12 @@ class TPS2(_tpsx.TPSX, \
     print("itnr = ", itnr, "; Err = ", itError)
 
 
+#------------------------------------------------------------------------------#
+
+
+  def _ComputeRHS(self): pass
+
+
 #------------------------------------------------------------------------------#
 
 

integrators/tps2ab.py

@@ -114,6 +114,9 @@ class TPS2AB(tps2.TPS2, _forms.lf_old.LinearForms_Old):
 
 #------------------------------------------------------------------------------#
 
+  def _ComputeRHS(self):
+    from . import _tpsx
+    _tpsx.TPSX._ComputeRHS(self)
 
   def _Compute_v(self):
     import scipy.sparse.linalg