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

Dario A. Bini

Dario A. Bini 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.OAmath.RA

Catalog footprint

What is connected

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

preprint2021arXiv

Geometric means of quasi-Toeplitz matrices

We study means of geometric type of quasi-Toeplitz matrices, that are semi-infinite matrices $A=(a_{i,j})_{i,j=1,2,\ldots}$ of the form $A=T(a)+E$, where $E$ represents a compact operator, and $T(a)$ is a semi-infinite Toeplitz matrix associated with the function $a$, with Fourier series $\sum_{\ell=-\infty}^{\infty} a_\ell e^{\mathfrak i \ell t}$, in the sense that $(T(a))_{i,j}=a_{j-i}$. If $a$ is \rv\ and essentially bounded, then these matrices represent bounded self-adjoint operators on $\ell^2$. We consider the case where $a$ is a continuous function, where quasi-Toeplitz matrices coincide with a classical Toeplitz algebra, and the case where $a$ is in the Wiener algebra, that is, has absolutely convergent Fourier series. We prove that if $a_1,\ldots,a_p$ are continuous and positive functions, or are in the Wiener algebra with some further conditions, then means of geometric type, such as the ALM, the NBMP and the Karcher mean of quasi-Toeplitz positive definite matrices associated with $a_1,\ldots,a_p$, are quasi-Toeplitz matrices associated with the geometric mean $(a_1\cdots a_p)^{1/p}$, which differ only by the compact correction. We show by numerical tests that these operator means can be practically approximated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A computational framework for two-dimensional random walks with restarts

The treatment of two-dimensional random walks in the quarter plane leads to Markov processes which involve semi-infinite matrices having Toeplitz or block Toeplitz structure plus a low-rank correction. Finding the steady state probability distribution of the process requires to perform operations involving these structured matrices. We propose an extension of the framework of [5] which allows to deal with more general situations such as processes involving restart events. This is motivated by the need for modeling processes that can incur in unexpected failures like computer system reboots. Algebraically, this gives rise to corrections with infinite support that cannot be treated using the tools currently available in the literature. We present a theoretical analysis of an enriched Banach algebra that, combined with appropriate algorithms, enables the numerical treatment of these problems. The results are applied to the solution of bidimensional Quasi-Birth-Death processes with infinitely many phases which model random walks in the quarter plane, relying on the matrix analytic approach. This methodology reduces the problem to solving a quadratic matrix equation with coefficients of infinite size. We provide conditions on the transition probabilities which ensure that the solution of interest of the matrix equation belongs to the enriched algebra. The reliability of our approach is confirmed by extensive numerical experimentation on some case studies.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Solving quadratic matrix equations arising in random walks in the quarter plane

Quadratic matrix equations of the kind $A_1X^2+A_0X+A_{-1}=X$ are encountered in the analysis of Quasi--Birth-Death stochastic processes where the solution of interest is the minimal nonnegative solution $G$. In many queueing models, described by random walks in the quarter plane, the coefficients $A_1,A_0,A_{-1}$ are infinite tridiagonal matrices with an almost Toeplitz structure. Here, we analyze some fixed point iterations, including Newton's iteration, for the computation of $G$ and introduce effective algorithms and acceleration strategies which fully exploit the Toeplitz structure of the matrix coefficients and of the current approximation. Moreover, we provide a structured perturbation analysis for the solution $G$. The results of some numerical experiments which demonstrate the effectiveness of our approach are reported.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficient cyclic reduction for QBDs with rank structured blocks

We provide effective algorithms for solving block tridiagonal block Toeplitz systems with $m\times m$ quasiseparable blocks, as well as quadratic matrix equations with $m\times m$ quasiseparable coefficients, based on cyclic reduction and on the technology of rank-structured matrices. The algorithms rely on the exponential decay of the singular values of the off-diagonal submatrices generated by cyclic reduction. We provide a formal proof of this decay in the Markovian framework. The results of the numerical experiments that we report confirm a significant speed up over the general algorithms, already starting with the moderately small size $m\approx 10^2$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

General solution of the Poisson equation for Quasi-Birth-and-Death processes

We consider the Poisson equation $(I-P)\boldsymbol{u}=\boldsymbol{g}$, where $P$ is the transition matrix of a Quasi-Birth-and-Death (QBD) process with infinitely many levels, $\bm g$ is a given infinite dimensional vector and $\bm u$ is the unknown. Our main result is to provide the general solution of this equation. To this purpose we use the block tridiagonal and block Toeplitz structure of the matrix $P$ to obtain a set of matrix difference equations, which are solved by constructing suitable resolvent triples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On Functions of quasi Toeplitz matrices

Let $a(z)=\sum_{i\in\mathbb Z}a_iz^i$ be a complex valued continuous function, defined for $|z|=1$, such that $\sum_{i=-\infty}^{+\infty}|ia_i|<\infty$. Consider the semi-infinite Toeplitz matrix $T(a)=(t_{i,j})_{i,j\in\mathbb Z^+}$ associated with the symbol $a(z)$ such that $t_{i,j}=a_{j-i}$. A quasi-Toeplitz matrix associated with the continuous symbol $a(z)$ is a matrix of the form $A=T(a)+E$ where $E=(e_{i,j})$, $\sum_{i,j\in\mathbb Z^+}|e_{i,j}|<\infty$, and is called a CQT-matrix. Given a function $f(x)$ and a CQT matrix $M$, we provide conditions under which $f(M)$ is well defined and is a CQT matrix. Moreover, we introduce a parametrization of CQT matrices and algorithms for the computation of $f(M)$. We treat the case where $f(x)$ is assigned in terms of power series and the case where $f(x)$ is defined in terms of a Cauchy integral. This analysis is applied also to finite matrices which can be written as the sum of a Toeplitz matrix and of a low rank correction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the exponential of semi-infinite quasi-Toeplitz matrices

Let $a(z)=\sum_{i\in\mathbb Z}a_iz^i$ be a complex valued function defined for $|z|=1$, such that $\sum_{i\in\mathbb Z}|ia_i|<\infty$, and let $E=(e_{i,j})_{i,j\in\mathbb {Z}^+}$ be such that $\sum_{i,j\in\mathbb{Z}^+}|e_{i,j}|<\infty$. A semi-infinite quasi-Toeplitz matrix is a matrix of the kind $A=T(a)+E$, where $T(a)=(t_{i,j})_{i,j\in\mathbb{Z}^+}$ is the semi-infinite Toeplitz matrix associated with the symbol $a(z)$, that is, $t_{i,j}=a_{j-i}$ for $i,j\in\mathbb Z^+$. We analyze theoretical and computational properties of the exponential of $A$. More specifically, it is shown that $\exp(A)=T(\exp(a))+F$ where $F=(f_{i,j})_{i,j\in\mathbb{Z}^+}$ is such that $\sum_{i,j\in\mathbb{Z}^+}|f_{i,j}|$ is finite, i.e., $\exp(A)$ is a semi-infinite quasi-Toeplitz matrix as well, and an effective algorithm for its computation is given. These results can be extended from the function $\exp(z)$ to any function $f(z)$ satisfying mild conditions, and can be applied to finite quasi-Toeplitz matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Semi-Infinite Quasi-Toeplitz Matrices with Applications to QBD Stochastic Processes

Denote by $\mathcal{W}_1$ the set of complex valued functions of the form $a(z)=\sum_{i=-\infty}^{+\infty}a_iz^i$ which are continuous on the unit circle, and such that $\sum_{i=-\infty}^{+\infty}|ia_i|<\infty$. We call CQT matrix a quasi-Toeplitz matrix $A$, associated with a continuous symbol $a(z)\in\mathcal W_1$, of the form $A=T(a)+E$, where $T(a)=(t_{i,j})_{i,j\in\mathbb{Z}^+}$ is the semi-infinite Toeplitz matrix such that $t_{i,j}=a_{j-i}$, for $i,j\in\mathbb Z^+$, and $E=(e_{i,j})_{i,j\in\mathbb{Z}^+}$ is a semi-infinite matrix such that $\sum_{i,j=1}^{+\infty}|e_{i,j}|$ is finite. We prove that the class of CQT matrices is a Banach algebra with a suitable sub-multiplicative matrix norm $\|\cdot\|$. We introduce a finite representation of CQT matrices together with algorithms which implement elementary matrix operations. An application to solving quadratic matrix equations of the kind $AX^2+BX+C=0$, encountered in the solution of Quasi-Birth and Death (QBD) stochastic processes with a denumerable set of phases, is presented where $A,B,C$ are CQT matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shift techniques for Quasi-Birth and Death processes: canonical factorizations and matrix equations

We revisit the shift technique applied to Quasi-Birth and Death (QBD) processes (He, Meini, Rhee, SIAM J. Matrix Anal. Appl., 2001) by bringing the attention to the existence and properties of canonical factorizations. To this regard, we prove new results concerning the solutions of the quadratic matrix equations associated with the QBD. These results find applications to the solution of the Poisson equation for QBDs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generalization of the Brauer Theorem to Matrix Polynomials and Matrix Laurent Series

Given a square matrix $A$, Brauer's theorem [Duke Math. J. 19 (1952), 75--91] shows how to modify one single eigenvalue of $A$ via a rank-one perturbation, without changing any of the remaining eigenvalues. We reformulate Brauer's theorem in functional form and provide extensions to matrix polynomials and to matrix Laurent series $A(z)$ together with generalizations to shifting a set of eigenvalues. We provide conditions under which the modified function $\widetilde A(z)$ has a canonical factorization $\widetilde A(z)=\widetilde U(z)\widetilde L(z^{-1})$ and we provide explicit expressions of the factors $\widetilde U(z)$ and $\widetilde L(z)$. Similar conditions and expressions are given for the factorization of $\widetilde A(z^{-1})$. Some applications are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On a Class of Matrix Pencils and $\ell$-ifications Equivalent to a Given Matrix Polynomial

A new class of linearizations and $\ell$-ifications for $m\times m$ matrix polynomials $P(x)$ of degree $n$ is proposed. The $\ell$-ifications in this class have the form $A(x) = D(x) + (e\otimes I_m) W(x)$ where $D$ is a block diagonal matrix polynomial with blocks $B_i(x)$ of size $m$, $W$ is an $m\times qm$ matrix polynomial and $e=(1,\ldots,1)^t\in\mathbb C^q$, for a suitable integer $q$. The blocks $B_i(x)$ can be chosen a priori, subjected to some restrictions. Under additional assumptions on the blocks $B_i(x)$ the matrix polynomial $A(x)$ is a strong $\ell$-ification, i.e., the reversed polynomial of $A(x)$ defined by $A^\#(x) := x^{\mathrm{deg} A(x)} A(x^{-1})$ is an $\ell$-ification of $P^\#(x)$. The eigenvectors of the matrix polynomials $P(x)$ and $A(x)$ are related by means of explicit formulas. Some practical examples of $\ell$-ifications are provided. A strategy for choosing $B_i(x)$ in such a way that $A(x)$ is a well conditioned linearization of $P(x)$ is proposed. Some numerical experiments that validate the theoretical results are reported

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quasiseparable Hessenberg reduction of real diagonal plus low rank matrices and applications

We present a novel algorithm to perform the Hessenberg reduction of an $n\times n$ matrix $A$ of the form $A = D + UV^*$ where $D$ is diagonal with real entries and $U$ and $V$ are $n\times k$ matrices with $k\le n$. The algorithm has a cost of $O(n^2k)$ arithmetic operations and is based on the quasiseparable matrix technology. Applications are shown to solving polynomial eigenvalue problems and some numerical experiments are reported in order to analyze the stability of the approach

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Locating the eigenvalues of matrix polynomials

Some known results for locating the roots of polynomials are extended to the case of matrix polynomials. In particular, a theorem by A.E. Pellet [Bulletin des Sciences Mathématiques, (2), vol 5 (1881), pp.393-395], some results of D.A. Bini [Numer. Algorithms 13:179-200, 1996] based on the Newton polygon technique, and recent results of M. Akian, S. Gaubert and M. Sharify (see in particular [LNCIS, 389, Springer p.p.291-303] and [M. Sharify, Ph.D. thesis, École Polytechnique, ParisTech, 2011]). These extensions are applied for determining effective initial approximations for the numerical computation of the eigenvalues of matrix polynomials by means of simultaneous iterations, like the Ehrlich-Aberth method. Numerical experiments that show the computational advantage of these results are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solving polynomial eigenvalue problems by means of the Ehrlich-Aberth method

Given the $n\times n$ matrix polynomial $P(x)=\sum_{i=0}^kP_i x^i$, we consider the associated polynomial eigenvalue problem. This problem, viewed in terms of computing the roots of the scalar polynomial $\det P(x)$, is treated in polynomial form rather than in matrix form by means of the Ehrlich-Aberth iteration. The main computational issues are discussed, namely, the choice of the starting approximations needed to start the Ehrlich-Aberth iteration, the computation of the Newton correction, the halting criterion, and the treatment of eigenvalues at infinity. We arrive at an effective implementation which provides more accurate approximations to the eigenvalues with respect to the methods based on the QZ algorithm. The case of polynomials having special structures, like palindromic, Hamiltonian, symplectic, etc., where the eigenvalues have special symmetries in the complex plane, is considered. A general way to adapt the Ehrlich-Aberth iteration to structured matrix polynomial is introduced. Numerical experiments which confirm the effectiveness of this approach are reported.

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