BZPEER
More

Discover

About BZPEERSearchGraph

Research tools

HomeCollectionsResearchCollaborationPublish

Network

ExpertsCommunities

Opportunities

Listings

Account

Log inCreate account
Log inCreate account

Source author record

Federico Poloni

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record
math.NANumerical Analysismath.OCmath.PRMethodologyphysics.soc-phPopulations and Evolution

Catalog footprint

Research overview

15works
7topics
4close collaborators

Actions

Is this your work?

Log in to claim
Open graphBrowse papers

Published work

15 publications

preprint2022arXiv

On the use of Markov chains for epidemic modeling on networks

We discuss various models for epidemics on networks that rely on Markov chains. Random walks on graphs are often used to predict epidemic spread and to investigate possible control actions to mitigate them. In this study, we demonstrate that they do not fully reflect the dynamics of epidemics, as they overestimate infection times. Accordingly, we explain how Markov chains may still be used to accurately model the virus spread, and to correctly predict infection times. We also provide an algorithm that efficiently estimates infection times via a sampling strategy. Finally, we present a novel indicator based on infection times, and we compare its node ranking properties with other centrality measures based on random walks.

Open access
preprint2021arXiv

A continuation method for computing the multilinear Pagerank

The multilinear Pagerank model [Gleich, Lim and Yu, 2015] is a tensor-based generalization of the Pagerank model. Its computation requires solving a system of polynomial equations that contains a parameter α∈[0,1)α\in [0,1). For α≈1α\approx 1, this computation remains a challenging problem, especially since the solution may be non-unique. Extrapolation strategies that start from smaller values of αα and `follow' the solution by slowly increasing this parameter have been suggested; however, there are known cases where these strategies fail, because a globally continuous solution curve cannot be defined as a function of αα. In this paper, we improve on this idea, by employing a predictor-corrector continuation algorithm based on a more general representation of the solutions as a curve in Rn+1\mathbb{R}^{n+1}. We prove several global properties of this curve that ensure the good behavior of the algorithm, and we show in our numerical experiments that this method is significantly more reliable than the existing alternatives.

Open access
preprint2021arXiv

Nearest ΩΩ-stable matrix via Riemannian optimization

We study the problem of finding the nearest ΩΩ-stable matrix to a certain matrix AA, i.e., the nearest matrix with all its eigenvalues in a prescribed closed set ΩΩ. Distances are measured in the Frobenius norm. An important special case is finding the nearest Hurwitz or Schur stable matrix, which has applications in systems theory. We describe a reformulation of the task as an optimization problem on the Riemannian manifold of orthogonal (or unitary) matrices. The problem can then be solved using standard methods from the theory of Riemannian optimization. The resulting algorithm is remarkably fast on small-scale and medium-scale matrices, and returns directly a Schur factorization of the minimizer, sidestepping the numerical difficulties associated with eigenvalues with high multiplicity.

Open access
preprint2020arXiv

Iterative and doubling algorithms for Riccati-type matrix equations: a comparative introduction

We review a family of algorithms for Lyapunov- and Riccati-type equations which are all related to each other by the idea of \emph{doubling}: they construct the iterate Qk=X2kQ_k = X_{2^k} of another naturally-arising fixed-point iteration (Xh)(X_h) via a sort of repeated squaring. The equations we consider are Stein equations X−A∗XA=QX - A^*XA=Q, Lyapunov equations A∗X+XA+Q=0A^*X+XA+Q=0, discrete-time algebraic Riccati equations X=Q+A∗X(I+GX)−1AX=Q+A^*X(I+GX)^{-1}A, continuous-time algebraic Riccati equations Q+A∗X+XA−XGX=0Q+A^*X+XA-XGX=0, palindromic quadratic matrix equations A+QY+A∗Y2=0A+QY+A^*Y^2=0, and nonlinear matrix equations X+A∗X−1A=QX+A^*X^{-1}A=Q. We draw comparisons among these algorithms, highlight the connections between them and to other algorithms such as subspace iteration, and discuss open issues in their theory.

Open access
preprint2016arXiv

Componentwise accurate Brownian motion computations using Cyclic Reduction

Markov-modulated Brownian motion is a popular tool to model continuous-time phenomena in a stochastic context. The main quantity of interest is the invariant density, which satisfies a differential equation associated with the quadratic matrix polynomial P(z)=Vz2−Dz+QP(z) = Vz^2-Dz +Q, where the matrices VV and DD are diagonal and QQ is the transition matrix of a discrete-time Markov chain. Its solution is typically constructed by computing an invariant pair of P(z)P(z) associated with its eigenvalues in the left half-plane, or by solving the matrix equation X2V−XD+Q=0X^2V-XD+Q=0. We show that these tasks can be solved using a componentwise accurate algorithm based on Cyclic Reduction, generalizing the recently appeared algorithms for the linear case (V=0V=0). We give a proof of the numerical stability of our algorithm in the componentwise sense; the same proof applies to Cyclic Reduction in a more general M-matrix setting which appears in other applications such as the modelling of QBD processes.

Open access
preprint2015arXiv

Multivariate trend-cycle extraction with the Hodrick-Prescott filter

The Hodrick-Prescott filter represents one of the most popular method for trend-cycle extraction in macroeconomic time series. In this paper we provide a multivariate generalization of the Hodrick-Prescott filter, based on the seemingly unrelated time series approach. We first derive closed-form expressions linking the signal-noise matrix ratio to the parameters of the VARMA representation of the model. We then show that the parameters can be estimated using a recently introduced method, called "Moment Estimation Through Aggregation (META)". This method replaces the traditional multivariate likelihood estimation with a procedure that requires estimating univariate processes only. This makes the estimation simpler, faster and better-behaved numerically. We prove that our estimation method is consistent and asymptotically normal distributed for the proposed framework. Finally, we present an empirical application focusing on the industrial production of several European countries.

Open access
preprint2011arXiv

A duality relation for matrix pencils with application to linearizations

The aim of this paper is twofold. First, we introduce a new class of linearizations, based on the generalization of a construction used in polynomial algebra to find the zeros of a system of (scalar) polynomial equations. We show that one specific linearization in this class, which is constructed naturally from the QR factorization of the matrix obtained by stacking the coefficients of A(x)A(x), has good conditioning and stability properties. Moreover, while analyzing this class, we introduce a general technique to derive new linearizations from existing ones. This technique generalizes some ad-hoc arguments used in dealing with the existing linearization classes, and can hopefully be used to derive a simpler and more general theory of linearizations. This technique relates linearizations to \emph{pencil arithmetic}, a technique used in solving matrix equations that allows to extend some algebraic operations from matrix to matrix pencils.

Open access
preprint2011arXiv

A Subspace Shift Technique for Nonsymmetric Algebraic Riccati Equations

The worst situation in computing the minimal nonnegative solution of a nonsymmetric algebraic Riccati equation associated with an M-matrix occurs when the corresponding linearizing matrix has two very small eigenvalues, one with positive and one with negative real part. When both these eigenvalues are exactly zero, the problem is called critical or null recurrent. While in this case the problem is ill-conditioned and the convergence of the algorithms based on matrix iterations is slow, there exist some techniques to remove the singularity and transform the problem to a well-behaved one. Ill-conditioning and slow convergence appear also in close-to-critical problems, but when none of the eigenvalues is exactly zero the techniques used for the critical case cannot be applied. In this paper, we introduce a new method to accelerate the convergence properties of the iterations also in close-to-critical cases, by working on the invariant subspace associated with the problematic eigenvalues as a whole. We present a theoretical analysis and several numerical experiments which confirm the efficiency of the new method.

Open access
preprint2011arXiv

Quadratic Vector Equations

We study in an unified fashion several quadratic vector and matrix equations with nonnegativity hypotheses. Specific cases of such problems (QBD equations, nonsymmetric algebraic Riccati equations, Lu's simple equation, Markovian binary trees equations) have been studied extensively in the past by several authors. Many of the results appearing here have already been proved for one or more of the single instances of the problems, resorting to specific characteristics of the problem. In some cases the proofs we present here are mere rewriting of the original proofs with a little change of notation to adapt them to our framework, but in some cases we are effectively able to remove some hypotheses and generalize the results by abstracting the specific aspects of each problem.

Open access
preprint2010arXiv

A Perron iteration for the solution of a quadratic vector equation arising in Markovian Binary Trees

We propose a novel numerical method for solving a quadratic vector equation arising in Markovian Binary Trees. The numerical method consists in a fixed point iteration, expressed by means of the Perron vectors of a sequence of nonnegative matrices. A theoretical convergence analysis is performed. The proposed method outperforms the existing methods for close-to-critical problems.

Open access
preprint2010arXiv

On the solution of a quadratic vector equation arising in Markovian Binary Trees

We present some advances, both from a theoretical and from a computational point of view, on a quadratic vector equation (QVE) arising in Markovian Binary Trees. Concerning the theoretical advances, some irreducibility assumptions are relaxed, and the minimality of the solution of the QVE is expressed in terms of properties of the Jacobian of a suitable function. From the computational point of view, we elaborate on the Perron vector-based iteration proposed in [http://arxiv.org/abs/1006.0577]. In particular we provide a condition which ensures that the Perron iteration converges to the sought solution of the QVE. Moreover we introduce a variant of the algorithm which consists in applying the Newton method instead of a fixed-point iteration. This method has the same convergence behaviour as the Perron iteration, since its convergence speed increases for close-to-critical problems. Moreover, unlike the Perron iteration, the method has a quadratic convergence. Finally, we show that it is possible to alter the bilinear form defining the QVE in several ways without changing the solution. This modification has an impact on convergence speed of the algorithms.

Open access
preprint2009arXiv

A note on the O(n)-storage implementation of the GKO algorithm

We propose a new O(n)-space implementation of the GKO-Cauchy algorithm for the solution of linear systems with Cauchy-like matrix. Despite its slightly higher computational cost, this new algorithm makes a more efficient use of the processor cache memory. Thus, for matrices of size larger than about 500-1000, it outperforms the existing algorithms. We present an applicative case of Cauchy-like matrices with non-reconstructible main diagonal. In this special instance, the O(n) space algorithms can be adapted nicely to provide an efficient implementation of basic linear algebra operations in terms of the low displacement-rank generators.

Open access

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.