Related papers: A study of defect-based error estimates for the Kr…
We establish quantitative bounds for H\"older exponents in the Krylov--Safonov and Evans--Krylov theories when the ellipticity ratio is close to one. Our analysis relies on the Ishii--Lions method for the Krylov--Safonov theory and a…
We present a framework that relates preconditioning with a posteriori error estimates in finite element methods. In particular, we use standard tools in subspace correction methods to obtain reliable and efficient error estimators. As a…
We present a class of algorithms based on rational Krylov methods to compute the action of a generalized matrix function on a vector. These algorithms incorporate existing methods based on the Golub-Kahan bidiagonalization as a special…
This paper introduces a new recovery-based a posteriori error estimator for the lowest order Nedelec finite element approximation to the H(curl) interface problem. The error estimator is analyzed by establishing both the reliability and the…
Adaptive quasicontinuum (QC) methods are important methodologies in molecular mechanics for the simulations of materials with defects, intending to achieve the optimal balance of accuracy and efficiency on the fly. In this study, we propose…
We derive a posteriori error estimates for the hybridizable discontinuous Galerkin (HDG) methods, including both the primal and mixed formulations, for the approximation of a linear second-order elliptic problem on conforming simplicial…
Estimators derived from a divergence criterion such as $\varphi-$divergences are generally more robust than the maximum likelihood ones. We are interested in particular in the so-called MD$\varphi$DE, an estimator built using a dual…
Numerical simulations of crystal defects are necessarily restricted to finite computational domains, supplying artificial boundary conditions that emulate the effect of embedding the defect in an effectively infinite crystalline…
The computation of f(A)b, the action of a matrix function on a vector, is a task arising in many areas of scientific computing. In many applications, the matrix A is sparse but so large that only a rather small number of Krylov basis…
Approximations of the Dirac delta distribution are commonly used to create sequences of smooth functions approximating nonsmooth (generalized) functions, via convolution. In this work, we show a priori rates of convergence of this…
In this paper we introduce and analyze the residual-based a posteriori error estimation of the partially penalized immersed finite element method for solving elliptic interface problems. The immersed finite element method can be naturally…
We propose a new practical adaptive refinement strategy for $hp$-finite element approximations of elliptic problems. Following recent theoretical developments in polynomial-degree-robust a posteriori error analysis, we solve two types of…
Scalable quantum Fisher information (QFI) estimation becomes actionable when the numerical estimate is paired with a trustworthy stopping decision. Krylov-shadow QFI estimation has two resource directions: the Krylov order sets the…
Rank-constrained matrix problems appear frequently across science and engineering. The convergence analysis of iterative algorithms developed for these problems often hinges on local error bounds, which correlate the distance to the…
We study structure-preserving Krylov subspace methods for approximating the matrix-vector products f(H)b, where H is a large Hamiltonian matrix and f denotes either the matrix exponential or the related phi-function. Such computations are…
As computational machines become larger and more complex, the probability of hardware failure rises. ``Silent errors'', or bit flips, may not be immediately apparent but can cause detrimental effects to algorithm behavior. In this work, we…
We introduced and analyzed robust recovery-based a posteriori error estimators for various lower order finite element approximations to interface problems in [9, 10], where the recoveries of the flux and/or gradient are implicit (i.e.,…
Building on the successes of local kernel methods for approximating the solutions to partial differential equations (PDE) and the evaluation of definite integrals (quadrature/cubature), a local estimate of the error in such approximations…
Mixed-dimensional elliptic equations exhibiting a hierarchical structure are commonly used to model problems with high aspect ratio inclusions, such as flow in fractured porous media. We derive general abstract estimates based on the theory…
We give a goal-oriented a posteriori error estimator for the atomistic-continuum modeling error in the quasicontinuum method, and we use this estimator to design an adaptive algorithm to compute a quantity of interest to a given tolerance…