Algorithms for Scientific Computing
|
|
- Charla Stephens
- 5 years ago
- Views:
Transcription
1 Algorithms for Scientific Computing Hierarchical Methods and Sparse Grids d-dimensional Hierarchical Basis Michael Bader Technical University of Munich Summer 208
2 Intermezzo/ Big Picture : Archimedes Quadrature Start with 2d example (compare tutorials): f := 6x (x )x 2 (x 2 ), Ω = [0, ] 2 f Ω = 0 Consider hierarchical surplus at grid points with n = 3, h 3 = Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 2
3 Big Picture : Archimedes Quadrature (2) Ω f d x = 4/9 = 0. 4 = = Consider volume of subvolumes (pagodas) for quadrature Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 3
4 Big Picture : Archimedes Quadrature (3) What, if we leave out (adaptively) all subvolumes with volume < ɛ =? 49 grid points (full grid) 7 grid points (sparse grid) 0 0 Approximation of volume: = = Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 4
5 Part I Hierarchical Decomposition, d-dimensional Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 5
6 Recall: Archimedes and Hierarchical Basis (in D) Archimedes Quadrature: Hierarchical Basis: t= t=2 ¼ 0 ½ ¼. x, Φ, Φ 2, Φ 2,3. x 2, x 2,3 0 ½ 0 ½ Φ 3, Φ 3,3 Φ 3,5 Φ 3,7. x 3, x 3,3 x 3,5 x 3,7 use nodal basis functions φ l,i with I l := {i : i < 2 l, i odd} hierarchical basis Ψ n := n l= {φ l,i : i I l } hierachical function spaces W l := span {φ l,i : i I l } and V l = V l W l unique hierachical representation u = n l= w l = n l= i I l v l,i φ l,i size of surpluses v l,i roughly decays with 4 n for smooth functions Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 6
7 Hierarchical Decomposition Step by Step Now (and more formally), starting with d-dimensional hierarchical decompositions... Transfer from d = to d > Functions in multiple variables x = (x,..., x d ) Domain Ω := [0, ] d We consider only functions u which are 0 on Ω (on the edges of the square, faces of the cube,... ) Each hierarchical grid described by multi-index l = (l,..., l d ) N d Grids have different mesh-widths in different dimensions: h l := (h,..., h d ) := (2 l,..., 2 l d ) =: 2 l Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 7
8 Hierarchical Decomposition, d > Introducing further notation (which we ll need later on): Grid points (for function evaluations): x l, i = (i h l,..., i d h ld ) Comparisons of multi-indices component-wise: Two norms for multi-indices l l i lk i k, k =,..., d index sum: l := l l d maximum index: l := max { l,..., l d } Note: taking the absolute values,, for l k N is not necessary, but is part of the usual definition of and Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 8
9 Practicing Identifiers l, h l, x l, i l = l =2 l =3 l l 2 = l 2 =2 l 2 =3 l 2 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 9
10 Piecewise d-linear Functions Suitable generalization of piecewise linear functions Piecewise d-linear functions w.r.t. h l grid If you fix d coordinates, they are linear in remaining x j V l : space of all functions for given l Alternative point of view: Define suitable basis Φ l Regard V l as span of Φ l d-dimensional basis functions: products of one-dimensional hat functions: φ l, i ( x) = d φ lj,i j (x j ) = φ l,i (x ) φ l2,i 2 (x 2 )... φ ld,i d (x d ) j= Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 0
11 d-dimensional Basis Functions Basis functions are pagoda functions (not pyramids!) Examples: φ (,),(,), and φ (2,3),(3,5) : x x x x Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208
12 Function Spaces V l and V n Basis for space of piecewise linear functions w.r.t. h l grid Φ l := {φ l, i, i < 2 l } Function space V l := span{φ l } with dim V l = (2 l )... (2 l d ) O(2 l ) Special case l =... = l d function space denoted as V n : V n := V (n,...,n) Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 2
13 Hierarchical Increments W l Analogous to d: Omit grid points with even index (exist on coarser grid) Now in all directions I l := { i : i < 2 l, all i j odd} Hierarchical increment spaces l 2 = l 2 =2 l = l =2 l =3 l W l := span{φ l, i } i I l contain all functions of V l that vanish at all grid points of all coarser grids l 2 =3 l 2 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 3
14 Hierarchical Increments W l vs. Nodal Basis Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 4
15 Hierarchical Subspace Decomposition For l N d we obtain a unique representation of each u V l as u = l l w l with w l W l Representation in the hierarchical basis u = w l = v l, i φ l, i l l l l i I l with d-dimensional hierarchical surpluses v l, i Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 5
16 Determining the Hierarchical Surpluses in 2D We now compute the hierarchical surpluses v l, i for some u V n = V (n,n) : u( x) = φ l, i Φ (n,n) u(x l, i ) φ l, i ( x) = 2 n i = 2 n i 2 = First step Hierarchization in x -direction (fix x 2 and employ d hierarchization in x -direction): u( x) = n 2 n l = i I l i 2 = u(x n, i ) φ n,i (x )φ n,i2 (x 2 ) v l,i (x n,i2 ) φ l,i (x ) φ n,i2 (x 2 ) with d surplus (still depending on x 2, evaluated at all x 2 = x n,i2 ) v l,i (x 2 ) = u(x l,i, x 2 ) u(x l,i, x 2 ) + u(x l,i +, x 2 ) 2 Note: the indices i ± of the grid points x l,i and x l,i + are even, such that the corresponding hierarchical basis functions belong to a parent/ancestor level. Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 6
17 Determining the Hierarchical Surpluses in 2D (2) A bit more intuitive: We mark the grid points of the corresponding ansatz functions we use (before and after) l = l =2 l =3 l l = l =2 l =3 l l 2 = l 2 = l 2 =2 l 2 =2 l 2 =3 l 2 =3 l 2 l 2 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 7
18 Determining the Hierarchical Surpluses in 2D (3) Second step Hierarchize every v l,i (x 2 ) (separately) in x 2 dimension:: n n u( x) = v (l,l 2 ),(i,i 2 ) φ l,i (x ) φ l2,i 2 (x 2 ) l = l 2 = i 2 I l2 i I l l = l =2 l =3 l l = l =2 l =3 l l 2 = l 2 = l 2 =2 l 2 =2 l 2 =3 l 2 =3 l 2 l 2 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 8
19 Determining the Hierarchical Surpluses (the general d-dimensional case) Now: compute the d-dim. hierarchical surpluses v l, i for some u( x) V n : u( x) = u(x l, i ) φ l, i ( x) = u(x l, i ) φ l,i (x )... φ ld,i d (x d ) φ l, i Φ (n,...,n) φ l, i Φ (n,...,n) First step Hierarchization in x d -direction (fix x,..., x d and employ d hierarchization): n u = v ld,i d (x n,(i,...,i d )) φ ld,i d (x d ) φ l,i (x )... φ ld,i d (x d ) φ l, i Φ (n,...,n) l d = i d I ld with d surplus evaluted at (x,..., x d ) = x n,(i,...,i d ): v ld,i d (x,..., x d ) = u(x,..., x d, x ld,i d ) u(x,..., x d, x ld,i d ) + u(x,..., x d, x ld,i d +) 2 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 9
20 Determining the Hierarchical Surpluses (the general d-dimensional case, second step) Second step Hierarchize every v ld,i d : R d R (separately) in its first argument: n n ( u( x) = l d = i d I ld l d = i d I ld φ l, Φ i (n,...,n) ) v ld,i d (x n,(i,...,i d 2 )) φ ld,i d (x d ) φ ld,i d (x d ) φ l,i (x )... φ ld 2,i d 2 (x d 2 ) Steps 3 to d All steps correspondingly for each remaining dimension Afterwards we have computed surpluses v l, i (functions in zero parameters / scalar values) u( x) = v l, i φ l,i (x )... φ ld,i d (x d ) = v l, i φ l, i ( x) = l i I l l i I l where l i I l is short for n l d = i d I ld n l = i I l w l ( x) Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer l
21 Comparison 2D with Wavelet Transform Apply Level-Wise Hierarchisation First level: First step: split into nodal basis and hierarchical surpluses (st argument) Second step: split into nodal basis and hierarchical surpluses (2nd argument) l = l =2 l =3 l l = l =2 l =3 l l 2 = l 2 = l 2 =2 l 2 =2 l 2 =3 l 2 =3 l 2 l 2 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 2
22 Comparison 2D with Wavelet Transform (2) Apply Level-Wise Hierarchisation Second level: First step: split into nodal basis and hierarchical surpluses (st argument) Second step: split into nodal basis and hierarchical surpluses (2nd argument) l = l =2 l =3 l l = l =2 l =3 l l 2 = l 2 = l 2 =2 l 2 =2 l 2 =3 l 2 =3 l 2 l 2 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
23 Part II Hierarchical Decomposition Outlook on Cost and Accuracy Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
24 Analysis of Hierarchical Decomposition Contribution of summands in hierarchical decomposition in D: n n u = w l = v l,i φ l,i l= l= i I l in dd: u = l w l = l start analysis in univariate setting i I l v l, i φ l, i ( x) and port to mulitvariate setting Cost/benefit analysis quantifies reduction of effort Need several norms to measure w l Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
25 Norms of Functions As always, we assume sufficiently smooth functions u : [0, ] R, then: Maximum norm L 2 norm u := max x [0,] u(x) for the L 2 scalar product u 2 := 0 u(x) 2 dx, (u, v) 2 := 0 u(x)v(x) dx Energy norm u E := u 2 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
26 Norms of Basis Functions For the basis functions φ l,i, we obtain φ l,i = φ l,i 2 = φ l,i E = 2hl 3 2 h l h -2 Φ Φ 2 (Φ')2... x l,i- x l,i x l,i+ x l,i- x l,i x l,i+ x l,i- x l,i x l,i+ Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
27 Estimation of Surpluses Consider surplus v l,i of basis function φ l,i : u two times differentiable v l,i := u(x l,i ) 2 (u(x l,i ) + u(x l,i+ )) We can then write v l,i as (see separate proof) v l,i = 0 ψ l,i (x)u (x) dx with ψ l,i := h l 2 φ l,i v l,i depends on u, thus we define for future use µ 2 (u) := u 2 and µ (u) := u. note: µ 2 (u) and µ (u) are properties of the function u Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
28 Estimation of Surplusses (2) With integral representation v l,i = 0 h l 2 φ l,i(x)u (x) dx, we can bound ( ) v l,i h l 2 φ l,i dx µ (u) = h2 l 2 µ (u) O(hl 2 ) 0 and, via Cauchy-Schwartz inequality (u, v) u v, v l,i h l 2 φ hl 3 l,i 2 µ 2 (u Ti ) = 6 µ 2(u Ti ), where u Ti restricts u to the support T i = [x l,i, x l,i+ ] of φ l,i Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
29 Estimation of w l Estimate contribution of entire level l in hierarchical decomposition of u, i.e. w l = i I l v l,i φ l,i Use that supports of φ l,i are pairwise disjoint Maximum norm w l h2 l 2 µ (u) O(h 2 l ) L 2 norm w l 2 2 = v l,i 2 φ l,i 2 2 h3 l 6 2h l 3 µ 2 (u Ti ) 2 = h4 l 9 µ 2(u) 2 i I l i I l w l 2 O(h 2 l ) Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
30 Estimation of w l (2) Energy norm w l 2 E = i I l v l,i 2 φ l,i 2 E = i I l v l,i 2 2 h l 2 h4 l h l 4 µ (u) 2 = h2 l 2h l 4 µ (u) 2 (2 l = /(2h l ) summands) w l E O(h l ) Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
31 Estimation of w l (3) We can write u (twice differentiable) as infinite series u = l= w l Convergent in all three norms Approximation error given as u u n := u in maximum and L 2 norm: O(h 2 n) in energy norm: O(h n ) n w l = l= l=n+ w l Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 3
32 Towards d Dimensions: Norms of φ l, i Estimating the w l will enable us to select those subspaces that contribute most to overall solution (best cost-benefit ratios) Same procedure as for d =, but slightly more complicated functions Start with norms Maximum norm: L 2 norm: φ l, i 2 := [0,] d φ l, i ( x) 2 d x = φ l, i := max φ l, i ( x) = x [0,] d d ) d d (2 ) d φ lj,i j 2 = ( 2 h lj = l j= j= Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
33 Norms of φ l, i (2) Energy norm (defined as L 2 norm of the Euclidean norm of the gradient φ l, i ): φ l, i E := φ l, i ( x) φ l, i ( x) d x =... = [0,] d ( ) d d = 2 h... h d 2 3 h 2 (here always: h j := h lj ) j= j ( ) d = 2 d l 2 2l j For the two-dimensional settings (d = 2), we obtain ( 4 h φ l, i E = + h ) 2 3 h 2 h j= Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
34 Estimation of Surpluses Hierarchical surpluses now depend on mixed 2nd derivatives If we define ψ l, i := 2d 2d u u := x 2... x d 2 d d ψ lj,i j = j= j= h j 2 φ l, i = ( ) d 2 l d φ l, i we can derive an integral representation similar to D: v l, i = ψ l, i 2d u d x [0,] d (Proof: Fubini s theorem and d integral representation) Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
35 Estimation of Surpluses (2) We define (correspondingly to d) We can thus bound v l, i as ( d v l, i j= µ 2 (u) := 2d u 2 and µ (u) := 2d u ) ( ) ( d h j h 2 ) j φ l, 2 i d x µ (u) = µ (u) = 2 2 l d µ (u) 2 [0,] d j= and v l, i ( d j= ) h j h 3 φ l, 2 i 2 µ 2 (u T i ) =... h3 d 6 d µ 2 (u T i ) = ( 6) d/2 2 3 l /2 µ 2 (u T i ) Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
36 Estimation of w l Obtain estimates for w l in subspace W l analogously as in d: Make use of the fact that supports of basis functions for a grid are disjoint (apart from the boundaries) Maximum norm L 2 norm d w l j= d w l 2 Energy norm w l E ( ) d d 4 2 j= j= h 2 j 2 h 2 j 3 h 4... h4 d h 2 j µ (u) = 2 2 l d µ (u), µ 2 (u) = 3 d 2 2 l µ 2 (u), µ (u) = ( ) d d l 2 2l j µ (u) j= Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
37 Analysis of Cost-Benefit Ratio Consider not individual basis functions, but whole hierarchical increments From the tableau of subspaces, select those subspaces that minimize the cost, or maximize the benefit respectively, for u : [0, ] d R (u sufficiently often differentiable) Cost Measure cost in number of grid points ( coefficients ) c( l) = I l = 2 l d Benefit How to measure benefit? interpolation error Let L N d be the set of indices of all selected grids, then u L := l L w l and u u L = l L w l Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
38 Analysis of Cost-Benefit Ratio (2) For each component w l, we have derived bounds of the type w l s( l) µ(u) with s( l) = 2 d 2 2 l or s( l) = 3 d 2 2 l and appropriate indices for norm and µ We obtain u u L ( ) w l s( l) µ(u) l L l L = ( l N d ) ( ) s( l) s( l) µ(u) l L st factor depends only on selected subspaces, 2nd factor only on u Justifies to interpret s( l) as benefit/contribution of subspace W l Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
39 Quality of Approximation of Full Grid V n Examine cost c( l) and benefit s( l) for full grid Regular grid with mesh-width h = 2 n in each direction (full grid) for function space V n total cost dim V n O(2 dn ), or dim V n O(h d ) Considered subset of hierarchical increments: L n := { l : l n}. Bounds in L 2 and maximum norm involve factor s( l) = C 2 2 l In the following estimation, leave out l-independent factor C can be appended to the estimate in the end Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
40 Quality of Approximation of Full Grid V n (2) We can estimate s( l) = 2 2 l = l Ln l Ln = ( 4 n 4 4 n l = ( n n ) d 2 2(l + +l n) = 2 2k l d = k= ) d ( = d ( 3) ) 2 2n d ( ) d ( ) 3 d 2 2n using ( ɛ) d dɛ for 0 ɛ and d N For n we obtain ( ) d s( l) = and thus s( 3 l) = l N d l Ln ( ) d d 2 2n 3 Leads to bounds for the approximation error in L 2 - and maximum norm u u Ln C s( l) C d 3 d 2 2n O(hn) 2 l Ln with constant C (independent of n) Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
41 Sparse Grids Final steps to high-dimensional numerics Consider sum of benefits/contributions (for L 2 and maximum norm) l Ln 2 2 l Equal benefit of hierarchical increments W l for constant l Same for cost c( l) = 2 l d (number of grid points of W l ) Constant cost-benefit ratio c( l)/s( l) for constant l Full grids? Quadratic extract of subspaces is not economical: We take large subgrids with low contribution We could have taken others with much higher contribution Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer 208 4
42 Part III Sparse Grids Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
43 Sparse Grids! cost-benefit analysis: equal contribution of hierarchical increments W l for constant l Best choice: Cut diagonally in tableau of subspaces: L n := { l : l n + d } Resulting sparse grid space Vn := l n+d W l l 2 = l 2 =2 l 2 =3 l 2 l = l =2 l =3 l Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
44 Sparse Grids! Diagonall cut in tableau of subspaces: L n := { l : l n + d } Resulting sparse grid space Vn := l n+d W l Sparse grid for d = 2 and overall level n = 5 Grid points x l, i of same cost/benefit ratio in same color Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
45 Sparse Grids Cost Number of grid points? For d = 2: dim V n = For d = 3: l n+ dim V n = dim W l = n k= l n+ k(k + ) 2 2 l 2 = n k 2 k = 2 n (n ) +, k= ( n 2 k = 2 n 2 2 n ) 2 +, Both in O(2 n n d ) Holds for general d as well (proof with some combinatorics) Expressed in terms of N = 2 n (max. points per dimension): O ( N(log N) d ) Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
46 Sparse Grids Cost (2) In numbers... Compare cost for full grid V n and sparse grid V n : d = 2: n dim V n = (2 n ) ,046,529 dim V n = 2 n (n ) ,27 Even more distinct for d = 3: n dim V n = (2 n ) ,375...,070,590,67 ( ) dim Vn = 2 n n ,03 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
47 Sparse Grids Cost (3)... and for overall level n = 5 in different dimensions d V 5 V , , ,629,5, ,503,68 2, ,52,64, 4, ,89,037,44 6, ,439,622,60,67 9, ,628,286,980,80 3,44 The higher the dimension, the higher the benefit of sparse grids! Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
48 Sparse Grids Examples Sparse Grids of overall level n = 6 in d = 2 and d = 3 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
49 Sparse Grids Accuracy Much fewer grid points much lower accuracy? Would force us to choose larger n to obtain similar accuracy (and spoil everything) Error in L 2 and maximum norm: Compute sum ( l = k + ): l L n s( l) = k=n+ And for d = 3 (with l = k + 2): l L n s( l) = k=n+3 k(k + ) 2 k 2 2(k+) = 2 2(k+2) = ( n 2 + ) 2 2n 9 ( n n ) 2 2n 27 Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
50 Sparse Grids Accuracy (2) In general, it can be shown Error of interpolation in L 2 and maximum norm is O(2 2n n d ) or, expressed in mesh size h := 2 n : O ( h ( ) 2 log d ) h Only polynomial (in n) factor worse than full grid with O(2 2n ) or, expressed in mesh size h := 2 n : O(h 2 ) Outlook on Energy norm: ( Algorithms for Scientific Computing II) Analysis is more complicated (lines through subspaces with similar s( l), and thus c( l)/s( l), are more complicated) Overall result even better: obtain accuracy of O(2 n ) with only O(2 n ) grid points no polynomial terms (of type n d ) left! Michael Bader Algorithms for Scientific Computing d-dim. Hierarchical Basis Summer
Algorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationMultilevel stochastic collocations with dimensionality reduction
Multilevel stochastic collocations with dimensionality reduction Ionut Farcas TUM, Chair of Scientific Computing in Computer Science (I5) 27.01.2017 Outline 1 Motivation 2 Theoretical background Uncertainty
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationRandomized Coordinate Descent Methods on Optimization Problems with Linearly Coupled Constraints
Randomized Coordinate Descent Methods on Optimization Problems with Linearly Coupled Constraints By I. Necoara, Y. Nesterov, and F. Glineur Lijun Xu Optimization Group Meeting November 27, 2012 Outline
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationLocally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem
56 Chapter 7 Locally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem Recall that C(X) is not a normed linear space when X is not compact. On the other hand we could use semi
More informationWe simply compute: for v = x i e i, bilinearity of B implies that Q B (v) = B(v, v) is given by xi x j B(e i, e j ) =
Math 395. Quadratic spaces over R 1. Algebraic preliminaries Let V be a vector space over a field F. Recall that a quadratic form on V is a map Q : V F such that Q(cv) = c 2 Q(v) for all v V and c F, and
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
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 information4th Preparation Sheet - Solutions
Prof. Dr. Rainer Dahlhaus Probability Theory Summer term 017 4th Preparation Sheet - Solutions Remark: Throughout the exercise sheet we use the two equivalent definitions of separability of a metric space
More informationANALYSIS QUALIFYING EXAM FALL 2017: SOLUTIONS. 1 cos(nx) lim. n 2 x 2. g n (x) = 1 cos(nx) n 2 x 2. x 2.
ANALYSIS QUALIFYING EXAM FALL 27: SOLUTIONS Problem. Determine, with justification, the it cos(nx) n 2 x 2 dx. Solution. For an integer n >, define g n : (, ) R by Also define g : (, ) R by g(x) = g n
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationLena Leitenmaier & Parikshit Upadhyaya KTH - Royal Institute of Technology
Data mining with sparse grids Lena Leitenmaier & Parikshit Upadhyaya KTH - Royal Institute of Technology Classification problem Given: training data d-dimensional data set (feature space); typically relatively
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationNotions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy
Banach Spaces These notes provide an introduction to Banach spaces, which are complete normed vector spaces. For the purposes of these notes, all vector spaces are assumed to be over the real numbers.
More informationMathematics Department Stanford University Math 61CM/DM Inner products
Mathematics Department Stanford University Math 61CM/DM Inner products Recall the definition of an inner product space; see Appendix A.8 of the textbook. Definition 1 An inner product space V is a vector
More information4. Numerical Quadrature. Where analytical abilities end... continued
4. Numerical Quadrature Where analytical abilities end... continued Where analytical abilities end... continued, November 30, 22 1 4.3. Extrapolation Increasing the Order Using Linear Combinations Once
More informationDefinition 1. A set V is a vector space over the scalar field F {R, C} iff. there are two operations defined on V, called vector addition
6 Vector Spaces with Inned Product Basis and Dimension Section Objective(s): Vector Spaces and Subspaces Linear (In)dependence Basis and Dimension Inner Product 6 Vector Spaces and Subspaces Definition
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationFunction Approximation
1 Function Approximation This is page i Printer: Opaque this 1.1 Introduction In this chapter we discuss approximating functional forms. Both in econometric and in numerical problems, the need for an approximating
More informationIowa State University. Instructor: Alex Roitershtein Summer Exam #1. Solutions. x u = 2 x v
Math 501 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 015 Exam #1 Solutions This is a take-home examination. The exam includes 8 questions.
More informationIntroduction to Signal Spaces
Introduction to Signal Spaces Selin Aviyente Department of Electrical and Computer Engineering Michigan State University January 12, 2010 Motivation Outline 1 Motivation 2 Vector Space 3 Inner Product
More informationFourier Transform & Sobolev Spaces
Fourier Transform & Sobolev Spaces Michael Reiter, Arthur Schuster Summer Term 2008 Abstract We introduce the concept of weak derivative that allows us to define new interesting Hilbert spaces the Sobolev
More informationMath 320-2: Midterm 2 Practice Solutions Northwestern University, Winter 2015
Math 30-: Midterm Practice Solutions Northwestern University, Winter 015 1. Give an example of each of the following. No justification is needed. (a) A metric on R with respect to which R is bounded. (b)
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationBest approximation in the 2-norm
Best approximation in the 2-norm Department of Mathematical Sciences, NTNU september 26th 2012 Vector space A real vector space V is a set with a 0 element and three operations: Addition: x, y V then x
More information(1)(a) V = 2n-dimensional vector space over a field F, (1)(b) B = non-degenerate symplectic form on V.
18.704 Supplementary Notes March 21, 2005 Maximal parabolic subgroups of symplectic groups These notes are intended as an outline for a long presentation to be made early in April. They describe certain
More informationFunctional Analysis Exercise Class
Functional Analysis Exercise Class Week 2 November 6 November Deadline to hand in the homeworks: your exercise class on week 9 November 13 November Exercises (1) Let X be the following space of piecewise
More informationDimension-adaptive Sparse Grids
Dimension-adaptive Sparse Grids Jörg Blank April 14, 2008 1 Data Mining Function Reconstruction 2 Sparse Grids Motivation Introduction 3 Combination technique Motivation Dimension-adaptive 4 Examples 5
More informationOverview. A Posteriori Error Estimates for the Biharmonic Equation. Variational Formulation and Discretization. The Biharmonic Equation
Overview A Posteriori rror stimates for the Biharmonic quation R Verfürth Fakultät für Mathematik Ruhr-Universität Bochum wwwruhr-uni-bochumde/num1 Milan / February 11th, 013 The Biharmonic quation Summary
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationMath 5520 Homework 2 Solutions
Math 552 Homework 2 Solutions March, 26. Consider the function fx) = 2x ) 3 if x, 3x ) 2 if < x 2. Determine for which k there holds f H k, 2). Find D α f for α k. Solution. We show that k = 2. The formulas
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
More informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More informationMath 205, Summer I, Week 3a (continued): Chapter 4, Sections 5 and 6. Week 3b. Chapter 4, [Sections 7], 8 and 9
Math 205, Summer I, 2016 Week 3a (continued): Chapter 4, Sections 5 and 6. Week 3b Chapter 4, [Sections 7], 8 and 9 4.5 Linear Dependence, Linear Independence 4.6 Bases and Dimension 4.7 Change of Basis,
More informationProblem Set 2: Solutions Math 201A: Fall 2016
Problem Set 2: s Math 201A: Fall 2016 Problem 1. (a) Prove that a closed subset of a complete metric space is complete. (b) Prove that a closed subset of a compact metric space is compact. (c) Prove that
More informationLecture 2: Linear Algebra Review
EE 227A: Convex Optimization and Applications January 19 Lecture 2: Linear Algebra Review Lecturer: Mert Pilanci Reading assignment: Appendix C of BV. Sections 2-6 of the web textbook 1 2.1 Vectors 2.1.1
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationLecture 6: The Pigeonhole Principle and Probability Spaces
Lecture 6: The Pigeonhole Principle and Probability Spaces Anup Rao January 17, 2018 We discuss the pigeonhole principle and probability spaces. Pigeonhole Principle The pigeonhole principle is an extremely
More informationCauchy s Theorem (rigorous) In this lecture, we will study a rigorous proof of Cauchy s Theorem. We start by considering the case of a triangle.
Cauchy s Theorem (rigorous) In this lecture, we will study a rigorous proof of Cauchy s Theorem. We start by considering the case of a triangle. Given a certain complex-valued analytic function f(z), for
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
More informationPart III. 10 Topological Space Basics. Topological Spaces
Part III 10 Topological Space Basics Topological Spaces Using the metric space results above as motivation we will axiomatize the notion of being an open set to more general settings. Definition 10.1.
More informationMid Term-1 : Practice problems
Mid Term-1 : Practice problems These problems are meant only to provide practice; they do not necessarily reflect the difficulty level of the problems in the exam. The actual exam problems are likely to
More informationMH 7500 THEOREMS. (iii) A = A; (iv) A B = A B. Theorem 5. If {A α : α Λ} is any collection of subsets of a space X, then
MH 7500 THEOREMS Definition. A topological space is an ordered pair (X, T ), where X is a set and T is a collection of subsets of X such that (i) T and X T ; (ii) U V T whenever U, V T ; (iii) U T whenever
More informationMath 328 Course Notes
Math 328 Course Notes Ian Robertson March 3, 2006 3 Properties of C[0, 1]: Sup-norm and Completeness In this chapter we are going to examine the vector space of all continuous functions defined on the
More informationLINEAR ALGEBRA: THEORY. Version: August 12,
LINEAR ALGEBRA: THEORY. Version: August 12, 2000 13 2 Basic concepts We will assume that the following concepts are known: Vector, column vector, row vector, transpose. Recall that x is a column vector,
More informationAugust 23, 2017 Let us measure everything that is measurable, and make measurable everything that is not yet so. Galileo Galilei. 1.
August 23, 2017 Let us measure everything that is measurable, and make measurable everything that is not yet so. Galileo Galilei 1. Vector spaces 1.1. Notations. x S denotes the fact that the element x
More informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More information1. Subspaces A subset M of Hilbert space H is a subspace of it is closed under the operation of forming linear combinations;i.e.,
Abstract Hilbert Space Results We have learned a little about the Hilbert spaces L U and and we have at least defined H 1 U and the scale of Hilbert spaces H p U. Now we are going to develop additional
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationSpanning and Independence Properties of Finite Frames
Chapter 1 Spanning and Independence Properties of Finite Frames Peter G. Casazza and Darrin Speegle Abstract The fundamental notion of frame theory is redundancy. It is this property which makes frames
More informationLecture 3: Haar MRA (Multi Resolution Analysis)
U U U WAVELETS AND MULTIRATE DIGITAL SIGNAL PROCESSING Lecture 3: Haar MRA (Multi Resolution Analysis) Prof.V.M.Gadre, EE, IIT Bombay 1 Introduction The underlying principle of wavelets is to capture incremental
More informationFunctional Analysis MATH and MATH M6202
Functional Analysis MATH 36202 and MATH M6202 1 Inner Product Spaces and Normed Spaces Inner Product Spaces Functional analysis involves studying vector spaces where we additionally have the notion of
More informationWarm-up Quantifiers and the harmonic series Sets Second warmup Induction Bijections. Writing more proofs. Misha Lavrov
Writing more proofs Misha Lavrov ARML Practice 3/16/2014 and 3/23/2014 Warm-up Using the quantifier notation on the reference sheet, and making any further definitions you need to, write the following:
More information3. Numerical Quadrature. Where analytical abilities end...
3. Numerical Quadrature Where analytical abilities end... Numerisches Programmieren, Hans-Joachim Bungartz page 1 of 32 3.1. Preliminary Remarks The Integration Problem Numerical quadrature denotes numerical
More informationDefining the Discrete Wavelet Transform (DWT)
Defining the Discrete Wavelet Transform (DWT) can formulate DWT via elegant pyramid algorithm defines W for non-haar wavelets (consistent with Haar) computes W = WX using O(N) multiplications brute force
More informationLecture Notes 1: Vector spaces
Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector
More information2 Metric Spaces Definitions Exotic Examples... 3
Contents 1 Vector Spaces and Norms 1 2 Metric Spaces 2 2.1 Definitions.......................................... 2 2.2 Exotic Examples...................................... 3 3 Topologies 4 3.1 Open Sets..........................................
More informationSome Background Material
Chapter 1 Some Background Material In the first chapter, we present a quick review of elementary - but important - material as a way of dipping our toes in the water. This chapter also introduces important
More informationLecture 9. d N(0, 1). Now we fix n and think of a SRW on [0,1]. We take the k th step at time k n. and our increments are ± 1
Random Walks and Brownian Motion Tel Aviv University Spring 011 Lecture date: May 0, 011 Lecture 9 Instructor: Ron Peled Scribe: Jonathan Hermon In today s lecture we present the Brownian motion (BM).
More informationVector spaces. DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis.
Vector spaces DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Vector space Consists of: A set V A scalar
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 (Many figures from C. M. Bishop, "Pattern Recognition and ") 1of 254 Part V
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationDO NOT OPEN THIS QUESTION BOOKLET UNTIL YOU ARE TOLD TO DO SO
QUESTION BOOKLET EECS 227A Fall 2009 Midterm Tuesday, Ocotober 20, 11:10-12:30pm DO NOT OPEN THIS QUESTION BOOKLET UNTIL YOU ARE TOLD TO DO SO You have 80 minutes to complete the midterm. The midterm consists
More informationLecture 2: Vector Spaces, Metric Spaces
CCS Discrete II Professor: Padraic Bartlett Lecture 2: Vector Spaces, Metric Spaces Week 2 UCSB 2015 1 Vector Spaces, Informally The two vector spaces 1 you re probably the most used to working with, from
More informationDecoupling course outline Decoupling theory is a recent development in Fourier analysis with applications in partial differential equations and
Decoupling course outline Decoupling theory is a recent development in Fourier analysis with applications in partial differential equations and analytic number theory. It studies the interference patterns
More informationFitting multidimensional data using gradient penalties and combination techniques
1 Fitting multidimensional data using gradient penalties and combination techniques Jochen Garcke and Markus Hegland Mathematical Sciences Institute Australian National University Canberra ACT 0200 jochen.garcke,markus.hegland@anu.edu.au
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationDot Products. K. Behrend. April 3, Abstract A short review of some basic facts on the dot product. Projections. The spectral theorem.
Dot Products K. Behrend April 3, 008 Abstract A short review of some basic facts on the dot product. Projections. The spectral theorem. Contents The dot product 3. Length of a vector........................
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationAccuracy, Precision and Efficiency in Sparse Grids
John, Information Technology Department, Virginia Tech.... http://people.sc.fsu.edu/ jburkardt/presentations/ sandia 2009.pdf... Computer Science Research Institute, Sandia National Laboratory, 23 July
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More informationThe Finite Element Method for the Wave Equation
The Finite Element Method for the Wave Equation 1 The Wave Equation We consider the scalar wave equation modelling acoustic wave propagation in a bounded domain 3, with boundary Γ : 1 2 u c(x) 2 u 0, in
More information(K + L)(c x) = K(c x) + L(c x) (def of K + L) = K( x) + K( y) + L( x) + L( y) (K, L are linear) = (K L)( x) + (K L)( y).
Exercise 71 We have L( x) = x 1 L( v 1 ) + x 2 L( v 2 ) + + x n L( v n ) n = x i (a 1i w 1 + a 2i w 2 + + a mi w m ) i=1 ( n ) ( n ) ( n ) = x i a 1i w 1 + x i a 2i w 2 + + x i a mi w m i=1 Therefore y
More informationFinite Element Methods
Solving Operator Equations Via Minimization We start with several definitions. Definition. Let V be an inner product space. A linear operator L: D V V is said to be positive definite if v, Lv > for every
More informationLECTURE NOTE #11 PROF. ALAN YUILLE
LECTURE NOTE #11 PROF. ALAN YUILLE 1. NonLinear Dimension Reduction Spectral Methods. The basic idea is to assume that the data lies on a manifold/surface in D-dimensional space, see figure (1) Perform
More informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More information3 Orthogonality and Fourier series
3 Orthogonality and Fourier series We now turn to the concept of orthogonality which is a key concept in inner product spaces and Hilbert spaces. We start with some basic definitions. Definition 3.1. Let
More informationA Generic Bound on Cycles in Two-Player Games
A Generic Bound on Cycles in Two-Player Games David S. Ahn February 006 Abstract We provide a bound on the size of simultaneous best response cycles for generic finite two-player games. The bound shows
More information1 Adjacency matrix and eigenvalues
CSC 5170: Theory of Computational Complexity Lecture 7 The Chinese University of Hong Kong 1 March 2010 Our objective of study today is the random walk algorithm for deciding if two vertices in an undirected
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationIn class midterm Exam - Answer key
Fall 2013 In class midterm Exam - Answer key ARE211 Problem 1 (20 points). Metrics: Let B be the set of all sequences x = (x 1,x 2,...). Define d(x,y) = sup{ x i y i : i = 1,2,...}. a) Prove that d is
More informationInner Product Spaces 6.1 Length and Dot Product in R n
Inner Product Spaces 6.1 Length and Dot Product in R n Summer 2017 Goals We imitate the concept of length and angle between two vectors in R 2, R 3 to define the same in the n space R n. Main topics are:
More informationMinimal Element Interpolation in Functions of High-Dimension
Minimal Element Interpolation in Functions of High-Dimension J. D. Jakeman, A. Narayan, D. Xiu Department of Mathematics Purdue University West Lafayette, Indiana Random Model as a function We can consider
More informationSPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationLecture 13: Lower Bounds using the Adversary Method. 2 The Super-Basic Adversary Method [Amb02]
Quantum Computation (CMU 18-859BB, Fall 015) Lecture 13: Lower Bounds using the Adversary Method October 1, 015 Lecturer: Ryan O Donnell Scribe: Kumail Jaffer 1 Introduction There are a number of known
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 information