“…At this point, we use ( 40) and obtain Unfortunately, in practice, this fixed point method converges only for very small values of ∆t. To enhance the stability region, and then to allow to take larger values of ∆t, we use the ∆ κ acceleration procedure introduced in [11] and applied in [1,2,19] for Allen-Cahn's, weakly damped Schrödinger and BBM equations respectively. In two words, the ∆ κ procedure consists in replacing the Picard iterates by…”