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

C. H. Jeffrey Pang

C. H. Jeffrey Pang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.OCmath.NAmath.APmath.CA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
4topics
3close 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

15 published item(s)

preprint2014arXiv

Accelerating the alternating projection algorithm for the case of affine subspaces using supporting hyperplanes

The von Neumann-Halperin method of alternating projections converges strongly to the projection of a given point onto the intersection of finitely many closed affine subspaces. We propose acceleration schemes making use of two ideas: Firstly, each projection onto an affine subspace identifies a hyperplane of codimension 1 containing the intersection, and secondly, it is easy to project onto a finite intersection of such hyperplanes. We give conditions for which our accelerations converge strongly. Finally, we perform numerical experiments to show that these accelerations perform well for a matrix model updating problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improved analysis of algorithms based on supporting halfspaces and quadratic programming for the convex intersection and feasibility problems

This paper improves the algorithms based on supporting halfspaces and quadratic programming for convex set intersection problems in our earlier paper in several directions. First, we give conditions so that much smaller quadratic programs (QPs) and approximate projections arising from partially solving the QPs are sufficient for multiple-term superlinear convergence for nonsmooth problems. Second, we identify additional regularity, which we call the second order supporting hyperplane property (SOSH), that gives multiple-term quadratic convergence. Third, we show that these fast convergence results carry over for the convex inequality problem. Fourth, we show that infeasibility can be detected in finitely many operations. Lastly, we explain how we can use the dual active set QP algorithm of Goldfarb and Idnani to get useful iterates by solving the QPs partially, overcoming the problem of solving large QPs in our algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improving the distance reduction step in the von Neumann algorithm

A known first order method to find a feasible solution to a conic problem is an adapted von Neumann algorithm. We improve the distance reduction step there by projecting onto the convex hull of previously generated points using a primal active set quadratic programming (QP) algorithm. The convergence theory is improved when the QPs are as large as possible. For problems in R^2, we analyze our algorithm by epigraphs and the monotonicity of subdifferentials. Logically, the larger the set to project onto, the better the performance per iteration, and this is indeed seen in our numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Characterizing generalized derivatives of set-valued maps: Extending the tangential and normal approaches

For a set-valued map, we characterize, in terms of its (unconvexified or convexified) graphical derivatives near the point of interest, positively homogeneous maps that are generalized derivatives in the sense of [20]. This result generalizes the Aubin criterion in [9]. A second characterization of these generalized derivatives is easier to check in practice, especially in the finite dimensional case. Finally, the third characterization in terms of limiting normal cones and coderivatives generalizes the Mordukhovich criterion in the finite dimensional case. The convexified coderivative has a bijective relationship with the set of possible generalized derivatives. We conclude by illustrating a few applications of our result.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finding saddle points of mountain pass type with quadratic models on affine spaces

The problem of computing saddle points is important in certain problems in numerical partial differential equations and computational chemistry, and is often solved numerically by a minimization problem over a set of mountain passes. We propose an algorithm to find saddle points of mountain pass type to find the bottlenecks of optimal mountain passes. The key step is to minimize the distance between level sets by using quadratic models on affine spaces similar to the strategy in the conjugate gradient algorithm. We discuss parameter choices, convergence results, and how to augment the algorithm to a path based method. Finally, we perform numerical experiments to test the convergence of our algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Set intersection problems: Supporting hyperplanes and quadratic programming

We study how the supporting hyperplanes produced by the projection process can complement the method of alternating projections and its variants for the convex set intersection problem. For the problem of finding the closest point in the intersection of closed convex sets, we propose an algorithm that, like Dykstra's algorithm, converges strongly in a Hilbert space. Moreover, this algorithm converges in finitely many iterations when the closed convex sets are cones in $\mathbb{R}^{n}$ satisfying an alignment condition. Next, we propose modifications of the alternating projection algorithm, and prove its convergence. The algorithm converges superlinearly in $\mathbb{R}^{n}$ under some nice conditions. Under a conical condition, the convergence can be finite. Lastly, we discuss the case where the intersection of the sets is empty.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Some principles for mountain pass algorithms, and the parallel distance

The problem of computing saddle points is important in certain problems in numerical partial differential equations and computational chemistry, and is often solved numerically by a minimization problem over a set of mountain passes. We point out that a good global mountain pass algorithm should have good local and global properties. Next, we define the parallel distance, and show that the square of the parallel distance has a quadratic property. We show how to design algorithms for the mountain pass problem based on perturbing parameters of the parallel distance, and that methods based on the parallel distance have midrange local and global properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Subdifferential analysis of differential inclusions via discretization

The framework of differential inclusions encompasses modern optimal control and the calculus of variations. Necessary optimality conditions in the literature identify potentially optimal paths, but do not show how to perturb paths to optimality. We first look at the corresponding discretized inclusions, estimating the subdifferential dependence of the optimal value in terms of the endpoints of the feasible paths. Our approach is to first estimate the coderivative of the reachable map. The discretized (nonsmooth) Euler-Lagrange and transversality conditions follow as a corollary. We obtain corresponding results for differential inclusions by passing discretized inclusions to the limit.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Generalized differentiation with positively homogeneous maps: Applications in set-valued analysis and metric regularity

We propose a new concept of generalized differentiation of set-valued maps that captures the first order information. This concept encompasses the standard notions of Frechet differentiability, strict differentiability, calmness and Lipschitz continuity in single-valued maps, and the Aubin property and Lipschitz continuity in set-valued maps. We present calculus rules, sharpen the relationship between the Aubin property and coderivatives, and study how metric regularity and open covering can be refined to have a directional property similar to our concept of generalized differentiation. Finally, we discuss the relationship between the robust form of generalization differentiation and its one sided counterpart.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Level set methods for finding critical points of mountain pass type

Computing mountain passes is a standard way of finding critical points. We describe a numerical method for finding critical points that is convergent in the nonsmooth case and locally superlinearly convergent in the smooth finite dimensional case. We apply these techniques to describe a strategy for the Wilkinson problem of calculating the distance of a matrix to a closest matrix with repeated eigenvalues. Finally, we relate critical points of mountain pass type to nonsmooth and metric critical point theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Level set methods for finding saddle points of general Morse index

For a real valued function, a point is critical if its derivatives are zero, and a critical point is a saddle point if it is not a local extrema. In this paper, we study algorithms to find saddle points of general Morse index. Our approach is motivated by the multidimensional mountain pass theorem, and extends our earlier work on methods (based on studying the level sets) to find saddle points of mountain pass type. We prove the convergence of our algorithms in the nonsmooth case, and the local superlinear convergence of another algorithm in the smooth finite dimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Continuity of set-valued maps revisited in the light of tame geometry

Continuity of set-valued maps is hereby revisited: after recalling some basic concepts of variational analysis and a short description of the State-of-the-Art, we obtain as by-product two Sard type results concerning local minima of scalar and vector valued functions. Our main result though, is inscribed in the framework of tame geometry, stating that a closed-valued semialgebraic set-valued map is almost everywhere continuous (in both topological and measure-theoretic sense). The result, depending on stratification techniques, holds true in a more general setting of o-minimal (or tame) set-valued maps. Some applications are briefly discussed at the end.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Lipschitz behavior of the robust regularization

To minimize or upper-bound the value of a function "robustly", we might instead minimize or upper-bound the "epsilon-robust regularization", defined as the map from a point to the maximum value of the function within an epsilon-radius. This regularization may be easy to compute: convex quadratics lead to semidefinite-representable regularizations, for example, and the spectral radius of a matrix leads to pseudospectral computations. For favorable classes of functions, we show that the robust regularization is Lipschitz around any given point, for all small epsilon > 0, even if the original function is nonlipschitz (like the spectral radius). One such favorable class consists of the semi-algebraic functions. Such functions have graphs that are finite unions of sets defined by finitely-many polynomial inequalities, and are commonly encountered in applications.

0 citations0 reviewsNo code linkedOpen access