Alternating Projection Methods
|
|
- Esther O’Connor’
- 5 years ago
- Views:
Transcription
1 Alternating Projection Methods Failure in the Absence of Convexity Matthew Tam AMSI Vacation Scholar with generous support from CSIRO Supervisor: Prof. Jon Borwein CSIRO Big Day In: 2nd & 3rd February 2012 Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
2 The Problem The setting: A Hilbert space, H. The players: r 2 sets S 1,..., S r with corresponding projections P Si. The problem: Given an initial point x 0 we seek a feasible point in r i=1 S i. Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
3 The Problem The setting: A Hilbert space, H. The players: r 2 sets S 1,..., S r with corresponding projections P Si. The problem: Given an initial point x 0 we seek a feasible point in r i=1 S i. von Neumann (1933) Suppose S 1, S 2 are closed subspaces, then x 0 H: (P S2 P S1 ) n x 0 P S1 S 2 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
4 The Problem The setting: A Hilbert space, H. The players: r 2 sets S 1,..., S r with corresponding projections P Si. The problem: Given an initial point x 0 we seek a feasible point in r i=1 S i. von Neumann (1933) Bregman (1965) Suppose S 1, S 2 are closed convex sets, then x 0 H: (P S2 P S1 ) n x 0 w. x where x S 1 S 2 Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
5 Why It Works? S 2 x 0 P S1 S 2 (x 0 ) S 1 Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
6 Why It Works? S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
7 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
8 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
9 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
10 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
11 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
12 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
13 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
14 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
15 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
16 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
17 Why It Works? P S2 P S1 (x 0 ) S 2 x 0 P S1 S 2 (x 0 ) S 1 P S1 (x 0 ) Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
18 Some Variants Douglas-Rachford (1959) Suppose S 1, S 2 are closed convex sets, then x 0 H: x n+1 := x n + R S2 R S1 (x n ) 2 where R Si (x) := 2P Si (x) x then x n w. x, a fixed point, with P S1 (x) S 1 S 2. Dysktra (1986) Suppose S 1, S 2 are closed convex sets, then x 0 H: x 1 n := x 2 n 1, x i n := P Si (x i 1 n I n 1 i ), In i := xn i (xn i 1 In 1) i with initial values x 2 0 := x 0, I i 0 := 0, then x n P S1 S 2 (x 0 ). Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
19 The Hubble Telescope 1 Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
20 The Hubble Telescope Before correction 1 Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
21 The Hubble Telescope Before correction After correction 1 Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
22 Project Aims Investigate behaviour of three alternating projection variants: von Neumann Dykstra Douglas-Rachford Particularly, cases when the underlying subsets are non-convex. 1 Partially answer the question: When does convergence fail? Develop some tools to help better understand behaviour, which are: Visual Interactive Hands-on (literally!) Even behaviour in R 2 is poorly understood. 2 Douglas Rachford Algorithm in the Absence of Convexity, Borwein, JM & Sims, B, (2011). Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
23 Cinderella in Action a 1 x + a 2 y = b x y 1 2 = 1 Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
24 A Line and 1/2-sphere We consider the case when where the sets are the 1/2-sphere and a line. Projection onto the line is simply the orthogonal projection. The 1/2-sphere is more difficult: 3 x y 1 2 = Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
25 Results Sufficient conditions for failure of von Neumann and Dysktra s methods. Behaviour of Douglas-Rachford in particular cases. Some examples of failure. Despite theoretical justification, these methods appear fairly robust. Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
26 Example: von Neumann Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
27 Example: von Neumann Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
28 Example: von Neumann Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
29 Example: Dykstra Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
30 Example: Dykstra Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
31 Example: Dykstra Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
32 Example: Dykstra Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
33 Example: Dykstra Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
34 Example: Dykstra Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
35 Example: Dykstra Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
36 Example: Douglas-Rachford Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
37 Example: Douglas-Rachford Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
38 Example: Douglas-Rachford Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
39 Where To Next? Explain algorithms, in terms of: Local convergence results. Behaviour of iterations, in the case of divergence. Behaviour for infeasible problems. (Can infeasibility be detected?) Convergence Rate. Acceleration schemes. More examples. Ultimately, a more complete theory for non-convex sets. Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
40 Questions? A big thanks the Big Day In organisers, AMSI and CSIRO; And of course to my supervisor Jon Borwein. Matthew Tam (UoN) Alternating Projection Methods CSIRO Big Day In / 14
The Method of Alternating Projections
Matthew Tam Variational Analysis Session 56th Annual AustMS Meeting 24th 27th September 2012 My Year So Far... Closed Subspaces Honours student supervised by Jon Borwein. Thesis topic: alternating projections.
More informationTHE CYCLIC DOUGLAS RACHFORD METHOD FOR INCONSISTENT FEASIBILITY PROBLEMS
THE CYCLIC DOUGLAS RACHFORD METHOD FOR INCONSISTENT FEASIBILITY PROBLEMS JONATHAN M. BORWEIN AND MATTHEW K. TAM Abstract. We analyse the behaviour of the newly introduced cyclic Douglas Rachford algorithm
More informationA Cyclic Douglas Rachford Iteration Scheme
oname manuscript o. (will be inserted by the editor) A Cyclic Douglas Rachford Iteration Scheme Jonathan M. Borwein Matthew K. Tam Received: date / Accepted: date Abstract In this paper we present two
More informationAnalysis of the convergence rate for the cyclic projection algorithm applied to semi-algebraic convex sets
Analysis of the convergence rate for the cyclic projection algorithm applied to semi-algebraic convex sets Liangjin Yao University of Newcastle June 3rd, 2013 The University of Newcastle liangjin.yao@newcastle.edu.au
More informationConvex Feasibility Problems
Laureate Prof. Jonathan Borwein with Matthew Tam http://carma.newcastle.edu.au/drmethods/paseky.html Spring School on Variational Analysis VI Paseky nad Jizerou, April 19 25, 2015 Last Revised: May 6,
More informationSparse and TV Kaczmarz solvers and the linearized Bregman method
Sparse and TV Kaczmarz solvers and the linearized Bregman method Dirk Lorenz, Frank Schöpfer, Stephan Wenger, Marcus Magnor, March, 2014 Sparse Tomo Days, DTU Motivation Split feasibility problems Sparse
More informationarxiv: v1 [math.oc] 15 Apr 2016
On the finite convergence of the Douglas Rachford algorithm for solving (not necessarily convex) feasibility problems in Euclidean spaces arxiv:1604.04657v1 [math.oc] 15 Apr 2016 Heinz H. Bauschke and
More informationThe Method of Alternating Projections
The Method of Alternating Projections by Matthew Tam c3090179 A Thesis submitted in partial fulfilment of the requirements for the degree of Bachelor of Mathematics (Honours) November, 2012. Supervisor:
More informationOn the order of the operators in the Douglas Rachford algorithm
On the order of the operators in the Douglas Rachford algorithm Heinz H. Bauschke and Walaa M. Moursi June 11, 2015 Abstract The Douglas Rachford algorithm is a popular method for finding zeros of sums
More informationA Dykstra-like algorithm for two monotone operators
A Dykstra-like algorithm for two monotone operators Heinz H. Bauschke and Patrick L. Combettes Abstract Dykstra s algorithm employs the projectors onto two closed convex sets in a Hilbert space to construct
More informationChapter 6. Metric projection methods. 6.1 Alternating projection method an introduction
Chapter 6 Metric projection methods The alternating projection method is an algorithm for computing a point in the intersection of some convex sets. The common feature of these methods is that they use
More informationVictoria Martín-Márquez
THE METHOD OF CYCLIC PROJECTIONS ON INCONSISTENT CONVEX FEASIBILITY PROBLEMS Victoria Martín-Márquez Dpto. de Análisis Matemático Universidad de Sevilla Spain V Workshop on Fixed Point Theory and Applications
More informationHybrid projection reflection method for phase retrieval
Bauschke et al. Vol. 20, No. 6/June 2003/J. Opt. Soc. Am. A 1025 Hybrid projection reflection method for phase retrieval Heinz H. Bauschke Department of Mathematics and Statistics, University of Guelph,
More informationDouglas-Rachford; a very versatile algorithm
Douglas-Rachford; a very versatile algorithm Brailey Sims CARMA, University of Newcastle, AUSTRALIA ACFPTO2018, Chiang Mai, 16 18 July 2018 https://carma.newcastle.edu.au/brailey/ Outline I Feasibility
More informationCARMA RETREAT. 19 July An overview of Pure mathematical research within CARMA. Brailey Sims
CARMA RETREAT 1 19 July 2011 An overview of Pure mathematical research within CARMA Brailey Sims Pure mathematics in CARMA is represented by: 2 Principal researchers - Laureate Professor Jonathan Borwein
More informationConvergence rate analysis for averaged fixed point iterations in the presence of Hölder regularity
Convergence rate analysis for averaged fixed point iterations in the presence of Hölder regularity Jonathan M. Borwein Guoyin Li Matthew K. Tam October 23, 205 Abstract In this paper, we establish sublinear
More information566 Y. CENSOR AND M. ZAKNOON
566 Y. CENSOR AND M. ZAKNOON Section 1.3]. In this approach, constraints of the real-world problem are represented by the demand that a solution should belong to sets C i, called constraint sets. If C
More information8 Queens, Sudoku, and Projection Methods
8 Queens, Sudoku, and Projection Methods Heinz Bauschke Mathematics, UBC Okanagan Research supported by NSERC and CRC Program heinz.bauschke@ubc.ca The Mathematical Interests of Peter Borwein The IRMACS
More informationConvergence of Fixed-Point Iterations
Convergence of Fixed-Point Iterations Instructor: Wotao Yin (UCLA Math) July 2016 1 / 30 Why study fixed-point iterations? Abstract many existing algorithms in optimization, numerical linear algebra, and
More informationOn a result of Pazy concerning the asymptotic behaviour of nonexpansive mappings
On a result of Pazy concerning the asymptotic behaviour of nonexpansive mappings arxiv:1505.04129v1 [math.oc] 15 May 2015 Heinz H. Bauschke, Graeme R. Douglas, and Walaa M. Moursi May 15, 2015 Abstract
More informationOn a result of Pazy concerning the asymptotic behaviour of nonexpansive mappings
On a result of Pazy concerning the asymptotic behaviour of nonexpansive mappings Heinz H. Bauschke, Graeme R. Douglas, and Walaa M. Moursi September 8, 2015 (final version) Abstract In 1971, Pazy presented
More information2.4 Hilbert Spaces. Outline
2.4 Hilbert Spaces Tom Lewis Spring Semester 2017 Outline Hilbert spaces L 2 ([a, b]) Orthogonality Approximations Definition A Hilbert space is an inner product space which is complete in the norm defined
More informationHeinz H. Bauschke and Walaa M. Moursi. December 1, Abstract
The magnitude of the minimal displacement vector for compositions and convex combinations of firmly nonexpansive mappings arxiv:1712.00487v1 [math.oc] 1 Dec 2017 Heinz H. Bauschke and Walaa M. Moursi December
More informationDouglas Rachford Feasibility Methods for Matrix Completion Problems
Douglas Rachford Feasibility Methods for Matrix Completion Problems Francisco J. Aragón Artacho 1 Jonathan M. Borwein 2 Matthew K. Tam 3 March 31, 2014 Abstract In this paper we give general recommendations
More informationAn Algorithm for Solving the Convex Feasibility Problem With Linear Matrix Inequality Constraints and an Implementation for Second-Order Cones
An Algorithm for Solving the Convex Feasibility Problem With Linear Matrix Inequality Constraints and an Implementation for Second-Order Cones Bryan Karlovitz July 19, 2012 West Chester University of Pennsylvania
More informationJon Borwein s Remembrance Day
Jon Borwein s Remembrance Day Institut Henri Poincaré - Paris February, 10, 2017 Speaker: Francisco Javier Aragón Artacho, Department of Mathematics University of Alicante Title: Walking guided by Dr Pi
More informationFINDING BEST APPROXIMATION PAIRS RELATIVE TO TWO CLOSED CONVEX SETS IN HILBERT SPACES
FINDING BEST APPROXIMATION PAIRS RELATIVE TO TWO CLOSED CONVEX SETS IN HILBERT SPACES Heinz H. Bauschke, Patrick L. Combettes, and D. Russell Luke Abstract We consider the problem of finding a best approximation
More informationChapter 2: Linear Programming Basics. (Bertsimas & Tsitsiklis, Chapter 1)
Chapter 2: Linear Programming Basics (Bertsimas & Tsitsiklis, Chapter 1) 33 Example of a Linear Program Remarks. minimize 2x 1 x 2 + 4x 3 subject to x 1 + x 2 + x 4 2 3x 2 x 3 = 5 x 3 + x 4 3 x 1 0 x 3
More informationSplitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches
Splitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches Patrick L. Combettes joint work with J.-C. Pesquet) Laboratoire Jacques-Louis Lions Faculté de Mathématiques
More informationProjection methods in geodesic metric spaces
Projection methods in geodesic metric spaces Computer Assisted Research Mathematics & Applications University of Newcastle http://carma.newcastle.edu.au/ 7th Asian Conference on Fixed Point Theory and
More informationThe Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment
he Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment William Glunt 1, homas L. Hayden 2 and Robert Reams 2 1 Department of Mathematics and Computer Science, Austin Peay State
More informationComputing Minmax; Dominance
Computing Minmax; Dominance CPSC 532A Lecture 5 Computing Minmax; Dominance CPSC 532A Lecture 5, Slide 1 Lecture Overview 1 Recap 2 Linear Programming 3 Computational Problems Involving Maxmin 4 Domination
More informationShock robustness of Mobile devices with HDD. Jan Ruigrok
Shock robustness of Mobile devices with HDD Jan Ruigrok 1 Content of this presentation 1. Introduction MP3-player 2. Development of new test method (pendulum) Requirements Solved problems 3. Robustness
More informationOn Nesterov s Random Coordinate Descent Algorithms - Continued
On Nesterov s Random Coordinate Descent Algorithms - Continued Zheng Xu University of Texas At Arlington February 20, 2015 1 Revisit Random Coordinate Descent The Random Coordinate Descent Upper and Lower
More informationABSTRACT CONDITIONAL EXPECTATION IN L 2
ABSTRACT CONDITIONAL EXPECTATION IN L 2 Abstract. We prove that conditional expecations exist in the L 2 case. The L 2 treatment also gives us a geometric interpretation for conditional expectation. 1.
More informationThe rate of linear convergence of the Douglas Rachford algorithm for subspaces is the cosine of the Friedrichs angle
The rate of linear convergence of the Douglas Rachford algorithm for subspaces is the cosine of the Friedrichs angle Heinz H. Bauschke, J.Y. Bello Cruz, Tran T.A. Nghia, Hung M. Phan, and Xianfu Wang April
More informationFunctional Analysis HW #5
Functional Analysis HW #5 Sangchul Lee October 29, 2015 Contents 1 Solutions........................................ 1 1 Solutions Exercise 3.4. Show that C([0, 1]) is not a Hilbert space, that is, there
More informationDistributed Kernel Regression: An Algorithm for Training Collaboratively
1 Distributed Kernel Regression: An Algorithm for Training Collaboratively J. B. Predd, Member, IEEE, S. R. Kulkarni, Fellow, IEEE, and H. V. Poor, Fellow, IEEE arxiv:cs/0601089v1 [cs.lg] 20 Jan 2006 Abstract
More informationA note on alternating projections for ill-posed semidefinite feasibility problems
Noname manuscript No. (will be inserted by the editor A note on alternating projections for ill-posed semidefinite feasibility problems Dmitriy Drusvyatskiy Guoyin Li Henry Wolkowicz Received: date / Accepted:
More informationCoordinate Update Algorithm Short Course Operator Splitting
Coordinate Update Algorithm Short Course Operator Splitting Instructor: Wotao Yin (UCLA Math) Summer 2016 1 / 25 Operator splitting pipeline 1. Formulate a problem as 0 A(x) + B(x) with monotone operators
More informationThe Direct Extension of ADMM for Multi-block Convex Minimization Problems is Not Necessarily Convergent
The Direct Extension of ADMM for Multi-block Convex Minimization Problems is Not Necessarily Convergent Yinyu Ye K. T. Li Professor of Engineering Department of Management Science and Engineering Stanford
More informationLecture 3: Semidefinite Programming
Lecture 3: Semidefinite Programming Lecture Outline Part I: Semidefinite programming, examples, canonical form, and duality Part II: Strong Duality Failure Examples Part III: Conditions for strong duality
More informationComputing Minmax; Dominance
Computing Minmax; Dominance CPSC 532A Lecture 5 Computing Minmax; Dominance CPSC 532A Lecture 5, Slide 1 Lecture Overview 1 Recap 2 Linear Programming 3 Computational Problems Involving Maxmin 4 Domination
More informationFIXED POINT ITERATION
FIXED POINT ITERATION The idea of the fixed point iteration methods is to first reformulate a equation to an equivalent fixed point problem: f (x) = 0 x = g(x) and then to use the iteration: with an initial
More informationNotes on Iterated Expectations Stephen Morris February 2002
Notes on Iterated Expectations Stephen Morris February 2002 1. Introduction Consider the following sequence of numbers. Individual 1's expectation of random variable X; individual 2's expectation of individual
More informationMath 21b: Linear Algebra Spring 2018
Math b: Linear Algebra Spring 08 Homework 8: Basis This homework is due on Wednesday, February 4, respectively on Thursday, February 5, 08. Which of the following sets are linear spaces? Check in each
More informationA NEW ITERATIVE METHOD FOR THE SPLIT COMMON FIXED POINT PROBLEM IN HILBERT SPACES. Fenghui Wang
A NEW ITERATIVE METHOD FOR THE SPLIT COMMON FIXED POINT PROBLEM IN HILBERT SPACES Fenghui Wang Department of Mathematics, Luoyang Normal University, Luoyang 470, P.R. China E-mail: wfenghui@63.com ABSTRACT.
More informationALGORITHMS AND CONVERGENCE RESULTS OF PROJECTION METHODS FOR INCONSISTENT FEASIBILITY PROBLEMS: A REVIEW
ALGORITHMS AND CONVERGENCE RESULTS OF PROJECTION METHODS FOR INCONSISTENT FEASIBILITY PROBLEMS: A REVIEW YAIR CENSOR AND MAROUN ZAKNOON Abstract. The convex feasibility problem (CFP) is to find a feasible
More informationORDERED INVOLUTIVE OPERATOR SPACES
ORDERED INVOLUTIVE OPERATOR SPACES DAVID P. BLECHER, KAY KIRKPATRICK, MATTHEW NEAL, AND WEND WERNER Abstract. This is a companion to recent papers of the authors; here we consider the selfadjoint operator
More informationNew Douglas-Rachford Algorithmic Structures and Their Convergence Analyses
New Douglas-Rachford Algorithmic Structures and Their Convergence Analyses Yair Censor and Raq Mansour Department of Mathematics University of Haifa Mt. Carmel, Haifa 3498838, Israel December 24, 2014.
More informationOptimisation and Operations Research
Optimisation and Operations Research Lecture 22: Linear Programming Revisited Matthew Roughan http://www.maths.adelaide.edu.au/matthew.roughan/ Lecture_notes/OORII/ School
More informationMassachusetts Institute of Technology 6.854J/18.415J: Advanced Algorithms Friday, March 18, 2016 Ankur Moitra. Problem Set 6
Massachusetts Institute of Technology 6.854J/18.415J: Advanced Algorithms Friday, March 18, 2016 Ankur Moitra Problem Set 6 Due: Wednesday, April 6, 2016 7 pm Dropbox Outside Stata G5 Collaboration policy:
More informationDual and primal-dual methods
ELE 538B: Large-Scale Optimization for Data Science Dual and primal-dual methods Yuxin Chen Princeton University, Spring 2018 Outline Dual proximal gradient method Primal-dual proximal gradient method
More informationStochastic Histories. Chapter Introduction
Chapter 8 Stochastic Histories 8.1 Introduction Despite the fact that classical mechanics employs deterministic dynamical laws, random dynamical processes often arise in classical physics, as well as in
More informationAnalysis of the convergence rate for the cyclic projection algorithm applied to basic semi-algebraic convex sets
Analysis of the convergence rate for the cyclic projection algorithm applied to basic semi-algebraic convex sets Jonathan M. Borwein, Guoyin Li, and Liangjin Yao Second revised version: January 04 Abstract
More informationConstants and Normal Structure in Banach Spaces
Constants and Normal Structure in Banach Spaces Satit Saejung Department of Mathematics, Khon Kaen University, Khon Kaen 4000, Thailand Franco-Thai Seminar in Pure and Applied Mathematics October 9 31,
More informationLinear and strong convergence of algorithms involving averaged nonexpansive operators
arxiv:1402.5460v1 [math.oc] 22 Feb 2014 Linear and strong convergence of algorithms involving averaged nonexpansive operators Heinz H. Bauschke, Dominikus Noll and Hung M. Phan February 21, 2014 Abstract
More informationSpooky Memory Models Jerome Busemeyer Indiana Univer Univ sity er
Spooky Memory Models Jerome Busemeyer Indiana University Associative Memory star t = planet c=universe space Mars Notice: No return from space or mars to planet Failure of Spreading Activation Model Why
More informationConvergence Analysis of Processes with Valiant Projection Operators in. Hilbert Space. Noname manuscript No. (will be inserted by the editor)
Noname manuscript No. will be inserted by the editor) Convergence Analysis of Processes with Valiant Projection Operators in Hilbert Space Yair Censor 1 and Raq Mansour 2 December 4, 2016. Revised: May
More informationMarkov decision processes and interval Markov chains: exploiting the connection
Markov decision processes and interval Markov chains: exploiting the connection Mingmei Teo Supervisors: Prof. Nigel Bean, Dr Joshua Ross University of Adelaide July 10, 2013 Intervals and interval arithmetic
More informationProjections, Radon-Nikodym, and conditioning
Projections, Radon-Nikodym, and conditioning 1 Projections in L 2......................... 1 2 Radon-Nikodym theorem.................... 3 3 Conditioning........................... 5 3.1 Properties of
More informationBelief Propagation, Information Projections, and Dykstra s Algorithm
Belief Propagation, Information Projections, and Dykstra s Algorithm John MacLaren Walsh, PhD Department of Electrical and Computer Engineering Drexel University Philadelphia, PA jwalsh@ece.drexel.edu
More informationMATHEMATICS. Course Syllabus. Section A: Linear Algebra. Subject Code: MA. Course Structure. Ordinary Differential Equations
MATHEMATICS Subject Code: MA Course Structure Sections/Units Section A Section B Section C Linear Algebra Complex Analysis Real Analysis Topics Section D Section E Section F Section G Section H Section
More informationVictoria Martín-Márquez
A NEW APPROACH FOR THE CONVEX FEASIBILITY PROBLEM VIA MONOTROPIC PROGRAMMING Victoria Martín-Márquez Dep. of Mathematical Analysis University of Seville Spain XIII Encuentro Red de Análisis Funcional y
More informationNear-Potential Games: Geometry and Dynamics
Near-Potential Games: Geometry and Dynamics Ozan Candogan, Asuman Ozdaglar and Pablo A. Parrilo September 6, 2011 Abstract Potential games are a special class of games for which many adaptive user dynamics
More informationShort proofs of the Quantum Substate Theorem
Short proofs of the Quantum Substate Theorem Rahul Jain CQT and NUS, Singapore Ashwin Nayak University of Waterloo Classic problem Given a biased coin c Pr(c = H) = 1 - e Pr(c = T) = e Can we generate
More informationLecture 5 : Projections
Lecture 5 : Projections EE227C. Lecturer: Professor Martin Wainwright. Scribe: Alvin Wan Up until now, we have seen convergence rates of unconstrained gradient descent. Now, we consider a constrained minimization
More informationConductivity imaging from one interior measurement
Conductivity imaging from one interior measurement Amir Moradifam (University of Toronto) Fields Institute, July 24, 2012 Amir Moradifam (University of Toronto) A convergent algorithm for CDII 1 / 16 A
More informationAbout Split Proximal Algorithms for the Q-Lasso
Thai Journal of Mathematics Volume 5 (207) Number : 7 http://thaijmath.in.cmu.ac.th ISSN 686-0209 About Split Proximal Algorithms for the Q-Lasso Abdellatif Moudafi Aix Marseille Université, CNRS-L.S.I.S
More informationA Unified Approach to Proximal Algorithms using Bregman Distance
A Unified Approach to Proximal Algorithms using Bregman Distance Yi Zhou a,, Yingbin Liang a, Lixin Shen b a Department of Electrical Engineering and Computer Science, Syracuse University b Department
More informationThe PANORAMA. Numerical Ranges in. Modern Times. Man-Duen Choi 蔡文端. For the 12 th WONRA.
The PANORAMA of Numerical Ranges in Modern Times Man-Duen Choi 蔡文端 choi@math.toronto.edu For the 12 th WONRA In the last century, Toeplitz-Hausdorff Theorem, Ando-Arveson Dilation Theorem, In the 21 st
More informationLinear programs Optimization Geoff Gordon Ryan Tibshirani
Linear programs 10-725 Optimization Geoff Gordon Ryan Tibshirani Review: LPs LPs: m constraints, n vars A: R m n b: R m c: R n x: R n ineq form [min or max] c T x s.t. Ax b m n std form [min or max] c
More informationA strongly polynomial algorithm for linear systems having a binary solution
A strongly polynomial algorithm for linear systems having a binary solution Sergei Chubanov Institute of Information Systems at the University of Siegen, Germany e-mail: sergei.chubanov@uni-siegen.de 7th
More informationDistributed Optimization via Alternating Direction Method of Multipliers
Distributed Optimization via Alternating Direction Method of Multipliers Stephen Boyd, Neal Parikh, Eric Chu, Borja Peleato Stanford University ITMANET, Stanford, January 2011 Outline precursors dual decomposition
More informationZero sum games Proving the vn theorem. Zero sum games. Roberto Lucchetti. Politecnico di Milano
Politecnico di Milano General form Definition A two player zero sum game in strategic form is the triplet (X, Y, f : X Y R) f (x, y) is what Pl1 gets from Pl2, when they play x, y respectively. Thus g
More informationMAT 578 FUNCTIONAL ANALYSIS EXERCISES
MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.
More informationREAL RENORMINGS ON COMPLEX BANACH SPACES
REAL RENORMINGS ON COMPLEX BANACH SPACES F. J. GARCÍA PACHECO AND A. MIRALLES Abstract. In this paper we provide two ways of obtaining real Banach spaces that cannot come from complex spaces. In concrete
More informationA GENERAL FRAMEWORK FOR A CLASS OF FIRST ORDER PRIMAL-DUAL ALGORITHMS FOR TV MINIMIZATION
A GENERAL FRAMEWORK FOR A CLASS OF FIRST ORDER PRIMAL-DUAL ALGORITHMS FOR TV MINIMIZATION ERNIE ESSER XIAOQUN ZHANG TONY CHAN Abstract. We generalize the primal-dual hybrid gradient (PDHG) algorithm proposed
More informationINERTIAL ACCELERATED ALGORITHMS FOR SOLVING SPLIT FEASIBILITY PROBLEMS. Yazheng Dang. Jie Sun. Honglei Xu
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX INERTIAL ACCELERATED ALGORITHMS FOR SOLVING SPLIT FEASIBILITY PROBLEMS Yazheng Dang School of Management
More informationFINDING BEST APPROXIMATION PAIRS RELATIVE TO A CONVEX AND A PROX-REGULAR SET IN A HILBERT SPACE
FINDING BEST APPROXIMATION PAIRS RELATIVE TO A CONVEX AND A PROX-REGULAR SET IN A HILBERT SPACE D. RUSSELL LUKE Abstract. We study the convergence of an iterative projection/reflection algorithm originally
More informationProjection and proximal point methods: convergence results and counterexamples
Projection and proximal point methods: convergence results and counterexamples Heinz H. Bauschke, Eva Matoušková, and Simeon Reich June 14, 2003 Version 1.15 Abstract Recently, Hundal has constructed a
More informationExperimental Computation and Visual Theorems
Experimental Computation and Visual Theorems Jonathan M. Borwein 1 CARMA, University of Newcastle, NSW jon.borwein@gmail.com, www.carma.newcastle.edu.au Abstract. Long before current graphic, visualisation
More informationL 2 -cohomology of hyperplane complements
(work with Tadeusz Januskiewicz and Ian Leary) Oxford, Ohio March 17, 2007 1 Introduction 2 The regular representation L 2 -(co)homology Idea of the proof 3 Open covers Proof of the Main Theorem Statement
More informationInverse Eigenvalue Problems in Wireless Communications
Inverse Eigenvalue Problems in Wireless Communications Inderjit S. Dhillon Robert W. Heath Jr. Mátyás Sustik Joel A. Tropp The University of Texas at Austin Thomas Strohmer The University of California
More informationarxiv: v4 [math.oc] 28 Aug 2017
GENERALIZED SINKHORN ITERATIONS FOR REGULARIZING INVERSE PROBLEMS USING OPTIMAL MASS TRANSPORT JOHAN KARLSSON AND AXEL RINGH arxiv:1612.02273v4 [math.oc] 28 Aug 2017 Abstract. The optimal mass transport
More informationVariational inequalities for fixed point problems of quasi-nonexpansive operators 1. Rafał Zalas 2
University of Zielona Góra Faculty of Mathematics, Computer Science and Econometrics Summary of the Ph.D. thesis Variational inequalities for fixed point problems of quasi-nonexpansive operators 1 by Rafał
More informationSummary of the simplex method
MVE165/MMG631,Linear and integer optimization with applications The simplex method: degeneracy; unbounded solutions; starting solutions; infeasibility; alternative optimal solutions Ann-Brith Strömberg
More informationMATH 650. THE RADON-NIKODYM THEOREM
MATH 650. THE RADON-NIKODYM THEOREM This note presents two important theorems in Measure Theory, the Lebesgue Decomposition and Radon-Nikodym Theorem. They are not treated in the textbook. 1. Closed subspaces
More informationMATH 12 CLASS 5 NOTES, SEP
MATH 12 CLASS 5 NOTES, SEP 30 2011 Contents 1. Vector-valued functions 1 2. Differentiating and integrating vector-valued functions 3 3. Velocity and Acceleration 4 Over the past two weeks we have developed
More informationKAKUTANI S FIXED POINT THEOREM AND THE MINIMAX THEOREM IN GAME THEORY
KAKUTANI S FIXED POINT THEOREM AND THE MINIMAX THEOREM IN GAME THEORY YOUNGGEUN YOO Abstract. The imax theorem is one of the most important results in game theory. It was first introduced by John von Neumann
More informationInformation Projection Algorithms and Belief Propagation
χ 1 π(χ 1 ) Information Projection Algorithms and Belief Propagation Phil Regalia Department of Electrical Engineering and Computer Science Catholic University of America Washington, DC 20064 with J. M.
More informationTelescopes. Astronomy 320 Wednesday, February 14, 2018
Telescopes Astronomy 320 Wednesday, February 14, 2018 Telescopes gather light and resolve detail A telescope is sometimes called a light bucket. Number of photons collected per second is proportional to
More informationLecture December 2009 Fall 2009 Scribe: R. Ring In this lecture we will talk about
0368.4170: Cryptography and Game Theory Ran Canetti and Alon Rosen Lecture 7 02 December 2009 Fall 2009 Scribe: R. Ring In this lecture we will talk about Two-Player zero-sum games (min-max theorem) Mixed
More informationDerivatives. Differentiability problems in Banach spaces. Existence of derivatives. Sharpness of Lebesgue s result
Differentiability problems in Banach spaces David Preiss 1 Expanded notes of a talk based on a nearly finished research monograph Fréchet differentiability of Lipschitz functions and porous sets in Banach
More informationHomework Assignment #5 Due Wednesday, March 3rd.
Homework Assignment #5 Due Wednesday, March 3rd. 1. In this problem, X will be a separable Banach space. Let B be the closed unit ball in X. We want to work out a solution to E 2.5.3 in the text. Work
More informationConvex Analysis and Monotone Operator Theory in Hilbert Spaces
Heinz H. Bauschke Patrick L. Combettes Convex Analysis and Monotone Operator Theory in Hilbert Spaces Springer Foreword This self-contained book offers a modern unifying presentation of three basic areas
More informationAssignment #9: Orthogonal Projections, Gram-Schmidt, and Least Squares. Name:
Assignment 9: Orthogonal Projections, Gram-Schmidt, and Least Squares Due date: Friday, April 0, 08 (:pm) Name: Section Number Assignment 9: Orthogonal Projections, Gram-Schmidt, and Least Squares Due
More informationPart II Probability and Measure
Part II Probability and Measure Theorems Based on lectures by J. Miller Notes taken by Dexter Chua Michaelmas 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly)
More informationProf. M. Saha Professor of Mathematics The University of Burdwan West Bengal, India
CHAPTER 9 BY Prof. M. Saha Professor of Mathematics The University of Burdwan West Bengal, India E-mail : mantusaha.bu@gmail.com Introduction and Objectives In the preceding chapters, we discussed normed
More informationNonconvex notions of regularity and convergence of fundamental algorithms for feasibility problems
Nonconvex notions of regularity and convergence of fundamental algorithms for feasibility problems Robert Hesse and D. Russell Luke December 12, 2012 Abstract We consider projection algorithms for solving
More information