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

Federico Poloni

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

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

Catalog footprint

What is connected

15works
7topics
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

15 published item(s)

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.

0 citations0 reviewsNo code linkedOpen 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 $α\in [0,1)$. For $α\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 $\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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nearest $Ω$-stable matrix via Riemannian optimization

We study the problem of finding the nearest $Ω$-stable matrix to a certain matrix $A$, 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.

0 citations0 reviewsNo code linkedOpen 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 $Q_k = X_{2^k}$ of another naturally-arising fixed-point iteration $(X_h)$ via a sort of repeated squaring. The equations we consider are Stein equations $X - A^*XA=Q$, Lyapunov equations $A^*X+XA+Q=0$, discrete-time algebraic Riccati equations $X=Q+A^*X(I+GX)^{-1}A$, continuous-time algebraic Riccati equations $Q+A^*X+XA-XGX=0$, palindromic quadratic matrix equations $A+QY+A^*Y^2=0$, and nonlinear matrix equations $X+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.

0 citations0 reviewsNo code linkedOpen 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) = Vz^2-Dz +Q$, where the matrices $V$ and $D$ are diagonal and $Q$ is the transition matrix of a discrete-time Markov chain. Its solution is typically constructed by computing an invariant pair of $P(z)$ associated with its eigenvalues in the left half-plane, or by solving the matrix equation $X^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=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.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen 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)$, 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.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The SDA Method for Numerical Solution of Lur'e Equations

We introduce a numerical method for the numerical solution of the so-called Lur'e matrix equations that arise in balancing-related model reduction and linear-quadratic infinite time horizon optimal control. Based on the fact that the set of solutions can be characterized in terms of deflating subspaces of even matrix pencils, an iterative scheme is derived that converges linearly to the maximal solution.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen access