BUG in KW1 fixed. (also influences KW2(AB))

f7447355 · Carl-Martin Pfeiler · af7edc51 · f7447355
Commit f7447355 authored 5 years ago by Carl-Martin Pfeiler
--- a/integrators/kw1.py
+++ b/integrators/kw1.py
@@ -199,7 +199,7 @@ class KW1(_integrator._Integrator, \

    ### lf_iter begin
    rhs_iter_cf = -self._DerivLinearTerms_pi()
-    rhs_iter_cf = -self._DerivNonLinearTerms_Pi()
+    rhs_iter_cf -= self._DerivNonLinearTerms_Pi()

    self._f_iter_parts = LinearForm(self._X)
    self._f_iter_parts += SymbolicLFI(rhs_iter_cf.Compile( \