updated readme and todos

Author: pancetta

README.md

@@ -10,6 +10,8 @@ iterators are implemented as well, giving full access to all values at any time.
 News
 ----
 
+* February 2, 2016: the fast-wave slow-wave release, containing all the code for the FWSW-SDC [paper](http://arxiv.org/abs/1602.01626)
+* October 27, 2015: the fault-tolerant release, containing all the code for the FT-PFASST [paper](http://arxiv.org/abs/1510.08334)
 * June 17, 2015: Gray-Scott example using FEniCS is up and running
 * May 7, 2015: added first [FEniCS](http://fenicsproject.org/) example
 * May 5, 2015: added [Clawpack](http://www.clawpack.org/) examples using Sharpclaw

ToDo.md

@@ -20,4 +20,10 @@ ToDo list
 
 * fix behavior of hooks during PFASST runs: initialization is called per level, but some things may need to be done only once (file openings, figures)
 
-* add temporal coarsening
+* add temporal coarsening
+
+* remove mass_matrix sweeper, since this can all be done within solve_system 
+
+* full weak formulation for finite element examples?
+
+* one matrix to rule them all: can we extract the very costly matrix buildups? 

examples/boussinesq_2d_imex/rungmrescounter.py

@@ -99,27 +99,27 @@
     dirkp = dirk(P, np.min([4,dirk_order]))
     u0 = uinit.values.flatten()
     udirk = np.copy(u0)
-    print "Running DIRK ...."
+    print("Running DIRK ....")
     for i in range(0,Nsteps):
       udirk = dirkp.timestep(udirk, dt)  
 
     rkimex = rk_imex(P, dirk_order)
     uimex  = np.copy(u0)
     dt_imex = dt
-    print "Running RK-IMEX ...."
+    print("Running RK-IMEX ....")
     for i in range(0,Nsteps):
       uimex = rkimex.timestep(uimex, dt_imex)
 
     # call main function to get things done...
-    print "Running SDC..."
+    print("Running SDC...")
     uend,stats = mp.run_pfasst(MS,u0=uinit,t0=t0,dt=dt,Tend=Tend)
 
     # For reference solution, increase GMRES tolerance
     P.gmres_tol_limit = 1e-10
     rkimexref = rk_imex(P, 5)
     uref      = np.copy(u0)
     dt_ref    = dt/10.0
-    print "Running RK-IMEX reference...."
+    print("Running RK-IMEX reference....")
     for i in range(0,10*Nsteps):
       uref = rkimexref.timestep(uref, dt_ref)
 

pySDC/sweeper_classes/mass_matrix_imex.py

@@ -95,7 +95,7 @@ def update_nodes(self):
             rhs = P.dtype_u(integral[m])
             for j in range(m+1):
                 rhs += L.dt*(self.QI[m+1,j]*L.f[j].impl + self.QE[m+1,j]*L.f[j].expl)
-            # apply mass matrix to rhs
+            # apply mass matrix to rhs #TODO: Hloy crap, I can put this into solve_system right???
             rhs = P.apply_mass_matrix(rhs)
 
             # implicit solve with prefactor stemming from QI