BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Shuhuang Xiang

Shuhuang Xiang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - Emerging
7works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 published item(s)

preprint2021arXiv

Fast linear barycentric rational interpolation for singular functions via scaled transformations

In this paper, applied strictly monotonic increasing scaled maps, a kind of well-conditioned linear barycentric rational interpolations are proposed to approximate functions of singularities at the origin, such as $x^α$ for $α\in (0,1)$ and $\log(x)$. It just takes $O(N)$ flops and can achieve fast convergence rates with the choice the scaled parameter, where $N$ is the maximum degree of the denominator and numerator. The construction of the rational interpolant couples rational polynomials in the barycentric form of second kind with the transformed Jacobi-Gauss-Lobatto points. Numerical experiments are considered which illustrate the accuracy and efficiency of the algorithms. The convergence of the rational interpolation is also considered.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence rates of spectral orthogonal projection approximation for functions of algebraic and logarithmatic regularities

Based on the Hilb type formula between Jacobi polynomials and Bessel functions, optimal decay rates on Jacobi expansion coefficients are derived, by applying van der Corput type lemmas, for functions of logarithmatic singularities, which leads to the optimal convergence rates on the Jacobi, Gegenbauer and Chebyshev orthogonal projections. It is interesting to see that for boundary singularities, one may get faster convergence rate on the Jacobi or Gegenbauer projection as $(α,β)$ and $λ$ increases. The larger values of parameter, the higher convergence rates can be achieved. In particular, the Legendre projection has one half order higher than Chebyshev. Moreover, if $\min\{α,β\}>0$ and $λ>\frac{1}{2}$, the Jacobi and Gegenbauer orthogonal projections have higher convergence orders compared with Legendre. While for interior singularity, the convergence order is independent of $(α,β)$ and $λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computation of highly oscillatory Bessel transforms with algebraic singularities

In this paper, we consider the Clenshaw-Curtis-Filon method for the highly oscillatory Bessel transform $\int_0^1x^α(1-x)^βf(x) J_ν(ωx)dx$, where $f$ is a smooth function on $[0, 1]$, and $ν\geq0.$ The method is based on Fast Fourier Transform (FFT) and fast computation of the modified moments. We give a recurrence relation for the modified moments and present an efficient method for the evaluation of modified moments by using recurrence relation. Moreover, the corresponding error bound in inverse powers of $N$ for this method for the integral is presented. Numerical examples are provided to support our analysis and show the efficiency and accuracy of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On Interpolation Approximation: Convergence rates for polynomial interpolation for functions of limited regularity

The convergence rates on polynomial interpolation in most cases are estimated by Lebesgue constants. These estimates may be overestimated for some special points of sets for functions of limited regularities. In this paper, by applying the Peano kernel theorem and Wainerman's lemma, new formulas on the convergence rates are considered. Based upon these new estimates, it shows that the interpolation at strongly normal pointsystems can achieve the optimal convergence rate, the same as the best polynomial approximation. Furthermore, by using the asymptotics on Jacobi polynomials, the convergence rates are established for Gauss-Jacobi, Jacobi-Gauss-Lobatto or Jacobi-Gauss-Radau pointsystems. From these results, we see that the interpolations at the Gauss-Legendre, Legendre-Gauss-Lobatto pointsystem, or at strongly normal pointsystems, has essentially the same approximation accuracy compared with those at the two Chebyshev piontsystems, which also illustrates the equally accuracy of the Gauss and Clenshaw-Curtis quadrature. In addition, numerical examples illustrate the perfect coincidence with the estimates, which means the convergence rates are optimal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Fast Implementation of Higher Order Hermite-Fejer Interpolation

The problem of barycentric Hermite interpolation is highly susceptible to overflows or underflows. In this paper, based on Sturm-Liouville equations for Jacobi orthogonal polynomials, we consider the fast implementation on the second barycentric formula for higher order Hermite-Fejér interpolation at Gauss-Jacobi or Jacobi-Gauss-Lobatto pointsystems, where the barycentric weights can be efficiently evaluated and cost linear operations corresponding to the number of grids totally. Furthermore, due to the division of the second barycentric form, the exponentially increasing common factor in the barycentric weights can be canceled, which yields a superiorly stable method for computing the simplified barycentric weights, and leads to a fast implementation of the higher order Hermite-Fejér interpolation with linear operations on the number of grids. In addition, the convergence rates are derived for Hermite-Fejér interpolation at Gauss-Jacobi pointsystems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Optimal Rates of Convergence for Quadratures Derived from Chebyshev Points

In this paper, we study the optimal general convergence rates for quadratures derived from Chebyshev points. By building on the aliasing errors on integration of Chebyshev polynomials, together with the asymptotic formulae on the coefficients of Chebyshev expansions, new and optimal convergence rates for $n$-point Clenshaw-Curtis, Fejér's first and second quadrature rules are established for Jacobi weights or Jacobi weights multiplied by $\ln((x+1)/2)$. The convergence orders are attainable for some functions of finite regularities. In addition, by using refined estimates on aliasing errors on integration of Chebyshev polynomials by Gauss-Legendre quadrature, an improved convergence rate for Gauss-Legendre is given too.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the convergence rates of Gauss and Clenshaw-Curtis quadrature for functions of limited regularity

We study the optimal general rate of convergence of the n-point quadrature rules of Gauss and Clenshaw-Curtis when applied to functions of limited regularity: if the Chebyshev coefficients decay at a rate O(n^{-s-1}) for some s > 0, Clenshaw-Curtis and Gauss quadrature inherit exactly this rate. The proof (for Gauss, if 0 < s < 2, there is numerical evidence only) is based on work of Curtis, Johnson, Riess, and Rabinowitz from the early 1970s and on a refined estimate for Gauss quadrature applied to Chebyshev polynomials due to Petras (1995). The convergence rate of both quadrature rules is up to one power of n better than polynomial best approximation; hence, the classical proof strategy that bounds the error of a quadrature rule with positive weights by polynomial best approximation is doomed to fail in establishing the optimal rate.

0 citations0 reviewsNo code linkedOpen access