We prove sharp wavenumber-explicit error bounds for first- or second-type-N\'ed\'elec-element (a.k.a. edge-element) conforming discretisations, of arbitrary (fixed) order, of the variable-coefficient time-harmonic Maxwell equations posed in a bounded domain with perfect electric conductor (PEC) boundary conditions. The PDE coefficients are allowed to be piecewise regular and complex-valued; this set-up therefore includes scattering from a PEC obstacle and/or variable real-valued coefficients, with the radiation condition approximated by a perfectly matched layer (PML). In the analysis of the $h$-version of the finite-element method, with fixed polynomial degree $p$, applied to the time-harmonic Maxwell equations, the $\textit{asymptotic regime}$ is when the meshwidth, $h$, is small enough (in a wavenumber-dependent way) that the Galerkin solution is quasioptimal independently of the wavenumber, while the $\textit{preasymptotic regime}$ is the complement of the asymptotic regime. The results of this paper are the first preasymptotic error bounds for the time-harmonic Maxwell equations using first-type N\'ed\'elec elements or higher-than-lowest-order second-type N\'ed\'elec elements. Furthermore, they are the first wavenumber-explicit results, even in the asymptotic regime, for Maxwell scattering problems with a non-empty scatterer.
相關內容
We discuss the second-order differential uniformity of vectorial Boolean functions. The closely related notion of second-order zero differential uniformity has recently been studied in connection to resistance to the boomerang attack. We prove that monomial functions with univariate form $x^d$ where $d=2^{2k}+2^k+1$ and $\gcd(k,n)=1$ have optimal second-order differential uniformity. Computational results suggest that, up to affine equivalence, these might be the only optimal cubic power functions. We begin work towards generalising such conditions to all monomial functions of algebraic degree 3. We also discuss further questions arising from computational results.
Finding vertex-to-vertex correspondences in real-world graphs is a challenging task with applications in a wide variety of domains. Structural matching based on graphs connectivities has attracted considerable attention, while the integration of all the other information stemming from vertices and edges attributes has been mostly left aside. Here we present the Graph Attributes and Structure Matching (GASM) algorithm, which provides high-quality solutions by integrating all the available information in a unified framework. Parameters quantifying the reliability of the attributes can tune how much the solutions should rely on the structure or on the attributes. We further show that even without attributes GASM consistently finds as-good-as or better solutions than state-of-the-art algorithms, with similar processing times.
We formulate a uniform tail bound for empirical processes indexed by a class of functions, in terms of the individual deviations of the functions rather than the worst-case deviation in the considered class. The tail bound is established by introducing an initial ``deflation'' step to the standard generic chaining argument. The resulting tail bound is the sum of the complexity of the ``deflated function class'' in terms of a generalization of Talagrand's $\gamma$ functional, and the deviation of the function instance, both of which are formulated based on the natural seminorm induced by the corresponding Cram\'{e}r functions. Leveraging another less demanding natural seminorm, we also show similar bounds, though with implicit dependence on the sample size, in the more general case where finite exponential moments cannot be assumed. We also provide approximations of the tail bounds in terms of the more prevalent Orlicz norms or their ``incomplete'' versions under suitable moment conditions.
In causal inference on directed acyclic graphs, the orientation of edges is in general only recovered up to Markov equivalence classes. We study Markov equivalence classes of uniformly random directed acyclic graphs. Using a tower decomposition, we show that the ratio between the number of Markov equivalence classes and directed acyclic graphs approaches a positive constant when the number of sites goes to infinity. For a typical directed acyclic graph, the expected number of elements in its Markov equivalence class remains bounded. More precisely, we prove that for a uniformly chosen directed acyclic graph, the size of its Markov equivalence class has super-polynomial tails.
We derive a robust error estimate for a recently proposed numerical method for $\alpha$-dissipative solutions of the Hunter-Saxton equation, where $\alpha \in [0, 1]$. In particular, if the following two conditions hold: i) there exist a constant $C > 0$ and $\beta \in (0, 1]$ such that the initial spatial derivative $\bar{u}_{x}$ satisfies $\|\bar{u}_x(\cdot + h) - \bar{u}_x(\cdot)\|_2 \leq Ch^{\beta}$ for all $h \in (0, 2]$, and ii), the singular continuous part of the initial energy measure is zero, then the numerical wave profile converges with order $O(\Delta x^{\frac{\beta}{8}})$ in $L^{\infty}(\mathbb{R})$. Moreover, if $\alpha=0$, then the rate improves to $O(\Delta x^{\frac{1}{4}})$ without the above assumptions, and we also obtain a convergence rate for the associated energy measure - it converges with order $O(\Delta x^{\frac{1}{2}})$ in the bounded Lipschitz metric. These convergence rates are illustrated by several examples.
We propose a randomized lattice algorithm for approximating multivariate periodic functions over the $d$-dimensional unit cube from the weighted Korobov space with mixed smoothness $\alpha > 1/2$ and product weights $\gamma_1,\gamma_2,\ldots\in [0,1]$. Building upon the deterministic lattice algorithm by Kuo, Sloan, and Wo\'{z}niakowski (2006), we incorporate a randomized quadrature rule by Dick, Goda, and Suzuki (2022) to accelerate the convergence rate. This randomization involves drawing the number of points for function evaluations randomly, and selecting a good generating vector for rank-1 lattice points using the randomized component-by-component algorithm. We prove that our randomized algorithm achieves a worst-case root mean squared $L_2$-approximation error of order $M^{-\alpha/2 - 1/8 + \varepsilon}$ for an arbitrarily small $\varepsilon > 0$, where $M$ denotes the maximum number of function evaluations, and that the error bound is independent of the dimension $d$ if the weights satisfy $\sum_{j=1}^\infty \gamma_j^{1/\alpha} < \infty$. Our upper bound converges faster than a lower bound on the worst-case $L_2$-approximation error for deterministic rank-1 lattice-based approximation proved by Byrenheid, K\"{a}mmerer, Ullrich, and Volkmer (2017). We also show a lower error bound of order $M^{-\alpha/2-1/2}$ for our randomized algorithm, leaving a slight gap between the upper and lower bounds open for future research.
The main purpose of this paper is to design a local discontinuous Galerkin (LDG) method for the Benjamin-Ono equation. We analyze the stability and error estimates for the semi-discrete LDG scheme. We prove that the scheme is $L^2$-stable and it converges at a rate $\mathcal{O}(h^{k+1/2})$ for general nonlinear flux. Furthermore, we develop a fully discrete LDG scheme using the four-stage fourth order Runge-Kutta method and ensure the devised scheme is strongly stable in case of linear flux using two-step and three-step stability approach under an appropriate time step constraint. Numerical examples are provided to validate the efficiency and accuracy of the method.
We present a fully discrete Crank-Nicolson Fourier-spectral-Galerkin (FSG) scheme for approximating solutions of the fractional Korteweg-de Vries (KdV) equation, which involves a fractional Laplacian with exponent $\alpha \in [1,2]$ and a small dispersion coefficient of order $\varepsilon^2$. The solution in the limit as $\varepsilon \to 0$ is known as the zero dispersion limit. We demonstrate that the semi-discrete FSG scheme conserves the first three integral invariants, thereby structure preserving, and that the fully discrete FSG scheme is $L^2$-conservative, ensuring stability. Using a compactness argument, we constructively prove the convergence of the approximate solution to the unique solution of the fractional KdV equation in $C([0,T]; H_p^{1+\alpha}(\mathbb{R}))$ for the periodic initial data in $H_p^{1+\alpha}(\mathbb{R})$. The devised scheme achieves spectral accuracy for the initial data in $H_p^r,$ $r \geq 1+\alpha$ and exponential accuracy for the analytic initial data. Additionally, we establish that the approximation of the zero dispersion limit obtained from the fully discrete FSG scheme converges to the solution of the Hopf equation in $L^2$ as $\varepsilon \to 0$, up to the gradient catastrophe time $t_c$. Beyond $t_c$, numerical investigations reveal that the approximation converges to the asymptotic solution, which is weakly described by the Whitham's averaged equation within the oscillatory zone for $\alpha = 2$. Numerical results are provided to demonstrate the convergence of the scheme and to validate the theoretical findings.
Oscillatory second order linear ordinary differential equations arise in many scientific calculations. Because the running times of standard solvers increase linearly with frequency when they are applied to such problems, a variety of specialized methods, most of them quite complicated, have been proposed. Here, we point out that one of the simplest approaches not only works, but yields a scheme for solving oscillatory second order linear ordinary differential equations which is significantly faster than current state-of-the-art techniques. Our method, which operates by constructing a slowly varying phase function representing a basis of solutions of the differential equation, runs in time independent of the frequency and can be applied to second order equations whose solutions are oscillatory in some regions and slowly varying in others. In the high-frequency regime, our algorithm discretizes the nonlinear Riccati equation satisfied by the derivative of the phase function via a Chebyshev spectral collocation method and applies the Newton-Kantorovich method to the resulting system of nonlinear algebraic equations. We prove that the iterates converge quadratically to a nonoscillatory solution of the Riccati equation. The quadratic convergence of the Newton-Kantorovich method and the simple form of the linearized equations ensure that this procedure is extremely efficient. Our algorithm then extends the slowly varying phase function calculated in the high-frequency regime throughout the solution domain by solving a certain third order linear ordinary differential equation related to the Riccati equation. We describe the results of numerical experiments showing that our algorithm is orders of magnitude faster than existing schemes, including the modified Magnus method [18], the current state-of-the-art approach [7] and the recently introduced ARDC method [1].
We consider non-linear Bayesian inverse problems of determining the parameter $f$. For the posterior distribution with a class of Gaussian process priors, we study the statistical performance of variational Bayesian inference to the posterior with variational sets consisting of Gaussian measures or a mean-field family. We propose certain conditions on the forward map $\mathcal{G}$, the variational set $\mathcal{Q}$ and the prior such that, as the number $N$ of measurements increases, the resulting variational posterior distributions contract to the ground truth $f_0$ generating the data, and derive a convergence rate with polynomial order or logarithmic order. As specific examples, we consider a collection of non-linear inverse problems, including the Darcy flow problem, the inverse potential problem for a subdiffusion equation, and the inverse medium scattering problem. Besides, we show that our convergence rates are minimax optimal for these inverse problems.
注冊地址： 北京市海淀區羊坊店路18號2幢3層301-191