Gianmarco Gurioli

Luigi Brugnano, Gianmarco Gurioli, Chengjian Zhang

In this paper we study the geometric solution of the so called "good" Boussinesq equation. This goal is achieved by using a convenient space semi-discretization, able to preserve the corresponding Hamiltonian structure, then using energy-conserving Runge-Kutta methods in the HBVM class for the time integration. Numerical tests are reported, confirming the effectiveness of the proposed method.

Luigi Brugnano, Gianmarco Gurioli, Felice Iavernaro

In this paper we are concerned with the analysis of a class of geometric integrators, at first devised in [14, 18], which can be regarded as an energy-conserving variant of Gauss collocation methods. With these latter they share the property of conserving quadratic first integrals but, in addition, they also conserve the Hamiltonian function itself. We here reformulate the methods in a more convenient way, and propose a more refined analysis than that given in [18] also providing, as a by-product, a practical procedure for their implementation. A thorough comparison with the original Gauss methods is carried out by means of a few numerical tests solving Hamiltonian and Poisson problems.

Luigi Brugnano, Gianmarco Gurioli, Felice Iavernaro et al.

In this paper, we propose a second-order energy-conserving approximation procedure for Hamiltonian systems with holonomic constraints. The derivation of the procedure relies on the use of the so-called line integral framework. We provide numerical experiments to illustrate theoretical findings.

Luigi Brugnano, Gianmarco Gurioli, Felice Iavernaro et al.

The efficient numerical solution of fractional differential equations has been recently tackled through the definition of Fractional HBVMs (FHBVMs), a class of Runge-Kutta type methods. Corresponding Matlab (c) codes have been also made available on the internet, proving to be very competitive w.r.t. existing ones. However, so far, FHBVMs have been given for solving systems of fractional differential equations with the same order of fractional derivative, whereas the numerical solution of multi-order problems (i.e., problems in which different orders of fractional derivatives occur) has not been handled, yet. Due to their relevance in applications, in this paper we propose an extension of FHBVMs for addressing fractional multi-order problems, providing full details for such an approach. A corresponding Matlab (c) code, handling the case of two different fractional orders, is also made available, proving very effective for numerically solving these problems.