On lower bounds for integration of multivariate permutation-invariant functions

Size: px
Start display at page:

Download "On lower bounds for integration of multivariate permutation-invariant functions"

Transcription

1 arxiv: v1 [math.na] 15 Oct 2013 On lower bouns for integration of multivariate permutation-invariant functions Markus Weimar October 16, 2013 Abstract In this note we stuy multivariate integration for permutation-invariant functions fromacertain Banach spacee,α of Korobov typeinthe worst case setting. We present a lower error boun which particularly implies that in imension every cubature rule which reuces the initial error necessarily uses at least + 1 function values. Since this hols inepenently of the number of permutation-invariant coorinates, this shows that the integration problem can never be strongly polynomially tractable in this setting. Our assertions generalize results ue to Sloan an Woźniakowski [SW97]. Moreover, for large smoothness parameters α our boun can not be improve. Finally, we exten our results to the case of permutation-invariant functions from Korobov-type spaces equippe with prouct weights. 1 Introuction an main result Consier the integration problem Int = (Int ) N, Int : E,α C, Int (f) = f(x)x, [0,1] for perioic, complex-value functions in the Korobov class { E,α := f L 1 ([0,1] ) f := f E,α := sup f(k) ( } ) α k 1... k < k Z Philipps-University Marburg, Faculty of Mathematics an Computer Science, Hans-Meerwein-Straße, Lahnberge, Marburg, Germany. weimar@mathematik.uni-marburg.e This work has been supporte by Deutsche Forschungsgemeinschaft DFG (DA 360/19-1). 1

2 where N an α > 1. Here Z enotes the set of integers, N := {1,2,...}, an we set k m := max{1, k m }. Moreover, for f L 1 ([0,1] ) f(k) := f,e 2πik L 2 := f(x)e 2πikx x, k = (k 1,...,k ) Z, [0,1] enotes its kth Fourier coefficient, where kx = m=1 k m x m, an i = 1. To approximate Int (f), withoutloss ofgenerality, weconsier algorithmsfromtheclassofalllinear cubature rules A(f) := A N, (f) := N w n f ( t (n)), N N 0 := {0} N, (1) n=1 that use at most N values of the input function f at some points t (n) [0,1], n = 1...,N. The weights w n can be arbitrary complex numbers. Clearly, every function f E,α has a 1-perioic extension since their Fourier series are absolutely convergent: 2πik ( ) α f(k)e f k1... k = f (1+2ζ(α)) <. k Z k Z As usual, ζ(s) = m=1 m s is the Riemann zeta function evaluate at s > 1. In [SW97] Sloan an Woźniakowski showe that for every N the Nth minimal worst case error of Int = (Int ) N, e(n,;int,e,α ) := inf A N, sup f E,α 1 Int (f) A N, (f), equals the initial error e(0,;int,e,α ) = 1 provie that N < 2. In other wors, the integration problem on the full spaces (E,α ) N suffers from the curse of imensionality, since for every fixe ε (0,1) its information complexity grows exponentially with the imension : n(ε,) := n(ε,;int,e,α ) := min{n N 0 e(n,;int,e,α ) ε} 2, N. We generalize this result to the case of permutation-invariant 1 subspaces in the sense of [Wei12]. To this en, for N let I {1,...,} be some subset of coorinates an consier the integration problem Int = (Int ) N restricte to the subspace S I (E,α ) of all 1 In [Wei12] we use the name symmetric what cause some confusion. 2

3 I -permutation-invariant functions f E,α. That is, in imension we restrict ourselves to functions f that satisfy f(x) = f(σ(x)) for all x [0,1] (2) an any permutation σ from { S I := σ: {1,...,} {1,...,} σ bijective an σ } {1,...,}\I = i (3) that leaves the elements in the complement of I fixe. For the ease of presentation we shall use the same notation for permutations σ S I an for the corresponing permutations σ : R R of -imensional vectors, given by x = (x 1,...,x ) σ (x) := (x σ(1),...,x σ() ). Observe that in the case I = we clearly have S I (E,α ) = E,α. One motivation to stuy the integration problem restricte to those subspaces is relate to approximate solutions of partial ifferential equations. Many approaches to obtain such solutions lea us to the problem of calculating high-imensional integrals, e.g., to obtain certain wavelet coefficients. Obviously, it is of interest whether this can be one efficiently since taking into account a large number of coefficients woul lea to better approximations to the exact solution. Therefore it is important to incorporate as many structural properties (such as permutation-invariance of the integrans uner consieration) as possible in orer to reuce the effort for every single calculation. In information-base complexity (IBC) this effort is measure by the behavior of the information complexity n(ε,) which can be formalize by several notions of tractability. We remin the reaer that a problem is calle polynomially tractable if its information complexity n(ε,) is boune from above by some polynomial in an ε 1, i.e., n(ε,) Cε p q for some C,p > 0, q 0 an all ε (0,1], N. If the last inequality remains vali even for q = 0 then we say that the problem is strongly polynomially tractable. Apart from that, several weaker notions of tractability were introuce recently; for etails see, e.g., Sielecki [Sie13]. If for some fixe s, t (0, 1] the information complexity satisfies lnn(ε,) lim = 0, (4) ε 1 + ε s + t then we have (s, t)-weak tractability. This generalizes the well-establishe notion of weakly tractable problems (which is inclue as the special case s = t = 1). It is use to escribe the 3

4 case of at most subexponentially growing information complexities such as, e.g., n(ε, ) = exp( ). Finally, a given problem is sai to be uniformly weakly tractable if the limit conition (4) hols 2 for every s,t (0,1]. Our main result states that the integration problem for permutation-invariant functions in the above sense can never be strongly polynomially tractable, inepenent of the size of the sets I which escribes the number of impose permutation-invariance conitions. The assertion reas as follows: Theorem 1. Let Then, for every N < N, an N := N (,I ) := (#I +1) 2 #I, N. (5) e(n,;int,s I (E,α )) = 1 (6) e(n,;int,s I (E,α )) for all N an α > 1. Consequently, for all N. lim α e(n,;int,s I (E,α )) = 0 ( 1+ ζ(α) ) 1 (7) 2 α 1 Before we give the proof of this theorem in Section 3 we a some comments on this result in the next section. 2 Discussion an further results Note that, in particular, (6) yiels that for every imension the initial error of the integration problem uner consieration equals 1. Thus the problem is well-scale an we o not nee to istinguish between the absolute an the normalize error criteria. We stress that ue to results of Smolyak an Bakhvalov the choice of linear, nonaaptive cubature rules A in (1) is inee without loss of generality. For etails an further references see, e.g., [SW97, Remark 1] an Novak an Woźniakowski [NW08, Section 4.2.2]. 2 Note that it clearly suffices to check (4) for every s = t (0,1] in orer to conclue uniform weak tractability. 4

5 Moreover, observe that ue to (7) the assertion state in (6) can not be extene to N N provie that the smoothness α is sufficiently large. That means that at least in this case our result is sharp. Note that (6) can be reformulate equivalently in terms of the information complexity: n(ε,;int,s I (E,α )) (#I +1) 2 #I for all N an every ε (0,1). Therefore our Theorem 1 inee generalizes [SW97, Theorem 1] since the latter is containe as the special case I = for every N. Furthermore, observe that in any case the right-han sie of the last inequality is lower boune by +1. Hence, even for the fully permutation-invariant problem, where we have I = {1,...,} for all N, the information complexity grows at least linearly with imension provie that ε < 1. Together with some obvious estimates this proves the following corollary. Corollary 1. Assume α > 1 an let b := #I. We stuy Int = (Int ) N for the sequence of I -permutation-invariant subspaces (S I (E,α )) N in the worst case setting: If the problem is polynomially tractable with the constants C, p, q then q 1 an (b ) N O(ln). In particular, the problem is never strongly polynomially tractable. If the problem is uniformly weakly tractable then (b ) N o( t ) for all t (0,1]. If the problem is (s,t)-weakly tractable for some s,t (0,1] then (b ) N o( t ). In particular, weak tractability implies (b ) N o(). If (b ) N / o() then we have the curse of imensionality. In turn, alreay the absence of the curse implies (b ) N o(). For some applications it might be useful to impose permutation-invariance conitions with respect to finitely many isjoint subsets I (1),...,I(R) {1,...,} of the coorinates, R > 1; see [Wei12] for etails. In this case the respective subspace S (1) I,...,I(R) (E,α ) := R r=1 S (r) I (E,α ) E,α consists of all f E,α which satisfy (2) for all σ R r=1 S, where S I (r) (r) I an (E,α ) are efine as before. It turns out that Theorem 1 remains vali for this S I (r) 5

6 case when we replace the efinition of N given in (5) by ( ( ) R ( ) ) 2 N := N,I (1),...,I(R) := #I (r) +1 R r=1 #I(r). Without going into etails we mention that our proof given in Section 3 below can be transferre almost literally to this case. Consequently, an analogue of Corollary 1 remains vali if we set b := R r=1 #I(r) for N. During the last two ecaes many numerical problems such as (high-imensional) integration were proven to be computationally har. Fortunately, it turne out that the introuction of weights to the norm of the unerlying source spaces can ramatically reuce the complexity of those problems such that they can be hanle efficiently; see, e.g., [NW08]. In [Woź09] Woźniakowski consiere the integration problem efine above for weighte Korobov-type spaces E γ,α := { f L 1 ([0,1] ) f γ := f E γ r=1 where N an α > 1, as well as γ := {γ,u(k) γ,u(k) := m u(k) γ,m are prouct weights with,α ( ) α } k1... k := sup f(k) <, k Z γ,u(k) (8) k Z }. Here the quantities u(k) := {m {1,...,} k m 0} an 1 = γ, γ,1... γ, 0 for N an k = (k 1,...,k ) Z. This generalizes [SW97] since if γ,m 1 for all m = 1,..., an N then E γ,α = E,α. Theorem 1 in [Woź09] states that η,n e(n,;int,e γ,α ) 1 for all α > 1, N, an N < 2, where η,n enotes the (N+1)st largest weight in the sequence (γ,u ) u {1,...,} [0,1]. So what about weights in our setting? First observe that the subset of I -permutation-invariant functions in E γ,α again forms a linear subspace. In accorance with our previous notation we enote it by S I (E γ,α ). Furthermore, we note that the Fourier coefficients of an I -permutation-invariant function f are I -permutation-invariant again, i.e., they satisfy f(k) = f(σ(k)) for every k Z an all σ S I. Since clearly also ( ) α k 1... k is I -permutation-invariant, f E γ,α belongs to the unit ball of the subspace S I (E γ,α ) if an only if γ,u(σ(k)) µ,u(k) f(k) min ( ) α =: ( ) α, k Z, k1... k k1... k σ S I 6

7 where the new weights µ,u are efine as a prouct of an orer-epenent part (w.r.t. coorinates from I ) an a part consisting of usual prouct weights. For the ease of notation, in what follows we assume permutation-invariance w.r.t. the first #I coorinates, i.e., I = {1,2,...,#I }. Then µ,u(k) = min σ S I m u I (σ(k)) γ,m m u I c (k) γ,m = #u I (k) i=1 γ,#i i+1 m u I c (k) γ,m. (9) Here we set u I (k) := u(k) I an u I c (k) := u(k)\i to enote the support of k Z w.r.t. the sets I an I c = {1,...,}\I, respectively. Thus we have S I (E γ,α ) = S I (E µ,α ), where E µ,α is efine by (8) with γ,u(k) replace by µ,u(k) an S I once more enotes the restriction to I -permutation-invariant functions. Again the lower boun for the error of numerical integration epens on the (N+1)th largest weight µ,u(k) that can appear. To formalize this point let ψ: {0,1,...,# 1} enote a rearrangement of the multi-inices k from the set := { k = (m,l) {0,1} #I {0,1} ( #I ) m 1... m #I } {0,1} Z (10) such that the corresponing weights µ,u(k) possess a non-increasing orering, i.e. ν,n := µ,u(ψ(n)) µ,u(ψ(n+1)) = ν,n+1 for all n = 0,1,...,# 2. (11) Then the weighte analogue of Theorem 1 reas as follows: Theorem 2. Assume α > 1, for every N let I {1,...,} be given, an consier N (,I ) efine as in (5). Then we have for all N an every N < N (,I ). We illustrate Theorem 2 by two examples. ν,n e ( N,;Int,S I (E µ,α )) 1 (12) 7

8 Example 1. Consier the fully permutation-invariant problem, i.e., assume I = {1,...,} for every N. For simplicity let γ,m := 1 m for every m {1,...,} an all N. Then = {k {0,1} k 1... k } an # = + 1. Moreover, we conclue µ,u(k) = #u(k) i=1 1/( i+1), k, an consequently ν,0 = 1, ν,1 = 1/, ν,2 = 1/( ( 1)),. ν, = 1/(!). Hence, for α > 1 an N we obtain e(0,;int ;S I (E µ,α )) = 1 an e(n,;int ;S I (E µ,α )) 1 ( 1)... ( N +1), if N {1,...,}. Our secon, more sophisticate example shows that permutation-invariance is not as powerful as aitional knowlege moele by weights. It generalizes an assertion ue to Woźniakowski [Woź09, p. 648]. Example 2. For N let #I {0,1,...,} be arbitrary an assume permutationinvariance w.r.t. the first #I coorinates in imension. Furthermore, let γ,1 :=... := γ,#i := 1 an 1 γ,#i γ, > 0. That is, we assume the weights only act on coorinates without permutation-invariance conitions an leave the remaining coorinates unweighte. Thus we ask how much permutation-invariance can relax (by now well-establishe) necessary conitions on the sequences (γ,m ) m=1, N, for (strong) polynomial tractability: Assuming polynomial tractability with constants C 1, p > 0 an q 0 it is easy to check that e(n,;int,s I (E µ,α )) (2C)1/p q/p N 1/p, 8

9 see [Woź09] for etails. For every N an all κ > p this implies e(n,;int,s I (E µ,α ))κ 1+(2C) κ/p qκ/p ζ(κ/p) C qκ/p N=0 with some finite constant C > 0 which only epens on C,κ, an p. On the other han, Theorem 2 provies the lower boun e(n,;int ;S I (E µ,α ))κ N=0 In conclusion lim sup = N (,I ) 1 N=0 m {0,1} #I m 1... m #I = (#I +1) ν κ,n = 1 # 1 N=0 l {0,1} #I #I j=1 µ κ,u(ψ(n)) = k µ κ,u(k) j u(l) ( ) 1+γ κ,#i +j. γ,#i +j (#I ( +1) #I j=1 1+γ κ,#i +j) < (13) qκ/p is a necessary conition for (strong) polynomial tractability. In particular, this implies that we can have strong polynomial tractability (q = 0) only if the number of permutation-invariant coorinates #I is uniformly boune. If so, then (13) yiels that lim sup j=1 γ κ,j < for all κ > p which resembles the assertion state in [Woź09, Cor. 1]. If we assume polynomial tractability with q > 0 then we fin that for all κ > p #I #I O( qκ/p j=1 γ,#i κ ) an limsup +j <. ln Note that if #I is boune then the latter conition again coincies with the known conition given by Woźniakowski. Otherwise, if #I grows like say β for some β (0,1], then our conition turns out to be less restrictive, i.e., we may have polynomial tractability although polynomially many coorinates are unweighte (but permutationinvariant). 9 κ

10 3 Proofs In orer to prove Theorem 1 we basically combine the ieas state in [SW97, Woź09] with the technique evelope in [Wei12]. Since the proof is a little bit technical we ivie it into three steps which are organize as follows. InafirststepweshowthatforanygivenintegrationrulethatusesN < N functionvalues there exists a certain fooling function which shows that e(n,;int,s I (E,α )) 1. Afterwars, in Step 2, we notice that this lower boun is sharp, because e(0,;int,s I (E,α )) 1 for every N an all α > 1. Finally, we present a cubature rule that uses at most N = N function values, whereas its worst case error is no larger than the boun state in (7). Proof (Theorem 1). Step 1. Following Sloan an Woźniakowski [SW97] we fix α > 1, as well as N, an consier an arbitrary linear cubature rule A := A N,, given by (1), with N < N (,I ) := (#I + 1) 2 #I. Without loss of generality let us assume that I = {1,...,#I }. We will show that there exists a function f N in the unit ball of S I (E,α ) such that A(f N ) = 0, whereas the integral of f N equals 1. Fork Z let e k := exp(2πik ). Following thelines of [Wei12] we efine alinear operator S I : E,α E,α calle symmetrizer by (S I e k )(x) := 1 #S I σ S I e k (σ(x)), x [0,1], ancontinuousextension. ThereinS I isgivenin(3). ThenS I e k = (1/#S I ) λ S I e λ 1 (k)( ) is I -permutation-invariant in the sense of (2), i.e., (S I e k )(σ(x)) = (S I e k )(x) for all x [0,1] an any σ S I. Moreover, for k Z let M I (k) be efine as in [Wei12]. That is, M I (k)! enotes the number of ifferent permutations σ S I such that σ(k) = λ(k) for any fixe λ S I. To prove the claim we consier the set introuce in (10). Observing that # = N (,I ) = (#I +1) 2 #I we choose a bijection ψ: {0,1,...,N (,I ) 1}. (14) Furthermore, we consier the homogeneous linear system N n=0 a n (S I e ψ(n) )(t (i) ) M I (ψ(n))! = 0, i = 1,...,N, which consists of N < N (,I ) linear equations in N + 1 complex variables a n, n = 0,1,...,N. Here the points t (i) [0,1], i = 1,...,N, enote the integration noes use by 10

11 thecubatureruleaapplietothefunction0 S I (E,α ). Clearly, wecanselect anon-trivial solution a = (a n ) N n=0 CN+1 of this system, scale such that for some n {0,1,...,N} a n = 1 max n=0,1,...,n a n. Next we efine the function f N : [0,1] C by ( N ) f N (x) := #S I (S I e ψ(n ))(x) a n (S I e ψ(n) )(x), M I (ψ(n))! x [0,1]. (15) n=0 Observe that then f N (t (i) ) = 0 for all i = 1,...,N an thus we have A(f N ) = 0. Moreover, weseethatf N isi -permutation-invariantasaprouctofi -permutation-invariantfunctions. It remains to show that f N is a suitable fooling function for our integration problem. That is, we show that f N is an element of the unit ball of S I (E,α ), i.e., its Fourier coefficients satisfy f N (k) 1 ( k1... k ) α for every k = (k 1,...,k ) Z, (16) an that its integral is as large as possible, i.e., Int (f N ) = f N (0) = 1. For k Z we calculate f N (k) = f N,e k L2 = #S I 1 #S I = 1 #S I σ S I λ S I n=0 σ S I e σ 1 ( ψ(n )) N N n=0 a n M I (ψ(n))! a n e σ(ψ(n M I (ψ(n))! ))e λ(ψ(n)),e k L 2 1 #S I λ S I e λ 1 (ψ(n)),e k L 2 an e σ(ψ(n ))e λ(ψ(n)),e k ( 2πi ( λ(ψ(n)) k σ(ψ(n )) ) ) x x L 2 = exp [0,1] = δ [λ(ψ(n)) k σ(ψ(n ))=0], whereδ [C] equals1iftheconitionc isfulfillean0otherwise. Inparticularfrom f N (k) 0 it follows that there exist λ an σ in S I, as well as n {0,1,...,N}, such that k = λ(ψ(n)) σ(ψ(n ) { 1,0,1}, (17) 11

12 since ψ(n) {0,1} for all n an for every permutation λ S I the multi-inex λ(ψ(n)) is an element of {0,1}, too. Hence we arrive at f N (k) = 0 if k Z \{ 1,0,1} an f N (k) = 1 #S I σ S I δ [ k+σ(ψ(n )) {0,1} ] λ S I ( N n=0 a n M I (ψ(n))! δ [λ(ψ(n))=k+σ(ψ(n ))] if k { 1,0,1}. We will show that the summation within the brackets can be reuce to at most one non-vanishing term. Therefore assume σ S I to be fixe an consier h := k+σ(ψ(n )) {0,1}. For this h there exists one (an only one) multi-inex j such that (j m ) m I c = (h m ) m I c an #{m I j m = 1} = #{m I h m = 1}, because λ an σ in S I leave the coorinates in I c = {#I + 1,...,} fixe. Since the mapping ψ was assume to be bijective, j uniquely efines n(k,n,σ) := ψ 1 (j) {0,1,...,# 1}. Hence, there can be at most one n {0,1,...,N} with λ(ψ(n)) = h. If so, then there exist exactly M I (ψ(n(k,n,σ)))! ifferent permutations λ S I such that λ(ψ(n(k,n,σ))) = h = k+σ(ψ(n )). If k = 0 Z then h = σ(ψ(n )) which implies j = ψ(n ) an thus n(0,n,σ) = ψ 1 (j) = n for all σ S I. Consequently, we obtain f N (k) = 1 #S I which particularly yiels σ S I δ [k+σ(ψ(n )) {0,1} ] δ [n(k,n,σ) N] a n(k,n,σ), k { 1,0,1}, (18) Int (f N ) = f N (0) = 1 #S I σ S I a n = a n = 1. ) 12

13 Furthermore, formula (18) implies f N (k) ( ) α k 1... k = fn (k) 1 #S I σ S I 1 #S I δ [ k+σ(ψ(n )) {0,1} ] σ S I a n = 1, δ[n(k,n an(k,n,σ) N],σ) for k { 1,0,1}. Together with f N (k) ( ) α k 1... k = 0 for k Z \ { 1,0,1} this finally proves (16) an completes this step. Step 2. Clearly, we have e(0,;int,s I (E,α )) 1 for all N an any I {1,...,} since the worst case error of the zero algorithm A 0, 0 is given by wor (A 0, ;Int,S I (E,α )) = sup f(0) 1. f S I (E,α ) 1 Step 3. To show (7) we once more follow the arguments given in [SW97]. There it has been shown that the 2 -point prouct-rectangle rule R 2,: E,α C, f R 2,(f) := 1 ( j1 f 2 2, j 2 2,..., j ), 2 j {0,1} is a suitable algorithm for Int on E,α. For N an α > 1 its worst case error was foun to be ( wor (R 2,;Int,E,α ) = 1+ ζ(α) ) 1. (19) 2 α 1 Since S I (E,α ) is a linear subspace of E,α, where S I (E,α ) = E,α, we can restrict the algorithm R 2, to I -permutation-invariant functions in E,α an its worst case error remains boune from above by (19). On the other han, it is easy to see that A N,(f) := k = 1 2 k #S I 2 M I (k)! f ( k 2 f(σ(k)/2) M I (k)! σ S I ) = 1 2 j {0,1} f ( ) j = R 2 2,(f) (20) 13

14 on S I (E,α ). Hence the restriction of R 2,, i.e., the algorithm A N, efine by (20), is a cubature rule in the sense of (1) that uses not more than # = N = N (,I ) function values in imension ; see (10). Consequently, e(n,;int,s I (E,α )) wor (A N,;Int,S I (E,α )) ( wor (R 2,;Int,E,α ) = 1+ ζ(α) ) 1 2 α 1 which completes the proof. The proof of Theorem 2 can be erive using only one aitional argument. As in the previous proof, we construct a suitable fooling function g N that lies in the unit ball of S I (E µ,α ) an for which Int (g N ) is as large as possible while A(g N ) = 0. Proof (Theorem 2). Since we now eal with the weighte case we nee to show that for every k = (k 1,...,k ) Z we have µ,u(k) ĝ N (k) ( ) α k1... k instea of (16). We start by selecting a bijection ψ: {0,1,...,N (,I ) 1} which provies a non-increasing orering of the weights; see (14) an (11). Setting g N (x) := µ,u(ψ(n)) µ,u(ψ(n )) f N (x), x [0,1], it now remains to show that for every k Z with ĝ N (k) 0 we can estimate µ,u(ψ(n)) µ,u(ψ(n )) µ,u(k). (21) Recall that here the relation of k an n is given by k = λ(ψ(n)) σ(ψ(n )) for some λ,σ S I an a certain n {0,1,...,N}; see (17). To prove (21) we first note that the latter representation of k yiels #u I (k) min{#u I (ψ(n))+#u I (ψ(n )),#I }, since λ an σ o not change the size of the support of the respective multi-inices ψ(n) an 14

15 ψ(n ) in. Consequently we obtain #u I (k) i=1 γ,#i i+1 = #u I (ψ(n)) i=1 #u I (ψ(n)) i=1 #u I (ψ(n)) i=1 γ,#i i+1 γ,#i i+1 γ,#i i+1 min{#u I (ψ(n))+#u I (ψ(n )),#I } i=#u I (ψ(n))+1 min{#u I (ψ(n )),#I #u I (ψ(n))} #u I (ψ(n )) i=1 i=1 γ,#i i+1 γ,#i #u I (ψ(n)) i+1 γ,#i i+1, (22) by exploiting the general assumption 1 γ,1... γ,#i 0. Moreover, for the coorinates relate to the prouct weight part we conclue which immeiately implies that γ,m m u I c (k) u I c (k) u I c (ψ(n)) u I c (ψ(n )) m u I c (ψ(n)) γ,m m u I c (ψ(n )) γ,m. (23) Together with the representation of µ state in (9) the estimates (22) an (23) imply µ,u(ψ(n)) µ,u(ψ(n )) µ,u(k). This leas us to the claime boun (21) using the impose orering µ,u(ψ(n)) µ,u(ψ(n)) for n = 0,1,...,N. Finally, since f N (0) = 1, the integral of g N is lower boune by Int (g N ) = ĝ N (0) = µ,u(ψ(n)) µ,u(ψ(n )) f N (0) µ,u(ψ(n)) = ν,n. The fact that A(g N ) = µ,u(ψ(n)) µ,u(ψ(n )) A(f N ) = 0 completes the proof. Acknowlegements I want to thank Henryk Woźniakowski for pointing out the existence of the reference [Woź09] to me. Moreover, I like to thank Dirk Nuyens an Gowri Suryanarayana for motivating me to stuy integration problems of permutation-invariant functions, as well as for their kin hospitality uring my stays in Leuven in October 12 an in March

16 References [NW08] E. Novak an H. Woźniakowski - Tractability of Multivariate Problems. Vol. I: Linear Information. EMS Tracts in Mathematics 6. European Mathematical Society (EMS), Zürich [Sie13] P. Sielecki - Uniform weak tractability. J. Complexity 29, 2013, pp [SW97] I. Sloan an H. Woźniakowski - An intractability result for multiple integration. Mathematics of Computation 66, 1997, pp [Wei12] M. Weimar - The complexity of linear tensor prouct problems in (anti)symmetric Hilbert spaces. J. Approx. Theory 164(10), 2012, pp [Woź09] H. Woźniakowski - Tractability of multivariate integration for weighte Korobov spaces: my 15 year partnership with Ian Sloan, in: P. L Ecuyer an A.B. Owen (Es.) - Monte Carlo an Quasi-Monte Carlo Methos Springer, Berlin. 2009, pp

On lower bounds for integration of multivariate permutation-invariant functions

On lower bounds for integration of multivariate permutation-invariant functions On lower bounds for of multivariate permutation-invariant functions Markus Weimar Philipps-University Marburg Oberwolfach October 2013 Research supported by Deutsche Forschungsgemeinschaft DFG (DA 360/19-1)

More information

Tractability results for weighted Banach spaces of smooth functions

Tractability results for weighted Banach spaces of smooth functions Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March

More information

Integration of permutation-invariant. functions

Integration of permutation-invariant. functions permutation-invariant Markus Weimar Philipps-University Marburg Joint work with Dirk Nuyens and Gowri Suryanarayana (KU Leuven, Belgium) MCQMC2014, Leuven April 06-11, 2014 un Outline Permutation-invariant

More information

Generalized Tractability for Multivariate Problems

Generalized Tractability for Multivariate Problems Generalize Tractability for Multivariate Problems Part II: Linear Tensor Prouct Problems, Linear Information, an Unrestricte Tractability Michael Gnewuch Department of Computer Science, University of Kiel,

More information

ON THE OPTIMAL CONVERGENCE RATE OF UNIVERSAL AND NON-UNIVERSAL ALGORITHMS FOR MULTIVARIATE INTEGRATION AND APPROXIMATION

ON THE OPTIMAL CONVERGENCE RATE OF UNIVERSAL AND NON-UNIVERSAL ALGORITHMS FOR MULTIVARIATE INTEGRATION AND APPROXIMATION ON THE OPTIMAL CONVERGENCE RATE OF UNIVERSAL AN NON-UNIVERSAL ALGORITHMS FOR MULTIVARIATE INTEGRATION AN APPROXIMATION MICHAEL GRIEBEL AN HENRYK WOŹNIAKOWSKI Abstract. We stuy the optimal rate of convergence

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

Discrete Mathematics

Discrete Mathematics Discrete Mathematics 309 (009) 86 869 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: wwwelseviercom/locate/isc Profile vectors in the lattice of subspaces Dániel Gerbner

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

Monte Carlo Methods with Reduced Error

Monte Carlo Methods with Reduced Error Monte Carlo Methos with Reuce Error As has been shown, the probable error in Monte Carlo algorithms when no information about the smoothness of the function is use is Dξ r N = c N. It is important for

More information

Lower bounds on Locality Sensitive Hashing

Lower bounds on Locality Sensitive Hashing Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

The Generalized Incompressible Navier-Stokes Equations in Besov Spaces

The Generalized Incompressible Navier-Stokes Equations in Besov Spaces Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne

More information

Witt#5: Around the integrality criterion 9.93 [version 1.1 (21 April 2013), not completed, not proofread]

Witt#5: Around the integrality criterion 9.93 [version 1.1 (21 April 2013), not completed, not proofread] Witt vectors. Part 1 Michiel Hazewinkel Sienotes by Darij Grinberg Witt#5: Aroun the integrality criterion 9.93 [version 1.1 21 April 2013, not complete, not proofrea In [1, section 9.93, Hazewinkel states

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Discrete Operators in Canonical Domains

Discrete Operators in Canonical Domains Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Convergence of Random Walks

Convergence of Random Walks Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

Acute sets in Euclidean spaces

Acute sets in Euclidean spaces Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of

More information

Step 1. Analytic Properties of the Riemann zeta function [2 lectures]

Step 1. Analytic Properties of the Riemann zeta function [2 lectures] Step. Analytic Properties of the Riemann zeta function [2 lectures] The Riemann zeta function is the infinite sum of terms /, n. For each n, the / is a continuous function of s, i.e. lim s s 0 n = s n,

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Abstract A nonlinear partial differential equation of the following form is considered:

Abstract A nonlinear partial differential equation of the following form is considered: M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent

More information

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial

More information

The effect of dissipation on solutions of the complex KdV equation

The effect of dissipation on solutions of the complex KdV equation Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS

SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS GEORGE A HAGEDORN AND CAROLINE LASSER Abstract We investigate the iterate Kronecker prouct of a square matrix with itself an prove an invariance

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Permanent vs. Determinant

Permanent vs. Determinant Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED

More information

arxiv: v1 [math-ph] 5 May 2014

arxiv: v1 [math-ph] 5 May 2014 DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

arxiv: v4 [cs.ds] 7 Mar 2014

arxiv: v4 [cs.ds] 7 Mar 2014 Analysis of Agglomerative Clustering Marcel R. Ackermann Johannes Blömer Daniel Kuntze Christian Sohler arxiv:101.697v [cs.ds] 7 Mar 01 Abstract The iameter k-clustering problem is the problem of partitioning

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

arxiv: v1 [math.co] 15 Sep 2015

arxiv: v1 [math.co] 15 Sep 2015 Circular coloring of signe graphs Yingli Kang, Eckhar Steffen arxiv:1509.04488v1 [math.co] 15 Sep 015 Abstract Let k, ( k) be two positive integers. We generalize the well stuie notions of (k, )-colorings

More information

A. Incorrect! The letter t does not appear in the expression of the given integral

A. Incorrect! The letter t does not appear in the expression of the given integral AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

A Sketch of Menshikov s Theorem

A Sketch of Menshikov s Theorem A Sketch of Menshikov s Theorem Thomas Bao March 14, 2010 Abstract Let Λ be an infinite, locally finite oriente multi-graph with C Λ finite an strongly connecte, an let p

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition

More information

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae

More information

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number

More information

REAL ANALYSIS I HOMEWORK 5

REAL ANALYSIS I HOMEWORK 5 REAL ANALYSIS I HOMEWORK 5 CİHAN BAHRAN The questions are from Stein an Shakarchi s text, Chapter 3. 1. Suppose ϕ is an integrable function on R with R ϕ(x)x = 1. Let K δ(x) = δ ϕ(x/δ), δ > 0. (a) Prove

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

LECTURE NOTES ON DVORETZKY S THEOREM

LECTURE NOTES ON DVORETZKY S THEOREM LECTURE NOTES ON DVORETZKY S THEOREM STEVEN HEILMAN Abstract. We present the first half of the paper [S]. In particular, the results below, unless otherwise state, shoul be attribute to G. Schechtman.

More information

All s Well That Ends Well: Supplementary Proofs

All s Well That Ends Well: Supplementary Proofs All s Well That Ens Well: Guarantee Resolution of Simultaneous Rigi Boy Impact 1:1 All s Well That Ens Well: Supplementary Proofs This ocument complements the paper All s Well That Ens Well: Guarantee

More information

Robustness and Perturbations of Minimal Bases

Robustness and Perturbations of Minimal Bases Robustness an Perturbations of Minimal Bases Paul Van Dooren an Froilán M Dopico December 9, 2016 Abstract Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

More information

Function Spaces. 1 Hilbert Spaces

Function Spaces. 1 Hilbert Spaces Function Spaces A function space is a set of functions F that has some structure. Often a nonparametric regression function or classifier is chosen to lie in some function space, where the assume structure

More information

arxiv: v4 [math.pr] 27 Jul 2016

arxiv: v4 [math.pr] 27 Jul 2016 The Asymptotic Distribution of the Determinant of a Ranom Correlation Matrix arxiv:309768v4 mathpr] 7 Jul 06 AM Hanea a, & GF Nane b a Centre of xcellence for Biosecurity Risk Analysis, University of Melbourne,

More information

VIBRATIONS OF A CIRCULAR MEMBRANE

VIBRATIONS OF A CIRCULAR MEMBRANE VIBRATIONS OF A CIRCULAR MEMBRANE RAM EKSTROM. Solving the wave equation on the isk The ynamics of vibrations of a two-imensional isk D are given by the wave equation..) c 2 u = u tt, together with the

More information

Noether s theorem applied to classical electrodynamics

Noether s theorem applied to classical electrodynamics Noether s theorem applie to classical electroynamics Thomas B. Mieling Faculty of Physics, University of ienna Boltzmanngasse 5, 090 ienna, Austria (Date: November 8, 207) The consequences of gauge invariance

More information

Proof of SPNs as Mixture of Trees

Proof of SPNs as Mixture of Trees A Proof of SPNs as Mixture of Trees Theorem 1. If T is an inuce SPN from a complete an ecomposable SPN S, then T is a tree that is complete an ecomposable. Proof. Argue by contraiction that T is not a

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable

More information

3 The variational formulation of elliptic PDEs

3 The variational formulation of elliptic PDEs Chapter 3 The variational formulation of elliptic PDEs We now begin the theoretical stuy of elliptic partial ifferential equations an bounary value problems. We will focus on one approach, which is calle

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

Analyzing the Structure of Multidimensional Compressed Sensing Problems through Coherence

Analyzing the Structure of Multidimensional Compressed Sensing Problems through Coherence Analyzing the Structure of Multiimensional Compresse Sensing Problems through Coherence A. Jones University of Cambrige B. Acock Purue Univesity A. Hansen University of Cambrige Abstract In previous work

More information

On the enumeration of partitions with summands in arithmetic progression

On the enumeration of partitions with summands in arithmetic progression AUSTRALASIAN JOURNAL OF COMBINATORICS Volume 8 (003), Pages 149 159 On the enumeration of partitions with summans in arithmetic progression M. A. Nyblom C. Evans Department of Mathematics an Statistics

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

The Ehrenfest Theorems

The Ehrenfest Theorems The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent

More information

Introduction to variational calculus: Lecture notes 1

Introduction to variational calculus: Lecture notes 1 October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1 Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

7.1 Support Vector Machine

7.1 Support Vector Machine 67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to

More information

Resistant Polynomials and Stronger Lower Bounds for Depth-Three Arithmetical Formulas

Resistant Polynomials and Stronger Lower Bounds for Depth-Three Arithmetical Formulas Resistant Polynomials an Stronger Lower Bouns for Depth-Three Arithmetical Formulas Maurice J. Jansen University at Buffalo Kenneth W.Regan University at Buffalo Abstract We erive quaratic lower bouns

More information

Analysis IV, Assignment 4

Analysis IV, Assignment 4 Analysis IV, Assignment 4 Prof. John Toth Winter 23 Exercise Let f C () an perioic with f(x+2) f(x). Let a n f(t)e int t an (S N f)(x) N n N then f(x ) lim (S Nf)(x ). N a n e inx. If f is continuously

More information

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the

More information

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything

More information

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

GCD of Random Linear Combinations

GCD of Random Linear Combinations JOACHIM VON ZUR GATHEN & IGOR E. SHPARLINSKI (2006). GCD of Ranom Linear Combinations. Algorithmica 46(1), 137 148. ISSN 0178-4617 (Print), 1432-0541 (Online). URL https://x.oi.org/10.1007/s00453-006-0072-1.

More information