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: 21 pages, 3 figures

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

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

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 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 exploit the idea to use the maximal-entropy method, successfully tested in information theory and statistical thermodynamics, to determine approximating function's coefficients and squared errors' weights simultaneously as output of one single problem in least-square approximation. We provide evidence of the method's capabilities and performance through its application to representative test cases by working with polynomials as a first step. We conclude by formulating suggestions for future work to improve the version of the method we present in this paper., Comment: 15 pages, 1 table, 10 figures. Accepted in Mathematical Methods in the Applied Sciences. Changes: fixed typo in Eq. (27)

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 the present paper we analyze and discuss some mathematical aspects of the fluid-static configurations of a self-gravitating perfect gas enclosed in a spherical solid shell. The mathematical model we consider is based on the well-known Lane-Emden equation, albeit under boundary conditions that differ from those usually assumed in the astrophysical literature. The existence of multiple solutions requires particular attention in devising appropriate numerical schemes apt to deal with and catch the solution multiplicity as efficiently and accurately as possible. In sequence, we describe some analytical properties of the model, the two algorithms used to obtain numerical solutions, and the numerical results for two selected cases., Comment: Revision accepted for publication in "Communications in Computational Physics"; 23 pages, 9 figures

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