The Method of Alternating Projections
|
|
- Brice Rich
- 5 years ago
- Views:
Transcription
1 Matthew Tam Variational Analysis Session 56th Annual AustMS Meeting 24th 27th September 2012
2 My Year So Far... Closed Subspaces Honours student supervised by Jon Borwein. Thesis topic: alternating projections. Over the past year I ve learnt about: Classical alternating projection results. Difficulties of nonconvex alternating projections (AMSI Vacation) including development of an interactive Cinderella interface. Alternating Bregman projection in Banach Spaces. Computational experiments including nonconvex instances.
3 My Year So Far... Closed Subspaces Honours student supervised by Jon Borwein. Thesis topic: alternating projections. Over the past year I ve learnt about: Classical alternating projection results. Difficulties of nonconvex alternating projections (AMSI Vacation) including development of an interactive Cinderella interface. Alternating Bregman projection in Banach Spaces. Computational experiments including nonconvex instances. Now routinely applied to nonconvex problems (convergence not guaranteed) including hard combinatorial problems with good results. Eg. Diophantine equations, protein folding, sphere-packing, 3SAT, Sudoku, image reconstruction
4 My Year So Far... Closed Subspaces Honours student supervised by Jon Borwein. Thesis topic: alternating projections. Over the past year I ve learnt about: Classical alternating projection results. Difficulties of nonconvex alternating projections (AMSI Vacation) including development of an interactive Cinderella interface. Alternating Bregman projection in Banach Spaces. Computational experiments including nonconvex instances. Now routinely applied to nonconvex problems (convergence not guaranteed) including hard combinatorial problems with good results. Eg. Diophantine equations, protein folding, sphere-packing, 3SAT, Sudoku, image reconstruction We are designing large scale experiments to understand this better.
5 Introduction Closed Subspaces Let H be a Hilbert space. The (metric) projection of x H onto the set M is a point p M such that p x m x for all m M We write P M (x) = p when p exists uniquely.
6 Introduction Closed Subspaces Let H be a Hilbert space. The (metric) projection of x H onto the set M is a point p M such that p x m x for all m M We write P M (x) = p when p exists uniquely. Given sets M, N such that M N can we: Compute P M N (x) given x H? (Best approximation) Find a point x M N? (Feasibility)
7 Introduction Closed Subspaces Let H be a Hilbert space. The (metric) projection of x H onto the set M is a point p M such that p x m x for all m M We write P M (x) = p when p exists uniquely. Given sets M, N such that M N can we: Compute P M N (x) given x H? (Best approximation) Find a point x M N? (Feasibility) We address the question: Can these problems be solved knowing only P M and P N?
8 Two Closed Subspaces Let M, N be closed subspaces. Then: Fact P M, P N commute if and only if their composition is equal to P M N. If the projections commute then their composition gives P M N.
9 Two Closed Subspaces Let M, N be closed subspaces. Then: Fact P M, P N commute if and only if their composition is equal to P M N. If the projections commute then their composition gives P M N. Otherwise, try projecting alternatively: What happens in the limit? P x M P 0 x N P 1 M P x2 x N P 3 M P x4 x N 5...
10 Why It Works? Closed Subspaces N x 0 P M N (x 0 ) M If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
11 Why It Works? Closed Subspaces N x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
12 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
13 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
14 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
15 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
16 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
17 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
18 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
19 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
20 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
21 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
22 Why It Works? Closed Subspaces N P N P M (x 0 ) x 0 P M N (x 0 ) M P M (x 0 ) If H is the hyperplane given by a, x = b then. P H (x) = x a, x a 2 a
23 von Neumann s Alternating Projections Theorem (von Neumann, 1933) Let M, N H be closed subspaces then x H: (P M P N ) n (x) P M N (x) Proof. To show that (x n ) is Cauchy: P N (... P M P N P M P N ) = (P N P N... P N P M P N ) }{{}}{{} k terms (k+1) terms or (P N P M... P N P M P N ) }{{} (k+1) terms
24 von Neumann s Alternating Projections Theorem (von Neumann, 1933) Let M, N H be closed subspaces then x H: (P M P N ) n (x) P M N (x) Proof. To show that (x n ) is Cauchy: P N (... P M P N P M P N ) }{{} k terms = (P N... P N P M P N ) }{{} k terms or (P N P M... P N P M P N ) }{{} (k+1) terms
25 Halperin s Extension Theorem (Halperin, 1962) Let S 1, S 2,..., S r H be closed subspaces then x H: (P Sr... P S2 P S1 ) n (x) (P r i=1 S i )(x) Proof. If T linear, nonexpansive then H = ker(i T ) cl(range(i T )).
26 Halperin s Extension Theorem (Halperin, 1962) Let S 1, S 2,..., S r H be closed subspaces then x H: (P Sr... P S2 P S1 ) n (x) (P r i=1 S i )(x) Proof. If T linear, nonexpansive then H = ker(i T ) cl(range(i T )). If T linear, idempotent, nonexpansive then T = P ker(i T ).
27 Halperin s Extension Theorem (Halperin, 1962) Let S 1, S 2,..., S r H be closed subspaces then x H: T n (x) = (P Sr... P S2 P S1 ) n (x) (P r i=1 S i )(x) Proof. If T linear, nonexpansive then H = ker(i T ) cl(range(i T )). If T linear, idempotent, nonexpansive then T = P ker(i T ). T n x T n+1 x 0 hence T n (I T )x 0,
28 Halperin s Extension Theorem (Halperin, 1962) Let S 1, S 2,..., S r H be closed subspaces then x H: T n (x) = (P Sr... P S2 P S1 ) n (x) (P r i=1 S i )(x) Proof. If T linear, nonexpansive then H = ker(i T ) cl(range(i T )). If T linear, idempotent, nonexpansive then T = P ker(i T ). T n x T n+1 x 0 hence T n y 0, y range(i T )
29 Halperin s Extension Theorem (Halperin, 1962) Let S 1, S 2,..., S r H be closed subspaces then x H: T n (x) = (P Sr... P S2 P S1 ) n (x) (P r i=1 S i )(x) Proof. If T linear, nonexpansive then H = ker(i T ) cl(range(i T )). If T linear, idempotent, nonexpansive then T = P ker(i T ). T n x T n+1 x 0 hence T n y 0, y cl range(i T )
30 Best Approximation for Half-Spaces? x 0 N M Beyond hyperplanes, half-spaces, spheres, balls,... projections are difficult to compute. Even for the ellipse in R 2, given by x 2 /a 2 + y 2 /b 2 = 1, we have: ( a 2 u P E (u, v) = a 2 t, b 2 ) v a 2 u 2 b 2 where t a 2 t + b2 v 2 b 2 t = 1
31 Best Approximation for Half-Spaces? x 0 N P N (x 0 ) M Beyond hyperplanes, half-spaces, spheres, balls,... projections are difficult to compute. Even for the ellipse in R 2, given by x 2 /a 2 + y 2 /b 2 = 1, we have: ( a 2 u P E (u, v) = a 2 t, b 2 ) v a 2 u 2 b 2 where t a 2 t + b2 v 2 b 2 t = 1
32 Best Approximation for Half-Spaces? x 0 N P N (x 0 ) M P M P N (x 0 ) Beyond hyperplanes, half-spaces, spheres, balls,... projections are difficult to compute. Even for the ellipse in R 2, given by x 2 /a 2 + y 2 /b 2 = 1, we have: ( a 2 u P E (u, v) = a 2 t, b 2 ) v a 2 u 2 b 2 where t a 2 t + b2 v 2 b 2 t = 1
33 Best Approximation for Half-Spaces? x 0 N P N (x 0 ) M P N P M P N (x 0 ) Beyond hyperplanes, half-spaces, spheres, balls,... projections are difficult to compute. Even for the ellipse in R 2, given by x 2 /a 2 + y 2 /b 2 = 1, we have: ( a 2 u P E (u, v) = a 2 t, b 2 ) v a 2 u 2 b 2 where t a 2 t + b2 v 2 b 2 t = 1
34 Best Approximation for Half-Spaces? x 0 N P N (x 0 ) M P N P M P M P N (x 0 ) Beyond hyperplanes, half-spaces, spheres, balls,... projections are difficult to compute. Even for the ellipse in R 2, given by x 2 /a 2 + y 2 /b 2 = 1, we have: ( a 2 u P E (u, v) = a 2 t, b 2 ) v a 2 u 2 b 2 where t a 2 t + b2 v 2 b 2 t = 1
35 Best Approximation for Half-Spaces? x 0 N P N (x 0 ) M P M N (x 0 ) P N P M P M P N (x 0 ) Beyond hyperplanes, half-spaces, spheres, balls,... projections are difficult to compute. Even for the ellipse in R 2, given by x 2 /a 2 + y 2 /b 2 = 1, we have: ( a 2 u P E (u, v) = a 2 t, b 2 ) v a 2 u 2 b 2 where t a 2 t + b2 v 2 b 2 t = 1
36 Bregman s Generalisation Theorem (Bregman, 1965) Let C 1, C 2,..., C r H be closed convex sets then x H: (P Cr... P C2 P C1 ) n (x) w. x r i=1 C i Proof. Use weak compactness to extract a weakly convergence subsequence.
37 Bregman s Generalisation Theorem (Bregman, 1965) Let C 1, C 2,..., C r H be closed convex sets then x H: (P Cr... P C2 P C1 ) n (x) w. x r i=1 C i Proof. Use weak compactness to extract a weakly convergence subsequence. Can MAP fail to converge in norm?
38 The Hundal Example (Revisited) Failure of Norm Convergence (Hundal, 2004) Let H = l 2 and {e i } an orthonormal basis. Take x 0 = e 3 and C 1 = ker(e 1 ) and C 2 = cl conv k=2 epi C 0,k then MAP fails to converge is norm. Note: C 1 is a hyperplane and C 2 a closed convex cone. h k (t) = exp( (t + kπ/2) 3 ) C 1 C 2 = {0}
39 The Hundal Example (Revisited) Failure of Norm Convergence (Hundal, 2004)(Matoušková & Reich, 2003) Let H = l 2 and {e i } an orthonormal basis. Take x 0 = e k0 and C 1 = ker(e 1 ) and C 2 = cl conv k=k 0 epi C 0,k then MAP fails to converge is norm. Note: C 1 is a hyperplane and C 2 a closed convex cone. h k (t) = exp( (t + kπ/2) 3 ) C 1 C 2 = {0}
40 The Hundal Example (Revisited) Failure of Norm Convergence (Hundal, 2004)(Matoušková & Reich, 2003) Let H = l 2 and {e i } an orthonormal basis. Take x 0 = e k0 and C 1 = ker(e 1 ) and C 2 = cl conv k=k 0 epi C 0,k then MAP fails to converge is norm. Note: C 1 is a hyperplane and C 2 a closed convex cone. h k (t) = exp( (t + kπ/2) 3 ) C 1 C 2 = {0}
41 The Hundal Example (Revisited) Failure of Norm Convergence (Hundal, 2004)(Matoušková & Reich, 2003) Let H = l 2 and {e i } an orthonormal basis. Take x 0 = e k0 and C 1 = ker(e 1 ) and C 2 = cl conv k=k 0 epi C 0,k then MAP fails to converge is norm. Note: C 1 is a hyperplane and C 2 a closed convex cone. h k (t) = exp( (t + kπ/2) 3 ) Final step: C 1 C 2 = {0} (P C2 P C1 ) Nm e k0 e m < 1/7 = (P C2 P C1 ) Nm e k0 > 6/7
42 The Hundal Example (Revisited) Can it fail to converge in norm on a real problem? Conjecture (Borwein & Bauschke, 1993) If C 1 is closed and affine with finite codimension, C 2 is the nonnegative cone in l 2 (N) then MAP is norm convergent. True when C 1 is a hyperplane (unlike Hundal). This captures most concrete applications.
43 Averaged Projections (The Crucial Product Space Trick) Given C 1, C 2,..., C r H consider H r = H H. Define: C = {(x 1, x 2,..., x r ) : x i C i }, D = {(x 1, x 2,..., x r ) : x 1 = x i }
44 Averaged Projections (The Crucial Product Space Trick) Given C 1, C 2,..., C r H consider H r = H H. Define: C = {(x 1, x 2,..., x r ) : x i C i }, D = {(x 1, x 2,..., x r ) : x 1 = x i } It is easily verified that: (P C x) i = P Ci x i, (P D x) i = 1 r r j=1 x j
45 Averaged Projections (The Crucial Product Space Trick) Given C 1, C 2,..., C r H consider H r = H H. Define: C = {(x 1, x 2,..., x r ) : x i C i }, D = {(x 1, x 2,..., x r ) : x 1 = x i } It is easily verified that: (P C x) i = P Ci x i, (P D x) i = 1 r r j=1 x j Each iteration, P D P C : H r H r, can be described by: (P D P C x) i = 1 r r P Cj x j j=1
46 Averaged Projections (The Crucial Product Space Trick) Given C 1, C 2,..., C r H consider H r = H H. Define: C = {(x 1, x 2,..., x r ) : x i C i }, D = {(x 1, x 2,..., x r ) : x 1 = x i } It is easily verified that: (P C x) i = P Ci x i, (P D x) i = 1 r r j=1 x j Each iteration, P D P C : H r H r, can be described by: (P D P C x) i = 1 r r P Cj x j=1
47 Averaged Projections (The Crucial Product Space Trick) Given C 1, C 2,..., C r H consider H r = H H. Define: C = {(x 1, x 2,..., x r ) : x i C i }, D = {(x 1, x 2,..., x r ) : x 1 = x i } It is easily verified that: (P C x) i = P Ci x i, (P D x) i = 1 r r j=1 x j Each iteration, T : H H, can be described by: T x = 1 r r P Cj x j=1
48 Douglas-Rachford and Dykstra Methods Theorem (Lions & Mercier, 1979) Let C 1, C 2 H be closed convex sets, x H iterate: x n+1 := x n + R C2 R C1 (x n ) 2 where R Ci (x) := 2P Ci (x) x then x n w. x, a fixed point, with P C1 (x) C 1 C 2. Theorem (Boyle & Dykstra, 1980) Let C 1,..., C r H be closed convex sets, x H iterate: x i n := P Ci (x i 1 n I i n 1), I i n := x i n (x i 1 n I i n 1), x 0 n := x r n 1 with initial values x 0 1 := x, I i 0 := 0 then x n (P r i=1 C i )(x).
49 Closed Subspaces Projections onto non-convex sets are no longer guaranteed to be: In R n : Nonexpansive/Firmly nonexpansive Unique (i.e. P C is set-valued). The method becomes: x 2n+1 P C1 (x 2n ), x 2n P C2 (x 2n 1 ) Local linear convergence for alternating and averaged nonconvex projections, Lewis, Luke & Malick (2009). Restricted normal cones and the method of alternating projections, Bauschke, Luke, Phan & Wang (2012). The Douglas-Rachford algorithm in the absence of convexity, Borwein & Sims (2011).
50 Sudoku: Modelling an NP-Complete Non-Convex Problem
51 Sudoku: Modelling an NP-Complete Non-Convex Problem
52 Sudoku: Modelling an NP-Complete Non-Convex Problem
53 Sudoku: Modelling an NP-Complete Non-Convex Problem Let A (R 9 ) 3 indexed by (i, j, k). Constraint types are: C 1 = {A ij is a standard unit vector} C 2 = {A ik is a standard unit vector} C 3 = {A jk is a standard unit vector} C 4 = {3 3 submatrix = standard unit vector} A solution is x C 1 C 2 C 3 C 4.
54 Sudoku: Modelling an NP-Complete Non-Convex Problem Let A (R 9 ) 3 indexed by (i, j, k). Constraint types are: C 1 = {A ij is a standard unit vector} C 2 = {A ik is a standard unit vector} C 3 = {A jk is a standard unit vector} C 4 = {3 3 submatrix = standard unit vector} A solution is x C 1 C 2 C 3 C 4.
55 Sudoku: Modelling an NP-Complete Non-Convex Problem Let A (R 9 ) 3 indexed by (i, j, k). Constraint types are: C 1 = {A ij is a standard unit vector} C 2 = {A ik is a standard unit vector} C 3 = {A jk is a standard unit vector} C 4 = {3 3 submatrix = standard unit vector} A solution is x C 1 C 2 C 3 C 4.
56 Sudoku: Modelling an NP-Complete Non-Convex Problem Let A (R 9 ) 3 indexed by (i, j, k). Constraint types are: C 1 = {A ij is a standard unit vector} C 2 = {A ik is a standard unit vector} C 3 = {A jk is a standard unit vector} C 4 = {3 3 submatrix = standard unit vector} A solution is x C 1 C 2 C 3 C 4.
57 Sudoku: Modelling an NP-Complete Non-Convex Problem Let A (R 9 ) 3 indexed by (i, j, k). Constraint types are: C 1 = {A ij is a standard unit vector} C 2 = {A ik is a standard unit vector} C 3 = {A jk is a standard unit vector} C 4 = {3 3 submatrix = standard unit vector} A solution is x C 1 C 2 C 3 C 4.
58 Sudoku: Modelling an NP-Complete Non-Convex Problem Let A (R 9 ) 3 indexed by (i, j, k). Constraint types are: C 1 = {A ij is a standard unit vector} C 2 = {A ik is a standard unit vector} C 3 = {A jk is a standard unit vector} C 4 = {3 3 submatrix = standard unit vector} A solution is x C 1 C 2 C 3 C 4. Similar modelling can be done for: N-queens 3-SAT (NP-Complete) TetraVex (NP-Complete) Solves large instances! (Sudoku = R 2916 )
59 J. von Neumann. Functional operators Vol. II. The geometry of orthogonal spaces. Annal. of Math. Studies, 22, I. Halperin. The product of projection operators. Acta. Sci. Math., 23:96 99, L.M. Bregman. The method of successive projection for finding a common point of convex sets. Soviet Mathematics, 6: , H.S. Hundal. An alternating projection that does not converge in norm. Nonlinear analysis, 57(1):35 61, E. Matoušková & S. Reich. The Hundal example revisited. J. Nonlinear Convex Anal, 4: , P.L. Lions & B. Mercier. Splitting algorithms for the sum of two nonlinear operators. SIAM Journal on Numerical Analysis, , J.P. Boyle & R.L. Dykstra. A method for finding projections onto the intersection of convex sets in Hilbert spaces. Advances in Order Restricted Statistical Inference, J. Schaad. Modelling the 8-queens problem and Sudoku using an algorithm based on projections onto nonconvex sets, Master s Thesis, 2010.
Alternating Projection Methods
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
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationarxiv: v1 [math.fa] 20 Aug 2015
STRANGE PRODUCTS OF PROJECTIONS EVA KOPECKÁ AND ADAM PASZKIEWICZ arxiv:1508.05029v1 [math.fa] 20 Aug 2015 Abstract. Let H be an infinite dimensional Hilbert space. We show that there exist three orthogonal
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 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 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 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 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 informationOn Slater s condition and finite convergence of the Douglas Rachford algorithm for solving convex feasibility problems in Euclidean spaces
On Slater s condition and finite convergence of the Douglas Rachford algorithm for solving convex feasibility problems in Euclidean spaces Heinz H. Bauschke, Minh N. Dao, Dominikus Noll and Hung M. Phan
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 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 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 informationAPPLICATIONS IN FIXED POINT THEORY. Matthew Ray Farmer. Thesis Prepared for the Degree of MASTER OF ARTS UNIVERSITY OF NORTH TEXAS.
APPLICATIONS IN FIXED POINT THEORY Matthew Ray Farmer Thesis Prepared for the Degree of MASTER OF ARTS UNIVERSITY OF NORTH TEXAS December 2005 APPROVED: Elizabeth M. Bator, Major Professor Paul Lewis,
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 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 informationGENERAL NONCONVEX SPLIT VARIATIONAL INEQUALITY PROBLEMS. Jong Kyu Kim, Salahuddin, and Won Hee Lim
Korean J. Math. 25 (2017), No. 4, pp. 469 481 https://doi.org/10.11568/kjm.2017.25.4.469 GENERAL NONCONVEX SPLIT VARIATIONAL INEQUALITY PROBLEMS Jong Kyu Kim, Salahuddin, and Won Hee Lim Abstract. In this
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 informationThere is no variational characterization of the cycles in the method of periodic projections
There is no variational characterization of the cycles in the method of periodic projections J.-B. Baillon, 1 P. L. Combettes, 2 and R. Cominetti 3 1 Université Paris 1 Panthéon-Sorbonne SAMM EA 4543 75013
More informationSTRONG CONVERGENCE OF AN IMPLICIT ITERATION PROCESS FOR ASYMPTOTICALLY NONEXPANSIVE IN THE INTERMEDIATE SENSE MAPPINGS IN BANACH SPACES
Kragujevac Journal of Mathematics Volume 36 Number 2 (2012), Pages 237 249. STRONG CONVERGENCE OF AN IMPLICIT ITERATION PROCESS FOR ASYMPTOTICALLY NONEXPANSIVE IN THE INTERMEDIATE SENSE MAPPINGS IN BANACH
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 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 informationA splitting minimization method on geodesic spaces
A splitting minimization method on geodesic spaces J.X. Cruz Neto DM, Universidade Federal do Piauí, Teresina, PI 64049-500, BR B.P. Lima DM, Universidade Federal do Piauí, Teresina, PI 64049-500, BR P.A.
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 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 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 informationConvergence Theorems for Bregman Strongly Nonexpansive Mappings in Reflexive Banach Spaces
Filomat 28:7 (2014), 1525 1536 DOI 10.2298/FIL1407525Z Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Convergence Theorems for
More informationPARALLEL SUBGRADIENT METHOD FOR NONSMOOTH CONVEX OPTIMIZATION WITH A SIMPLE CONSTRAINT
Linear and Nonlinear Analysis Volume 1, Number 1, 2015, 1 PARALLEL SUBGRADIENT METHOD FOR NONSMOOTH CONVEX OPTIMIZATION WITH A SIMPLE CONSTRAINT KAZUHIRO HISHINUMA AND HIDEAKI IIDUKA Abstract. In this
More informationSTRONG CONVERGENCE THEOREMS FOR COMMUTATIVE FAMILIES OF LINEAR CONTRACTIVE OPERATORS IN BANACH SPACES
STRONG CONVERGENCE THEOREMS FOR COMMUTATIVE FAMILIES OF LINEAR CONTRACTIVE OPERATORS IN BANACH SPACES WATARU TAKAHASHI, NGAI-CHING WONG, AND JEN-CHIH YAO Abstract. In this paper, we study nonlinear analytic
More informationSubgradient Method. Ryan Tibshirani Convex Optimization
Subgradient Method Ryan Tibshirani Convex Optimization 10-725 Consider the problem Last last time: gradient descent min x f(x) for f convex and differentiable, dom(f) = R n. Gradient descent: choose initial
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 informationChapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.
Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space
More informationFinding a best approximation pair of points for two polyhedra
Finding a best approximation pair of points for two polyhedra Ron Aharoni, Yair Censor Zilin Jiang, July 30, 2017 Abstract Given two disjoint convex polyhedra, we look for a best approximation pair relative
More informationShih-sen Chang, Yeol Je Cho, and Haiyun Zhou
J. Korean Math. Soc. 38 (2001), No. 6, pp. 1245 1260 DEMI-CLOSED PRINCIPLE AND WEAK CONVERGENCE PROBLEMS FOR ASYMPTOTICALLY NONEXPANSIVE MAPPINGS Shih-sen Chang, Yeol Je Cho, and Haiyun Zhou Abstract.
More informationReferences 1. Aleyner, A., & Reich, S. (2008). Block-iterative algorithms for solving convex feasibility problems in Hilbert and Banach spaces. Journa
References 1. Aleyner, A., & Reich, S. (2008). Block-iterative algorithms for solving convex feasibility problems in Hilbert and Banach spaces. Journal of Mathematical Analysis and Applications, 343, 427
More informationAuerbach bases and minimal volume sufficient enlargements
Auerbach bases and minimal volume sufficient enlargements M. I. Ostrovskii January, 2009 Abstract. Let B Y denote the unit ball of a normed linear space Y. A symmetric, bounded, closed, convex set A in
More informationProjection and proximal point methods: convergence results and counterexamples
Nonlinear Analysis 56 (004) 75 738 www.elsevier.com/locate/na Projection and proximal point methods: convergence results and counterexamples Heinz H. Bauschke a, Eva Matouskova b, Simeon Reich c; a Department
More informationConvergence theorems for mixed type asymptotically nonexpansive mappings in the intermediate sense
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 2016), 5119 5135 Research Article Convergence theorems for mixed type asymptotically nonexpansive mappings in the intermediate sense Gurucharan
More information1 The Projection Theorem
Several Important Theorems by Francis J. Narcowich November, 14 1 The Projection Theorem Let H be a Hilbert space. When V is a finite dimensional subspace of H and f H, we can always find a unique p V
More informationCompressed Sensing: a Subgradient Descent Method for Missing Data Problems
Compressed Sensing: a Subgradient Descent Method for Missing Data Problems ANZIAM, Jan 30 Feb 3, 2011 Jonathan M. Borwein Jointly with D. Russell Luke, University of Goettingen FRSC FAAAS FBAS FAA Director,
More informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
More informationConvergence in shape of Steiner symmetrized line segments. Arthur Korneychuk
Convergence in shape of Steiner symmetrized line segments by Arthur Korneychuk A thesis submitted in conformity with the requirements for the degree of Master of Science Graduate Department of Mathematics
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 information1991 Mathematics Subject Classification: 47H09, 47H10.
æ THE LERAY-SCHAUDER ALTERNATIVE FOR NONEXPANSIVE MAPS FROM THE BALL CHARACTERIZE HILBERT SPACE Michael Ireland Department of Mathematics The University of Newcastle Newcastle 2308, NSW, Australia William
More informationIterative Convex Optimization Algorithms; Part One: Using the Baillon Haddad Theorem
Iterative Convex Optimization Algorithms; Part One: Using the Baillon Haddad Theorem Charles Byrne (Charles Byrne@uml.edu) http://faculty.uml.edu/cbyrne/cbyrne.html Department of Mathematical Sciences
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 informationThe small ball property in Banach spaces (quantitative results)
The small ball property in Banach spaces (quantitative results) Ehrhard Behrends Abstract A metric space (M, d) is said to have the small ball property (sbp) if for every ε 0 > 0 there exists a sequence
More informationResearch Article Some Inequalities Concerning the Weakly Convergent Sequence Coefficient in Banach Spaces
Hindawi Publishing Corporation Abstract and Applied Analysis Volume 008, Article ID 80387, 8 pages doi:0.55/008/80387 Research Article Some Inequalities Concerning the Weakly Convergent Sequence Coefficient
More informationSpectral theory for compact operators on Banach spaces
68 Chapter 9 Spectral theory for compact operators on Banach spaces Recall that a subset S of a metric space X is precompact if its closure is compact, or equivalently every sequence contains a Cauchy
More informationIterative algorithms based on the hybrid steepest descent method for the split feasibility problem
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (206), 424 4225 Research Article Iterative algorithms based on the hybrid steepest descent method for the split feasibility problem Jong Soo
More informationThe method of cyclic intrepid projections: convergence analysis and numerical experiments
Chapter 1 The method of cyclic intrepid projections: convergence analysis and numerical experiments Heinz H. Bauschke, Francesco Iorio, and Valentin R. Koch Abstract The convex feasibility problem asks
More informationA NEW GEOMETRIC ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD
Electronic Transactions on Numerical Analysis. Volume 45,. 330 341, 2016. Coyright c 2016,. ISSN 1068 9613. ETNA A NEW GEOMETRIC ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD WILLIAMS LÓPEZ
More informationRenormings of c 0 and the minimal displacement problem
doi: 0.55/umcsmath-205-0008 ANNALES UNIVERSITATIS MARIAE CURIE-SKŁODOWSKA LUBLIN POLONIA VOL. LXVIII, NO. 2, 204 SECTIO A 85 9 ŁUKASZ PIASECKI Renormings of c 0 and the minimal displacement problem Abstract.
More informationvorgelegt von Robert Hesse aus Northeim
Fixed Point Algorithms for Nonconvex Feasibility with Applications Dissertation zur Erlangung des mathematisch-naturwissenschaftlichen Doktorgrades Doctor rerum naturalium im Promotionsprogramm PhD School
More informationSplitting methods for decomposing separable convex programs
Splitting methods for decomposing separable convex programs Philippe Mahey LIMOS - ISIMA - Université Blaise Pascal PGMO, ENSTA 2013 October 4, 2013 1 / 30 Plan 1 Max Monotone Operators Proximal techniques
More informationExistence and Approximation of Fixed Points of. Bregman Nonexpansive Operators. Banach Spaces
Existence and Approximation of Fixed Points of in Reflexive Banach Spaces Department of Mathematics The Technion Israel Institute of Technology Haifa 22.07.2010 Joint work with Prof. Simeon Reich General
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 informationExistence of Solutions to Split Variational Inequality Problems and Split Minimization Problems in Banach Spaces
Existence of Solutions to Split Variational Inequality Problems and Split Minimization Problems in Banach Spaces Jinlu Li Department of Mathematical Sciences Shawnee State University Portsmouth, Ohio 45662
More informationITERATIVE SCHEMES FOR APPROXIMATING SOLUTIONS OF ACCRETIVE OPERATORS IN BANACH SPACES SHOJI KAMIMURA AND WATARU TAKAHASHI. Received December 14, 1999
Scientiae Mathematicae Vol. 3, No. 1(2000), 107 115 107 ITERATIVE SCHEMES FOR APPROXIMATING SOLUTIONS OF ACCRETIVE OPERATORS IN BANACH SPACES SHOJI KAMIMURA AND WATARU TAKAHASHI Received December 14, 1999
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More information08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms
(February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops
More informationA New Modified Gradient-Projection Algorithm for Solution of Constrained Convex Minimization Problem in Hilbert Spaces
A New Modified Gradient-Projection Algorithm for Solution of Constrained Convex Minimization Problem in Hilbert Spaces Cyril Dennis Enyi and Mukiawa Edwin Soh Abstract In this paper, we present a new iterative
More informationOn Total Convexity, Bregman Projections and Stability in Banach Spaces
Journal of Convex Analysis Volume 11 (2004), No. 1, 1 16 On Total Convexity, Bregman Projections and Stability in Banach Spaces Elena Resmerita Department of Mathematics, University of Haifa, 31905 Haifa,
More informationA New Use of Douglas-Rachford Splitting and ADMM for Identifying Infeasible, Unbounded, and Pathological Conic Programs
A New Use of Douglas-Rachford Splitting and ADMM for Identifying Infeasible, Unbounded, and Pathological Conic Programs Yanli Liu Ernest K. Ryu Wotao Yin October 14, 2017 Abstract In this paper, we present
More informationOptimal rates of linear convergence of relaxed alternating projections and generalized Douglas-Rachford methods for two subspaces
Optimal rates of linear convergence of relaxed alternating projections and generalized Douglas-Rachford methods for two subspaces Heinz H. Bauschke, J.Y. Bello Cruz, Tran T.A. Nghia, Hung M. Phan, and
More informationOn the split equality common fixed point problem for quasi-nonexpansive multi-valued mappings in Banach spaces
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (06), 5536 5543 Research Article On the split equality common fixed point problem for quasi-nonexpansive multi-valued mappings in Banach spaces
More informationProjection Methods for Quantum Channel Construction
Noname manuscript No. (will be inserted by the editor) Projection Methods for Quantum Channel Construction Dmitriy Drusvyatskiy Chi-Kwong Li Diane Christine Pelejo Yuen-Lam Voronin Henry Wolkowicz Received:
More informationWeak and strong convergence theorems of modified SP-iterations for generalized asymptotically quasi-nonexpansive mappings
Mathematica Moravica Vol. 20:1 (2016), 125 144 Weak and strong convergence theorems of modified SP-iterations for generalized asymptotically quasi-nonexpansive mappings G.S. Saluja Abstract. The aim of
More informationStrong Convergence Theorems for Nonself I-Asymptotically Quasi-Nonexpansive Mappings 1
Applied Mathematical Sciences, Vol. 2, 2008, no. 19, 919-928 Strong Convergence Theorems for Nonself I-Asymptotically Quasi-Nonexpansive Mappings 1 Si-Sheng Yao Department of Mathematics, Kunming Teachers
More informationAsynchronous Algorithms for Conic Programs, including Optimal, Infeasible, and Unbounded Ones
Asynchronous Algorithms for Conic Programs, including Optimal, Infeasible, and Unbounded Ones Wotao Yin joint: Fei Feng, Robert Hannah, Yanli Liu, Ernest Ryu (UCLA, Math) DIMACS: Distributed Optimization,
More informationIn English, this means that if we travel on a straight line between any two points in C, then we never leave C.
Convex sets In this section, we will be introduced to some of the mathematical fundamentals of convex sets. In order to motivate some of the definitions, we will look at the closest point problem from
More informationSPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint
More informationLocal convergence of the method of alternating projections
Local convergence of the method of alternating projections Dominikus Noll Université Paul Sabatier Alternating projections invented by Hermann Amandus Schwarz in 1869 H. A. Schwarz. Über einen Grenzübergang
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 informationCONVERGENCE OF HYBRID FIXED POINT FOR A PAIR OF NONLINEAR MAPPINGS IN BANACH SPACES
International Journal of Analysis and Applications ISSN 2291-8639 Volume 8, Number 1 2015), 69-78 http://www.etamaths.com CONVERGENCE OF HYBRID FIXED POINT FOR A PAIR OF NONLINEAR MAPPINGS IN BANACH SPACES
More informationThe Split Common Fixed Point Problem for Asymptotically Quasi-Nonexpansive Mappings in the Intermediate Sense
International Mathematical Forum, Vol. 8, 2013, no. 25, 1233-1241 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2013.3599 The Split Common Fixed Point Problem for Asymptotically Quasi-Nonexpansive
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 informationThe Geometric Approach for Computing the Joint Spectral Radius
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 12-15, 2005 TuB08.2 The Geometric Approach for Computing the Joint Spectral
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 informationIdentifying Active Constraints via Partial Smoothness and Prox-Regularity
Journal of Convex Analysis Volume 11 (2004), No. 2, 251 266 Identifying Active Constraints via Partial Smoothness and Prox-Regularity W. L. Hare Department of Mathematics, Simon Fraser University, Burnaby,
More informationCONVERGENCE TO FIXED POINTS OF INEXACT ORBITS FOR BREGMAN-MONOTONE AND FOR NONEXPANSIVE OPERATORS IN BANACH SPACES
CONVERGENCE TO FIXED POINTS OF INEXACT ORBITS FOR BREGMAN-MONOTONE AND FOR NONEXPANSIVE OPERATORS IN BANACH SPACES Dan Butnariu, Simeon Reich, and Alexander J. Zaslavski Abstract. Fixed points of Bregman-monotone
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 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 information2nd Symposium on System, Structure and Control, Oaxaca, 2004
263 2nd Symposium on System, Structure and Control, Oaxaca, 2004 A PROJECTIVE ALGORITHM FOR STATIC OUTPUT FEEDBACK STABILIZATION Kaiyang Yang, Robert Orsi and John B. Moore Department of Systems Engineering,
More informationECON 4117/5111 Mathematical Economics
Test 1 September 29, 2006 1. Use a truth table to show the following equivalence statement: (p q) (p q) 2. Consider the statement: A function f : X Y is continuous on X if for every open set V Y, the pre-image
More information