BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Haiyong Wang

Haiyong Wang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.CANumerical Analysis

Catalog footprint

What is connected

9works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2023arXiv

Convergence analysis of Laguerre approximations for analytic functions

Laguerre spectral approximations play an important role in the development of efficient algorithms for problems in unbounded domains. In this paper, we present a comprehensive convergence rate analysis of Laguerre spectral approximations for analytic functions. By exploiting contour integral techniques from complex analysis, we prove that Laguerre projection and interpolation methods of degree $n$ converge at the root-exponential rate $O(\exp(-2ρ\sqrt{n}))$ with $ρ>0$ when the underlying function is analytic inside and on a parabola with focus at the origin and vertex at $z=-ρ^2$. As far as we know, this is the first rigorous proof of root-exponential convergence of Laguerre approximations for analytic functions. Several important applications of our analysis are also discussed, including Laguerre spectral differentiations, Gauss-Laguerre quadrature rules, the scaling factor and the Weeks method for the inversion of Laplace transform, and some sharp convergence rate estimates are derived. Numerical experiments are presented to verify the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analysis of multivariate Gegenbauer approximation in the hypercube

In this paper, we are concerned with multivariate Gegenbauer approximation of functions defined in the $d$-dimensional hypercube. Two new and sharper bounds for the coefficients of multivariate Gegenbauer expansion of analytic functions are presented based on two different extensions of the Bernstein ellipse. We then establish an explicit error bound for the multivariate Gegenbauer approximation associated with an $\ell^q$ ball index set in the uniform norm. We also consider the multivariate approximation of functions with finite regularity and derive the associated error bound on the full grid in the uniform norm. As an application, we extend our arguments to obtain some new tight bounds for the coefficients of tensorized Legendre expansions in the context of polynomial approximation of parameterized PDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A unified framework for asymptotic analysis and computation of finite Hankel transform

In this paper we present a unified framework for asymptotic analysis and computation of the finite Hankel transform. This framework enables us to derive asymptotic expansions of the transform, including the cases where the oscillator has zeros and stationary points. As a consequence, two efficient and affordable methods for computing the transform numerically are developed and a detailed analysis of their asymptotic error estimate is carried out. Numerical examples are provided to confirm our analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast and highly accurate computation of Chebyshev expansion coefficients of analytic functions

Chebyshev expansion coefficients can be computed efficiently by using the FFT, and for smooth functions the resulting approximation is close to optimal, with computations that are numerically stable. Given sufficiently accurate function samples, the Chebyshev expansion coefficients can be computed to machine precision accuracy. However, the accuracy is only with respect to absolute error, and this implies that very small expansion coefficients typically have very large relative error. Upon differentiating a Chebyshev expansion, this relative error in the small coefficients is magnified and accuracy may be lost, especially after repeated differentiation. At first sight, this seems unavoidable. Yet, in this paper, we focus on an alternative computation of Chebyshev expansion coefficients using contour integrals in the complex plane. The main result is that the coefficients can be computed with machine precision relative error, rather than absolute error. This implies that even very small coefficients can be computed with full floating point accuracy, even when they are themselves much smaller than machine precision. As a result, no accuracy is lost after differentiating the expansion, and even the 100th derivative of an analytic function can be computed with near machine precision accuracy using standard floating point arithmetic. In some cases, the contour integrals can be evaluated using the FFT, making the approach both highly accurate and fast.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence rate and acceleration of Clenshaw-Curtis quadrature for functions with endpoint singularities

In this paper, we investigate the rate of convergence of Clenshaw-Curtis quadrature and its acceleration for functions with endpoint singularities in X^s, where X^s denotes the space of functions whose Chebyshev coefficients decay asymptotically as a_k = O(k^{-s-1}) for some positive s. For such unctions, we show that the convergence rate of (n + 1)-point Clenshaw-Curtis quadrature is O(n^{-s-2}). Furthermore, an asymptotic error expansion for Clenshaw-Curtis quadrature is presented which enables us to employ some extrapolation techniques to accelerate its convergence. Numerical examples are provided to confirm our analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On Fast Implementation of Clenshaw-Curtis and Fejér-type Quadrature Rules

Based upon the fast computation of the coefficients of the interpolation polynomials at Chebyshev-type points by FFT, DCT and IDST, respectively, together with the efficient evaluation of the modified moments by forwards recursions or by the Oliver's algorithm, this paper presents interpolating integration algorithms, by using the coefficients and modified moments, for Clenshaw-Curtis, Fejér's first and second-type rules for Jacobi or Jacobi weights multiplied by a logarithmic function. The corresponding {\sc Matlab} codes are included. Numerical examples illustrate the stability, accuracy of the Clenshaw-Curtis, Fejér's first and second rules, and show that the three quadratures have nearly the same convergence rates as Gauss-Jacobi quadrature for functions of finite regularities for Jacobi weights, and are more efficient upon the cpu time than the Gauss evaluated by fast computation of the weights and nodes by {\sc Chebfun}.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Gaussian quadrature rule for oscillatory integrals on a bounded interval

We investigate a Gaussian quadrature rule and the corresponding orthogonal polynomials for the oscillatory weight function $e^{iωx}$ on the interval $[-1,1]$. We show that such a rule attains high asymptotic order, in the sense that the quadrature error quickly decreases as a function of the frequency $ω$. However, accuracy is maintained for all values of $ω$ and in particular the rule elegantly reduces to the classical Gauss-Legendre rule as $ω\to 0$. The construction of such rules is briefly discussed, and though not all orthogonal polynomials exist, it is demonstrated numerically that rules with an even number of points are always well defined. We show that these rules are optimal both in terms of asymptotic order as well as in terms of polynomial order.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Explicit barycentric weights for polynomial interpolation in the roots or extrema of classical orthogonal polynomials

Barycentric interpolation is arguably the method of choice for numerical polynomial interpolation. The polynomial interpolant is expressed in terms of function values using the so-called barycentric weights, which depend on the interpolation points. Few explicit formulae for these barycentric weights are known. In [H. Wang and S. Xiang, Math. Comp., 81 (2012), 861--877], the authors have shown that the barycentric weights of the roots of Legendre polynomials can be expressed explicitly in terms of the weights of the corresponding Gaussian quadrature rule. This idea was subsequently implemented in the Chebfun package [L. N. Trefethen and others, The Chebfun Development Team, 2011] and in the process generalized by the Chebfun authors to the roots of Jacobi, Laguerre and Hermite polynomials. In this paper, we explore the generality of the link between barycentric weights and Gaussian quadrature and show that such relationships are related to the existence of lowering operators for orthogonal polynomials. We supply an exhaustive list of cases, in which all known formulae are recovered and also some new formulae are derived, including the barycentric weights for Gauss-Radau and Gauss-Lobatto points. Based on a fast ${\mathcal O}(n)$ algorithm for the computation of Gaussian quadrature, due to Hale and Townsend, this leads to an ${\mathcal O}(n)$ computational scheme for barycentric weights.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Asymptotic expansions and fast computation of oscillatory Hilbert transforms

In this paper, we study the asymptotics and fast computation of the one-sided oscillatory Hilbert transforms of the form $$H^{+}(f(t)e^{iωt})(x)=-int_{0}^{\infty}e^{iωt}\frac{f(t)}{t-x}dt,\qquad ω>0,\qquad x\geq 0,$$ where the bar indicates the Cauchy principal value and $f$ is a real-valued function with analytic continuation in the first quadrant, except possibly a branch point of algebraic type at the origin. When $x=0$, the integral is interpreted as a Hadamard finite-part integral, provided it is divergent. Asymptotic expansions in inverse powers of $ω$ are derived for each fixed $x\geq 0$, which clarify the large $ω$ behavior of this transform. We then present efficient and affordable approaches for numerical evaluation of such oscillatory transforms. Depending on the position of $x$, we classify our discussion into three regimes, namely, $x=\mathcal{O}(1)$ or $x\gg1$, $0<x\ll 1$ and $x=0$. Numerical experiments show that the convergence of the proposed methods greatly improve when the frequency $ω$ increases. Some extensions to oscillatory Hilbert transforms with Bessel oscillators are briefly discussed as well.

0 citations0 reviewsNo code linkedOpen access