We analyze different measures for the backward error of a set of numerical approximations for the roots of a polynomial. We focus mainly on the element-wise mixed backward error introduced by Mastronardi and Van Dooren, and the tropical backward error introduced by Tisseur and Van Barel. We show that these measures are equivalent under suitable assumptions. We also show relations between these measures and the classical element-wise and norm-wise backward error measures.
In two dimensions, we propose and analyze an a posteriori error estimator for finite element approximations of the stationary Navier Stokes equations with singular sources on Lipschitz, but not necessarily convex, polygonal domains. Under a smallness assumption on the continuous and discrete solutions, we prove that the devised error estimator is reliable and locally efficient. We illustrate the theory with numerical examples.
Various methods can obtain certified estimates for roots of polynomials. Many applications in science and engineering additionally utilize the value of functions evaluated at roots. For example, critical values are obtained by evaluating an objective function at critical points. For analytic evaluation functions, Newtons method naturally applies to yield certified estimates. These estimates no longer apply, however, for Holder continuous functions, which are a generalization of Lipschitz continuous functions where continuous derivatives need not exist. This work develops and analyzes an alternative approach for certified estimates of evaluating locally Holder continuous functions at roots of polynomials. An implementation of the method in Maple demonstrates efficacy and efficiency.
Let $G$ be a finite union of disjoint and bounded Jordan domains in the complex plane, let $mathcal{K}$ be a compact subset of $G$ and consider the set $G^star$ obtained from $G$ by removing $mathcal{K}$; i.e., $G^star:=Gsetminus mathcal{K}$. We refer to $G$ as an archipelago and $G^star$ as an archipelago with lakes. Denote by ${p_n(G,z)}_{n=0}^infty$ and ${p_n(G^star,z)}_{n=0}^infty$, the sequences of the Bergman polynomials associated with $G$ and $G^star$, respectively; that is, the orthonormal polynomials with respect to the area measure on $G$ and $G^star$. The purpose of the paper is to show that $p_n(G,z)$ and $p_n(G^star,z)$ have comparable asymptotic properties, thereby demonstrating that the asymptotic properties of the Bergman polynomials for $G^star$ are determined by the boundary of $G$. As a consequence we can analyze certain asymptotic properties of $p_n(G^star,z)$ by using the corresponding results for $p_n(G,z)$, which were obtained in a recent work by B. Gustafsson, M. Putinar, and two of the present authors. The results lead to a reconstruction algorithm for recovering the shape of an archipelago with lakes from a partial set of its complex moments.
A thorough backward stability analysis of Hotellings deflation, an explicit external deflation procedure through low-rank updates for computing many eigenpairs of a symmetric matrix, is presented. Computable upper bounds of the loss of the orthogonality of the computed eigenvectors and the symmetric backward error norm of the computed eigenpairs are derived. Sufficient conditions for the backward stability of the explicit external deflation procedure are revealed. Based on these theoretical results, the strategy for achieving numerical backward stability by dynamically selecting the shifts is proposed. Numerical results are presented to corroborate the theoretical analysis and to demonstrate the stability of the procedure for computing many eigenpairs of large symmetric matrices arising from applications.
In this work, we study the numerical approximation of a class of singular fully coupled forward backward stochastic differential equations. These equations have a degenerate forward component and non-smooth terminal condition. They are used, for example, in the modeling of carbon market[9] and are linked to scalar conservation law perturbed by a diffusion. Classical FBSDEs methods fail to capture the correct entropy solution to the associated quasi-linear PDE. We introduce a splitting approach that circumvent this difficulty by treating differently the numerical approximation of the diffusion part and the non-linear transport part. Under the structural condition guaranteeing the well-posedness of the singular FBSDEs [8], we show that the splitting method is convergent with a rate $1/2$. We implement the splitting scheme combining non-linear regression based on deep neural networks and conservative finite difference schemes. The numerical tests show very good results in possibly high dimensional framework.