Related papers: Constructing Sobolev orthonormal rational function…
In this paper we develop randomized Krylov subspace methods for efficiently computing regularized solutions to large-scale linear inverse problems. Building on the recently developed randomized Gram-Schmidt process, where sketched inner…
Matrices whose adjoint is a low rank perturbation of a rational function of the matrix naturally arise when trying to extend the well known Faber-Manteuffel theorem, which provides necessary and sufficient conditions for the existence of a…
This work develops novel rational Krylov methods for updating a large-scale matrix function f(A) when A is subject to low-rank modifications. It extends our previous work in this context on polynomial Krylov methods, for which we present a…
In this paper, we study the sequence of orthogonal polynomials $\{S_n\}_{n=0}^{\infty}$ with respect to the Sobolev-type inner product $$\langle f,g \rangle= \int_{-1}^{1} f(x) g(x) \,d\mu(x) +\sum_{j=1}^{N} \eta_{j} \,f^{(d_j)}(c_{j})…
Our main objective in this work is to show how Sobolev orthogonal polynomials emerge as a useful tool within the framework of spectral methods for boundary-value problems. The solution of a boundary-value problem for a stationary…
The randomized Arnoldi process has been used in large-scale scientific computing because it produces a well-conditioned basis for the Krylov subspace more quickly than the standard Arnoldi process. However, the resulting Hessenberg matrix…
This contribution aims to obtain several connection formulae for the polynomial sequence, which is orthogonal with respect to the discrete Sobolev inner product \[ \langle f, g\rangle_n=\langle {\bf u}, fg\rangle+ \sum_{j=1}^M \mu_{j}…
We study the connection between block Krylov subspaces and matrix orthogonal functions. Under a no-deflation assumption, we show that polynomial block Krylov subspaces are isometrically isomorphic to spaces of matrix polynomials of bounded…
We introduce an extended Sobolev scale on a smooth compact manifold with boundary. The scale is formed by inner-product H\"ormander spaces for which an RO-varying radial function serves as a regularity index. These spaces do not depend on a…
We glue two families of Bernstein-Szego polynomials to construct the eigenbasis of an associated finite-dimensional Jacobi matrix. This gives rise to finite orthogonality relations for this composite eigenbasis of Bernstein-Szego…
In this paper we deal with a family of non--standard polynomials orthogonal with respect to an inner product involving differences. This type of inner product is the so--called $\Delta$--Sobolev inner product. Concretely, we consider the…
A fast and weakly stable method for computing the zeros of a particular class of hypergeometric polynomials is presented. The studied hypergeometric polynomials satisfy a higher order differential equation and generalize Laguerre…
Flexible Krylov methods are a common standpoint for inverse problems. In particular, they are used to address the challenges associated with explicit variational regularization when it goes beyond the two-norm, for example involving an…
Orthogonal polynomials with respect to the weight function $w_{\beta,\gamma}(t) = t^\beta (1-t)^\gamma$, $\gamma > -1$, on the conic surface $\{(x,t): \|x\| = t, \, x \in \mathbb{R}^d, \, t \le 1\}$ are studied recently, and are shown to be…
We introduce and develop a theory of orthogonality with respect to Sobolev inner products on the real line for sequences of functions with a tridiagonal, skew-Hermitian differentiation matrix. While a theory of such L2-orthogonal systems is…
By closely following a construction by Ganelius, we construct Faber rational functions that allow us to derive tight and explicit bounds on Zolotarev numbers. We use our results to bound the singular values of matrices, including…
In this paper we give a geometric condition which ensures that $(q,p)$-Poincar\'e-Sobolev inequalities are implied from generalized $(1,1)$-Poincar\'e inequalities related to $L^1$ norms in the context of product spaces. The concept of…
The parallel orbital-updating approach is an orbital/eigenfunction iteration based approach for solving eigenvalue problems when many eigenpairs are required. It has been proven to be efficient, for instance, in electronic structure…
When integrating functions that have poles outside the interval of integration, but are regular otherwise, it is suggested that the quadrature rule in question ought to integrate exactly not only polynomials (if any), but also suitable…
A new representation is proposed for functions in a Sobolev space with dominating mixed smoothness on an $N$-dimensional hyperrectangle. In particular, it is shown that these functions can be expressed in terms of their highest-order mixed…