Low-discrepancy point sets lifted to the unit sphere
|
|
- Alyson Tucker
- 6 years ago
- Views:
Transcription
1 Low-discrepancy point sets lifted to the unit sphere Johann S. Brauchart School of Mathematics and Statistics, UNSW MCQMC 2012 (UNSW) February 17 based on joint work with Ed Saff (Vanderbilt) and Josef Dick, Ian Sloan, Rob Womersley (UNSW) and Christoph Aistleitner (TU-Graz)
2 Rob Womersley ( Mathematics & Statistics in all Spheres of Life
3 How do you uniformly arrange N (repulsive) particles on a sphere?
4 How do you quasi-uniformly arrange N (repulsive) particles on a sphere? pre-scribed density, well-separated, small holes
5 Preliminaries UNIFORM DISTRIBUTION ON THE UNIT SPHERE WORST-CASE ERROR (WCE) CUI AND FREEDEN S GENERALIZED DISCREPANCY STOLARSKY S INVARIANCE PRINCIPLE
6 Uniform Distribution on the Unit Sphere S d in R d+1 Definition A sequence {X N } is asymptotically uniformly distributed over S d if # {k : x k,n B} lim = σ d (B) N N for every σ d -measurable clopen set B in S d. Informally: A reasonable set gets a fair share of points as N becomes large. I[f ]:= f d σ d, S d Q N [f ]:= 1 N N f (x k ) Equiv.: {X N } is asymptotically uniformly distributed over S d if a k=1 lim Q N[f ] = I[f ] for every f C(S d ). N a I.e., ν(x (s) N ) σ d as N (in the weak- limit).
7 Sobolev Space over S d Sequence of positive real numbers a = (a (s) l ) l 0 satisfying lim l a (s) l = a(d, s) for some a(d, s) > 0, 2s [l + (d 1)/2] Define an inner product for functions f and g in C (S d ) by means of (f, g) H s a := l=0 Z (d,l) k=1 1 a (s) l ˆf (d) l,k ĝ(d) l,k, where the Laplace-Fourier coefficients are given by (d) (d) ˆf l,k = (f, Y l,k ) L 2 (S d ) = f (x)y (d) l,k (x) d σ d(x). S d The Sobolov space H s a(s d ) then is the completion of the space { } Z (d,l) f C (S d 1 ] 2 (d) ) : [ˆf l,k < l=0 k=1 a (s) l with respect to the Sobolev norm H s a := (, ) H s a.
8 Cui and Freeden s Generalized Discrepancy [SIAM J. SCI. COMPUT. 18:2 (1997)]... based on pseudodifferential operators yielding A Koksma-Hlawka like inequality QN [f ] I[f ] 6 DCF (X N ) f H 3/2 (S 2 ), where f is from the Sobolev space H 3/2 (S 2 ). D CF (X N ) in terms of elementary function expressible: 4π [D CF (X N )] 2 = 1 1 N 2 N j,k=1 log (1 + x j x k /2) 2. Definition (Equidistribution in H 3/2 (S 2 )) A sequence {X N } is equidistributed in H 3/2 (S 2 ) if lim N D CF(X N ) = an asymptotically uniformly distributed {X N } is equidistributed in C(S 2 ).
9 Interpretation as Worst Case Error (WCE) Sloan and Womersley [Adv. Comput. Math. 21 (2004)]... show that D CF (X N ) has a natural interpretation, as the worst-case error (apart from a constant factor) for the equal-weight cubature rule E n (f ) = 4π d n d n k=1 f (x k ) applied to a function f B(H), where B(H) is the unit ball in a certain Hilbert space H.
10 Worst Case Error (WCE) in a Sobolev Space H s a(s d ) Sequence of positive real numbers a = (a (s) l ) l 0 satisfying lim l a (s) l = a(d, s) for some a(d, s) > 0, 2s [l + (d 1)/2] The worst-case error of Q N [f ] wce(q N ; H s a(s d )) = sup { QN [f ] I[f ] : f H s (S d ), f d,s 1 }. By Cauchy-Schwarz inequality [ wce(qn ; H s a(s d )) ] 2 = 1 N 2 S d N K (a; x j x k ) j,k=1 S d K (a; x y) d σ d (x) d σ d (y). Example A Example B
11 Stolarsky s Invariance Principle for Euclidean Distance Spherical Cap centered at z with height t C(z; t):= { x S d : x, z t }, z S d, 1 t 1. Spherical cap L 2 -discrepancy: { 1 2 1/2 DL C 2 (X N ):= X N C(x; t) σ d (C(x; t)) N d σ d (x) d t}. 1 S d Theorem (Stolarsky [Proc. of the AMS 41:2 (1973)]) 1 N 2 N j,k=1 x j x k + H d(s d ) [ ] 2 H d (B d DL C ) 2 (X m ) = x y d σ d (x) d σ d (y). S d S d Stolarsky s proof makes use of Haar integrals over SO(d + 1). Reproved using reproducing kernel Hilbert space techniques (B-Dick [Proc. of the AMS, in press])
12 Sum of distance integral: V 1 (S d ) = S d x y d σ d (x) d σ d (y). S d Theorem (B. Womersley) Let d 2 and H s (S d ), s = (d + 1)/2, with reproducing kernel Then K d,s (x, y) = 2V 1 (S d ) x y, x, y S d. [ wce(qn ; H s (S d )) ] 2 = V 1 (S d ) 1 N N 2 x i x j i,j=1 = H d(s d ) [ 2 H d (B d D L 2 C N)] ) (X.
13 Point Constructions SPHERICAL DIGITAL NETS SPHERICAL FIBONACCI LATTICE APPROXIMATIVE SPHERICAL DESIGNS (ersatz FOR SPHERICAL DESIGNS?)
14 Area preserving map to S 2 & Consequences Lambert cylindrical equal-area projection Map to S d x:=φ s (y):= Φ s (φ 1 (y 1 ), φ 2 (y 2 )) S 2, where φ 1 (y 1 ):=2πy 1, φ 2 (y 2 ):=1 2y 2, and ( ) Φ s (φ, t):= 1 t 2 cos φ, 1 t 2 sin φ, t. Lemma (B-Dick, Numerische Mathematik (2011)) D ( ) N (S2, Ω ( ) ; Z N ) = D ( ) N ([0, 1)2, R ( ) ; Z N ). Theorem (B-Dick, Numerische Mathematik (2011)) Sequence (Z N ) N 2 s.t. Z N = Φ 2 (Z N ), Z N [0, 1) 2. [ ] ( 2 wce(q N ; H 3/2 24 a (S 2 )) + 2 ) 2 DN([0, 1) 2, R ( ) ; Z N ). 3
15 Definition (Isotropic Discrepancy) J N (Z N ):=D N ([0, 1] 2, A; Z N ), where A is the family of convex subsets of [0, 1] 2 Theorem (Aistleitner-B-Dick, submitted) D C (Φ(Z N )) 11J N (Z N ), Z N [0, 1) 2.
16 DIGITAL NETS AND SPHERICAL DIGITAL NETS
17 Digital Nets Definition (Star Discrepancy and local Discrepancy) D N(X N ) = sup δ XN (x), x [0,1] s δ XN (x) = X N [0, x) N VOL([0, x)) Digital Net (Characterization) A (t, m, s)-net in base b with N = b m points has the property that each elementary interval contains exactly b t points. s i=1 [ ai b d i, a ) 0 a i < b d i i + 1 d b d d s = m t i 0 d 1,..., d s Z (elementary interval) Example Construction Theorem X N (t, m, s)-net in base b: D N(X N ) C (m t) s 1 /b m t.
18 Low Discrepancy Point Sets and Sequences D N(X N ) = sup δ XN (x), x [0,1] s δ XN (x) = X N [0, x) N VOL([0, x)) THM.: X N (0, m, s)-net in base b: D N(X N ) C m s 1 /b m 1. HAMMERSLEY [Ann. New York Acad. Sci ] and HALTON [Numer. Math ] (i) for any N 2 there exists x 1,..., x N [0, 1) s s.t. D N((x 1,..., x N )) = O( (log N)s 1 ) N (ii) there exists a sequence (x n ) n 1 in [0, 1) s s.t. D N((x n ) n 1 ) = O( (log N)s ) N
19 Spherical Rectangle discrepancy Families of half-open axis-parallel rectangles lifted to S 2 Ω = {Φ s (R x,y ) : 0 x y 1}, Ω = {Φ s (R 0,y ) : 0 y 1}. Definition (Spherical Rectangle (Star) Discrepancy) D N (Ω; X N ):= sup δ N (R; X N ), D N(Ω ; X N ):= sup δ N (R; X N ). R Ω R Ω Theorem (B-Dick, Numerische Mathematik) For Z N (a (0, m, 2)-net in base b lifted to S 2 ) D b m(ω, Z N ) b2 m b + 1 b m + 1 b m Db m(ω, Z N ) b2 /4 m b + 1 b m + 1 b m By Roth [Mathematika (1954)] ( ) + 1 b b 2m ( ) + 1 b b 2m Recall ( 2b 1 4b + 3 (b + 1) 2 ), ( b b + 3 4(b + 1) 2 D N (Ω, X N ) D N(Ω, X N ) ( log 2 N + 3)/(2 8 N). ).
20 Isotropic & Spherical Cap Discrepancy Theorem (Aistleitner-B-Dick, submitted) Z N... (0, m, 2)-net in base b (N = b m ) J N (Z N ) 4 2b N. Z N... first N points of (0, 2)-sequence J N (Z N ) 4 2(b 2 + b 3/2 ) N. Remark Optimal rate of convergence for J N (cf. Beck and Chen, 1987 & 2008): N 2/3 (log N) c for some 0 c 4!
21 WCE of Numerical Integration Rules based on Digital Nets (Sobol points) on Spheres m swce N 3/2 swce e e e e e e e e e e e e e e e e e e e e Conjecture (B-Dick, Numerische Mathematik) d = 2 d = 200 A sequence of Sobol points in [0, 1) s (s 2) lifted to the s-sphere achieves optimal convergence rate of the sum of distances.
22 FIBONACCI LATTICE AND SPHERICAL FIBONACCI LATTICES
23 Fibonacci lattice points in [0, 1] 2 Fibonacci sequence: F 1 = 1, F 2 = 1, and F n+1 = F n + F n 1, n 1. Fibonacci lattice in [0, 1] 2 ( k f k :=, F n { kfn 1 F n }), 0 k < F n, where {x} = x x for nonnegative real numbers x. The set F n := {f 0,..., f Fn 1} is called a Fibonacci lattice point set.
24 Fibonacci lattice points on the unit sphere S 2 Lambert transformation Φ : [0, 1] 2 S 2 Φ(x, y) = ( ) 2 cos(2πy) x x 2, 2 sin(2πy) x x 2, 1 2x. Note: Area preserving map. Fibonacci lattice on S 2 The spherical Fibonacci lattice points are then given by and the point set z k = Φ(f k ), 0 k < F n, Z n = {z 0,..., z Fn 1} is the spherical Fibonacci lattice point set.
25 Separation results (B Dick, work in progress) Minimum distance: d min (P) = min 0 k<l<n y k y l, P = {y 0,..., y N 1 } Rd+1. Theorem (Fibonacci points in [0, 1] 2 ) d min (F 2m+1 ) = 1 F 2m+1, m 1, d min (F 2m ) = 2Fm F 2m, m 2. Theorem (Fibonacci points in S 2 ) d min (Z 2m+1 ) 1 F 2m+1, m 1. Conjecture (Fibonacci points in S 2 ) d min (Z n) = 2 Fn, n 2. n Fn Zndmin (Zn) ± ± ± ± ± ± ± ± ± ±
26 Isotropic & Spherical Cap Discrepancy Corollary (Aistleitner-B-Dick, submitted) { D C 44 2/F m if m is odd, (Z Fm ) 44 8/F m if m is even. Table: D(Z Fm )... maximum of the absolute values of the local discrepancies of Z Fm at the spherical caps centered at the spherical Fibonacci points Z Fm. m Fm D(Z Fm ) Fm 3/4 log Fm D(Z Fm ) Fm 3/4 log Fm m Fm D(Z Fm ) Fm 3/4 log Fm D(Z Fm ) F 3/4 m log Fm
27 Spherical Cap L 2 -discrepancy (B Dick, work in progress) [ 4 D L C ] 2 4 (Zn) = Fn 2 Fn 1 z j z k j,k=0 [ [ ] 2 n Fn 4 D L C ] 2 (Zn) 2 Fn 3/2 4Fn 3/2 D L C (Zn) e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e cf. B [Uniform Distribution Theory 6:2 (2011)]
28 SEQUENCES OF APPROXIMATE SPHERICAL DESIGNS
29 Motivation Lower bound in terms of number of nodes: Theorem (He-Sl 2005 (d = 2); He-Sl 2006 (d 2)) wce(q N ; H s (S d )) c (d, s) N s/d. Upper bound in terms of exactness degree: Theorem (He-Sl 2005 (H 3/2 (S 2 )); He-Sl 2006 (H s (S 2 ), s > 1); B.-He 2007 (H s (S d ), s > d/2)) (i) Q N(t) is exact for all polynomials of degree t, and (ii) (Q N(t) ) t N satisfies a local regularity condition (Property (R) ). a wce(q N(t) ; H s (S d )) C (d, s) t s. a Automatically fulfilled for positive weight numerical integration rules (Reimer 2000). Optimality: There are always Q N(t) with N(t) C t d.
30 Spherical Designs Definition (Delsarte, Goethals and Seidel 1977 ) Spherical t-designs x 1,..., x N(t) with N(t) points have the property S d P(x) d σ d (x) = 1 N(t) for all polynomials P with deg P t. N(t) P(x j ) j=1
31 Theorem (Upper Bound) A sequence of N(t)-point spherical designs (X N(t) ) t 1 on S d has the property that there exists C(s, d) > 0, independent of t and N(t), s. t. 1 C(s, d) sup f (x) f (x) d σ d (x) N(t) x X S t s as t, d N(t) f H s (S d ), f H s 1 where H s denotes the norm in the Sobolev space H s. H 3/2 (S 2 ): K. Hesse & I. H. Sloan [Bull. Austral. Math. Soc. 71 (2005)] H s (S 2 ), s > 1: K. Hesse & I. H. Sloan [J. Approx. Theory 141 (2006)] H s (S d ), s > d/2: B. & K. Hesse [Constr. Approx. 25 (2007)]
32 Approximate Spherical Designs Definition (B-Saff-Sloan-Womersley (work in progress)) A sequence (X N(t) ) t 1 on S d, with t and N(t), is said to be a sequence of approximate spherical t-designs for H s (S d ) for some s > d/2 if there exists c(s, d) > 0, independent of t and N(t), s. t. sup f H s (S d ), f H s 1 1 f (x) f (x) d σ d (x) N(t) x X S d N(t) c(s, d) t s as t. Proposition Spherical designs are approximative spherical designs for all s > d/2.
33 Approximate Spherical Designs are better than average Let Set Definition A(x y) = l=1 A N [X N ]:= 1 N 2 a l }{{} 0 N Z (d, l)p (d) l (x y). j=1 k=1 N A(x j x k ). Spherical average of A N [X N ] over all N-point config. on S d EA N := S d A N [{x 1,..., x N }] d σ d (x 1 ) d σ d (x N ). S d Proposition EA N = A(1) N.
34 Average WCE E [ wce(q N ; H s a(s d )) ] 2 K (a; 1) a (s) = 0. N Average Approximate Integration E 1 N N N 2 K t (a; x j x k ) a (s) 0 = K t(a; 1) a (s) 0. N j=1 k=1 Compare with optimal rate N 2s/d if N = O(t d ).
35 Characterization Theorem (B-Saff-Sloan-Womersley (work in progress)) Let s > 1. a = (a (s) l ) l 0 monotonically decreasing positive numbers satisfying lim l a (s) l = a(d, s) for some a(d, s) > 0, 2s [l + (d 1)/2] a (s) l (l + 1/2) 2s a(2; s)= O(l 2 ) a (s) l+1 (l + 3/2) 2s a (s) l (l + 1/2) 2s = O(l 3 ). Then a sequence (X N(t) ) t 0 of N(t)-point configurations is a sequence of approximate spherical t-designs for H s a(s 2 ) iff 1 N 2 N N j=1 k=1 K t (a; x j x k ) a (s) 0 = O(t 2s ) as t.
36 Variational characterization of sph. designs; Gr-Ti 1993, Sl-Wo 2009 An N(t)-point config. which minimizes the non-negative functional 1 N 2 N N j=1 k=1 K t (a; x j x k ) a (s) 0 and, in addition, evaluates it to zero is already a spherical t-design.
37 Families of Approximate Spherical Designs Theorem (B-Saff-Sloan-Womersley (work in progress)) Let d 2, d/2 < s < d/ The maximizer of the generalized distance functional N x j x k 2s d. j,k=1 form a sequence of approximate spherical t-designs for H s (S d ), provided N(t) t d.
38 The End Thank You!
39 Numerical Integration in the Unit Cube Koksma-Hlawka Inequality N 1 1 f (x n ) f (x) dx N [0,1] D N(X N )V HK (f ). s n=0 Return
40 Example (Digital (0, 3, 2)-net in base 2 with 8 points)
41 Example (Digital (0, 3, 2)-net in base 2 with 8 points)
42 Example (Digital (0, 3, 2)-net in base 2 with 8 points)
43 Example (Digital (0, 3, 2)-net in base 2 with 8 points)
44 Example (Digital (0, 3, 2)-net in base 2 with 8 points)
45 Example (Digital (0, 3, 2)-net in base 2 with 8 points) Return
46 Construction of (0, m, s)-nets by means of Sobol pts Let l = b 1 + b b n 2 n 2 n. The j-th component of the l-th point in the Sobol sequence x l,j = 2 n [b 1 v 1,j b n v n,j ], v k,j = m k,j 2 n k. The direction numbers {m 1,j, m 2,j, } are recursively defined by s j 1 m k,j = 2 ν a ν,j m k ν,j 2 s j m k sj,j m k sj,j, k s j + 1, ν=1 in terms of the coefficients of the primitive polynomial x s j + a 1,j x s j a sj 1,jx + 1 over Z 2 and the initial direction numbers m 1,j,..., m sj,j. (For the first component one may chose m k,1 = 1 for all k.) decoded as the integer a 1,j 2 s j a sj 1,j 2 is Return
47 88 x x x x2 86 Example Figure: 2-dim. projections of Sobol points in [0, 1)200 with N = 256, 512. Return
48 Proof. Substitution of d σ d (x d ) = ω d 1 ω d ( 1 τ 2 d ) d/2 1 d τd d σ d 1 (x d 1 ) d = 2, 3,..., yields σ 1 (Φ 1 (R y )) = 2πy 1 0 d φ/(2π) = y 1 and σ d (Φ d (R y )) = ω d 1 ω d 1 t d ( 1 τ 2 d ) d/2 1 d τd = T d (t d ) σ d 1 (Φ d 1 (R y )), where (since 2 d 1 ω d 1 /ω d = 1/B(d/2, d/2)) T d (t) = ω d 1 ω d 1 t = I (1 t)/2 (d/2, d/2). ( 1 τ 2 ) d/2 1 d τ = 2 d 1 ω d 1 ω d Φ d 1 (R y) (1 t)/2 Note that T 2 (t) = (1 t)/2 and t d = 1 2y d (2 d s). 0 d σ d 1 (x d 1 ) u d/2 1 d u Return
49 A(J; X N )... number of points of X N in J Definition (Local Discrepancy) δ N (J; X N ):= A(J; X N) N VOL(J). Return
50 Beck (1984) To any N-point cfg. {x 1,..., x N } on S d there exists a spherical cap C s.t. c 1 N 1/2 1/(2d) < # {k : x k C} N σ d (C). There exists an N-point cfg. {x 1,..., x N } on S d s.t. for any spherical cap C # {k : x k C} N σ d (C) < c 2 N 1/2 1/(2d) log N. Return
51 d = 2; N = 2,..., Return
52 d = 200; N = 2,..., Comparison, pseudo-random points Return 10 5
53 N Return
54 Spherical Designs Much progress has been made since their introduction by Delsarte, Goethals and Seidel [Geometriae Dedicata 6 (1977)]: existence of spherical designs in general and existence of so-called tight spherical designs in particular, their construction by algebraic and by variational means culminating in a recently proposed proof of the N 2 -conjecture (or N d -conjecture) of Korevaar and Meyers [Integral Transform. Spec. Funct. 1 (1993)] that N(t) = O(t d ) points are not only necessary but also sufficient to form a spherical t-design by Andriy Bondarenko, Danylo Radchenko and Maryna Viazovska (arxiv: v3 [math.mg]). Return
55 Property (R) Return Definition A sequence (Q N(t) ) t 1 of weighted numerical integration rules Q N(t) on S d is said to have the property (R) (or to be quadrature regular ), if there exist positive numbers c 0 and c 1 independent of n with c 1 π/2, such that for all t 1 the weights w 1,..., w N(t) associated with the nodes x 1,..., x N(t) satisfy N(t) j=1, x j C(x;φ) w j c 0 σ d (C(x; c 1 /t)), where C(x; φ):={y S d : x y cos φ} denotes the spherical cap.
56 Lifting Digital Nets to the Unit Sphere in R d+1, d 2 Cylindrical Coordinates on the d-sphere x = x s = ( 1 t 2 s x s 1, t s ), x 2 = ( 1 t 2 2 x 1, t 2 ), x 1 = ( cos φ, sin φ ) 1 t d 1, x d S d (d = 1,..., s) 0 φ 2π. Arbitrary x S s represented through angle φ and heights t 2,..., t s x = x(φ, t 2,..., t s ), 0 φ 2π, 1 t 2,..., t s 1. A point y = (y 1,..., y s ) [0, 1) s mapped to x S s using ϕ 1 (y 1 ) = 2πy 1, ϕ d (y d ) = 1 2y d (d = 2, 3,..., s) and cylinder coordinates x = Φ s (y) = x ( ϕ 1 (y 1 ), ϕ 2 (y 2 ),..., ϕ s (y s ) ).
57 Area preserving map Normalized surface area measure σ d on S d d σ 1 (x 1 ) = d φ/(2π), d σ d (x d ) = ω d 1 ω d ( 1 t 2 d ) d/2 1 d td d σ d 1 (x d 1 ), d = 2, 3,.... Let R y = [0, y 1 ) [0, y s ) [0, 1) s. Proposition s σ s (Φ s (R y )) = y 1 I yd (d/2, d/2). d=2 Proof Regularized Beta function I z (a, b) = B z (a, b)/ B(a, b), B z (a, b) = z 0 u a 1 (1 u) b 1 d u.
58 Let VOL(R z ) = z 1 z 2 z s for z [0, 1) s. Chose y 1 = z 1, y 2 = z 2 and y d such that I yd (d/2, d/2) = z d. Then again σ s (Φ s (R y )) = z 1 z 2 z s. Area preserving map Ψ s : ( z 1,..., z s ) Φs ( ( h 1 d (z 1), h 1 d (z 2),..., h 1 d (z s) ) ). h 1 d (y) is the inverse of h 1(y) = y and h d (y) = I y (d/2, d/2), d 2. I z 1 d d, d = 2,..., 10 and d = z Return
59 Example A Return From the well-known identity l=1 1 ( l (l + 1) P l(t) = 1 2 log 1 + one derives the sequence a = (a (s) l ) l 0 with 1 t 2 ), t [ 1, 1], a (3/2) 0 = 1, a (3/2) 1 l = l (l + 1) (2l + 1), l 1. The Sobolev space H (3/2) a (S 2 ) is a reproducing kernel Hilbert space with reproducing kernel K (a; x y) = K (3/2) H (x, y) = a l=0 ( a (3/2) x y ) l (2l + 1) P l (x y) = 2 2 log Considered by Cui and Freeden (1997); generalized discrepancy. a (3/2) l (l + 1/2)2 l + 1/2 = (l + 1/2) 3 l (l + 1) 2l + 1 = /8 + O(l 3 ) as l. l 2
60 Example B Return Let d/2 < s < d/ The weights a (s) 0 = V d 2s (S d ) = 2 2s 1 Γ((d + 1)/2)Γ(s) πγ(d/2 + s), a (s) l = V d 2s (S d ) (d/2 s) l (d/2 + s) l, in the sequence a = (a (s) l ) l 0 define the reproducing kernel K (a; x y) = l=0 a (s) l Z (d, l) P (d) l (x y) = 2V d 2s (S d ) x y 2s d, for the Sobolev space H (s) a (S d ), d/2 < s < d/ a (s) l Γ((d + 1)/2)Γ(s) = + a (d; s) + O(l 3 ) as l. (l + 1/2) 2s π[ Γ(d/2 s)] l 2
Hyperuniformity on the Sphere
Hyperuniformity on the Sphere Johann S. Brauchart j.brauchart@tugraz.at 14. June 2018 Week 2 Workshop MATRIX Research Program "On the Frontiers of High Dimensional Computation" MATRIX, Creswick [ 1 / 56
More informationQMC designs and covering of spheres by spherical caps
QMC designs and covering of spheres by spherical caps Ian H. Sloan University of New South Wales, Sydney, Australia Joint work with Johann Brauchart, Josef Dick, Ed Saff (Vanderbilt), Rob Womersley and
More informationSpherical designs with close to the minimal number of points 1
Spherical designs with close to the minimal number of points Robert S. Womersley School of Mathematics and Statistics, University of New South Wales, Sydney, NSW, 05, Australia Abstract In this paper we
More informationDISCRETE MINIMAL ENERGY PROBLEMS
DISCRETE MINIMAL ENERGY PROBLEMS Lecture III E. B. Saff Center for Constructive Approximation Vanderbilt University University of Crete, Heraklion May, 2017 Random configurations ΩN = {X1, X2,..., XN }:
More informationTutorial on quasi-monte Carlo methods
Tutorial on quasi-monte Carlo methods Josef Dick School of Mathematics and Statistics, UNSW, Sydney, Australia josef.dick@unsw.edu.au Comparison: MCMC, MC, QMC Roughly speaking: Markov chain Monte Carlo
More informationSpherical t ɛ -Design for Numerical Approximations on the Sphere, I
Spherical t ɛ -Design for Numerical Approximations on the Sphere, I Xiaojun Chen Department of Applied Mathematics The Hong Kong Polytechnic University April 2015, Shanghai Jiaotong University X. Chen
More informationEfficient Spherical Designs with Good Geometric Properties
Efficient Spherical Designs with Good Geometric Properties Rob Womersley, R.Womersley@unsw.edu.au School of Mathematics and Statistics, University of New South Wales 45-design with N = 1059 and symmetrtic
More informationDiscrete Minimum Energy Problems
Discrete Minimum Energy Problems D.P. Hardin Department of Mathematics Vanderbilt University TU Graz, July 3, 2009 Reisz s energy The Riesz s-energy of ω N = {x 1, x 2,..., x N } R p is, for s > 0, E s(ω
More informationQuadrature using sparse grids on products of spheres
Quadrature using sparse grids on products of spheres Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at 3rd Workshop on High-Dimensional Approximation UNSW,
More informationQuadrature using sparse grids on products of spheres
Quadrature using sparse grids on products of spheres Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ANZIAM NSW/ACT Annual Meeting, Batemans Bay. Joint
More informationThe Discrepancy Function and the Small Ball Inequality in Higher Dimensions
The Discrepancy Function and the Small Ball Inequality in Higher Dimensions Dmitriy Georgia Institute of Technology (joint work with M. Lacey and A. Vagharshakyan) 2007 Fall AMS Western Section Meeting
More informationPolynomial interpolation on the sphere, reproducing kernels and random matrices
Polynomial interpolation on the sphere, reproducing kernels and random matrices Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at MASCOS Workshop on Stochastics
More informationDiscrepancies of Point Sequences on the Sphere and Numerical Integration
Discrepancies of Point Sequences on the Sphere and umerical Integration Peter J. Grabner, Bernhard Klinger and Robert F. Tichy Abstract Several concepts of discrepancy for point sequences on the d-dimensional
More informationAPPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION. I.H. Sloan, X. Wang and H.
APPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION I.H. Sloan, X. Wang and H. Wozniakowski Published in Journal of Complexity, Volume 20, Number 1,
More informationQuasi-Monte Carlo Methods for Applications in Statistics
Quasi-Monte Carlo Methods for Applications in Statistics Weights for QMC in Statistics Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February 2012 1 / 24 Quasi-Monte Carlo Methods for Applications
More informationSTRONG CONVERGENCE OF AN IMPLICIT ITERATION PROCESS FOR ASYMPTOTICALLY NONEXPANSIVE IN THE INTERMEDIATE SENSE MAPPINGS IN BANACH SPACES
Kragujevac Journal of Mathematics Volume 36 Number 2 (2012), Pages 237 249. STRONG CONVERGENCE OF AN IMPLICIT ITERATION PROCESS FOR ASYMPTOTICALLY NONEXPANSIVE IN THE INTERMEDIATE SENSE MAPPINGS IN BANACH
More informationAPPLIED MATHEMATICS REPORT AMR04/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR. J. Dick and F. Pillichshammer
APPLIED MATHEMATICS REPORT AMR4/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR J. Dick and F. Pillichshammer January, 24 Diaphony, discrepancy, spectral test and worst-case error Josef Dick
More informationDigital lattice rules
Digital lattice rules Multivariate integration and discrepancy estimates A thesis presented to The University of New South Wales in fulfilment of the thesis requirement for the degree of Doctor of Philosophy
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca October 22nd, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationDelsarte-Yudin LP method and Universal Lower Bound on energy
Peter Dragnev, IPFW Delsarte-Yudin LP method and Universal Lower Bound on energy Peter Dragnev Indiana University-Purdue University Fort Wayne Joint work with: P. Boyvalenkov (BAS); D. Hardin, E. Saff
More informationNumerical Integration over unit sphere by using spherical t-designs
Numerical Integration over unit sphere by using spherical t-designs Congpei An 1 1,Institute of Computational Sciences, Department of Mathematics, Jinan University Spherical Design and Numerical Analysis
More informationProblem Set 6: Solutions Math 201A: Fall a n x n,
Problem Set 6: Solutions Math 201A: Fall 2016 Problem 1. Is (x n ) n=0 a Schauder basis of C([0, 1])? No. If f(x) = a n x n, n=0 where the series converges uniformly on [0, 1], then f has a power series
More informationREAL AND COMPLEX ANALYSIS
REAL AND COMPLE ANALYSIS Third Edition Walter Rudin Professor of Mathematics University of Wisconsin, Madison Version 1.1 No rights reserved. Any part of this work can be reproduced or transmitted in any
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationMulti-normed spaces and multi-banach algebras. H. G. Dales. Leeds Semester
Multi-normed spaces and multi-banach algebras H. G. Dales Leeds Semester Leeds, 2 June 2010 1 Motivating problem Let G be a locally compact group, with group algebra L 1 (G). Theorem - B. E. Johnson, 1972
More informationA trigonometric orthogonality with respect to a nonnegative Borel measure
Filomat 6:4 01), 689 696 DOI 10.98/FIL104689M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A trigonometric orthogonality with
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationWeak and strong convergence theorems of modified SP-iterations for generalized asymptotically quasi-nonexpansive mappings
Mathematica Moravica Vol. 20:1 (2016), 125 144 Weak and strong convergence theorems of modified SP-iterations for generalized asymptotically quasi-nonexpansive mappings G.S. Saluja Abstract. The aim of
More informationFinite Frames and Sigma Delta Quantization
Finite Frames and Sigma Delta Quantization with Matt Fickus Alex Powell and Özgür Yilmaz Aram Tangboondouangjit Finite Frames and Sigma Delta Quantization p.1/57 Frames Frames F = {e n } N n=1 for d-dimensional
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationNEEDLET APPROXIMATION FOR ISOTROPIC RANDOM FIELDS ON THE SPHERE
NEEDLET APPROXIMATION FOR ISOTROPIC RANDOM FIELDS ON THE SPHERE Quoc T. Le Gia Ian H. Sloan Yu Guang Wang Robert S. Womersley School of Mathematics and Statistics University of New South Wales, Australia
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More informationConstruction algorithms for plane nets in base b
Construction algorithms for plane nets in base b Gunther Leobacher, Friedrich Pillichshammer and Thomas Schell arxiv:1506.03201v1 [math.nt] 10 Jun 2015 Abstract The class of (0, m, s)-nets in base b has
More informationLow Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed?
Low Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed? Xiaoqun Wang 1,2 and Ian H. Sloan 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084,
More informationSTRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE CARLO ALGORITHMS
MATHEMATICS OF COMPUTATION Volume 72, Number 242, Pages 823 838 S 0025-5718(02)01440-0 Article electronically published on June 25, 2002 STRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE
More informationSparse grid quadrature on products of spheres
Sparse grid quadrature on products of spheres Markus Hegland a, Paul Leopardi a, a Centre for Mathematics and its Applications, Australian National University. Abstract This paper examines sparse grid
More informationSHRINKING PROJECTION METHOD FOR A SEQUENCE OF RELATIVELY QUASI-NONEXPANSIVE MULTIVALUED MAPPINGS AND EQUILIBRIUM PROBLEM IN BANACH SPACES
U.P.B. Sci. Bull., Series A, Vol. 76, Iss. 2, 2014 ISSN 1223-7027 SHRINKING PROJECTION METHOD FOR A SEQUENCE OF RELATIVELY QUASI-NONEXPANSIVE MULTIVALUED MAPPINGS AND EQUILIBRIUM PROBLEM IN BANACH SPACES
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationTHROUGHOUT this paper, we let C be a nonempty
Strong Convergence Theorems of Multivalued Nonexpansive Mappings and Maximal Monotone Operators in Banach Spaces Kriengsak Wattanawitoon, Uamporn Witthayarat and Poom Kumam Abstract In this paper, we prove
More informationViscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces
Viscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces YUAN-HENG WANG Zhejiang Normal University Department of Mathematics Yingbing Road 688, 321004 Jinhua
More informationConvergence of greedy approximation I. General systems
STUDIA MATHEMATICA 159 (1) (2003) Convergence of greedy approximation I. General systems by S. V. Konyagin (Moscow) and V. N. Temlyakov (Columbia, SC) Abstract. We consider convergence of thresholding
More informationLattice rules and sequences for non periodic smooth integrands
Lattice rules and sequences for non periodic smooth integrands dirk.nuyens@cs.kuleuven.be Department of Computer Science KU Leuven, Belgium Joint work with Josef Dick (UNSW) and Friedrich Pillichshammer
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationLocal radial basis function approximation on the sphere
Local radial basis function approximation on the sphere Kerstin Hesse and Q. T. Le Gia School of Mathematics and Statistics, The University of New South Wales, Sydney NSW 05, Australia April 30, 007 Abstract
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationConvergence to Common Fixed Point for Two Asymptotically Quasi-nonexpansive Mappings in the Intermediate Sense in Banach Spaces
Mathematica Moravica Vol. 19-1 2015, 33 48 Convergence to Common Fixed Point for Two Asymptotically Quasi-nonexpansive Mappings in the Intermediate Sense in Banach Spaces Gurucharan Singh Saluja Abstract.
More informationTHEOREMS, ETC., FOR MATH 515
THEOREMS, ETC., FOR MATH 515 Proposition 1 (=comment on page 17). If A is an algebra, then any finite union or finite intersection of sets in A is also in A. Proposition 2 (=Proposition 1.1). For every
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT
ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT Received: 31 July, 2008 Accepted: 06 February, 2009 Communicated by: SIMON J SMITH Department of Mathematics and Statistics La Trobe University,
More informationChapter One. The Calderón-Zygmund Theory I: Ellipticity
Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere
More informationMATH MEASURE THEORY AND FOURIER ANALYSIS. Contents
MATH 3969 - MEASURE THEORY AND FOURIER ANALYSIS ANDREW TULLOCH Contents 1. Measure Theory 2 1.1. Properties of Measures 3 1.2. Constructing σ-algebras and measures 3 1.3. Properties of the Lebesgue measure
More informationinto regions of equal measure and small diameter p. 1/2
into regions of equal measure and small diameter Paul Leopardi paulleopardi@unsweduau School of Mathematics, University of New South Wales For presentation at Vanderbilt University, Nashville, November
More information2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due 9/5). Prove that every countable set A is measurable and µ(a) = 0. 2 (Bonus). Let A consist of points (x, y) such that either x or y is
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationProgress in high-dimensional numerical integration and its application to stochastic optimization
Progress in high-dimensional numerical integration and its application to stochastic optimization W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Short course,
More informationStrong Convergence Theorems for Nonself I-Asymptotically Quasi-Nonexpansive Mappings 1
Applied Mathematical Sciences, Vol. 2, 2008, no. 19, 919-928 Strong Convergence Theorems for Nonself I-Asymptotically Quasi-Nonexpansive Mappings 1 Si-Sheng Yao Department of Mathematics, Kunming Teachers
More informationUniversal convex coverings
Bull. London Math. Soc. 41 (2009) 987 992 C 2009 London Mathematical Society doi:10.1112/blms/bdp076 Universal convex coverings Roland Bacher Abstract In every dimension d 1, we establish the existence
More informationInverse Function Theorem writeup for MATH 4604 (Advanced Calculus II) Spring 2015
Inverse Function Theorem writeup for MATH 4604 (Advanced Calculus II) Spring 2015 Definition: N := {1, 2, 3,...}. Definition: For all n N, let i n : R n R n denote the identity map, defined by i n (x)
More informationQMC methods in quantitative finance. and perspectives
, tradition and perspectives Johannes Kepler University Linz (JKU) WU Research Seminar What is the content of the talk Valuation of financial derivatives in stochastic market models using (QMC-)simulation
More informationL p Functions. Given a measure space (X, µ) and a real number p [1, ), recall that the L p -norm of a measurable function f : X R is defined by
L p Functions Given a measure space (, µ) and a real number p [, ), recall that the L p -norm of a measurable function f : R is defined by f p = ( ) /p f p dµ Note that the L p -norm of a function f may
More informationLecture 5. Ch. 5, Norms for vectors and matrices. Norms for vectors and matrices Why?
KTH ROYAL INSTITUTE OF TECHNOLOGY Norms for vectors and matrices Why? Lecture 5 Ch. 5, Norms for vectors and matrices Emil Björnson/Magnus Jansson/Mats Bengtsson April 27, 2016 Problem: Measure size of
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationKähler configurations of points
Kähler configurations of points Simon Salamon Oxford, 22 May 2017 The Hesse configuration 1/24 Let ω = e 2πi/3. Consider the nine points [0, 1, 1] [0, 1, ω] [0, 1, ω 2 ] [1, 0, 1] [1, 0, ω] [1, 0, ω 2
More informationShih-sen Chang, Yeol Je Cho, and Haiyun Zhou
J. Korean Math. Soc. 38 (2001), No. 6, pp. 1245 1260 DEMI-CLOSED PRINCIPLE AND WEAK CONVERGENCE PROBLEMS FOR ASYMPTOTICALLY NONEXPANSIVE MAPPINGS Shih-sen Chang, Yeol Je Cho, and Haiyun Zhou Abstract.
More informationDiscrepancy Bounds for Nets and Sequences
Discrepancy Bounds for Nets and Sequences Peter Kritzer Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences Linz, Austria Joint work with Henri Faure (Marseille)
More informationf(x n ) [0,1[ s f(x) dx.
ACTA ARITHMETICA LXXX.2 (1997) Dyadic diaphony by Peter Hellekalek and Hannes Leeb (Salzburg) 1. Introduction. Diaphony (see Zinterhof [13] and Kuipers and Niederreiter [6, Exercise 5.27, p. 162]) is a
More informationLARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS. S. G. Bobkov and F. L. Nazarov. September 25, 2011
LARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS S. G. Bobkov and F. L. Nazarov September 25, 20 Abstract We study large deviations of linear functionals on an isotropic
More informationLax Solution Part 4. October 27, 2016
Lax Solution Part 4 www.mathtuition88.com October 27, 2016 Textbook: Functional Analysis by Peter D. Lax Exercises: Ch 16: Q2 4. Ch 21: Q1, 2, 9, 10. Ch 28: 1, 5, 9, 10. 1 Chapter 16 Exercise 2 Let h =
More informationERRATA: Probabilistic Techniques in Analysis
ERRATA: Probabilistic Techniques in Analysis ERRATA 1 Updated April 25, 26 Page 3, line 13. A 1,..., A n are independent if P(A i1 A ij ) = P(A 1 ) P(A ij ) for every subset {i 1,..., i j } of {1,...,
More informationDeviation Measures and Normals of Convex Bodies
Beiträge zur Algebra und Geometrie Contributions to Algebra Geometry Volume 45 (2004), No. 1, 155-167. Deviation Measures Normals of Convex Bodies Dedicated to Professor August Florian on the occasion
More informationChapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.
Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space
More informationMath 321 Final Examination April 1995 Notation used in this exam: N. (1) S N (f,x) = f(t)e int dt e inx.
Math 321 Final Examination April 1995 Notation used in this exam: N 1 π (1) S N (f,x) = f(t)e int dt e inx. 2π n= N π (2) C(X, R) is the space of bounded real-valued functions on the metric space X, equipped
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationApproximate Fekete points and discrete Leja points based on equal area partitions of the unit sphere
Approximate Fekete points and discrete Leja points based on equal area partitions of the unit sphere Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at Computational
More informationEhrhart polynomial for lattice squares, cubes, and hypercubes
Ehrhart polynomial for lattice squares, cubes, and hypercubes Eugen J. Ionascu UWG, REU, July 10th, 2015 math@ejionascu.ro, www.ejionascu.ro 1 Abstract We are investigating the problem of constructing
More informationFunctional Analysis I
Functional Analysis I Course Notes by Stefan Richter Transcribed and Annotated by Gregory Zitelli Polar Decomposition Definition. An operator W B(H) is called a partial isometry if W x = X for all x (ker
More informationHAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM
Georgian Mathematical Journal Volume 9 (2002), Number 3, 591 600 NONEXPANSIVE MAPPINGS AND ITERATIVE METHODS IN UNIFORMLY CONVEX BANACH SPACES HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM
More informationOn the existence of higher order polynomial. lattices with large figure of merit ϱ α.
On the existence of higher order polynomial lattices with large figure of merit Josef Dick a, Peter Kritzer b, Friedrich Pillichshammer c, Wolfgang Ch. Schmid b, a School of Mathematics, University of
More informationAnother Low-Technology Estimate in Convex Geometry
Convex Geometric Analysis MSRI Publications Volume 34, 1998 Another Low-Technology Estimate in Convex Geometry GREG KUPERBERG Abstract. We give a short argument that for some C > 0, every n- dimensional
More informationASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE VOLTERRA EQUATIONS. Janusz Migda and Małgorzata Migda
Opuscula Math. 36, no. 2 (2016), 265 278 http://dx.doi.org/10.7494/opmath.2016.36.2.265 Opuscula Mathematica ASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE VOLTERRA EQUATIONS Janusz Migda and Małgorzata
More informationDuality of multiparameter Hardy spaces H p on spaces of homogeneous type
Duality of multiparameter Hardy spaces H p on spaces of homogeneous type Yongsheng Han, Ji Li, and Guozhen Lu Department of Mathematics Vanderbilt University Nashville, TN Internet Analysis Seminar 2012
More informationCHAPTER VIII HILBERT SPACES
CHAPTER VIII HILBERT SPACES DEFINITION Let X and Y be two complex vector spaces. A map T : X Y is called a conjugate-linear transformation if it is a reallinear transformation from X into Y, and if T (λx)
More informationGiovanni Migliorati. MATHICSE-CSQI, École Polytechnique Fédérale de Lausanne
Analysis of the stability and accuracy of multivariate polynomial approximation by discrete least squares with evaluations in random or low-discrepancy point sets Giovanni Migliorati MATHICSE-CSQI, École
More informationDyadic diaphony of digital sequences
Dyadic diaphony of digital sequences Friedrich Pillichshammer Abstract The dyadic diaphony is a quantitative measure for the irregularity of distribution of a sequence in the unit cube In this paper we
More informationApplication of QMC methods to PDEs with random coefficients
Application of QMC methods to PDEs with random coefficients a survey of analysis and implementation Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with Ivan Graham
More informationAnalysis-3 lecture schemes
Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space
More informationl(y j ) = 0 for all y j (1)
Problem 1. The closed linear span of a subset {y j } of a normed vector space is defined as the intersection of all closed subspaces containing all y j and thus the smallest such subspace. 1 Show that
More informationThe Heine-Borel and Arzela-Ascoli Theorems
The Heine-Borel and Arzela-Ascoli Theorems David Jekel October 29, 2016 This paper explains two important results about compactness, the Heine- Borel theorem and the Arzela-Ascoli theorem. We prove them
More information12.0 Properties of orthogonal polynomials
12.0 Properties of orthogonal polynomials In this section we study orthogonal polynomials to use them for the construction of quadrature formulas investigate projections on polynomial spaces and their
More informationStable Process. 2. Multivariate Stable Distributions. July, 2006
Stable Process 2. Multivariate Stable Distributions July, 2006 1. Stable random vectors. 2. Characteristic functions. 3. Strictly stable and symmetric stable random vectors. 4. Sub-Gaussian random vectors.
More informationMonte Carlo Method for Numerical Integration based on Sobol s Sequences
Monte Carlo Method for Numerical Integration based on Sobol s Sequences Ivan Dimov and Rayna Georgieva Institute of Information and Communication Technologies, Bulgarian Academy of Sciences, 1113 Sofia,
More informationp-adic Analysis Compared to Real Lecture 1
p-adic Analysis Compared to Real Lecture 1 Felix Hensel, Waltraud Lederle, Simone Montemezzani October 12, 2011 1 Normed Fields & non-archimedean Norms Definition 1.1. A metric on a non-empty set X is
More informationDISCRETE MINIMAL ENERGY PROBLEMS
DISCRETE MINIMAL ENERGY PROBLEMS Lecture II E. B. Saff Center for Constructive Approximation Vanderbilt University University of Crete, Heraklion May, 2017 Recall Notation ω N = {x 1,..., x N } A, A compact,
More informationExercise Solutions to Functional Analysis
Exercise Solutions to Functional Analysis Note: References refer to M. Schechter, Principles of Functional Analysis Exersize that. Let φ,..., φ n be an orthonormal set in a Hilbert space H. Show n f n
More informationDecomposition of Riesz frames and wavelets into a finite union of linearly independent sets
Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets Ole Christensen, Alexander M. Lindner Abstract We characterize Riesz frames and prove that every Riesz frame
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationPart III. Quasi Monte Carlo methods 146/349
Part III Quasi Monte Carlo methods 46/349 Outline Quasi Monte Carlo methods 47/349 Quasi Monte Carlo methods Let y = (y,...,y N ) be a vector of independent uniform random variables in Γ = [0,] N, u(y)
More informationA Comparison of Popular Point Configurations on 2
Volume 9 016 Pages 16 49 A Comparison of Popular Point Configurations on D. P. Hardin a T. Michaels ab E.B. Saff a Abstract There are many ways to generate a set of nodes on the sphere for use in a variety
More informationReal Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi
Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.
More information