#759 add some more comments

Author: martinjrobins

pybamm/solvers/base_solver.py

@@ -459,6 +459,9 @@ def solve(self, model, t_eval, external_variables=None, inputs=None):
                                     dtime + sys.float_info.epsilon])
         end_indices.append(len(t_eval))
 
+        # integrate separatly over each time segment and accumulate into the solution
+        # object, restarting the solver at each discontinuity (and recalculating a
+        # consistent state afterwards if a dae)
         old_y0 = model.y0
         solution = None
         for start_index, end_index in zip(start_indices, end_indices):