Search

In this paper we describe the efficient numerical implementation of Fractional HBVMs, a class of methods recently introduced for solving systems of fractional differential equations. The reported arguments are implemented in the Matlab code fhbvm, which is made available on the web. An extensive experimentation of the code is reported, to give evidence of its effectiveness., Comment: 26 pages, 6 figures

In this paper we consider the numerical solution of fractional terminal value problems (FDE-TVPs). In particular, the proposed procedure uses a Newton-type iteration which is particularly efficient when coupled with a recently-introduced step-by-step procedure for solving fractional initial value problems (FDE-IVPs), able to produce spectrally accurate solutions of FDE problems. Some numerical tests are reported to make evidence of its effectiveness., Comment: 23 pages, 4 figures, 7 tables, one typo fixed

We consider the weighted least squares spline approximation of a noisy dataset. By interpreting the weights as a probability distribution, we maximize the associated entropy subject to the constraint that the mean squared error is prescribed to a desired (small) value. Acting on this error yields a robust regression method that automatically detects and removes outliers from the data during the fitting procedure, by assigning them a very small weight. We discuss the use of both spline functions and spline curves. A number of numerical illustrations have been included to disclose the potentialities of the maximal-entropy approach in different application fields., Comment: 11 color pictures

Recently, the efficient numerical solution of Hamiltonian problems has been tackled by defining the class of energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs). Their derivation relies on the expansion of the vector field along the Legendre orthonormal basis. Interestingly, this approach can be extended to cope with other orthonormal bases and, in particular, we here consider the case of the Chebyshev polynomial basis. The corresponding Runge-Kutta methods were previously obtained by Costabile and Napoli [33]. In this paper, the use of a different framework allows us to carry out a novel analysis of the methods also when they are used as spectral formulae in time, along with some generalizations of the methods., Comment: 25 pages, 2 figures, 2 tables

In recent years, the efficient numerical solution of Hamiltonian problems has led to the definition of a class of energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs). Such methods admit an interesting interpretation in terms of continuous-stage Runge-Kutta methods, which is here recalled and revisited for general differential problems., Comment: 20 pages

In this paper we are concerned with numerical methods for the one-sided event location in discontinuous differential problems, whose event function is nonlinear (in particular, of polynomial type). The original problem is transformed into an equivalent Poisson problem, which is effectively solved by suitably adapting a recently devised class of energy-conserving methods for Poisson systems. The actual implementation of the methods is fully discussed, with a particular emphasis to the problem at hand. Some numerical tests are reported, to assess the theoretical findings., Comment: 14 pages, 1 figure, 2 tables

These are the lecture notes (in Italian) of a course held in Perugia, Italy, during the summer 2002. They concern the basic facts on the iterative solution of linear systems. The course is self-contained and requires only basic knowledge of numerical linear algebra. This is a corrected version of the original lecture notes, and are dedicated to the memory of Cecilia Magherini., Comment: 150 pagine. This is the corrected version of the lecture notes (in Italian). In memory of Cecilia Magherini

In this paper we are concerned with energy-conserving methods for Poisson problems, which are effectively solved by defining a suitable generalization of HBVMs, a class of energy-conserving methods for Hamiltonian problems. The actual implementation of the methods is fully discussed, with a particular emphasis on the conservation of Casimirs. Some numerical tests are reported, in order to assess the theoretical findings., Comment: 30 pages, 5 figures, 3 tables

In this paper we discuss a framework for the polynomial approximation to the solution of initial value problems for differential equations. The framework, initially devised for the approximation of ordinary differential equations, is further extended to cope with constant delay differential equations. Relevant classes of Runge-Kutta methods can be derived within this framework., Comment: 35 pages, 4 figures, 2 tables

We implement the PaperRank and AuthorRank indices introduced in [Amodio & Brugnano, 2014] in the Scopus database, in order to highlight quantitative and qualitative information that the bare number of citations and/or the h-index of an author are unable to provide. In addition to this, the new indices can be cheaply updated in Scopus, since this has a cost comparable to that of updating the number of citations. Some examples are reported to provide insight in their potentialities, as well as possible extensions., Comment: 20 pages, 3 figures, 4 tables. Accepted version on Journal of Infometrics

Background. The paper concerns the SARS-CoV2 (COVID-19) pandemic that, starting from the end of February 2020, began spreading along the Italian peninsula, by first attacking small communities in north regions, and then extending to the center and south of Italy, including the two main islands. Objective. The creation of a forecast model that manages to alert the decision-making bodies and, in particular, the healthcare system, to hinder the emergence of any other pandemic outbreaks, or the arrival of subsequent pandemic waves. Methods. A new mathematical model to describe the pandemic is given. The model includes the class of undiagnosed infected people, and has a multi-region extension, to cope with the in-time and in-space heterogeneity of the epidemic. Results. We obtain a robust and reliable tool for the forecast of the total and active cases, which can be also used to simulate different scenarios. Conclusions. We are able to address a number of issues, such as assessing the adoption of the lockdown in Italy, started from 11 March 2020, and how to employ a rapid screening test campaign for containing the epidemic., Comment: 21 pages, 10 figures

In this note, we describe simple generalizations of the basic SIR model for epidemic, in case of a multi-region scenario, to be used for predicting the COVID-19 epidemic spread in Italy., Comment: 7 pages (1 typo fixed)

We devise a variable precision floating-point arithmetic by exploiting the framework provided by the Infinity Computer. This is a computational platform implementing the Infinity Arithmetic system, a positional numeral system which can handle both infinite and infinitesimal quantities symbolized by the positive and negative finite powers of the radix grossone. The computational features offered by the Infinity Computer allows us to dynamically change the accuracy of representation and floating-point operations during the flow of a computation. When suitably implemented, this possibility turns out to be particularly advantageous when solving ill-conditioned problems. In fact, compared with a standard multi-precision arithmetic, here the accuracy is improved only when needed, thus not affecting that much the overall computational effort. An illustrative example about the solution of a nonlinear equation is also presented., Comment: 11 pages, 2 figures, 6 tables

Gyrocenter dynamics of charged particles plays a fundamental role in plasma physics. In particular, accuracy and conservation of energy are important features for correctly performing long-time simulations. For this purpose, we here propose arbitrarily high-order energy conserving methods for its simulation. The analysis and the efficient implementation of the methods are fully described, and some numerical tests are reported., Comment: 23 pages, 4 figures

In this paper, we study the numerical solution of Manakov systems by using a spectrally accurate Fourier decomposition in space, coupled with a spectrally accurate time integration. This latter relies on the use of spectral Hamiltonian boundary Value Methods. The used approach allows to conserve all the physical invariants of the systems. Some numerical tests are reported., Comment: 29 pages (typo in equation (20) fixed)

In recent years, the class of energy-conserving methods named Hamiltonian Boundary Value Methods (HBVMs) has been devised for numerically solving Hamiltonian problems. In this short note, we study their natural formulation as continuous-stage Runge-Kutta(-Nystr\"om) methods, which allows a deeper insight in the methods., Comment: 12 pages, accepted manuscript

In this paper, we report about recent findings in the numerical solution of Hamiltonian Partial Differential Equations (PDEs), by using energy-conserving line integral methods in the Hamiltonian Boundary Value Methods (HBVMs) class. In particular, we consider the semilinear wave equation, the nonlinear Schr\"odinger equation, and the Korteweg-de Vries equation, to illustrate the main features of this novel approach., Comment: 33 pages, 3 figures, 3 tables

Recently, the numerical solution of stiffly/highly-oscillatory Hamiltonian problems has been attacked by using Hamiltonian Boundary Value Methods (HBVMs) as spectral methods in time. While a theoretical analysis of this spectral approach has been only partially addressed, there is enough numerical evidence that it turns out to be very effective even when applied to a wider range of problems. Here we fill this gap by providing a thorough convergence analysis of the methods and confirm the theoretical results with the aid of a few numerical tests., Comment: 18 pages, 1 figure

Recently, the numerical solution of multi-frequency, highly-oscillatory Hamiltonian problems has been attacked by using Hamiltonian Boundary Value Methods (HBVMs) as spectral methods in time. When the problem derives from the space semi- discretization of (possibly Hamiltonian) partial differential equations (PDEs), the resulting problem may be stiffly-oscillatory, rather than highly-oscillatory. In such a case, a different implementation of the methods is needed, in order to gain the maximum efficiency., Comment: 17 pages, 3 figures

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., Comment: 25 pages, 4 figures

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., Comment: 28 pages, 2 figures, 4 tables

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., Comment: 30 pages, 3 figures, 4 tables

In this paper we discuss energy conservation issues related to the numerical solution of the nonlinear wave equation, when a Fourier expansion is considered for the space discretization. The obtained semi-discrete problem is then solved in time by means of energy-conserving Runge-Kutta methods in the HBVMs class., Comment: 25 pages, 7 figure

In this paper we discuss energy conservation issues related to the numerical solution of the nonlinear wave equation. As is well known, this problem can be cast as a Hamiltonian system that may be autonomous or not, depending on the specific boundary conditions at hand. We relate the conservation properties of the original problem to those of its semi-discrete version obtained by the method of lines. Subsequently, we show that the very same properties can be transferred to the solutions of the fully discretized problem, obtained by using energy-conserving methods in the HBVMs (Hamiltonian Boundary Value Methods) class. Similar arguments hold true for different types of Hamiltonian Partial Differential Equations, e.g., the nonlinear Schr\"odinger equation., Comment: 41 pages, 11 figure

Recently, the class of energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs), has been proposed for the efficient solution of Hamiltonian problems, as well as for other types of conservative problems. In this paper, we report further advances concerning such methods, resulting in their enhanced version (Enhanced HBVMs, or EHBVMs). The basic theoretical results are sketched, along with a few numerical tests on a Hamiltonian problem, taken from the literature, possessing multiple invariants., Comment: 4 pages

We here investigate the efficient implementation of the energy-conserving methods named Hamiltonian Boundary Value Methods (HBVMs) recently introduced for the numerical solution of Hamiltonian problems. In this note, we describe an iterative procedure, based on a triangular splitting, for solving the generated discrete problems, when the problem at hand is separable., Comment: 4 pages

In this paper we define an efficient implementation for the family of low-rank energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs), recently defined in the last years. The proposed implementation relies on the particular structure of the Butcher matrix defining such methods, for which we can derive an efficient splitting procedure. The very same procedure turns out to be automatically suited for the efficient implementation of Gauss-Legendre collocation methods, since these methods are a special instance of HBVMs. The linear convergence analysis of the splitting procedure exhibits excellent properties, which are confirmed by a few numerical tests., Comment: 19 pages, 7 tables, 2 figures

In a recent series of papers, the class of energy-conserving Runge-Kutta methods named Hamiltonian BVMs (HBVMs) has been defined and studied. Such methods have been further generalized for the efficient solution of general conservative problems, thus providing the class of Line Integral Methods (LIMs). In this paper we derive a further extension, which we name Enhanced Line Integral Methods (ELIMs), more tailored for Hamiltonian problems, allowing for the conservation of multiple invariants of the continuous dynamical system. The analysis of the methods is fully carried out and some numerical tests are reported, in order to confirm the theoretical achievements., Comment: 21 pages, 4 figures, 3 tables

These are the lecture notes of a course given by the first author on December 27, 2012 - January 4, 2013, held at the Academy of Mathematics and Systems Science Chinese Academy of Sciences in Beijing., Comment: 68 pages