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

Xiangmin Jiao

Xiangmin Jiao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.AP

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

5 published item(s)

preprint2022arXiv

ARPIST: Provably Accurate and Stable Numerical Integration over Spherical Triangles

Numerical integration on spheres, including the computation of the areas of spherical triangles, is a core computation in geomathematics. The commonly used techniques sometimes suffer from instabilities and significant loss of accuracy. We describe a new algorithm, called ARPIST, for accurate and stable integration of functions on spherical triangles. ARPIST is based on an easy-to-implement transformation to the spherical triangle from its corresponding linear triangle via radial projection to achieve high accuracy and efficiency. More importantly, ARPIST overcomes potential instabilities in computing the Jacobian of the transformation, even for poorly shaped triangles that may occur at poles in regular longitude-latitude meshes, by avoiding potential catastrophic rounding errors. We compare our proposed technique with L'Huilier's Theorem for computing the area of spherical triangles, and also compare it with the recently developed LSQST method (J. Beckmann, H.N. Mhaskar, and J. Prestin, GEM - Int. J. Geomath., 5:143-162, 2014) and a radial-basis-function-based technique (J. A. Reeger and B. Fornberg, Stud. Appl. Math., 137:174-188, 2015) for integration of smooth functions on spherical triangulations. Our results show that ARPIST enables superior accuracy and stability over previous methods while being orders of magnitude faster and significantly easier to implement.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Approximate Generalized Inverses with Iterative Refinement for $ε$-Accurate Preconditioning of Singular Systems

We introduce a new class of preconditioners to enable flexible GMRES to find a least-squares solution, and potentially the pseudoinverse solution, of large-scale sparse, asymmetric, singular, and potentially inconsistent systems. We develop the preconditioners based on a new observation that generalized inverses (i.e., $\boldsymbol{A}^{g}\in\{\boldsymbol{G}\mid\boldsymbol{A}\boldsymbol{G}\boldsymbol{A}=\boldsymbol{A}\}$) enable the preconditioned Krylov subspaces to converge in a single step. We then compute an approximate generalized inverse (AGI) efficiently using a hybrid incomplete factorization (HIF), which combines multilevel incomplete LU with rank-revealing QR on its final Schur complement. We define the criteria of $ε$-accuracy and stability of AGI to guarantee the convergence of preconditioned GMRES for consistent systems. For inconsistent systems, we fortify HIF with iterative refinement to obtain HIFIR, which allows accurate computations of the null-space vectors. By combining the two techniques, we then obtain a new solver, called PIPIT, for obtaining the pseudoinverse solutions for systems with low-dimensional null spaces. We demonstrate the robustness of HIF and HIFIR and show that they improve both accuracy and efficiency of the prior state of the art by orders of magnitude for systems with up to a million unknowns.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimal and Low-Memory Near-Optimal Preconditioning of Fully Implicit Runge-Kutta Schemes for Parabolic PDEs

Runge-Kutta (RK) schemes, especially Gauss-Legendre and some other fully implicit RK (FIRK) schemes, are desirable for the time integration of parabolic partial differential equations due to their A-stability and high-order accuracy. However, it is significantly more challenging to construct optimal preconditioners for them compared to diagonally implicit RK (or DIRK) schemes. To address this challenge, we first introduce mathematically optimal preconditioners called block complex Schur decomposition (BCSD), block real Schur decomposition (BRSD), and block Jordan form (BJF), motivated by block-circulant preconditioners and Jordan form solution techniques for IRK. We then derive an efficient, near-optimal singly-diagonal approximate BRSD (SABRSD) by approximating the quasi-triangular matrix in real Schur decomposition using an optimized upper-triangular matrix with a single diagonal value. A desirable feature of SABRSD is that it has comparable memory requirements and factorization (or setup) cost as singly DIRK (SDIRK). We approximate the diagonal blocks in these preconditioning techniques using an incomplete factorization with (near) linear complexity, such as multilevel ILU, ILU(0), or a multigrid method with an ILU-based smoother. We apply the block preconditioners in right-preconditioned GMRES to solve the advection-diffusion equation in 3D using finite element and finite difference methods. We show that BCSD, BRSD, and BJF significantly outperform other preconditioners in terms of GMRES iterations, and SABRSD is competitive with them and the prior state of the art in terms of computational cost while requiring the least amount of memory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Hybrid Method and Unified Analysis of Generalized Finite Differences and Lagrange Finite Elements

Finite differences, finite elements, and their generalizations are widely used for solving partial differential equations, and their high-order variants have respective advantages and disadvantages. Traditionally, these methods are treated as different (strong vs. weak) formulations and are analyzed using different techniques (Fourier analysis or Green's functions vs. functional analysis), except for some special cases on regular grids. Recently, the authors introduced a hybrid method, called Adaptive Extended Stencil FEM or AES-FEM (Int. J. Num. Meth. Engrg., 2016, DOI:10.1002/nme.5246), which combines features of generalized finite differences and Lagrange finite elements to achieve second-order accuracy over unstructured meshes. However, its analysis was incomplete due to the lack of existing mathematical theory that unifies the formulations and analysis of these different methods. In this work, we introduce the framework of generalized weighted residuals to unify the formulation of finite differences, finite elements, and AES-FEM. In addition, we propose a unified analysis of the well-posedness, convergence, and mesh-quality dependency of these different methods. We also report numerical results with AES-FEM to verify our analysis. We show that AES-FEM improves the accuracy of generalized finite differences while reducing the mesh-quality dependency and simplifying the implementation of high-order finite elements.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

WLS-ENO Remap: Superconvergent and Non-Oscillatory Weighted Least Squares Data Transfer on Surfaces

Data remap between non-matching meshes is a critical step in multiphysics coupling using a partitioned approach. The data fields being transferred often have jumps in function values or derivatives. It is important but very challenging to avoid spurious oscillations (a.k.a. the Gibbs Phenomenon) near discontinuities and at the same time to achieve high-order accuracy away from discontinuities. In this work, we introduce a new approach, called WLS-ENOR, or Weighted-Least-Squares-based Essentially Non-Oscillatory Remap, to address this challenge. Based on the WLS-ENO reconstruction technique proposed by Liu and Jiao (J. Comput. Phys. vol 314, pp 749--773, 2016), WLS-ENOR differs from WLS-ENO and other WENO schemes in that it resolves not only the O(1) oscillations due to C 0 discontinuities, but also the accumulated effect of O(h) oscillations due to C 1 discontinuities. To this end, WLS-ENOR introduces a robust detector of discontinuities and a new weighting scheme for WLS-ENO near discontinuities. We also optimize the weights at smooth regions to achieve superconvergence. As a result, WLS-ENOR is more than fifth-order accurate and highly conservative in smooth regions, while being non-oscillatory and minimally diffusive near discontinuities. We also compare WLS-ENOR with some commonly used methods based on L 2 projection, moving least squares, and radial basis functions.

0 citations0 reviewsNo code linkedOpen access