Source author record

Tsung-Ming Huang

Tsung-Ming Huang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

5works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Convergence Analysis of Dirichlet Energy Minimization for Spherical Conformal Parameterizations

In this paper, we first derive a theoretical basis for spherical conformal parameterizations between a simply connected closed surface $\mathcal{S}$ and a unit sphere $\mathbb{S}^2$ by minimizing the Dirichlet energy on $\overline{\mathbb{C}}$ by stereographic projection. The Dirichlet energy can be rewritten as the sum of the energies associated with the southern and northern hemispheres and can be decreased under an equivalence relation by alternatingly solving the corresponding Laplacian equations. Based on this theoretical foundation, we develop a modified Dirichlet energy minimization with nonequivalence deflation for the computation of the spherical conformal parameterization between $\mathcal{S}$ and $\mathbb{S}^2$. In addition, under some mild conditions, we verify the asymptotically R-linear convergence of the proposed algorithm. Numerical experiments on various benchmarks confirm that the assumptions for convergence always hold and indicate the efficiency, reliability and robustness of the developed modified Dirichlet energy minimization.

preprint2022arXiv

Fundamental Theory and R-linear Convergence of Stretch Energy Minimization for Equiareal parameterizations

In this paper, we first extend the finite distortion problem from the bounded domains in $\mathbb{R}^2$ to the closed genus-zero surfaces in $\mathbb{R}^3$ by the stereographic projection. Then we derive a theoretical foundation for spherical equiareal parameterizations between a simply connected closed surface $\mathcal{M}$ and a unit sphere $\mathbb{S}^2$ via minimizing the total area distortion energy on $\overline{\mathbb{C}}$. Provided we determine the minimizer of the total area distortion energy, the minimizer composed with the initial conformal map determines the equiareal map between the extended planes. Taking the inverse stereographic projection, we can derive the equiareal map between $\mathcal{M}$ and $\mathbb{S}^2$. The total area distortion energy can be rewritten into the sum of Dirichlet energies associated with the southern and northern hemispheres, respectively, and can be decreased by alternatingly solving the corresponding Laplacian equations. Based on this foundational theory, we develop a modified stretch energy minimization for the computation of the spherical equiareal parameterization between $\mathcal{M}$ and $\mathbb{S}^2$. In addition, under some mild conditions, we verify that our proposed algorithm has asymptotically R-linear convergence or forms a quasi-periodic solution. Numerical experiments on various benchmarks validate the assumptions for convergence always hold and indicate the efficiency, reliability and robustness of the developed modified stretch energy minimization.

preprint2016arXiv

Solving large-scale nonlinear eigenvalue problems by rational interpolation approach and resolvent sampling based Rayleigh-Ritz method

Numerical solution of nonlinear eigenvalue problems (NEPs) is frequently encountered in computational science and engineering. The applicability of most existing methods is limited by matrix structures, property of eigen-solutions, size of the problem, etc. This paper aims to break those limitations and to develop robust and universal NEP solvers for large-scale engineering applications. The novelty lies in two aspects. First, a rational interpolation approach (RIA) is proposed based on the Keldysh theorem for holomorphic matrix functions. Comparing with the existing contour integral approach (CIA), the RIA provides the possibility to select sampling points in more general regions and has advantages in improving accuracy and reducing computational cost. Second, a resolvent sampling scheme using the RIA is proposed for constructing reliable search spaces for the Rayleigh-Ritz procedure, based on which a robust eigen-solver, denoted by RSRR, is developed for solving general NEPs. RSRR can be easily implemented and parallelized. The advantages of the RIA and the performance of RSRR are demonstrated by a variety of benchmark and practical problems.

preprint2015arXiv

On the Convergence of Ritz Pairs and Refined Ritz Vectors for Quadratic Eigenvalue Problems

For a given subspace, the Rayleigh-Ritz method projects the large quadratic eigenvalue problem (QEP) onto it and produces a small sized dense QEP. Similar to the Rayleigh-Ritz method for the linear eigenvalue problem, the Rayleigh-Ritz method defines the Ritz values and the Ritz vectors of the QEP with respect to the projection subspace. We analyze the convergence of the method when the angle between the subspace and the desired eigenvector converges to zero. We prove that there is a Ritz value that converges to the desired eigenvalue unconditionally but the Ritz vector converges conditionally and may fail to converge. To remedy the drawback of possible non-convergence of the Ritz vector, we propose a refined Ritz vector that is mathematically different from the Ritz vector and is proved to converge unconditionally. We construct examples to illustrate our theory.

preprint2014arXiv

Singular Value Decompositions for Single-Curl Operators in Three-Dimensional Maxwell's Equations for Complex Media

This article focuses on solving the generalized eigenvalue problems (GEP) arising in the source-free Maxwell equation with magnetoelectric coupling effects that models three-dimensional complex media. The goal is to compute the smallest positive eigenvalues, and the main challenge is that the coefficient matrix in the discrete Maxwell equation is indefinite and degenerate. To overcome this difficulty, we derive a singular value decomposition (SVD) of the discrete single-curl operator and then explicitly express the basis of the invariant subspace corresponding to the nonzero eigenvalues of the GEP. Consequently, we reduce the GEP to a null space free standard eigenvalue problem (NFSEP) that contains only the nonzero (complex) eigenvalues of the GEP and can be solved by the shift-and-invert Arnoldi method without being disturbed by the null space. Furthermore, the basis of the eigendecomposition is chosen carefully so that we can apply fast Fourier transformation (FFT)-based matrix vector multiplication to solve the embedded linear systems efficiently by an iterative method. For chiral and pseudochiral complex media, which are of great interest in magnetoelectric applications, the NFSEP can be further transformed to a null space free generalized eigenvalue problem whose coefficient matrices are Hermitian and Hermitian positive definite (HHPD-NFGEP). This HHPD-NFGEP can be solved by using the invert Lanczos method without shifting. Furthermore, the embedded linear system can be solved efficiently by using the conjugate gradient method without preconditioning and the FFT-based matrix vector multiplications. Numerical results are presented to demonstrate the efficiency of the proposed methods.