Martín Maas

Gabriel Acosta, Juan Pablo Borthagaray, Oscar Bruno et al.

This paper presents regularity results and associated high-order numerical methods for one-dimensional Fractional-Laplacian boundary-value problems. On the basis of a factorization of solutions as a product of a certain edge-singular weight $ω$ times a "regular" unknown, a characterization of the regularity of solutions is obtained in terms of the smoothness of the corresponding right-hand sides. In particular, for right-hand sides which are analytic in a Bernstein Ellipse, analyticity in the same Bernstein Ellipse is obtained for the "regular" unknown. Moreover, a sharp Sobolev regularity result is presented which completely characterizes the co-domain of the Fractional-Laplacian operator in terms of certain weighted Sobolev spaces introduced in (Babuška and Guo, SIAM J. Numer. Anal. 2002). The present theoretical treatment relies on a full eigendecomposition for a certain weighted integral operator in terms of the Gegenbauer polynomial basis. The proposed Gegenbauer-based Nyström numerical method for the Fractional-Laplacian Dirichlet problem, further, is significantly more accurate and efficient than other algorithms considered previously. The sharp error estimates presented in this paper indicate that the proposed algorithm is spectrally accurate, with convergence rates that only depend on the smoothness of the right-hand side. In particular, convergence is exponentially fast (resp. faster than any power of the mesh-size) for analytic (resp. in nitely smooth) right-hand sides. The properties of the algorithm are illustrated with a variety of numerical results.

Oscar Bruno, Martín Maas

This paper introduces a fast algorithm, applicable throughout the electromagnetic spectrum, for the numerical solution of problems of scattering by periodic surfaces in two-dimensional space. The proposed algorithm remains highly accurate and efficient for challenging configurations including randomly rough surfaces, deep corrugations, large periods, near grazing incidences, and, importantly, Wood-anomaly resonant frequencies. The proposed approach is based on use of certain `shifted equivalent sources' which enable FFT acceleration of a Wood-anomaly-capable quasi-periodic Green function introduced recently (Bruno and Delourme, Jour. Computat. Phys., 262--290, 2014). The Green-function strategy additionally incorporates an exponentially convergent shifted version of the classical spectral series for the Green function. While the computing-cost asymptotics depend on the asymptotic configuration assumed, the computing costs rise at most linearly with the size of the problem for a number of important rough-surface cases we consider. In practice, single-core runs in computing times ranging from a fraction of a second to a few seconds suffice for the proposed algorithm to produce highly-accurate solutions in some of the most challenging contexts arising in applications.