Related papers: Multigrid methods for Toeplitz linear systems with…

We consider the Gabor system $\mathcal{G}(g,a\mathbb{Z}\times b\mathbb{Z})$ generated by a continuous, compactly supported function $g$ over the time-frequency lattice generated by the parameters $a$ and $b$. We show that, under an…

In this work, we propose a novel method for quantifying distances between Toeplitz structured covariance matrices. By exploiting the spectral representation of Toeplitz matrices, the proposed distance measure is defined based on an optimal…

In this manuscript, we present a collective multigrid algorithm to solve efficiently the large saddle-point systems of equations that typically arise in PDE-constrained optimization under uncertainty, and develop a novel convergence…

We address the estimation of the scatter matrix of a scale mixture of Gaussian stationary autoregressive vectors. This is equivalent to consider the estimation of a structured scatter matrix of a Spherically Invariant Random Vector (SIRV)…

We compute the minimal and the maximal bound on the number of generators of a minimal smooth monomial Togliatti system of forms of degree $d$ in $n+1$ variables, for any $d\ge 2$ and $n\geq 2$. We classify the Togliatti systems with number…

In the past decades, a remarkable amount of research has been carried out regarding fast solvers for large linear systems resulting from various discretizations of fractional differential equations (FDEs). In the current work, we focus on…

Motivated by applications in single-cell biology and metagenomics, we investigate the problem of matrix reordering based on a noisy disordered monotone Toeplitz matrix model. We establish the fundamental statistical limit for this problem…

The nearest circulant approximation of a real Toeplitz matrix in the Frobenius norm is derived. This matrix is symmetric. It is proven that symmetric circulant matrices are the only real circulant matrices with all real eigenvalues. The…

We introduce a framework and early results for massively scalable Gaussian processes (MSGP), significantly extending the KISS-GP approach of Wilson and Nickisch (2015). The MSGP framework enables the use of Gaussian processes (GPs) on…

This paper proposes two linear projection methods for supervised dimension reduction using only the first and second-order statistics. The methods, each catering to a different parameter regime, are derived under the general Gaussian model…

Machine learning-based unfolding has enabled unbinned and high-dimensional differential cross section measurements. Two main approaches have emerged in this research area: one based on discriminative models and one based on generative…

A new proof is presented of a theorem of L.~Gurvits, which states that the cone of positive block-Toeplitz matrices with matrix entries has no entangled elements. The proof of the Gurvits separation theorem is achieved by making use of the…

The theory of generalized locally Toeplitz (GLT) sequences is an apparatus for computing the spectral and singular value distribution of sequences of matrices that possess a (possibly hidden) Toeplitz-like structure. Sequences of this kind,…

With the emergence of mixed precision capabilities in hardware, iterative refinement schemes for solving linear systems $Ax=b$ have recently been revisited and reanalyzed in the context of three or more precisions. These new analyses show…

In this paper, we develop multigrid solvers for the biharmonic problem in the framework of isogeometric analysis (IgA). In this framework, one typically sets up B-splines on the unit square or cube and transforms them to the domain of…

In this paper two types of multgrid methods, i.e., the Rayleigh quotient iteration and the inverse iteration with fixed shift, are developed for solving the Maxwell eigenvalue problem with discontinuous relative magnetic permeability and…

We introduce two iterative methods, GPBiLQ and GPQMR, for solving unsymmetric partitioned linear systems. The basic mechanism underlying GPBiLQ and GPQMR is a novel simultaneous tridiagonalization via biorthogonality that allows for…

A new integration technique for multi-loop Feynman integrals, called the matrix method, is developed and then applied to the divergent part of the overlapping two-loop quark self-energy function $\,i\Sigma\,$ in the light- cone gauge. It is…

We consider the linear system Ax=b arising from one-dimensional Poisson's equation with Dirichlet boundary conditions, where A is the square matrix with the stencil form [-1 2 -1]. Here we show that a pairwise aggregation-based algebraic…

It is known that the generating function $f$ of a sequence of Toeplitz matrices $\{T_n(f)\}_n$ may not describe the asymptotic distribution of the eigenvalues of $T_n(f)$ if $f$ is not real. In a recent paper, we assume as a working…