BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Luigi Brugnano

Luigi Brugnano appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisDigital Librariesmath.DSPopulations and Evolution

Catalog footprint

What is connected

33works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2022arXiv

(Spectral) Chebyshev collocation methods for solving differential equations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Arbitrarily high-order energy-conserving methods for Poisson problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Arbitrary high-order methods for one-sided direct event location in discontinuous differential problems with nonlinear event function

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Continuous-Stage Runge-Kutta approximation to Differential Problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Arbitrarily high-order energy-preserving methods for simulating the gyrocenter dynamics of charged particles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the use of the Infinity Computer architecture to set up a dynamic precision floating-point arithmetic

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient implementation of Gauss collocation and Hamiltonian Boundary Value Methods

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient implementation of geometric integrators for separable Hamiltonian problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Enhanced HBVMs for the numerical solution of Hamiltonian problems with multiple invariants

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multiple invariants conserving Runge-Kutta type methods for Hamiltonian problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Recent advances in bibliometric indexes and the PaperRank problem

Bibliometric indexes are customary used in evaluating the impact of scientific research, even though it is very well known that in different research areas they may range in very different intervals. Sometimes, this is evident even within a single given field of investigation making very difficult (and inaccurate) the assessment of scientific papers. On the other hand, the problem can be recast in the same framework which has allowed to efficiently cope with the ordering of web-pages, i.e., to formulate the PageRank of Google. For this reason, we call such problem the PaperRank problem, here solved by using a similar approach to that employed by PageRank. The obtained solution, which is mathematically grounded, will be used to compare the usual heuristics of the number of citations with a new one here proposed. Some numerical tests show that the new heuristics is much more reliable than the currently used ones, based on the bare number of citations. Moreover, we show that our model improves on recently proposed ones.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Two Step, Fourth Order, Nearly-Linear Method with Energy Preserving Properties

We introduce a family of fourth order two-step methods that preserve the energy function of canonical polynomial Hamiltonian systems. Each method in the family may be viewed as a correction of a linear two-step method, where the correction term is O(h^5) (h is the stepsize of integration). The key tools the new methods are based upon are the line integral associated with a conservative vector field (such as the one defined by a Hamiltonian dynamical system) and its discretization obtained by the aid of a quadrature formula. Energy conservation is equivalent to the requirement that the quadrature is exact, which turns out to be always the case in the event that the Hamiltonian function is a polynomial and the degree of precision of the quadrature formula is high enough. The non-polynomial case is also discussed and a number of test problems are finally presented in order to compare the behavior of the new methods to the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A unifying framework for the derivation and analysis of effective classes of one-step methods for ODEs

In this paper, we provide a simple framework to derive and analyse several classes of effective one-step methods. The framework consists in the discretization of a local Fourier expansion of the continuous problem. Different choices of the basis lead to different classes of methods, even though we shall here consider only the case of an orthonormal polynomial basis, from which a large subclass of Runge-Kutta methods is derived. The obtained results are then applied to prove, in a simplified way, the order and stability properties of Hamiltonian BVMs (HBVMs), a recently introduced class of energy preserving methods for canonical Hamiltonian systems. A few numerical tests with such methods are also included, in order to confirm the effectiveness of the methods.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Blended General Linear Methods based on Boundary Value Methods in the GBDF family

Among the methods for solving ODE-IVPs, the class of General Linear Methods (GLMs) is able to encompass most of them, ranging from Linear Multistep Formulae (LMF) to RK formulae. Moreover, it is possible to obtain methods able to overcome typical drawbacks of the previous classes of methods. For example, order barriers for stable LMF and the problem of order reduction for RK methods. Nevertheless, these goals are usually achieved at the price of a higher computational cost. Consequently, many efforts have been made in order to derive GLMs with particular features, to be exploited for their efficient implementation. In recent years, the derivation of GLMs from particular Boundary Value Methods (BVMs), namely the family of Generalized BDF (GBDF), has been proposed for the numerical solution of stiff ODE-IVPs. In particular, this approach has been recently developed, resulting in a new family of L-stable GLMs of arbitrarily high order, whose theory is here completed and fully worked-out. Moreover, for each one of such methods, it is possible to define a corresponding Blended GLM which is equivalent to it from the point of view of the stability and order properties. These blended methods, in turn, allow the definition of efficient nonlinear splittings for solving the generated discrete problems. A few numerical tests, confirming the excellent potential of such blended methods, are also reported.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Energy and quadratic invariants preserving integrators of Gaussian type

Recently, a whole class of evergy-preserving integrators has been derived for the numerical solution of Hamiltonian problems. In the mainstream of this research, we have defined a new family of symplectic integrators depending on a real parameter. When it is zero, the corresponding method in the family becomes the classical Gauss collocation formula of order 2s, where s denotes the number of the internal stages. For any given non-null value of the parameter, the corresponding method remains symplectic and has order 2s-2: hence it may be interpreted as a symplectic perturbation of the Gauss method. Under suitable assumptions, it can be shown that the parameter a may be properly tuned, at each step of the integration procedure, so as to guarantee energy conservation in the numerical solution. The resulting method shares the same order 2s as the generating Gauss formula, and is able to preserve both energy and quadratic invariants.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hamiltonian Boundary Value Methods (Energy Conserving Discrete Line Integral Methods)

Recently, a new family of integrators (Hamiltonian Boundary ValueMethods) has been introduced, which is able to precisely conserve the energy function of polynomial Hamiltonian systems and to provide a practical conservation of the energy in the non-polynomial case. We settle the definition and the theory of such methods in a more general framework. Our aim is on the one hand to give account of their good behavior when applied to general Hamiltonian systems and, on the other hand, to find out what are the optimal formulae, in relation to the choice of the polynomial basis and of the distribution of the nodes. Such analysis is based upon the notion of extended collocation conditions and the definition of discrete line integral, and is carried out by looking at the limit of such family of methods as the number of the so called silent stages tends to infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Isospectral Property of Hamiltonian Boundary Value Methods (HBVMs) and their blended implementation

One main issue, when numerically integrating autonomous Hamiltonian systems, is the long-term conservation of some of its invariants, among which the Hamiltonian function itself. Recently, a new class of methods, named "Hamiltonian Boundary Value Methods (HBVMs)" has been introduced and analysed, which are able to exactly preserve polynomial Hamiltonians of arbitrarily high degree. We here study a further property of such methods, namely that of having, when cast as Runge-Kutta methods, a matrix of the Butcher tableau with the same spectrum (apart the zero eigenvalues) as that of the corresponding Gauss-Legendre method, independently of the considered abscissae. Consequently, HBVMs are always perfectly A-stable methods. Moreover, this allows their efficient "blended" implementation, for solving the generated discrete problems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Isospectral Property of Hamiltonian Boundary Value Methods (HBVMs) and their connections with Runge-Kutta collocation methods

One main issue, when numerically integrating autonomous Hamiltonian systems, is the long-term conservation of some of its invariants, among which the Hamiltonian function itself. Recently, a new class of methods, named Hamiltonian Boundary Value Methods (HBVMs) has been introduced and analysed, which are able to exactly preserve polynomial Hamiltonians of arbitrarily high degree. We here study a further property of such methods, namely that of having, when cast as a Runge-Kutta method, a matrix of the Butcher tableau with the same spectrum (apart from the zero eigenvalues) as that of the corresponding Gauss-Legendre method, independently of the considered abscissae. Consequently, HBVMs are always perfectly A-stable methods. This, in turn, allows to elucidate the existing connections with classical Runge-Kutta collocation methods.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical comparisons between Gauss-Legendre methods and Hamiltonian BVMs defined over Gauss points

Hamiltonian Boundary Value Methods are a new class of energy preserving one step methods for the solution of polynomial Hamiltonian dynamical systems. They can be thought of as a generalization of collocation methods in that they may be defined by imposing a suitable set of extended collocation conditions. In particular, in the way they are described in this note, they are related to Gauss collocation methods with the difference that they are able to precisely conserve the Hamiltonian function in the case where this is a polynomial of any high degree in the momenta and in the generalized coordinates. A description of these new formulas is followed by a few test problems showing how, in many relevant situations, the precise conservation of the Hamiltonian is crucial to simulate on a computer the correct behavior of the theoretical solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Solution of ODEs and the Columbus' Egg: Three Simple Ideas for Three Difficult Problems

On computers, discrete problems are solved instead of continuous ones. One must be sure that the solutions of the former problems, obtained in real time (i.e., when the stepsize h is not infinitesimal) are good approximations of the solutions of the latter ones. However, since the discrete world is much richer than the continuous one (the latter being a limit case of the former), the classical definitions and techniques, devised to analyze the behaviors of continuous problems, are often insufficient to handle the discrete case, and new specific tools are needed. Often, the insistence in following a path already traced in the continuous setting, has caused waste of time and efforts, whereas new specific tools have solved the problems both more easily and elegantly. In this paper we survey three of the main difficulties encountered in the numerical solutions of ODEs, along with the novel solutions proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the Existence of Energy-Preserving Symplectic Integrators Based upon Gauss Collocation Formulae

We introduce a new family of symplectic integrators depending on a real parameter. When the paramer is zero, the corresponding method in the family becomes the classical Gauss collocation formula of order 2s, where s denotes the number of the internal stages. For any given non-null value of the parameter, the corresponding method remains symplectic and has order 2s-2: hence it may be interpreted as an order 2s-2 (symplectic) perturbation of the Gauss method. Under suitable assumptions, we show that the free parameter may be properly tuned, at each step of the integration procedure, so as to guarantee energy conservation in the numerical solution. The resulting symplectic, energy conserving method shares the same order 2s as the generating Gauss formula.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Hamiltonian BVMs (HBVMs) Homepage

Hamiltonian Boundary Value Methods (in short, HBVMs) is a new class of numerical methods for the efficient numerical solution of canonical Hamiltonian systems. In particular, their main feature is that of exactly preserving, for the numerical solution, the value of the Hamiltonian function, when the latter is a polynomial of arbitrarily high degree. Clearly, this fact implies a practical conservation of any analytical Hamiltonian function. In this notes, we collect the introductory material on HBVMs contained in the HBVMs Homepage, available at http://web.math.unifi.it/users/brugnano/HBVM/index.html

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Lack of Continuity and the Role of Infinite and Infinitesimal in Numerical Methods for ODEs: the Case of Symplecticity

When numerically integrating canonical Hamiltonian systems, the long-term conservation of some of its invariants, among which the Hamiltonian function itself, assumes a central role. The classical approach to this problem has led to the definition of symplectic methods, among which we mention Gauss-Legendre collocation formulae. Indeed, in the continuous setting, energy conservation is derived from symplecticity via an infinite number of infinitesimal contact transformations. However, this infinite process cannot be directly transferred to the discrete setting. By following a different approach, in this paper we describe a sequence of methods, sharing the same essential spectrum (and, then, the same essential properties), which are energy preserving starting from a certain element of the sequence on, i.e., after a finite number of steps.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Analisys of Hamiltonian Boundary Value Methods (HBVMs): a class of energy-preserving Runge-Kutta methods for the numerical solution of polynomial Hamiltonian systems

One main issue, when numerically integrating autonomous Hamiltonian systems, is the long-term conservation of some of its invariants, among which the Hamiltonian function itself. For example, it is well known that classical symplectic methods can only exactly preserve, at most, quadratic Hamiltonians. In this paper, a new family of methods, called "Hamiltonian Boundary Value Methods (HBVMs)", is introduced and analyzed. HBVMs are able to exactly preserve, in the discrete solution, Hamiltonian functions of polynomial type of arbitrarily high degree. These methods turn out to be symmetric, precisely A-stable, and can have arbitrarily high order. A few numerical tests confirm the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fifty Years of Stiffness

The notion of stiffness, which originated in several applications of a different nature, has dominated the activities related to the numerical treatment of differential problems for the last fifty years. Contrary to what usually happens in Mathematics, its definition has been, for a long time, not formally precise (actually, there are too many of them). Again, the needs of applications, especially those arising in the construction of robust and general purpose codes, require nowadays a formally precise definition. In this paper, we review the evolution of such a notion and we also provide a precise definition which encompasses all the previous ones.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Iterative solution of piecewise linear systems for the numerical solution of obstacle problems

We investigate the use of piecewise linear systems, whose coefficient matrix is a piecewise constant function of the solution itself. Such systems arise, for example, from the numerical solution of linear complementarity problems and in the numerical solution of free-surface problems. In particular, we here study their application to the numerical solution of both the (linear) parabolic obstacle problem and the obstacle problem. We propose a class of effective semi-iterative Newton-type methods to find the exact solution of such piecewise linear systems. We prove that the semiiterative Newton-type methods have a global monotonic convergence property, i.e., the iterates converge monotonically to the exact solution in a finite number of steps. Numerical examples are presented to demonstrate the effectiveness of the proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Parallel Factorizations in Numerical Analysis

In this paper we review the parallel solution of sparse linear systems, usually deriving by the discretization of ODE-IVPs or ODE-BVPs. The approach is based on the concept of parallel factorization of a (block) tridiagonal matrix. This allows to obtain efficient parallel extensions of many known matrix factorizations, and to derive, as a by-product, a unifying approach to the parallel solution of ODEs.

0 citations0 reviewsNo code linkedOpen access