We consider the convergence of adaptive BEM for weakly-singular and hypersingular integral equations associated with the Laplacian and the Helmholtz operator in 2D and 3D. The local mesh-refinement is driven by some two-level error estimator. We show that the adaptive algorithm drives the underlying error estimates to zero. Moreover, we prove that the saturation assumption already implies linear convergence of the error with optimal algebraic rates.
In this paper, we study an adaptive planewave method for multiple eigenvalues of second-order elliptic partial equations. Inspired by the technique for the adaptive finite element analysis, we prove that the adaptive planewave method has the linear convergence rate and optimal complexity.
The paper considers a class of parametric elliptic partial differential equations (PDEs), where the coefficients and the right-hand side function depend on infinitely many (uncertain) parameters. We introduce a two-level a posteriori estimator to control the energy error in multilevel stochastic Galerkin approximations for this class of PDE problems. We prove that the two-level estimator always provides a lower bound for the unknown approximation error, while the upper bound is equivalent to a saturation assumption. We propose and empirically compare three adaptive algorithms, where the structure of the estimator is exploited to perform spatial refinement as well as parametric enrichment. The paper also discusses implementation aspects of computing multilevel stochastic Galerkin approximations.
In recent years, contour-based eigensolvers have emerged as a standard approach for the solution of large and sparse eigenvalue problems. Building upon recent performance improvements through non-linear least square optimization of so-called rational filters, we introduce a systematic method to design these filters by minimizing the worst-case convergence ratio and eliminate the parametric dependence on weight functions. Further, we provide an efficient way to deal with the box-constraints which play a central role for the use of iterative linear solvers in contour-based eigensolvers. Indeed, these parameter-free filters consistently minimize the number of iterations and the number of FLOPs to reach convergence in the eigensolver. As a byproduct, our rational filters allow for a simple solution to load balancing when the solution of an interior eigenproblem is approached by the slicing of the sought after spectral interval.
This paper presents an enhanced version of our previous work, hybrid non-uniform subdivision surfaces [19], to achieve optimal convergence rates in isogeometric analysis. We introduce a parameter $lambda$ ($frac{1}{4}<lambda<1$) to control the rate of shrinkage of irregular regions, so the method is called tuned hybrid non-uniform subdivision (tHNUS). Our previous work corresponds to the case when $lambda=frac{1}{2}$. While introducing $lambda$ in hybrid subdivision significantly complicates the theoretical proof of $G^1$ continuity around extraordinary vertices, reducing $lambda$ can recover the optimal convergence rates when tuned hybrid subdivision functions are used as a basis in isogeometric analysis. From the geometric point of view, the tHNUS retains comparable shape quality as [19] under non-uniform parameterization. Its basis functions are refinable and the geometric mapping stays invariant during refinement. Moreover, we prove that a tuned hybrid subdivision surface is globally $G^1$-continuous. From the analysis point of view, tHNUS basis functions form a non-negative partition of unity, are globally linearly independent, and their spline spaces are nested. We numerically demonstrate that tHNUS basis functions can achieve optimal convergence rates for the Poissons problem with non-uniform parameterization around extraordinary vertices.
Convergence of an adaptive collocation method for the stationary parametric diffusion equation with finite-dimensional affine coefficient is shown. The adaptive algorithm relies on a recently introduced residual-based reliable a posteriori error estimator. For the convergence proof, a strategy recently used for a stochastic Galerkin method with an hierarchical error estimator is transferred to the collocation setting. Extensions to other variants of adaptive collocation methods (including the classical one proposed in the paper Dimension-adaptive tensor-product quadratuture Computing (2003) by T. Gerstner and M. Griebel) is explored.