BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Donato Trigiante

Donato Trigiante contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
1topics
3close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

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

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 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

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

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