J.B. LASSERRE AND E.S. ZERON

Size: px
Start display at page:

Download "J.B. LASSERRE AND E.S. ZERON"

Transcription

1 L -NORMS, LOG-BARRIERS AND CRAMER TRANSFORM IN OPTIMIZATION J.B. LASSERRE AND E.S. ZERON Abstract. We show that the Lalace aroxiation of a sureu by L -nors has interesting consequences in otiization. For instance, the logarithic barrier functions (LBF) of a rial convex roble P and its dual P aear naturally when using this sile aroxiation technique for the value function g of P or its Legendre-Fenchel conjugate g. In addition, iniizing the LBF of the dual P is just evaluating the Craer transfor of the Lalace aroxiation of g. Finally, this technique erits to soeties define an exlicit dual roble P in cases when the Legendre-Fenchel conjugate g cannot be derived exlicitly fro its definition. 1. Introduction Let f : R and ω : R be a air of continuous aings defined on the convex cone R n. Consider the function g : R R given by the forula: (1.1) y g(y) := su f(x) : ω(x) y, x. x For each fixed y R, couting g(y) is solving the otiization roble (1.2) P : su f(x) : ω(x) y, x, x and g is called the value function associated with P. The value function g rovides a systeatic way to generate a dual roble P via its Legendre- Fenchel conjugate denoted g : R R. The concave version of g is defined by (1.3) g () := inf y R y g(y), and is finite on soe doain D R. Then the dual roble reduces to (1.4) (1.5) P : g(y) = inf y g () = inf y g () : D = inf R + su f(x) + (y ω(x) Matheatics Subject Classification. 90C25 90C26. Key words and hrases. otiization; Logarithic Barrier Function; Legendre-Fenchel and Craer transfors. 1

2 2 J.B. LASSERRE AND E.S. ZERON Weak duality ilies g(y) g(y), that is, (1.6) (1.7) g(y) = su inf R + inf f(x) + (y ω(x) R + su f(x) + (y ω(x) = g(y), and the equality g(y) = g(y) holds true under soe convexity assution. However, in general g cannot be obtained exlicitly fro its definition (1.3), and for dual ethods to solve P, the inner axiization in (1.7) ust be done nuerically for each fixed. A notable excetion is the conic otiization roble where f and ω are both linear aings, for which one ay obtain an exlicit dual (1.4). Of course, alternative exlicit duals have been roosed but they involve both rial (x) and dual () variables. In articular, the Wolfe [14] and Mond-Weir [11] duals even allow to consider weakened notions of convexity like e.g. seudo- or quasi-convexity. For a nice exosition and related references on this toic, the interested reader is referred to Mond [12] and the references therein. Contribution. Our contribution is to show that the sile and wellknown Lalace aroxiation of a sureu via a converging sequence of L -nors has interesting consequences in otiization, for both rial and dual robles P and P. Recall that the celebrated logarithic barrier function (LBF in short) associated with a convex otiization roble P as in (1.2), or with its dual P, is well-known for its good technical roerties like e.g. the selfconcordance (which holds in soe cases) which exlains its good nuerical efficiency; see e.g. [5, 6]. But the LBF is only one articular choice aong any other interior enalty functions! We first show that the LBF of the rial roble P (with araeter ) aears naturally by using the sile and well-known Lalace aroxiation of a sureu via L -nors, alied to the inner infiu in (1.6). It is a bit surising to obtain an efficient ethod in this way. Indeed, the inner infiu in (1.6) (which is exactly equal to zero for a feasible x) is relaced with its naive Lalace aroxiation by L -nors, and to the best of our knowledge, the efficiency of this aroxiation has not been roved or even tested nuerically! Siilarly, when using the sae Lalace L -nor aroxiation technique for the infiu in the definition (1.3) of the conjugate function g, we obtain a function φ : R R which: (a) deends on an integer araeter and (b), is valid on the relative interior ri D of soe doain D R. In doing so for conic otiization robles, the set D is just the feasible set of the (known) exlicit dual roble P, and φ is (u to a constant) the LBF with araeter, associated with P. So again, for conic rogras, the sile Lalace aroxiation of a sureu by L -nors erits to retrieve the LBF of the dual roble P! Interestingly, the function y in φ (; y) is nothing less than the Craer transfor of the Lalace aroxiation e f L (Ω(y)), where Ω(y)

3 L -NORMS AND CRAMER TRANSFORM IN OPTIMIZATION 3 is the feasible set of roble P and L is the usual nor associated with the Lebesgue sace L. Analogies between the Lalace and Fenchel transfors via exonentials and logariths in the Craer transfor have been already exlored in other contexts, in order to establish nice arallels between otiization and robability via a change of algebra; see e.g. Bacelli et al. [1], Maslov [10], Lasserre [9], and the any references therein. In robability, the Craer transfor of a robability easure has also been used to rovide exact asytotics of soe integrals as well as to derive large deviation rinciles. For a nice survey on this toic the interested reader is referred to Piterbarg and Falatov [13]. In addition, an interesting feature of this Lalace aroxiation technique is to rovides us with a systeatic way to obtain a dual roble (1.4) in cases when g cannot be obtained exlicitly fro its definition (1.3). Naely, in a nuber of cases and in contrast with g, the function φ (; y) obtained by using the Lalace aroxiation of the conjugate function g by L -nors, can be couted in closed-for exlicitly. Exales of such situations are briefly discussed. In the general case, φ is of the for φ(; y) + ψ(, ) where: for every ri D fixed, the function ψ(, ) 0 as ; and for each fixed, the function ψ(, ) is a barrier for the doain D. This yields to consider the otiization roble P : in φ(; y) : D as a natural dual of P, and for which φ is an associated barrier function with araeter. If g is concave then strong duality holds. 2. Main result We need soe interediary helful results before stating our ain result Soe reliinary results. Let L q () be the usual Lebesgue sace of integrable functions defined on a Borel-easurable set R n, and h L q () (or soeties h q ) be the associated nor h L q () = h q := ( 1/q h(x) dx) q. To ake the aer self-contained we rove the following known result. Lea 1. Let R n be any Borel-easurable set, and h L q () for soe given q 1, so that h L q () <. Then: li h L () = h := ess su h(x). Proof. Notice that ay be an unbounded set. Suose that h q < for soe given q > 1, and define Λ to be the essential sureu of h in

4 4 J.B. LASSERRE AND E.S. ZERON. The result is trivial when Λ = 0, so we assue that Λ (0, ). Then ( ess su h(x) = Λ = li h/λ L q () [ = li [ ( h(x) Λ Λ ) q/λ ) q dx] 1/ 1/ (2.1) li h(x) dx] = li h L (). It is also obvious that Λ li h when Λ =. On the other hand, suose that the essential sureu Λ of h in is finite. Given an arbitrary araeter ɛ > 0, there exists a bounded subset B with ositive finite Lebesgue easure (B) (0, ) such that h(x) > Λ ɛ for every x B. Then li h L () li h L (B) li (B)1/ (Λ ɛ) = Λ ɛ. Therefore, since ɛ is arbitrary, cobining the revius identity with (2.1) yields the desired result li h L () = Λ. In the sae way, assue that the essential sureu of h in is infinite. Given an arbitrary natural nuber N N, there exists a bounded subset B with ositive finite Lebesgue easure (B) (0, ) such that h(x) > N for every x B. Then li h L () li h L (B) li (B)1/ N = N. Therefore, since N is arbitrary, cobining the revius identity with (2.1) yields the desired result li h L () = Λ =. Next we also need the following interediate result. Lea 2. For every N let U R n be soe oen subset, and let h : U R be a sequence of functions indexed by the araeter N. Suose that h converges ointwise to a function h defined on an oen subset U of R n. Then: li inf h (x) inf h(x), x U x U rovided that the liit in the left side of the equation exists in the extended interval [, ). Proof. Suose that the infiniu of h on U is equal to. For every N R there is a oint x U such that h(x) < N, and so there is also an index 0 such that x U and h (x) < N for every > 0. Hence the infiniu of h on U is strictly less than N, and so li inf h (x) = = inf h(x), x U x U because N R is arbitrary. On the other hand, assue that the infiniu of h on U is equal to R. For every ɛ > 0 there is a oint x U such that h(x) < +ɛ, and so there is also an index 0 such that x U and

5 L -NORMS AND CRAMER TRANSFORM IN OPTIMIZATION 5 h (x) < +ɛ for every > 0. Since the infiniu of h on U is strictly less than +ɛ and ɛ > 0 is arbitrary, li inf h (x) = inf h(x). x U x U 2.2. L -nor aroxiations for the rial. Let us go back to roble P in (1.1) where R n is a convex cone, and let Z := R +. Let R n be the dual cone associated with, and let x (x) be the universal logarithic barrier function associated with the convex cone, that is, ( ) x (x) := ln e x y dy, x int, where int denotes the interior of. See e.g. Güller [4] and Güler and Tuncel [5]. Next, let H R n be the set H := x R n : ω(x) < y; x int. Recalling that P is a axiization roble, the LBF associated with the (rial) roble P, and with araeter N, is the function ψ : H R defined by: (2.2) x ψ (x) := f(x) + 1 (x) + ln(y ω(x)) j, The LBF in convex rograing dates back to Frisch [3] and becae widely known later in Fiacco and McCorick [2]. For ore details and a discussion, see e.g. den Hertog [6, Chater 2]. It is well-known that under soe convexity assutions, and if g(y) <, (2.3) g(y) = li su x ψ (x) : x H, and the sequence of iniizers (x ) N H of ψ converges to a iniizer of P. We next rovide a sile rationale that exlains how the LBF naturally aears to solve roble P. Observe that (1.6) can be rewritten as (2.4) g(y) = su x R n inf f(x) + (,µ) Z (y ω(x)) + µ x, and in the above equation, the inner iniization inf f(x) + (,µ) Z (y ω(x)) + µ x, (whose value is 0 if ω(x) y and x ), ay be rewritten as f(x) su (ω(x) y) µ x. (,µ) Z

6 6 J.B. LASSERRE AND E.S. ZERON But for fixed x R n and y R, one ay aroxiate the above sureu via L -nors. Indeed, if y R and x H, then for every N, e (ω(x) y) µ x dµ d <, Z and so e (ω(x) y) µ x L (Z ) <, N. Therefore, by Lea 1, su (ω(x) y) µ x = li ln (ω(x) y) µ x (,µ) Z e L (Z ) = li ( ) 1 ln e (ω(x) y) µ x dµ d Z 1 (x) 1 ln(y ω(x)) j li if x H = + otherwise. Next, observe that for each N, (x) = n (x) because is a cone, and so (2.4) now reads 1 (2.5) g(y) = su f(x) + li (x) + ln(y ω(x)) j x H. A direct alication of Lea 2 to (2.5) yields g(y) li su 1 f(x) + li (x) + x H = li su x ψ (x) : x H,, ln(y ω(x)) j and so (2.3) states that in fact one also has the reverse inequality, that is, in (2.5) one ay interchange the su and li oerators. In other words, the LBF ψ aears naturally when one aroxiates inf f(x) + (,µ) Z (y ω(x)) + µ x, (whose value is exactly zero if ω(x) y and x ), by the quantity 1 (x) + ln(y ω(x)) j, which coes fro the Lalace aroxiation of a su by L -nors. For instance, in Linear rograing where = R n +, x c x and x ω(x) = Ax for soe vector c R n and soe atrix A R n, g(y) = li su c x + 1 n ln(y Ax) j + ln(x i ) : x H. x i=1

7 L -NORMS AND CRAMER TRANSFORM IN OPTIMIZATION L -nor aroxiations for the dual. We now use the sae aroxiation technique via L -nors to either retrieve the known dual P when it is exlicit, or to rovide an exlicit dual roble P in cases where g cannot be obtained exlicitly fro its definition (1.3). Recall that if g is concave, uer sei-continuous, and bounded fro above by soe linear function, then by Legendre-Fenchel duality, (2.6) (2.7) g(y) = inf g () := inf y y g (), where y g(y). One can exress g in ters of the definition (1.1) of g and the involved continuous transforations f and ω. Naely, g () = su g(y) y y (2.8) = su su f(x) y y, ω(x) y su f(x) ω(x) if 0, (2.9) = + otherwise. Therefore the doain of definition D R of g is given by: (2.10) D := R : 0, su f(x) ω(x) <, with relative interior denoted by ri D. Observe that D is convex because g is convex and roer on D. Theore 3. Let g and g be as in (1.1) and (2.7), resectively. Assue that g is concave, uer sei-continuous, and bounded fro above by soe linear function. Suose that the relative interior ri D is not ety and for every ri D there exists an exonent q 1 such that (2.11) e f(x) ω(x) <. L q () Then: (2.12) g(y) = li inf y + ln ri D Proof. In view of (2.8) (2.13) g () = [ ln su y e f(x) ω(x) L () ] su e y+f(x) if D ω(x) y + otherwise. ln( j )

8 8 J.B. LASSERRE AND E.S. ZERON Hyothesis (2.11) and Lea 1 allow us to relace the sureu in (2.13) by the liit of the L -nors as. Naely, ( ) 1/ g () = li ln e y+f(x) dy dx (2.14) = li ln ( = li ln ω(x) y 1/ e dx) ω(x)+f(x) e f(x) ω(x) L () ln( j ) Hence fro (2.14), equation (2.6) can be rewritten as follows: g(y) = inf y g () = (2.15) = inf ri D li ri D li inf ri D y + ln y + ln e f(x) ω(x) e f(x) ω(x) L () L (). ln( j ) ln( j ) ln( j ) where we have alied Lea 2 in order to interchange the inf and li oerators. Notice that the ters between the brackets are the functions h () of Lea 2. On the other hand, given y R, let Θ(y) := (x, z) R + : ω(x) + z y R n+, so that whenever ri D, e f(x) L (Θ(y)) e f(x)+ (y ω(x) z) L (Θ(y)) e y = e y e f(x) ω(x) z L ( R + ) e f(x) ω(x) L () ( j ) 1/. By hyothesis (2.11), given ri D fixed, the L -nor ter in the last above identity is finite for soe large enough. Therefore, by definition (1.1) and Lea 1, one obtains g(y) = su ln e f(x) e = li ln f(x) li (x,z) Θ(y) inf ri D y + ln e f(x) ω(x) L () L (Θ(y)) which cobined with (2.15) yields the desired result (2.12)., ln( j ),

9 L -NORMS AND CRAMER TRANSFORM IN OPTIMIZATION 9 Reark 1. Condition (2.11) can be easily checked in articular cases. For instance let x ω(x) := Ax for soe atrix A R n. If := R n + and x f(x) := c x + ln x d for soe c and d in R n with d 0. The notation x d stands for the onoial n xd k k. Then f is concave and e f(x) ω(x) dx = e (c A ) x x d dx R n + = R n + n Γ(1 + d k ) (A k c k) 1+ d k <, whenever ri D := R : > 0, A > c and N. If = R and x f(x) := x Qx + c x for soe c R n and a syetric (strictly) ositive definite atrix Q R n n, then e f(x) ω(x) dx = e x Qx e (c A ) x dx <, R n R n whenever ri D := : > 0 and N. Consider next the following functions for every N : (2.16) φ (; y) := y + ln e f(x) ω(x) L () defined on soe doain of R, and (2.17) y g (y) := inf φ (; y) : ri D ln( j ) defined on R. Recall that the Craer transfor (denoted C) alied to an integrable function u : R R, is the Legendre-Fenchel transfor (denoted F) of the logarith of the Lalace transfor (denoted L) of u, i.e., u C(u) = F ln L (u). The Craer transfor is natural in the sense that the logarith of the Lalace transfor is always a convex function. For our urose, we will consider the concave version of the Fenchel transfor (2.18) û [F(û)]() = inf y y + û(y), for û : R R convex, so that û is concave. We clai that: Theore 4. The function y g (y) defined in (2.17) is the Craer transfor of the function (2.19) y g (y) := e f(x) dx = e f, L (Ω(y)) Ω(y) where Ω(y) := x : ω(x) y R n.

10 10 J.B. LASSERRE AND E.S. ZERON Proof. The result follows fro the definition of the Craer transfor C. Hence Therefore, g C( g ) := F ln L ( g ) y C( g )(y) = inf y + [ln L( g )](). [L( g )]() = e y g (y) dy = y R [ ] = e y e f(x) dx dy y R, ω(x) y [ ] = e f(x) e y dy dx y ω(x) [ ] = e f(x) ω(x) 1 dx j = e f(x) ω(x) 1. L () j [ln L( g )]() = ln e f(x) ω(x) ln( j ). L () On the other hand, recall the definition of g (y) given in (2.17)-(2.16), g (y) = inf y + ln e f(x) ω(x) ln( j ). ri D L () Thus, with F as in (2.18) and D := z : z D, we obtain the desired result : g (y) = inf y + ln e f(x) ω(x) ln( j ) ri D L () = inf y + [ln L( g )]() ri D = inf z y + [ln L( g )](z) z ri D = [F ln L( g )](y) = [C( g )](y), For linear rograing, this result was already obatined in [8, 9]. Exale 1. (Linear Prograing) In this case set the cone = R n + and the functions f(x) := c x and ω(x) = Ax for soe vector c R n and atrix

11 L -NORMS AND CRAMER TRANSFORM IN OPTIMIZATION 11 A R n. We easily have that = L () e f(x) ω(x) e (c A ) x dx = n 1 A k c k for every N and each in the relative interior ri D of the set (2.20) D = R : A c, 0. Hence fro (2.16) φ (; y) = y + 1 e ln f(x) ω(x) L () = y n ln(a k c k) ln( j ) ln( j ) +n = ln, One easily recognizes (u to the constant (+n)[ln ]/) the LBF with araeter, of the dual roble: P : in y : A c, 0. Exale 2. (The general conic roble) Consider the conic otiization roble in c x : ω x y, x, x for soe convex cone R n, soe vector c R n, and soe linear aing ω : R n R with adjoint aing ω : R R n. We easily have that e c x ω x e = (c ω ) x = L () L () = e (c ω ) x dx = n e (c ω ) x dx, because is a cone. Clai (2.16) reads φ (; y) = y + 1 e ln c x ω x L () (2.21) = y + ψ(ω c) ln j +n ln( j ) ln, where ψ : R n R is the so-called universal LBF associated with the dual cone, and with doain ri D, where (2.22) D = R : ω c, 0. See e.g. Güller [4] and Güler and Tuncel [5]. In φ (and u to a constant), one easily recognizes the LBF with araeter, of the dual roble: P : in y : ω c, 0.

12 12 J.B. LASSERRE AND E.S. ZERON Exale 3. (Quadratic rograing: non conic forulation) Consider syetric ositive seidefinite atrixes Q [j] R n n and vectors c [j] R n for j = 0, 1,...,. The notation Q 0 (res. Q 0) stands for Q is ositive seidefinite (res. strictly ositive definite). Let := R n, f(x) := x Q [0] x 2c [0] x and let ω : Rn R have entries ω j (x) := x Q [j] x+2c [j] x for every j = 1,...,. For R with > 0, define the real syetric atrix Q [] R n n and vector c [] R n : so that Q [] := Q [0] + j Q [j] and c [] := c [0] + e f(x) ω(x) L () j c [j], = ex ( x Q [] x 2 c [] x) dx = π ex ( c n/2 [] Q 1 [] c ) [] det ( ) <, Q [] whenever N and Q [] 0. Therefore (2.23) in 0, Q 0 φ (; y) = y + 1 e ln c x ω x L () ax = y + c [] Q 1 [] c [] ln ( ) det Q [] 2 ln j ln π ln + n ln 2, ln( j ) on the doain of definition ri D := : > 0, Q 0. Again, in equation (2.23) one easily recognizes (u to a constant) the LBF with araeter, of the dual roble P : x Q [0] x 2c [0] x ( j x ) Q [j] x+2c [j] x y j = in 0,Q 0 y + ax x Q [] x 2c [] x) = in y + c [] Q 1 [] c [] : 0, Q [] 0 where we have used the fact that x = Q 1 [] c [] R n is the unique otial solution to the inner axiization roble in the second equation above. If Q 0 0 and Q j 0, j = 1,...,, then ri D := : > 0 because Q 0 whenever > 0; in this case P is a convex otiization roble and there is no duality ga between P and P.,

13 L -NORMS AND CRAMER TRANSFORM IN OPTIMIZATION An exlicit dual. In Exales 1, 2, and 3, the function φ defined in (2.16) can be decoosed into a su of the for: (2.24) φ (; y) = h 1 (; y) + h 2 (; ) where h 1 is indeendent of the araeter. Moreover, if h 2 (; ) < for soe > 0 fixed, the ter h 2 (; ) converges to zero when. One ay also verify that h 1 () = y g (), D, where g is the Legendre-Fenchel conjugate in (2.7). In fact, the above revious decoosition is ore general and can be deduced fro soe sile facts. Recall the roble P given in (1.1), y g(y) := su f(x) : ω(x) y, x, x for a convex cone R n and continuous aings f and ω. Assue that g is concave, uer sei-continuous, and bounded fro above by soe affine function, so that Legendre-Fenchel duality yields g(y) = inf g () := inf y y g (), where y g(y), and where the doain of g is the set D R in (2.10). Lea 5. Suose that ri D is not ety and for every ri D there exists an exonent q 1 such that (2.11) holds. Then: (2.25) h 1 (; y) := li φ (; y) = y g (), ri D. Proof. Let ri D be fixed. The given hyothesis and Lea 1 ily that e li ln f(x) ω(x) = su f(x) ω(x) = L () = su f(x) y = su g(y) y = g (), y y su, ω(x) y and so (2.25) follows fro the definition (2.16) of φ. As a consequence we obtain: Corollary 6. Let D be as in (2.10), φ as in (2.16) and let h 1 (; y) be as in (2.25). Then the otiization roble (2.26) P : in h 1 (; y) : ri D. is a dual of P. Moreover, if g is concave, uer-seicontinuous and bounded above by soe affine function, then strong duality holds.

14 14 J.B. LASSERRE AND E.S. ZERON Proof. By Lea 5, h 1 (; y) = y g () for all ri D. And so if in P (res. ax P) denotes the otial value of P (res. P), one has in P = in y g () : ri D g(y) = ax P, where we have used that g () = su z g(z) z g(y) y. Finally, if g is concave, uer sei-continuous and bounded above by soe affine function, then g is convex. Therefore, as a convex function is continuous on its doain D (which is convex) in P = in h 1 (; y) : ri D that is, strong duality holds. = in y g () : ri D = in y g () : D = g(y), In a nuber of cases, the L -nor aroxiation of g can be obtained exlicitly as a function of, whereas g itself cannot be obtained exlicitly fro (1.3). In this situation one obtains an exlicit LBF φ with araeter, for soe dual P of P, and soeties an exlicit dual roble P. Indeed if φ is known exlicitly, one ay soeties get its ointwise liit h 1 (, y) in (2.25), in closed for, and so P is defined exlicitly by (2.26). With fixed, couting φ (; y) reduces to coute the integral over a convex cone of an exonential of soe function araetrized by and. Soeties this can be done with the hel of soe known transfors like e.g. the Lalace or Weierstrass transfors, as illustrated below. Linear aings and Lalace transfor. Let ω : R n R be a linear aing, with ω(x) = Ax for soe real atrix A R n, and let = R n +. Then ln e f(x) ω(x) L () = 1 ( ) ln e (A ) x e f(x) dx = 1 ln (L[e f ](A ) That is, the L -nor aroxiation is the logarith of the Lalace transfor of the function e f, evaluated at the oint A R n. So if in roble P, the objective function f is such that e f has an exlicit Lalace transfor, then one obtains an exlicit exression for the LBF φ (; y) defined in (2.16). For instance if f(x) = c x + ln q(x) for soe vector c R n and soe olynoial q R[x], ositive on the feasible set of P, write q(x) = α N n q α x α, ).

15 L -NORMS AND CRAMER TRANSFORM IN OPTIMIZATION 15 for finitely any non zero coefficients (q α ), and where the notation x α stand for the onoial x α 1 1 xαn n. Then ln ( L[e f ](A ) ) can be couted in closed-for since we have: ( ) ln L[e f ](A ) where α x α = n i=1 α i x α i i = ln ( α N n q α e (c A ) x x α dx ( ) = n ln + ln q α 1 α x α n, α N n i=1 (A c) i. Of course the above exression can becoe quite colicated, esecially for large values of. But it is exlicit in the variables ( i ). If the function x ln q(x) is concave, then Corollary 6 alies. Siilarly if f is linear, i.e. x f(x) = c x for soe vector c R n, then ( ) ln e f(x) ω(x) L () = 1 ln L[e ω(x) ](c) ), and so if the function x e ω(x) has an exlicit Lalace transfor then so does the L -nor aroxiation, and again, φ is obtained in closed for. Exale 4. As a sile illustrative exale, consider the otiization roble: P : su x c x + n b k ln x k : Ax y, x 0, x R n, for soe given atrix A R n and vectors b, c R n and y R. We suose that b 0, so that c x + ln(x b ) is concave, and in which case P is a convex rogra. Notice that with = R n +, su c x + ln(x b ) Ax : x R n, x x whenever lies in D = R : A < c, 0. We have < e c x Ax x b L () = = n e (c A ) x x b dx Γ(1+ b k ) (A k c k) 1+ b k <,

16 16 J.B. LASSERRE AND E.S. ZERON whenever N and > 0 in R satisfies A > c. Next, φ in (2.16) reads φ (; y) = y 1 e ln c x Ax x b n ln( k ) L () n [ ] ln Γ(1+ = bk ) y + b k ln(a k c k) n ln(a k c k) ln j + n ln. Stirling s aroxiation Γ(1 + t) (t/e) t 2πt for real nubers t 1 allow us to calculate the liit when goes to infinity [ ] ln Γ(1+ b k ) li b k ln = li b bk k ln b k ln e Lea 5 ilies that y g () = li φ (; y) = y = b k ln(e 1 b k ). n [ A ] b k ln k c k e 1. b k And so the function φ is the LBF with araeter, of the dual roble n [ A P : inf ] y b k ln k c k e 1 : A c, 0. b k In articular, by Corollary 6, strong duality holds, i.e., ax P = in P. References [1] F. Bacelli, G. Cohen, J. Olsder and J.P Quadrat, Syncronization and Linearity, Wiley, New York, [2] A.V. Fiaco, G.P. McCorick, Nonlinear Prograing, Sequential Unconstrained Miniization Techniques, John Wiley & Sons, New York, [3] K.R. Frisch, The logarithic otential ethod of convex rograing, Meorandu, Institute of Econoics, Oslo, Norway, [4] O. Güller, Barrier functions in interior oint ethods, Math. Oer. Res. 21 (1996), [5] O. Güller and L. Tuncel, Characterization of the barrier araeter of hoogeneous convex cones, Math. Progr. 81 (1998), [6] D. den Hertog, Interior Point Aroach to Linear, Quadratic and Convex Prograing, Kluwer, Dordrecht, [7] J.B. Hiriart-Urruty, Conditions for global otiality, in: Handbook of Global Otiization, R. Horst and P. Pardalos (Eds.), Kluwer, Dordrecht (1995), [8] J.B. Lasserre, Why the logarithic barrier function in convex and linear rograing, Oer. Res. Letters 27 (2000), [9] J.B. Lasserre, Linear and Integer Prograing Versus Linear Integration and Counting, Sringer, New York, [10] V.P. Maslov, Méthodes Oératorielles, Editions Mir, Moscou 1973, Traduction Fran caise, 1987.

17 L -NORMS AND CRAMER TRANSFORM IN OPTIMIZATION 17 [11] B. Mond and T. Weir, Generalized convexity and higher order duality, J. Math. Sci ( ), [12] B. Mond, Mond-Weir duality, in C.E.M. Pearce and E. Hunt (Eds.), Otiization: Structures and Alications, Sringer, Dordrecht (2009), [13] V.I. Piterbarg, V.R. Fatalov, The Lalace ethod for robability easures on Banach saces, Russian Math. Surveys 50 (1995), [14] F. Wolfe, A duality theore in nonlinear rograing, Quart. Al. Math. 19 (1961), LAAS-CNRS and Institute of Matheatics, University of Toulouse, LAAS, 7 avenue du Colonel Roche, Toulouse Cédex 4,France E-ail address: lasserre@laas.fr Deto. Mateaticas, CINVESTAV-IPN, Ado. Postal , Mexico, D.F , Mexico E-ail address: eszeron@ath.cinvestav.x

Edinburgh Research Explorer

Edinburgh Research Explorer Edinburgh Research Exlorer ALMOST-ORTHOGONALITY IN THE SCHATTEN-VON NEUMANN CLASSES Citation for ublished version: Carbery, A 2009, 'ALMOST-ORTHOGONALITY IN THE SCHATTEN-VON NEUMANN CLASSES' Journal of

More information

Some simple continued fraction expansions for an in nite product Part 1. Peter Bala, January ax 4n+3 1 ax 4n+1. (a; x) =

Some simple continued fraction expansions for an in nite product Part 1. Peter Bala, January ax 4n+3 1 ax 4n+1. (a; x) = Soe sile continued fraction exansions for an in nite roduct Part. Introduction The in nite roduct Peter Bala, January 3 (a; x) = Y ax 4n+3 ax 4n+ converges for arbitrary colex a rovided jxj

More information

Approximation by Piecewise Constants on Convex Partitions

Approximation by Piecewise Constants on Convex Partitions Aroxiation by Piecewise Constants on Convex Partitions Oleg Davydov Noveber 4, 2011 Abstract We show that the saturation order of iecewise constant aroxiation in L nor on convex artitions with N cells

More information

PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY

PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 13,,. 8 14 M a t h e a t i c s ON BOUNDEDNESS OF A CLASS OF FIRST ORDER LINEAR DIFFERENTIAL OPERATORS IN THE SPACE OF n 1)-DIMENSIONALLY

More information

Ayşe Alaca, Şaban Alaca and Kenneth S. Williams School of Mathematics and Statistics, Carleton University, Ottawa, Ontario, Canada. Abstract.

Ayşe Alaca, Şaban Alaca and Kenneth S. Williams School of Mathematics and Statistics, Carleton University, Ottawa, Ontario, Canada. Abstract. Journal of Cobinatorics and Nuber Theory Volue 6, Nuber,. 17 15 ISSN: 194-5600 c Nova Science Publishers, Inc. DOUBLE GAUSS SUMS Ayşe Alaca, Şaban Alaca and Kenneth S. Willias School of Matheatics and

More information

Scaled Enflo type is equivalent to Rademacher type

Scaled Enflo type is equivalent to Rademacher type Scaled Enflo tye is equivalent to Radeacher tye Manor Mendel California Institute of Technology Assaf Naor Microsoft Research Abstract We introduce the notion of scaled Enflo tye of a etric sace, and show

More information

Real Analysis 1 Fall Homework 3. a n.

Real Analysis 1 Fall Homework 3. a n. eal Analysis Fall 06 Homework 3. Let and consider the measure sace N, P, µ, where µ is counting measure. That is, if N, then µ equals the number of elements in if is finite; µ = otherwise. One usually

More information

Asymptotic Behaviour Of Solutions For Some Weakly Dissipative Wave Equations Of p-laplacian Type

Asymptotic Behaviour Of Solutions For Some Weakly Dissipative Wave Equations Of p-laplacian Type Alied Matheatics E-Notes, (), 75 83 c IN 67-5 Available free at irror sites of htt://www.ath.nthu.edu.tw/aen/ Asytotic Behaviour Of olutions For oe Weakly Dissiative Wave Equations Of -Lalacian Tye Nour-Eddine

More information

Lecture 3: October 2, 2017

Lecture 3: October 2, 2017 Inforation and Coding Theory Autun 2017 Lecturer: Madhur Tulsiani Lecture 3: October 2, 2017 1 Shearer s lea and alications In the revious lecture, we saw the following stateent of Shearer s lea. Lea 1.1

More information

Numerical Method for Obtaining a Predictive Estimator for the Geometric Distribution

Numerical Method for Obtaining a Predictive Estimator for the Geometric Distribution British Journal of Matheatics & Couter Science 19(5): 1-13, 2016; Article no.bjmcs.29941 ISSN: 2231-0851 SCIENCEDOMAIN international www.sciencedoain.org Nuerical Method for Obtaining a Predictive Estiator

More information

Control and Stability of the Time-delay Linear Systems

Control and Stability of the Time-delay Linear Systems ISSN 746-7659, England, UK Journal of Inforation and Couting Science Vol., No. 4, 206,.29-297 Control and Stability of the Tie-delay Linear Systes Negras Tahasbi *, Hojjat Ahsani Tehrani Deartent of Matheatics,

More information

DISCRETE DUALITY FINITE VOLUME SCHEMES FOR LERAY-LIONS TYPE ELLIPTIC PROBLEMS ON GENERAL 2D MESHES

DISCRETE DUALITY FINITE VOLUME SCHEMES FOR LERAY-LIONS TYPE ELLIPTIC PROBLEMS ON GENERAL 2D MESHES ISCRETE UALITY FINITE VOLUME SCHEMES FOR LERAY-LIONS TYPE ELLIPTIC PROBLEMS ON GENERAL 2 MESHES BORIS ANREIANOV, FRANCK BOYER AN FLORENCE HUBERT Abstract. iscrete duality finite volue schees on general

More information

ON THE INTEGER PART OF A POSITIVE INTEGER S K-TH ROOT

ON THE INTEGER PART OF A POSITIVE INTEGER S K-TH ROOT ON THE INTEGER PART OF A POSITIVE INTEGER S K-TH ROOT Yang Hai Research Center for Basic Science, Xi an Jiaotong University, Xi an, Shaanxi, P.R.China Fu Ruiqin School of Science, Xi an Shiyou University,

More information

EXACT BOUNDS FOR JUDICIOUS PARTITIONS OF GRAPHS

EXACT BOUNDS FOR JUDICIOUS PARTITIONS OF GRAPHS EXACT BOUNDS FOR JUDICIOUS PARTITIONS OF GRAPHS B. BOLLOBÁS1,3 AND A.D. SCOTT,3 Abstract. Edwards showed that every grah of size 1 has a biartite subgrah of size at least / + /8 + 1/64 1/8. We show that

More information

Convex Analysis and Economic Theory Winter 2018

Convex Analysis and Economic Theory Winter 2018 Division of the Humanities and Social Sciences Ec 181 KC Border Conve Analysis and Economic Theory Winter 2018 Toic 16: Fenchel conjugates 16.1 Conjugate functions Recall from Proosition 14.1.1 that is

More information

Congruences involving Bernoulli and Euler numbers Zhi-Hong Sun

Congruences involving Bernoulli and Euler numbers Zhi-Hong Sun The aer will aear in Journal of Nuber Theory. Congruences involving Bernoulli Euler nubers Zhi-Hong Sun Deartent of Matheatics, Huaiyin Teachers College, Huaian, Jiangsu 300, PR China Received January

More information

Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems

Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various

More information

AN EXPLICIT METHOD FOR NUMERICAL SIMULATION OF WAVE EQUATIONS

AN EXPLICIT METHOD FOR NUMERICAL SIMULATION OF WAVE EQUATIONS The 4 th World Conference on Earthquake Engineering October -7, 8, Beiing, China AN EXPLICIT ETHOD FOR NUERICAL SIULATION OF WAVE EQUATIONS Liu Heng and Liao Zheneng Doctoral Candidate, Det. of Structural

More information

Parallelizing Spectrally Regularized Kernel Algorithms

Parallelizing Spectrally Regularized Kernel Algorithms Journal of Machine Learning Research 19 (2018) 1-29 Subitted 11/16; Revised 8/18; Published 8/18 Parallelizing Sectrally Regularized Kernel Algoriths Nicole Mücke nicole.uecke@atheatik.uni-stuttgart.de

More information

On Maximizing the Convergence Rate for Linear Systems With Input Saturation

On Maximizing the Convergence Rate for Linear Systems With Input Saturation IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 7, JULY 2003 1249 On Maxiizing the Convergence Rate for Linear Systes With Inut Saturation Tingshu Hu, Zongli Lin, Yacov Shaash Abstract In this note,

More information

1 Bounding the Margin

1 Bounding the Margin COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #12 Scribe: Jian Min Si March 14, 2013 1 Bounding the Margin We are continuing the proof of a bound on the generalization error of AdaBoost

More information

#A62 INTEGERS 16 (2016) REPRESENTATION OF INTEGERS BY TERNARY QUADRATIC FORMS: A GEOMETRIC APPROACH

#A62 INTEGERS 16 (2016) REPRESENTATION OF INTEGERS BY TERNARY QUADRATIC FORMS: A GEOMETRIC APPROACH #A6 INTEGERS 16 (016) REPRESENTATION OF INTEGERS BY TERNARY QUADRATIC FORMS: A GEOMETRIC APPROACH Gabriel Durha Deartent of Matheatics, University of Georgia, Athens, Georgia gjdurha@ugaedu Received: 9/11/15,

More information

Exploiting Matrix Symmetries and Physical Symmetries in Matrix Product States and Tensor Trains

Exploiting Matrix Symmetries and Physical Symmetries in Matrix Product States and Tensor Trains Exloiting Matrix Syetries and Physical Syetries in Matrix Product States and Tensor Trains Thoas K Huckle a and Konrad Waldherr a and Thoas Schulte-Herbrüggen b a Technische Universität München, Boltzannstr

More information

Quadratic Reciprocity. As in the previous notes, we consider the Legendre Symbol, defined by

Quadratic Reciprocity. As in the previous notes, we consider the Legendre Symbol, defined by Math 0 Sring 01 Quadratic Recirocity As in the revious notes we consider the Legendre Sybol defined by $ ˆa & 0 if a 1 if a is a quadratic residue odulo. % 1 if a is a quadratic non residue We also had

More information

MULTIPLIER IDEALS OF SUMS VIA CELLULAR RESOLUTIONS

MULTIPLIER IDEALS OF SUMS VIA CELLULAR RESOLUTIONS MULTIPLIER IDEALS OF SUMS VIA CELLULAR RESOLUTIONS SHIN-YAO JOW AND EZRA MILLER Abstract. Fix nonzero ideal sheaves a 1,..., a r and b on a noral Q-Gorenstein colex variety X. For any ositive real nubers

More information

Elementary Analysis in Q p

Elementary Analysis in Q p Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some

More information

Solutions 1. Introduction to Coding Theory - Spring 2010 Solutions 1. Exercise 1.1. See Examples 1.2 and 1.11 in the course notes.

Solutions 1. Introduction to Coding Theory - Spring 2010 Solutions 1. Exercise 1.1. See Examples 1.2 and 1.11 in the course notes. Solutions 1 Exercise 1.1. See Exaples 1.2 and 1.11 in the course notes. Exercise 1.2. Observe that the Haing distance of two vectors is the iniu nuber of bit flips required to transfor one into the other.

More information

Suppress Parameter Cross-talk for Elastic Full-waveform Inversion: Parameterization and Acquisition Geometry

Suppress Parameter Cross-talk for Elastic Full-waveform Inversion: Parameterization and Acquisition Geometry Suress Paraeter Cross-talk for Elastic Full-wavefor Inversion: Paraeterization and Acquisition Geoetry Wenyong Pan and Kris Innanen CREWES Project, Deartent of Geoscience, University of Calgary Suary Full-wavefor

More information

Cone-Constrained Linear Equations in Banach Spaces 1

Cone-Constrained Linear Equations in Banach Spaces 1 Journal of Convex Analysis Volume 4 (1997), No. 1, 149 164 Cone-Constrained Linear Equations in Banach Spaces 1 O. Hernandez-Lerma Departamento de Matemáticas, CINVESTAV-IPN, A. Postal 14-740, México D.F.

More information

Lecture 21. Interior Point Methods Setup and Algorithm

Lecture 21. Interior Point Methods Setup and Algorithm Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and

More information

5. Dimensional Analysis. 5.1 Dimensions and units

5. Dimensional Analysis. 5.1 Dimensions and units 5. Diensional Analysis In engineering the alication of fluid echanics in designs ake uch of the use of eirical results fro a lot of exerients. This data is often difficult to resent in a readable for.

More information

One- and multidimensional Fibonacci search very easy!

One- and multidimensional Fibonacci search very easy! One and ultidiensional ibonacci search One and ultidiensional ibonacci search very easy!. Content. Introduction / Preliinary rearks...page. Short descrition of the ibonacci nubers...page 3. Descrition

More information

Alireza Kamel Mirmostafaee

Alireza Kamel Mirmostafaee Bull. Korean Math. Soc. 47 (2010), No. 4, pp. 777 785 DOI 10.4134/BKMS.2010.47.4.777 STABILITY OF THE MONOMIAL FUNCTIONAL EQUATION IN QUASI NORMED SPACES Alireza Kael Mirostafaee Abstract. Let X be a linear

More information

Input-Output (I/O) Stability. -Stability of a System

Input-Output (I/O) Stability. -Stability of a System Inut-Outut (I/O) Stability -Stability of a Syste Outline: Introduction White Boxes and Black Boxes Inut-Outut Descrition Foralization of the Inut-Outut View Signals and Signal Saces he Notions of Gain

More information

ON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS

ON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS #A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,

More information

arxiv: v4 [math.st] 9 Aug 2017

arxiv: v4 [math.st] 9 Aug 2017 PARALLELIZING SPECTRAL ALGORITHMS FOR KERNEL LEARNING GILLES BLANCHARD AND NICOLE MÜCKE arxiv:161007487v4 [athst] 9 Aug 2017 Abstract We consider a distributed learning aroach in suervised learning for

More information

SHOUYU DU AND ZHANLE DU

SHOUYU DU AND ZHANLE DU THERE ARE INFINITELY MANY COUSIN PRIMES arxiv:ath/009v athgm 4 Oct 00 SHOUYU DU AND ZHANLE DU Abstract We roved that there are infinitely any cousin ries Introduction If c and c + 4 are both ries, then

More information

NONNEGATIVE matrix factorization finds its application

NONNEGATIVE matrix factorization finds its application Multilicative Udates for Convolutional NMF Under -Divergence Pedro J. Villasana T., Stanislaw Gorlow, Meber, IEEE and Arvind T. Hariraan arxiv:803.0559v2 [cs.lg 5 May 208 Abstract In this letter, we generalize

More information

Sums of independent random variables

Sums of independent random variables 3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where

More information

Elementary theory of L p spaces

Elementary theory of L p spaces CHAPTER 3 Elementary theory of L saces 3.1 Convexity. Jensen, Hölder, Minkowski inequality. We begin with two definitions. A set A R d is said to be convex if, for any x 0, x 1 2 A x = x 0 + (x 1 x 0 )

More information

Answers to Econ 210A Midterm, October A. The function f is homogeneous of degree 1/2. To see this, note that for all t > 0 and all (x 1, x 2 )

Answers to Econ 210A Midterm, October A. The function f is homogeneous of degree 1/2. To see this, note that for all t > 0 and all (x 1, x 2 ) Question. Answers to Econ 20A Midter, October 200 f(x, x 2 ) = ax {x, x 2 } A. The function f is hoogeneous of degree /2. To see this, note that for all t > 0 and all (x, x 2 ) f(tx, x 2 ) = ax {tx, tx

More information

Introduction to Banach Spaces

Introduction to Banach Spaces CHAPTER 8 Introduction to Banach Saces 1. Uniform and Absolute Convergence As a rearation we begin by reviewing some familiar roerties of Cauchy sequences and uniform limits in the setting of metric saces.

More information

On Wald-Type Optimal Stopping for Brownian Motion

On Wald-Type Optimal Stopping for Brownian Motion J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of

More information

ON THE NORM OF AN IDEMPOTENT SCHUR MULTIPLIER ON THE SCHATTEN CLASS

ON THE NORM OF AN IDEMPOTENT SCHUR MULTIPLIER ON THE SCHATTEN CLASS PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 000-9939XX)0000-0 ON THE NORM OF AN IDEMPOTENT SCHUR MULTIPLIER ON THE SCHATTEN CLASS WILLIAM D. BANKS AND ASMA HARCHARRAS

More information

Modi ed Local Whittle Estimator for Long Memory Processes in the Presence of Low Frequency (and Other) Contaminations

Modi ed Local Whittle Estimator for Long Memory Processes in the Presence of Low Frequency (and Other) Contaminations Modi ed Local Whittle Estiator for Long Meory Processes in the Presence of Low Frequency (and Other Containations Jie Hou y Boston University Pierre Perron z Boston University March 5, 203; Revised: January

More information

The Generalized Integer Gamma DistributionA Basis for Distributions in Multivariate Statistics

The Generalized Integer Gamma DistributionA Basis for Distributions in Multivariate Statistics Journal of Multivariate Analysis 64, 8610 (1998) Article No. MV971710 The Generalized Inteer Gaa DistributionA Basis for Distributions in Multivariate Statistics Carlos A. Coelho Universidade Te cnica

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Notes for EE7C (Spring 018: Convex Optiization and Approxiation Instructor: Moritz Hardt Eail: hardt+ee7c@berkeley.edu Graduate Instructor: Max Sichowitz Eail: sichow+ee7c@berkeley.edu October 15,

More information

3.8. External Source in Quantum-Statistical Path Integral

3.8. External Source in Quantum-Statistical Path Integral 3.8. External Source in Quantu-Statistical Path Integral This section studies the quantu statistical evolution alitude ( x b, x a ) ω x ex - ħ d M x + ω x - () x (3.) which will be evaluated in two different

More information

Minimizing Machinery Vibration Transmission in a Lightweight Building using Topology Optimization

Minimizing Machinery Vibration Transmission in a Lightweight Building using Topology Optimization 1 th World Congress on Structural and Multidiscilinary Otiization May 19-4, 13, Orlando, Florida, USA Miniizing Machinery Vibration ransission in a Lightweight Building using oology Otiization Niels Olhoff,

More information

On Doob s Maximal Inequality for Brownian Motion

On Doob s Maximal Inequality for Brownian Motion Stochastic Process. Al. Vol. 69, No., 997, (-5) Research Reort No. 337, 995, Det. Theoret. Statist. Aarhus On Doob s Maximal Inequality for Brownian Motion S. E. GRAVERSEN and G. PESKIR If B = (B t ) t

More information

Convex Optimization methods for Computing Channel Capacity

Convex Optimization methods for Computing Channel Capacity Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem

More information

PREPRINT 2006:17. Inequalities of the Brunn-Minkowski Type for Gaussian Measures CHRISTER BORELL

PREPRINT 2006:17. Inequalities of the Brunn-Minkowski Type for Gaussian Measures CHRISTER BORELL PREPRINT 2006:7 Inequalities of the Brunn-Minkowski Type for Gaussian Measures CHRISTER BORELL Departent of Matheatical Sciences Division of Matheatics CHALMERS UNIVERSITY OF TECHNOLOGY GÖTEBORG UNIVERSITY

More information

TRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES

TRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES Khayyam J. Math. DOI:10.22034/kjm.2019.84207 TRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES ISMAEL GARCÍA-BAYONA Communicated by A.M. Peralta Abstract. In this aer, we define two new Schur and

More information

Analysis of low rank matrix recovery via Mendelson s small ball method

Analysis of low rank matrix recovery via Mendelson s small ball method Analysis of low rank atrix recovery via Mendelson s sall ball ethod Maryia Kabanava Chair for Matheatics C (Analysis) ontdriesch 0 kabanava@athc.rwth-aachen.de Holger Rauhut Chair for Matheatics C (Analysis)

More information

arxiv: v1 [math.ds] 19 Jun 2012

arxiv: v1 [math.ds] 19 Jun 2012 Rates in the strong invariance rincile for ergodic autoorhiss of the torus Jérôe Dedecker a, Florence Merlevède b and Françoise Pène c 1 a Université Paris Descartes, Sorbonne Paris Cité, Laboratoire MAP5

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I Contents 1. Preliinaries 2. The ain result 3. The Rieann integral 4. The integral of a nonnegative

More information

The Methods of Solution for Constrained Nonlinear Programming

The Methods of Solution for Constrained Nonlinear Programming Research Inventy: International Journal Of Engineering And Science Vol.4, Issue 3(March 2014), PP 01-06 Issn (e): 2278-4721, Issn (p):2319-6483, www.researchinventy.co The Methods of Solution for Constrained

More information

Handout 6 Solutions to Problems from Homework 2

Handout 6 Solutions to Problems from Homework 2 CS 85/185 Fall 2003 Lower Bounds Handout 6 Solutions to Probles fro Hoewor 2 Ait Charabarti Couter Science Dartouth College Solution to Proble 1 1.2: Let f n stand for A 111 n. To decide the roerty f 3

More information

On Isoperimetric Functions of Probability Measures Having Log-Concave Densities with Respect to the Standard Normal Law

On Isoperimetric Functions of Probability Measures Having Log-Concave Densities with Respect to the Standard Normal Law On Isoerimetric Functions of Probability Measures Having Log-Concave Densities with Resect to the Standard Normal Law Sergey G. Bobkov Abstract Isoerimetric inequalities are discussed for one-dimensional

More information

Optimal Adaptive Computations in the Jaffard Algebra and Localized Frames

Optimal Adaptive Computations in the Jaffard Algebra and Localized Frames www.oeaw.ac.at Otial Adative Coutations in the Jaffard Algebra and Localized Fraes M. Fornasier, K. Gröchenig RICAM-Reort 2006-28 www.rica.oeaw.ac.at Otial Adative Coutations in the Jaffard Algebra and

More information

Metric Cotype. 1 Introduction. Manor Mendel California Institute of Technology. Assaf Naor Microsoft Research

Metric Cotype. 1 Introduction. Manor Mendel California Institute of Technology. Assaf Naor Microsoft Research Metric Cotye Manor Mendel California Institute of Technology Assaf Naor Microsoft Research Abstract We introduce the notion of cotye of a etric sace, and rove that for Banach saces it coincides with the

More information

Sharp gradient estimate and spectral rigidity for p-laplacian

Sharp gradient estimate and spectral rigidity for p-laplacian Shar gradient estimate and sectral rigidity for -Lalacian Chiung-Jue Anna Sung and Jiaing Wang To aear in ath. Research Letters. Abstract We derive a shar gradient estimate for ositive eigenfunctions of

More information

Chordal Sparsity, Decomposing SDPs and the Lyapunov Equation

Chordal Sparsity, Decomposing SDPs and the Lyapunov Equation Chordal Sarsity, Decoosing SDPs and the Lyaunov Equation Richard P. Mason and Antonis Paachristodoulou Abstract Analysis questions in control theory are often forulated as Linear Matrix Inequalities and

More information

1. Introduction In this note we prove the following result which seems to have been informally conjectured by Semmes [Sem01, p. 17].

1. Introduction In this note we prove the following result which seems to have been informally conjectured by Semmes [Sem01, p. 17]. A REMARK ON POINCARÉ INEQUALITIES ON METRIC MEASURE SPACES STEPHEN KEITH AND KAI RAJALA Abstract. We show that, in a comlete metric measure sace equied with a doubling Borel regular measure, the Poincaré

More information

ON THE SET a x + b g x (mod p) 1 Introduction

ON THE SET a x + b g x (mod p) 1 Introduction PORTUGALIAE MATHEMATICA Vol 59 Fasc 00 Nova Série ON THE SET a x + b g x (mod ) Cristian Cobeli, Marian Vâjâitu and Alexandru Zaharescu Abstract: Given nonzero integers a, b we rove an asymtotic result

More information

RESOLVENT ESTIMATES FOR ELLIPTIC SYSTEMS IN FUNCTION SPACES OF HIGHER REGULARITY

RESOLVENT ESTIMATES FOR ELLIPTIC SYSTEMS IN FUNCTION SPACES OF HIGHER REGULARITY Electronic Journal of Differential Equations, Vol. 2011 2011, No. 109,. 1 12. ISSN: 1072-6691. URL: htt://ejde.ath.txstate.edu or htt://ejde.ath.unt.edu ft ejde.ath.txstate.edu RESOLVENT ESTIMATES FOR

More information

B8.1 Martingales Through Measure Theory. Concept of independence

B8.1 Martingales Through Measure Theory. Concept of independence B8.1 Martingales Through Measure Theory Concet of indeendence Motivated by the notion of indeendent events in relims robability, we have generalized the concet of indeendence to families of σ-algebras.

More information

Sufficient and necessary conditions for generalized Ho lder s inequality in p-summable sequence spaces

Sufficient and necessary conditions for generalized Ho lder s inequality in p-summable sequence spaces Sufficient and necessary conditions for generalized Ho lder s inequality in -suable sequence saces l Masta*, S Fatiah 2, rsisari 3, Y Y Putra 4 and F riani 5,2 Deartent of Matheatics Education, Universitas

More information

Monoidal categories for the combinatorics of group representations

Monoidal categories for the combinatorics of group representations Monoidal categories for the cobinatorics of grou reresentations Ross Street February 1998 Categories are known to be useful for organiing structural asects of atheatics. However, they are also useful in

More information

3.8 Three Types of Convergence

3.8 Three Types of Convergence 3.8 Three Types of Convergence 3.8 Three Types of Convergence 93 Suppose that we are given a sequence functions {f k } k N on a set X and another function f on X. What does it ean for f k to converge to

More information

An Estimate For Heilbronn s Exponential Sum

An Estimate For Heilbronn s Exponential Sum An Estimate For Heilbronn s Exonential Sum D.R. Heath-Brown Magdalen College, Oxford For Heini Halberstam, on his retirement Let be a rime, and set e(x) = ex(2πix). Heilbronn s exonential sum is defined

More information

Sobolev Spaces with Weights in Domains and Boundary Value Problems for Degenerate Elliptic Equations

Sobolev Spaces with Weights in Domains and Boundary Value Problems for Degenerate Elliptic Equations Sobolev Saces with Weights in Domains and Boundary Value Problems for Degenerate Ellitic Equations S. V. Lototsky Deartment of Mathematics, M.I.T., Room 2-267, 77 Massachusetts Avenue, Cambridge, MA 02139-4307,

More information

Uniform Deviation Bounds for k-means Clustering

Uniform Deviation Bounds for k-means Clustering Unifor Deviation Bounds for k-means Clustering Olivier Bache Mario Lucic S. Haed Hassani Andreas Krause Abstract Unifor deviation bounds liit the difference between a odel s exected loss and its loss on

More information

ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE

ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE CHRISTOPHER J. HILLAR Abstract. A long-standing conjecture asserts that the polynoial p(t = Tr(A + tb ] has nonnegative coefficients whenever is

More information

4 A Survey of Congruent Results 12

4 A Survey of Congruent Results 12 4 A urvey of Congruent Results 1 ECTION 4.5 Perfect Nubers and the iga Function By the end of this section you will be able to test whether a given Mersenne nuber is rie understand what is eant be a erfect

More information

A GENERAL THEORY OF PARTICLE FILTERS IN HIDDEN MARKOV MODELS AND SOME APPLICATIONS. By Hock Peng Chan National University of Singapore and

A GENERAL THEORY OF PARTICLE FILTERS IN HIDDEN MARKOV MODELS AND SOME APPLICATIONS. By Hock Peng Chan National University of Singapore and Subitted to the Annals of Statistics A GENERAL THEORY OF PARTICLE FILTERS IN HIDDEN MARKOV MODELS AND SOME APPLICATIONS By Hock Peng Chan National University of Singaore and By Tze Leung Lai Stanford University

More information

On Conditions for Linearity of Optimal Estimation

On Conditions for Linearity of Optimal Estimation On Conditions for Linearity of Optial Estiation Erah Akyol, Kuar Viswanatha and Kenneth Rose {eakyol, kuar, rose}@ece.ucsb.edu Departent of Electrical and Coputer Engineering University of California at

More information

HEAT AND LAPLACE TYPE EQUATIONS WITH COMPLEX SPATIAL VARIABLES IN WEIGHTED BERGMAN SPACES

HEAT AND LAPLACE TYPE EQUATIONS WITH COMPLEX SPATIAL VARIABLES IN WEIGHTED BERGMAN SPACES Electronic Journal of ifferential Equations, Vol. 207 (207), No. 236,. 8. ISSN: 072-669. URL: htt://ejde.math.txstate.edu or htt://ejde.math.unt.edu HEAT AN LAPLACE TYPE EQUATIONS WITH COMPLEX SPATIAL

More information

Products of Composition, Multiplication and Differentiation between Hardy Spaces and Weighted Growth Spaces of the Upper-Half Plane

Products of Composition, Multiplication and Differentiation between Hardy Spaces and Weighted Growth Spaces of the Upper-Half Plane Global Journal of Pure and Alied Mathematics. ISSN 0973-768 Volume 3, Number 9 (207),. 6303-636 Research India Publications htt://www.riublication.com Products of Comosition, Multilication and Differentiation

More information

Uniform Deviation Bounds for k-means Clustering

Uniform Deviation Bounds for k-means Clustering Unifor Deviation Bounds for k-means Clustering Olivier Bache Mario Lucic S Haed Hassani Andreas Krause Abstract Unifor deviation bounds liit the difference between a odel s exected loss and its loss on

More information

STRONG TYPE INEQUALITIES AND AN ALMOST-ORTHOGONALITY PRINCIPLE FOR FAMILIES OF MAXIMAL OPERATORS ALONG DIRECTIONS IN R 2

STRONG TYPE INEQUALITIES AND AN ALMOST-ORTHOGONALITY PRINCIPLE FOR FAMILIES OF MAXIMAL OPERATORS ALONG DIRECTIONS IN R 2 STRONG TYPE INEQUALITIES AND AN ALMOST-ORTHOGONALITY PRINCIPLE FOR FAMILIES OF MAXIMAL OPERATORS ALONG DIRECTIONS IN R 2 ANGELES ALFONSECA Abstract In this aer we rove an almost-orthogonality rincile for

More information

A := A i : {A i } S. is an algebra. The same object is obtained when the union in required to be disjoint.

A := A i : {A i } S. is an algebra. The same object is obtained when the union in required to be disjoint. 59 6. ABSTRACT MEASURE THEORY Having developed the Lebesgue integral with respect to the general easures, we now have a general concept with few specific exaples to actually test it on. Indeed, so far

More information

Some results of convex programming complexity

Some results of convex programming complexity 2012c12 $ Ê Æ Æ 116ò 14Ï Dec., 2012 Oerations Research Transactions Vol.16 No.4 Some results of convex rogramming comlexity LOU Ye 1,2 GAO Yuetian 1 Abstract Recently a number of aers were written that

More information

ON JOINT CONVEXITY AND CONCAVITY OF SOME KNOWN TRACE FUNCTIONS

ON JOINT CONVEXITY AND CONCAVITY OF SOME KNOWN TRACE FUNCTIONS ON JOINT CONVEXITY ND CONCVITY OF SOME KNOWN TRCE FUNCTIONS MOHMMD GHER GHEMI, NHID GHRKHNLU and YOEL JE CHO Communicated by Dan Timotin In this aer, we rovide a new and simle roof for joint convexity

More information

Small Zeros of Quadratic Forms Mod P m

Small Zeros of Quadratic Forms Mod P m International Mathematical Forum, Vol. 8, 2013, no. 8, 357-367 Small Zeros of Quadratic Forms Mod P m Ali H. Hakami Deartment of Mathematics, Faculty of Science, Jazan University P.O. Box 277, Jazan, Postal

More information

[95/95] APPROACH FOR DESIGN LIMITS ANALYSIS IN VVER. Shishkov L., Tsyganov S. Russian Research Centre Kurchatov Institute Russian Federation, Moscow

[95/95] APPROACH FOR DESIGN LIMITS ANALYSIS IN VVER. Shishkov L., Tsyganov S. Russian Research Centre Kurchatov Institute Russian Federation, Moscow [95/95] APPROACH FOR DESIGN LIMITS ANALYSIS IN VVER Shishkov L., Tsyganov S. Russian Research Centre Kurchatov Institute Russian Federation, Moscow ABSTRACT The aer discusses a well-known condition [95%/95%],

More information

A PRIORI ESTIMATES AND APPLICATION TO THE SYMMETRY OF SOLUTIONS FOR CRITICAL

A PRIORI ESTIMATES AND APPLICATION TO THE SYMMETRY OF SOLUTIONS FOR CRITICAL A PRIORI ESTIMATES AND APPLICATION TO THE SYMMETRY OF SOLUTIONS FOR CRITICAL LAPLACE EQUATIONS Abstract. We establish ointwise a riori estimates for solutions in D 1, of equations of tye u = f x, u, where

More information

New Set of Rotationally Legendre Moment Invariants

New Set of Rotationally Legendre Moment Invariants New Set of Rotationally Legendre Moent Invariants Khalid M. Hosny Abstract Orthogonal Legendre oents are used in several attern recognition and iage rocessing alications. Translation and scale Legendre

More information

Notes on duality in second order and -order cone otimization E. D. Andersen Λ, C. Roos y, and T. Terlaky z Aril 6, 000 Abstract Recently, the so-calle

Notes on duality in second order and -order cone otimization E. D. Andersen Λ, C. Roos y, and T. Terlaky z Aril 6, 000 Abstract Recently, the so-calle McMaster University Advanced Otimization Laboratory Title: Notes on duality in second order and -order cone otimization Author: Erling D. Andersen, Cornelis Roos and Tamás Terlaky AdvOl-Reort No. 000/8

More information

Mistiming Performance Analysis of the Energy Detection Based ToA Estimator for MB-OFDM

Mistiming Performance Analysis of the Energy Detection Based ToA Estimator for MB-OFDM IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS Mistiing Perforance Analysis of the Energy Detection Based ToA Estiator for MB-OFDM Huilin Xu, Liuqing Yang contact author, Y T Jade Morton and Mikel M Miller

More information

Introduction to Discrete Optimization

Introduction to Discrete Optimization Prof. Friedrich Eisenbrand Martin Nieeier Due Date: March 9 9 Discussions: March 9 Introduction to Discrete Optiization Spring 9 s Exercise Consider a school district with I neighborhoods J schools and

More information

THE WEIGHTING METHOD AND MULTIOBJECTIVE PROGRAMMING UNDER NEW CONCEPTS OF GENERALIZED (, )-INVEXITY

THE WEIGHTING METHOD AND MULTIOBJECTIVE PROGRAMMING UNDER NEW CONCEPTS OF GENERALIZED (, )-INVEXITY U.P.B. Sci. Bull., Series A, Vol. 80, Iss. 2, 2018 ISSN 1223-7027 THE WEIGHTING METHOD AND MULTIOBJECTIVE PROGRAMMING UNDER NEW CONCEPTS OF GENERALIZED (, )-INVEXITY Tadeusz ANTCZAK 1, Manuel ARANA-JIMÉNEZ

More information

The Semantics of Data Flow Diagrams. P.D. Bruza. Th.P. van der Weide. Dept. of Information Systems, University of Nijmegen

The Semantics of Data Flow Diagrams. P.D. Bruza. Th.P. van der Weide. Dept. of Information Systems, University of Nijmegen The Seantics of Data Flow Diagras P.D. Bruza Th.P. van der Weide Det. of Inforation Systes, University of Nijegen Toernooiveld, NL-6525 ED Nijegen, The Netherlands July 26, 1993 Abstract In this article

More information

New upper bound for the B-spline basis condition number II. K. Scherer. Institut fur Angewandte Mathematik, Universitat Bonn, Bonn, Germany.

New upper bound for the B-spline basis condition number II. K. Scherer. Institut fur Angewandte Mathematik, Universitat Bonn, Bonn, Germany. New upper bound for the B-spline basis condition nuber II. A proof of de Boor's 2 -conjecture K. Scherer Institut fur Angewandte Matheati, Universitat Bonn, 535 Bonn, Gerany and A. Yu. Shadrin Coputing

More information

On Z p -norms of random vectors

On Z p -norms of random vectors On Z -norms of random vectors Rafa l Lata la Abstract To any n-dimensional random vector X we may associate its L -centroid body Z X and the corresonding norm. We formulate a conjecture concerning the

More information

GENERALIZED NORMS INEQUALITIES FOR ABSOLUTE VALUE OPERATORS

GENERALIZED NORMS INEQUALITIES FOR ABSOLUTE VALUE OPERATORS International Journal of Analysis Alications ISSN 9-8639 Volume 5, Number (04), -9 htt://www.etamaths.com GENERALIZED NORMS INEQUALITIES FOR ABSOLUTE VALUE OPERATORS ILYAS ALI, HU YANG, ABDUL SHAKOOR Abstract.

More information

Design of Linear-Phase Two-Channel FIR Filter Banks with Rational Sampling Factors

Design of Linear-Phase Two-Channel FIR Filter Banks with Rational Sampling Factors R. Bregović and. Saraäi, Design of linear hase two-channel FIR filter bans with rational saling factors, Proc. 3 rd Int. Sy. on Iage and Signal Processing and Analysis, Roe, Italy, Set. 3,. 749 754. Design

More information

arxiv:math/ v1 [math.fa] 5 Dec 2003

arxiv:math/ v1 [math.fa] 5 Dec 2003 arxiv:math/0323v [math.fa] 5 Dec 2003 WEAK CLUSTER POINTS OF A SEQUENCE AND COVERINGS BY CYLINDERS VLADIMIR KADETS Abstract. Let H be a Hilbert sace. Using Ball s solution of the comlex lank roblem we

More information

On the minimax inequality and its application to existence of three solutions for elliptic equations with Dirichlet boundary condition

On the minimax inequality and its application to existence of three solutions for elliptic equations with Dirichlet boundary condition ISSN 1 746-7233 England UK World Journal of Modelling and Simulation Vol. 3 (2007) No. 2. 83-89 On the minimax inequality and its alication to existence of three solutions for ellitic equations with Dirichlet

More information

SUPPORTING INFORMATION FOR. Mass Spectrometrically-Detected Statistical Aspects of Ligand Populations in Mixed Monolayer Au 25 L 18 Nanoparticles

SUPPORTING INFORMATION FOR. Mass Spectrometrically-Detected Statistical Aspects of Ligand Populations in Mixed Monolayer Au 25 L 18 Nanoparticles SUPPORTIG IFORMATIO FOR Mass Sectroetrically-Detected Statistical Asects of Lig Poulations in Mixed Monolayer Au 25 L 8 anoarticles Aala Dass,,a Kennedy Holt, Joseh F. Parer, Stehen W. Feldberg, Royce

More information