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

Long Lee

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

ResearcherUnclaimed source record
math.NAnlin.SImath-phmath.MPmath.SPNumerical Analysis

Catalog footprint

What is connected

11works
6topics
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

11 published item(s)

preprint2022arXiv

A FEAST variant incorporated with a power iteration

We present a variant of the FEAST matrix eigensolver for solving restricted real and symmetric eigenvalue problems. The method is derived from a combination of a variant of the FEAST method, which employs two contour integrals per iteration, and a power subspace iteration process. Compared with the original FEAST method, our new method does not require that the search subspace dimension must be greater than or equal to the number of eigenvalues inside a search interval, and can deal with narrow search intervals more effectively. Empirically, the FEAST iteration and the power subspace iteration are in a mutually beneficial collaboration to make the new method stable and robust.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Spectral Projection Preconditioner for Solving Ill Conditioned Linear Systems

We present a preconditioner based on spectral projection that is combined with a deflated Krylov subspace method for solving ill conditioned linear systems of equations. Our results show that the proposed algorithm requires many fewer iterations to achieve the convergence criterion for solving an ill conditioned problem than a Krylov subspace solver. In our numerical experiments, the solution obtained by the proposed algorithm is more accurate in terms of the norm of the distance to the exact solution of the linear system of equations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A class of fast geodesic shooting algorithms for template matching and its applications via the $N$-particle system of the Euler-Poincaré equations

The Euler-Poincaré (EP) equations describe the geodesic motion on the diffeomorphism group. For template matching (template deformation), the Euler-Lagrangian equation, arising from minimizing an energy function, falls into the Euler-Poincaré theory and can be recast into the EP equations. By casting the EP equations in the Lagrangian (or characteristics) form, we formulate the equations as a finite dimensional particle system. The evolution of this particle system describes the geodesic motion of landmark points on a Riemann manifold. In this paper we present a class of novel algorithms that take advantage of the structure of the particle system to achieve a fast matching process between the reference and the target templates. The strong suit of the proposed algorithms includes (1) the efficient feedback control iteration, which allows one to find the initial velocity field for driving the deformation from the reference template to the target one, (2) the use of the conical kernel in the particle system, which limits the interaction between particles and thus accelerates the convergence, and (3) the availability of the implementation of fast-multipole method for solving the particle system, which could reduce the computational cost from $O(N^2)$ to $O(N\log N)$, where $N$ is the number of particles. The convergence properties of the proposed algorithms are analyzed. Finally, we present several examples for both exact and inexact matchings, and numerically analyze the iterative process to illustrate the efficiency and the robustness of the proposed algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solitary waves and $N$-particle algorithms for a class of Euler-Poincaré equations

We study a class of partial differential equations (PDEs) in the family of the so-called Euler-Poincaré differential systems, with the aim of developing a foundation for numerical algorithms of their solutions. This requires particular attention to the mathematical properties of this system when the associated class of elliptic operators possesses non-smooth kernels. By casting the system in its Lagrangian (or characteristics) form, we first formulate a particles system algorithm in free space with homogeneous Dirichlet boundary conditions for the evolving fields. We next examine the deformation of the system when non-homogeneous "constant stream" boundary conditions are assumed. We show how this simple change at the boundary deeply affects the nature of the evolution, from hyperbolic-like to dispersive with a non-trivial dispersion relation, and examine the potentially regularizing properties of singular kernels offered by this deformation. From the particle algorithm viewpoint, kernel singularities affect the existence and uniqueness of solutions to the corresponding ordinary differential equations systems. We illustrate this with the case when the operator kernel assumes a conical shape over the spatial variables, and examine in detail two-particle dynamics under the resulting lack of Lipschitz-continuity. Curiously, we find that for the conically-shaped kernels the motion of the related two-dimensional waves can become completely integrable under appropriate initial data. This reduction projects the two-dimensional system to the one-dimensional completely integrable Shallow-Water equation [Camassa, R. and Holm, D. D., Phys. Rev. Lett., 71, 1961-1964, 1993], while retaining the full dependence on two spatial dimensions for the single channel solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A connection between the shallow-water equations and the Euler-Poincaré equations

The Euler-Poincaré differential (EPDiff) equations and the shallow water (SW) equations share similar wave characteristics. Using the Hamiltonian structure of the SW equations with flat bottom topography, we establish a connection between the EPDiff equations and the SW equations in one and multi-dimensions. Additionally, we show that the EPDiff equations can be recast in a curl formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A conservation formulation and a numerical algorithm for the double-gyre nonlinear shallow-water model

We present a conservation formulation and a numerical algorithm for the reduced-gravity shallow-water equations on a beta plane, subjected to a constant wind forcing that leads to the formation of double-gyre circulation in a closed ocean basin. The novelty of the paper is that we reformulate the governing equations into a nonlinear hyperbolic conservation law plus source terms. A second-order fractional-step algorithm is used to solve the reformulated equations. In the first step of the fractional-step algorithm, we solve the homogeneous hyperbolic shallow-water equations by the wave-propagation finite volume method. The resulting intermediate solution is then used as the initial condition for the initial-boundary value problem in the second step. As a result, the proposed method is not sensitive to the choice of viscosity and gives high-resolution results for coarse grids, as long as the Rossby deformation radius is resolved. We discuss the boundary conditions in each step, when no-slip boundary conditions are imposed to the problem. We validate the algorithm by a periodic flow on an f-plane with exact solutions. The order-of-accuracy for the proposed algorithm is tested numerically. We illustrate a quasi-steady-state solution of the double-gyre model via the height anomaly and the contour of stream function for the formation of double-gyre circulation in a closed basin. Our calculations are highly consistent with the results reported in the literature. Finally, we present an application, in which the double-gyre model is coupled with the advection equation for modeling transport of a pollutant in a closed ocean basin.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis and development of compact finite difference schemes with optimized numerical dispersion relations

Finite difference approximation, in addition to Taylor truncation errors, introduces numerical dispersion-and-dissipation errors into numerical solutions of partial differential equations. We analyze a class of finite difference schemes which are designed to minimize these errors (at the expense of formal order of accuracy), and we analyze the interplay between the Taylor truncation errors and the dispersion-and-dissipation errors during mesh refinement. In particular, we study the numerical dispersion relation of the fully discretized non-dispersive transport equation in one and two space dimensions. We derive the numerical phase error and the $L^2$-norm error of the solution in terms of the dispersion-and-dissipation error. Based on our analysis, we investigate the error dynamics among various optimized compact schemes and the unoptimized higher-order generalized Padé compact schemes. The dynamics shed light on the principles of designing suitable optimized compact schemes for a given problem. Using these principles as guidelines, we then propose an optimized scheme that prescribes the numerical dispersion relation before finding the corresponding discretization. This approach produces smaller numerical dispersion-and-dissipation errors for linear and nonlinear problems, compared with the unoptimized higher-order compact schemes and other optimized schemes developed in the literature. Finally, we discuss the difficulty of developing an optimized composite boundary scheme for problems with non-trivial boundary conditions. We propose a composite scheme that introduces a buffer zone to connect an optimized interior scheme and an unoptimized boundary scheme. Our numerical experiments show that this strategy produces small $L^2$-norm error when a wave packet passes through the non-periodic boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The WKB approximation of semiclassical eigenvalues of the Zakharov-Shabat problem

We numerically compute eigenvalues of the non-self-adjoint Zakharov--Shabat problem in the semiclassical regime. In particular, we compute the eigenvalues for a Gaussian potential and compare the results to the corresponding (formal) WKB approximations used in the approach to the semiclassical or zero-dispersion limit of the focusing nonlinear Schroedinger equation via semiclassical soliton ensembles. This numerical experiment, taken together with recent numerical experiments [17,18], speaks directly to the viability of this approach; in particular, our experiment suggests a value for the rate of convergence of the WKB eigenvalues to the true eigenvalues in the semiclassical limit. This information provides some hint as to how these approximations might be rigorously incorporated into the asymptotic analysis of the singular limit for the associated nonlinear partial differential equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A second look at the Gaussian semiclassical soliton ensemble for the focusing nonlinear Schrödinger equation

We present the results of a numerical experiment inspired by the semiclassical (zero-dispersion) limit of the focusing nonlinear Schroedinger (NLS) equation. In particular, we focus on the Gaussian semiclassical soliton ensemble, a family of exact multisoliton solutions obtained by repeatedly solving the initial-value problem for a particular sequence of initial data. The sequence of data is generated by adding an asymptotically vanishing sequence of perturbations to pure Gaussian initial data. These perturbations are obtained by applying the inverse-scattering transform to formal WKB approximations of eigenvalues of the associated spectral problem with a Gaussian potential. Recent results [Lee, Lyng, & Vankova, Physica D 24 (2012):1767--1781] suggest that, remarkably, these perturbations---interlaced as they are with the integrable structure of the equation---do not excite the acute modulational instabilities that are known to be present in the semiclassical regime. Here, we provide additional evidence to support the claim that these WKB-induced perturbations indeed have a very special structure. In particular, as a control experiment, we examine the evolution from a family of initial data created by an asymptotically vanishing family of analytic perturbations which are qualitatively indistinguishable from the WKB-induced perturbations that generate the Gaussian semiclassical soliton ensemble. We then compare this evolution to the (numerically computed) true evolution of the Gaussian and also to the evolution of the corresponding members of the semiclassical soliton ensemble. Our results both highlight the exceptional nature of the WKB-induced perturbations used to generate the semiclassical soliton ensemble and provide new insight into the sensitivity properties of the semiclassical limit problem for the focusing NLS equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Gaussian semiclassical soliton ensemble and numerical methods for the focusing nonlinear Schrödinger equation

We report on a number of careful numerical experiments motivated by the semiclassical (zero-dispersion, ε\downarrow 0) limit of the focusing nonlinear Schrödinger equation. Our experiments are designed to study the evolution of a particular family of perturbations of the initial data. These asymptotically small perturbations are precisely those that result from modifying the initial-data by using formal approximations to the spectrum of the associated spectral problem; such modified data has always been a standard part of the analysis of zero-dispersion limits of integrable systems. However, in the context of the focusing nonlinear Schrödinger equation, the ellipticity of the Whitham equations casts some doubt on the validity of this procedure. To carry out our experiments, we introduce an implicit finite difference scheme for the partial differential equation, and we validate both the proposed scheme and the standard split-step scheme against a numerical implementation of the inverse scattering transform for a special case in which the scattering data is known exactly. As part of this validation, we also investigate the use of the Krasny filter which is sometimes suggested as appropriate for nearly ill-posed problems such as we consider here. Our experiments show that that the O(ε) rate of convergence of the modified data to the true data is propagated to positive times including times after wave breaking.

0 citations0 reviewsNo code linkedOpen access