BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Wen-Wei Lin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We propose a new modification of Newton iteration for finding some nonnegative Z-eigenpairs of a nonnegative tensor. The method has local quadratic convergence to a nonnegative eigenpair of a nonnegative tensor, under the usual assumption guaranteeing the local quadratic convergence of the original Newton iteration.
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.
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.
The structure-preserving doubling algorithm (SDA) is a fairly efficient method for solving problems closely related to Hamiltonian (or Hamiltonian-like) matrices, such as computing the required solutions to algebraic Riccati equations. However, for large-scale problems in $\mathbb{C}^n$ (also $\mathbb{R}^n$), the SDA with an $O(n^3)$ computational complexity does not work well. In this paper, we propose a new decoupled form of the SDA (we name it as dSDA), building on the associated Krylov subspaces thus leading to the inherent low-rank structures. Importantly, the approach decouples the original two to four iteration formulae. The resulting dSDA is much more efficient since only one quantity (instead of the original two to four) is computed iteratively. For large-scale problems, further efficiency is gained from the low-rank structures. This paper presents the theoretical aspects of the dSDA. A practical algorithm dSDA t with truncation and many illustrative numerical results will appear in a second paper.
Morphing is the process of changing one figure into another. Some numerical methods of 3D surface morphing by deformable modeling and conformal mapping are shown in this study. It is well known that there exists a unique Riemann conformal mapping from a simply connected surface into a unit disk by the Riemann mapping theorem. The dilation and relative orientations of the 3D surfaces can be linked through the Möbius transformation due to the conformal characteristic of the Riemann mapping. On the other hand, a 3D surface deformable model can be built via various approaches such as mutual parameterization from direct interpolation or surface matching using landmarks. In this paper, we take the advantage of the unique representation of 3D surfaces by the mean curvatures and the conformal factors associated with the Riemann mapping. By registering the landmarks on the conformal parametric domains, the correspondence of the mean curvatures and the conformal factors for each surfaces can be obtained. As a result, we can construct the 3D deformation field from the surface reconstruction algorithm proposed by Gu and Yau. Furthermore, by composition of the Möbius transformation and the 3D deformation field, the morphing sequence can be generated from the mean curvatures and the conformal factors on a unified mesh structure by using the cubic spline homotopy. Several numerical experiments of the face morphing are presented to demonstrate the robustness of our approach.
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.
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.
We construct a nonlinear differential equation of matrix pairs $(\mathcal{M}(t),\mathcal{L}(t))$ that is invariant (the \textbf{Structure-Preserving Property}) in the class of symplectic matrix pairs \begin{align*} \mathbb{S}_{\mathcal{S}_1,\mathcal{S}_2}=\left\{\left(\mathcal{M},\mathcal{L}\right)| \ \mathcal{M}=\left[% \begin{array}{cc} X_{12} & 0 X_{22} & I \end{array}% \right]\mathcal{S}_2, \mathcal{L}=\left[% \begin{array}{cc} I & X_{11} 0 & X_{21} \end{array}% \right]\mathcal{S}_1\right.\nonumber \left. \text{ and }X=\left[% \begin{array}{cc} X_{11} & X_{12} X_{21} & X_{22} \end{array}% \right]\text{is Hermitian}\right\} \end{align*} for certain fixed symplectic matrices $\mathcal{S}_1$ and $\mathcal{S}_2$. Its solution also preserves invariant subspaces on the whole orbit (the \textbf{Eigenvector-Preserving Property}). Such a flow is called a \textit{structure-preserving flow} and is governed by a Riccati differential equation (RDE). In addition, Radon's lemma leads to an explicit form. Therefore, blow-ups for the structure-preserving flows may happen at a finite $t$. To continue, we then utilize the Grassmann manifolds to extend the domain of the structure-preserving flow to the whole $\mathbb{R}$ subtracting some isolated points.
Numerous methods for computing conformal mesh paramterizations has been developed due to the vast applications in the field of geometry processing. Spectral conformal parameterization (SCP) is one of these methods to computing a quality conformal parameterization based on the spectral technique. SCP focus on a generalized eigenvalue problem (GEP) $L_{C}\mathbf{f} = λB\mathbf{f}$ whose eigenvector(s) associated with the smallest positive eigenvalue(s) will provide the parameterization result. This paper devotes to study a novel eigensolver for this GEP. Based on the structures of matrix pair $(L_{C},B)$, we show that this GEP can be transformed into a small-scaled compressed deating standard eigenvalue problem with a symmetric positive definite skew-Hamiltonian operator. We then propose a skew-Hamiltonian isotropic Lanczos algorithm (SHILA) to solve the reducing problem. Numerical experiments show that our compressed deating skill remove the inuence of the kernel of $L_{C}$ and transform the original problem to a more robust system. The novel SHILA method can effective avoid the disturbance of duplicate eigenvalues. As a result, our numerical eigensolver can accurately and efficiently compute the conformal parameterization based on the spectral model of SCP.