“…Hereafter, we shall assume that u 0 (x) and v 0 (x) are such that the solution of problem (4) and (5) is regular enough, as a periodic function on [a, b], for all t ≥ 0. The numerical solution of (1), (3) or (4) has been developed along different directions, ranging from the pseudo-spectral or splitting approach [13][14][15][16][17][18][19]46], up to finite-difference and finite-element schemes [20][21][22][23][24]47], as well as structure-preserving methods [10,25,26] and energy-preserving methods [27,28]. In particular, [11,12] consider an energy-conserving strategy based on the Hamiltonian boundary value methods (HBVMs) for the "good" Boussinesq and the improved Boussinesq equation, respectively, while a second-order symplectic method preserving the energy and the momentum is considered in [29].…”