Scientific Computing II
|
|
- Sophia Thomas
- 5 years ago
- Views:
Transcription
1 Scientific Computing II Conjugate Gadient Methods Michael Bade Summe 2014 Conjugate Gadient Methods, Summe
2 Families of Iteative Solves elaxation methods: Jacobi-, Gauss-Seidel-Relaxation,... Ove-Relaxation-Methods Kylov methods: Steepest Descent, Conjugate Gadient,... GMRES,... Multilevel/Multigid methods, Domain Decomposition,... Conjugate Gadient Methods, Summe
3 Remembe: The Residual Equation fo Ax = b, we defined the esidual as: and the eo: e i) := x x i) leads to the esidual equation: i) = b Ax i) Ae i) = i) elaxation methods: solve a modified easie) SLE: B ê i) = i) whee B A multigid methods: coase-gid coection on esidual equation A H e i) H = i) H and x i+1) := x i) + IH h ei) H Conjugate Gadient Methods, Summe
4 Pat I Quadatic Foms and Steepest Descent Quadatic Foms Diection of Steepest Descent Steepest Descent Conjugate Gadient Methods, Summe
5 Quadatic Foms A quadatic fom is a scala, quadatic function of a vecto of the fom: f x) = 1 2 x T Ax b T x + c, whee A = A T y 2 x x y Conjugate Gadient Methods, Summe
6 Quadatic Foms 2) The gadient of a quadatic fom is defined as f x) = Ax b f x) = x 1 f x). x n f x) f x) = 0 Ax b = 0 Ax = b Ax = b equivalent to a minimisation poblem pope minimum only if A positive definite Conjugate Gadient Methods, Summe
7 Diection of Steepest Descent gadient f x): diection of steepest ascent f x) = Ax b = with esidual = b Ax) esidual : diection of steepest descent 4 y x Conjugate Gadient Methods, Summe
8 Solving SLE via Minimum Seach basic idea to find minimum: move into diection of steepest descent most simple scheme: x i+1) = x i) + α i) α constant Richadson iteation often consideed as a elaxation method) bette choice of α: move to lowest point in that diection Steepest Descent Conjugate Gadient Methods, Summe
9 Steepest Descent find an optimal α task: line seach along the line x 1) = x 0) + α 0) choose α such that f x 1) ) is minimal: use chain ule: α f x 1) ) = 0 α f x 1) ) = f x 1) ) T α x 1) = f x 1) ) T 0) emembe f x 1) ) = 1), thus: 1)) T 0)! = 0 hence, f x 1) ) = 1) should be othogonal to 0) Conjugate Gadient Methods, Summe
10 Steepest Descent find α 2) 1)) T 0) = b Ax 1)) T 0) ) T b Ax 0) + α 0) ) 0) b Ax 0)) T 0) α A 0)) T 0) 0)) T 0) α 0)) T A 0) = 0 = 0 = 0 = 0 Solve fo α: ) 0) T 0) α = ) 0) T A 0) Conjugate Gadient Methods, Summe
11 Steepest Descent Algoithm 4 1. i) = b Ax i) ) i) T i) 2. α i = ) i) T A i) 3. x i+1) = x i) + α i i) y 2 x Obsevations: slow convegence sim. to Jacobi elaxation) ) i e i) κ 1 e A 0) κ+1 A fo positive definite A: κ = λ max /λ min lagest/smallest eigenvalues of A) many steps in the same diection -6 Conjugate Gadient Methods, Summe
12 Pat II Conjugate Gadients Conjugate Diections A-Othogonality Conjugate Gadients A Miacle Occus... CG Algoithm Conjugate Gadient Methods, Summe
13 Conjugate Diections obsevation: Steepest Descent takes epeated steps in the same diection obvious idea: ty to do only one step in each diection possible appoach: choose othogonal seach diections d 0) d 1) d 2)... notice: eos othogonal to pevious diections: e 1) d 0), e 2) d 1) d 0),... Conjugate Gadient Methods, Summe
14 Conjugate Diections 2) compute α fom d 0)) T e 1) = d 0)) T e 0) αd 0)) = 0 equies popagation of the eo e 1) = x x 1) x 1) = x 0) + α i d 0) x x 1) = x x 0) α i d 0) e 1) = e 0) α i d 0) fomula fo α: ) d 0) T e 0) α = ) d 0) T d 0) but: we don t know e 0) Conjugate Gadient Methods, Summe
15 A-Othogonality make the seach diections A-othogonal: d i)) T Ad j) = 0 again: eos A-othogonal to pevious diections: e i+1)) T Ad i)! = 0 equiv. to minimisation in seach diection d i) : f x i+1)))t α f x i+1)) = α x i+1) = 0 i+1)) T d i) d i)) T Ae i+1) = 0 = 0 Conjugate Gadient Methods, Summe
16 A-Conjugate Diections emembe the fomula fo conjugate diections: ) d 0) T e 0) α = ) d 0) T d 0) same computation, but with A-othogonality: ) d i) T Ae i) α i = ) d i) T Ad i) ) d i) T i) = ) d i) T Ad i) fo the i-th iteation) these α i can be computed! still to do: find A-othogonal seach diections Conjugate Gadient Methods, Summe
17 A-Conjugate Diections 2) classical appoach to find othogonal diections conjugate Gam-Schmidt pocess: fom linealy independent vectos u 0), u 1),..., u i 1) constuct othogonal diections d 0), d 1),..., d i 1) i 1 d i) = u i) + β ik d k) k=0 β ik = ui) ) T Ad k) d k) ) T Ad k) needs to keep all old seach vectos in memoy On 3 ) computational complexity infeasible Conjugate Gadient Methods, Summe
18 Conjugate Gadients use esiduals i.e., u i) := i) ) to constuct conjugate diections: i 1 d i) = i) + β ik d k) k=0 new diection d i) should be A-othogonal to all d j) : 0 =! d i)) T Ad j) = i)) T i 1 Ad j) + β ik d k) ) T Ad j) all diections d k) fo k = 0,..., i 1) ae aleady A-othogonal and j < i), hence: 0 = i)) T Ad j) + β ij d j) ) ) T Ad j) i) T Ad j) β ij = ) d j) T Ad j) k=0 Conjugate Gadient Methods, Summe
19 Conjugate Gadients Status 1. conjugate diections and computation of α i : α i = d i) ) T i) d i) ) T Ad i) x i+1) = x i) + α i d i) 2. use esiduals to compute seach diections: i 1 d i) = i) + β ik d k) k=0 ) i) T Ad k) β ik = ) d k) T Ad k) still infeasible, as we need to stoe all vectos d k) Conjugate Gadient Methods, Summe
20 A Miacle Occus Pat 1 Two small contibutions: 1. popagation of the eo e i) = x x i) x i+1) = x i) + α i d i) x x i+1) = x x i) α i d i) e i+1) = e i) α i d i) we have used this once, aleady) 2. popagation of esiduals i+1) = Ae i+1) = A e i) α i d i)) i+1) = i) α i Ad i) Conjugate Gadient Methods, Summe
21 A Miacle Occus Pat 2 Othogonality of the esiduals: seach diections ae A-othogonal only one step in each diections hence: eo is A-othogonal to pevious seach diections: d i) ) T Ae j) = 0, fo i < j esiduals ae othogonal to pevious seach diections: d i) ) T j) = 0, fo i < j seach diections ae built fom esiduals: span { d 0),..., d i 1)} = span { 0),..., i 1)} hence: esiduals ae othogonal i) ) T j) = 0, i < j Conjugate Gadient Methods, Summe
22 A Miacle Occus Pat 3 combine othogonality and ecuence fo esiduals: i) ) T j+1) α j i) ) T Ad j) = i)) T j) α j i) ) T Ad j) = i)) T j) i)) T j+1) i)) T j) = 0, if i j: i) ) ) 1 T Ad j) α i i) T i), i = j = ) 1 α i 1 i) T i), i = j othewise. Conjugate Gadient Methods, Summe
23 A Miacle Occus Pat 4 computation of β ik fo k = 0,..., i 1): ) i) T Ad k) β ik = ) = d k) T Ad k) thus: seach diections ) i) T i) α ), if i = k + 1 i 1 d i 1) T Ad i 1) 0, if i > k + 1 i 1 d i) = i) + β ik d k) = i) + β i,i 1 d i 1) k=0 ) i) T i) β i := β i,i 1 = α ) i 1 d i 1) T Ad i 1) educes to a simple iteative scheme fo β i Conjugate Gadient Methods, Summe
24 A Miacle Occus Pat 5 build seach diections emembe: α i = d i) ) T i) d i) ) T Ad i) thus: α i d i) ) T Ad i) = d i) ) T i) d i+1) = i+1) + β i d i) β i+1 = ) i+1) T i+1) α ) i d i) T Ad i) β i+1 = i+1) ) T i+1) d i) ) T i) = i+1) ) T i+1) i) ) T i) last step: d i) ) T i) = i) + β i 1 d i 1)) T i) = i)) T i) + β i 1 d i 1) ) T i) = i)) T i) esidual i) othogonal to pevious seach diection d i 1) ) Conjugate Gadient Methods, Summe
25 Conjugate Gadients Algoithm Stat with d 0) = 0) = b Ax 0) While i) > ɛ iteate ove: 1. α i = i) ) T i) d i) ) T Ad i) 2. x i+1) = x i) + α i d i) 3. i+1) = i) α i Ad i) 4. β i+1 = i+1) ) T i+1) i) ) T i) 5. d i+1) = i+1) + β i+1 d i) Conjugate Gadient Methods, Summe
26 Pat III Peconditioning CG Convegence Peconditioning CG with Change-of-Basis Peconditioning CG with Matix Peconditione Peconditiones Examples ILU and Incomplete Cholesky Conjugate Gadient Methods, Summe
27 Conjugate Gadients Convegence Convegence Analysis: uses Kylow subspace: span { 0), A 0), A 2 0),..., A i 1 0)} Kylow subspace method Convegence Results: in pinciple: diect method n steps) howeve: othogonality lost due to ound-off eos exact solution not found) in pactice: iteative scheme ) i e i) A κ 1 e 0) 2, κ + 1 A κ = λmax/λ min Conjugate Gadient Methods, Summe
28 Peconditioning convegence depends on matix A idea: modify linea system Ax = b M 1 Ax = M 1 b, then: convegence depends on matix M 1 A optimal peconditione: M 1 = A 1 : in pactice: A 1 Ax = A 1 b x = A 1 b. avoid explicit computation of M 1 A find an M simila to A, compute effect of M 1 i.e., appoximate solution of SLE) o: find an M 1 simila to A 1 Conjugate Gadient Methods, Summe
29 CG and Peconditioning just eplace A by M 1 A in the algoithm?? poblem: M 1 A not necessaily symmetic even if M and A both ae) we will ty an altenative fist: symmetic peconditioning Ax = b L T ALˆx = L T b, x = Lˆx Remembe: fo Finite Element discetization, this coesponds to a change of basis functions! equies some e-computations in the CG algoithm see following slides) Conjugate Gadient Methods, Summe
30 Change-of-Basis Peconditioning peconditioned system of equations: Ax = b L } T {{ AL } )ˆx = }{{} L T b, x = Lˆx =:  =: ˆb computation of esidual: ˆ = ˆb  ˆx = LT b L T ALˆx = L T b Ax) = L T computation of α fo new system: i) i) i) i) i) i) ˆ )Tˆ ˆ )Tˆ ˆ )Tˆ α i := ˆd i) )T  ˆd = ˆd ) i) i) T LT AL ˆd = d ) i) i) T A d i) d i) whee we defined Lˆd i) =: update of solution: ˆx i+1) = ˆx i) i) + α i ˆd x i+1) = Lˆx i+1) = Lˆx i) + Lα i ˆd i) = x i) + α i d i) Conjugate Gadient Methods, Summe
31 Change-of-Basis Peconditioning 2) update esiduals ˆ: ˆ i+1) = ˆ i) α i ˆd i) = ˆ i) α i L T i) AL ˆd = ˆ i) α i L T A d i) computation of β i : β i+1 = i+1) i+1) ˆ )Tˆ ˆ i) )Tˆ i) update of seach diections: ˆd i+1) = ˆ i+1) i) + β i ˆd d i+1) = Lˆd i+1) = Lˆ i+1) i) + Lβ i ˆd = Lˆ i+1) i) + β i d Conjugate Gadient Methods, Summe
32 CG with Change-of-Basis Peconditioning Stat with ˆ 0) = L T b Ax 0) ) and d 0) = L ˆ 0) ; While ˆ i) > ɛ iteate ove: i) i) ˆ )Tˆ 1. α i = d ) i) T A d i) 2. x i+1) = x i) + α i d i) 3. ˆ i+1) = ˆ i) α i L T i) A d 4. β i+1 = i+1) i+1) ˆ )Tˆ ˆ i) )Tˆ i) 5. d i+1) = L ˆ i+1) + β i d i) Conjugate Gadient Methods, Summe
33 Hieachical Basis Peconditioning Some specifics fo the CG implementation: L tansfoms coefficient vecto fom hieachical basis to nodal basis, fo example ˆx = L x o d = L ˆd L T tansfoms the vecto of basis functions fom nodal basis to hieachical basis cmp. FEM), thus ˆ = L T effect of L and L T can be computed in ON) opeations HB-CG fo the Poisson poblem: in 1D: convegence afte a log N iteations! in this case: L T AL diagonal matix with log N diffeent eigenvalues) in 2D and 3D vey fast convegence! futhe impoved by additional diagonal peconditioning so-called hieachical geneating systems change to a multigid basis) achieve multigid-like pefomance Conjugate Gadient Methods, Summe
34 CG with Hieachical Geneating Systems Recall: system of linea equations A GS v GS = b GS given as A h A h P2h h A h P h 4h Rh 2hA h A 2h A 2h P4h 2h v h b h v 2h = Rh 2hb h Rh 4hA h R2h 4hA 2h A 4h v 4h Rh 4hb h Peconditioning fo CG? system A GS v GS = b GS is singula subspace of solutions v GS minima of the quadatic fom!) howeve: any of the solutions will do! convegence esult: i e i) A ˆκGS 1 e 0) 2, ˆκ GS = ˆκGS + 1) ˆλ max /ˆλ min A whee ˆκ GS is the atio of lagest vs. smallest non-zeo eigenvalue fo Poisson eq.: ˆκ GS independent of h multigid convegence Conjugate Gadient Methods, Summe
35 CG and Peconditioning evisited) peconditioning: eplace A by M 1 A poblem: M 1 A not necessaily symmetic compae symmetic peconditioning Ax = b L T ALˆx = Lb, x = Lˆx wokaound: find E T E = M Cholesky fact.), then Ax = b E T AE 1ˆx = E T b, ˆx = Ex what if E cannot be computed efficiently)? neithe M no M 1 might be known explicitly!) E, E T, E 1 can be eliminated fom algoithm again equies some e-computations): set ˆd = Ed and use ˆ = E T, ˆx = Ex, E 1 E T = M 1 Conjugate Gadient Methods, Summe
36 CG with Peconditione Stat: 0) = b Ax 0) ; d 0) = M 1 0) 1. α i = i) ) T M 1 i) d i) ) T Ad i) 2. x i+1) = x i) + α i d i) 3. i+1) = i) α i Ad i) 4. β i+1 = i+1) ) T M 1 i+1) i) ) T M 1 i) 5. d i+1) = M 1 i+1) + β i+1 d i) fo detailed deivation, see Shewchuck) Conjugate Gadient Methods, Summe
37 Implementation Peconditioning steps: M 1 i), M 1 i+1) M 1 known then multiply M 1 i) M known? Then solve My = i) to obtain y = M 1 i) neithe M, no M 1 ae known explicitly: algoithm to solve My = i) is sufficient! any appoximate solve fo Ae = i) algoithm to compute M 1 is sufficient! compute spase) appoximate invese SPAI) Examples: Multigid, Jacobi, ILU, SPAI,... Conjugate Gadient Methods, Summe
38 Peconditiones fo CG Examples find M A and compute effect of M 1 : Jacobi peconditioning: M := D A Symmetic) Gauss-Seidel peconditioning: M := L A o M = D A + L A )D 1 A D A + L A )T ), etc. just compute effect of M 1 : any appoximate solve might do incl. multigid methods incomplete LU-decomposition ILU) should be symmetic incomplete Cholesky factoization use a multigid method as peconditione?) wothwhile only) in situations whee multigid does not wok well) as stand-alone solve find an M 1 simila to A 1 spase appoximate invese SPAI) ties to minimise I MA F, whee M is a matix with given) spase non-zeo patten Conjugate Gadient Methods, Summe
39 Peconditiones ILU and Incomplete Cholesky Recall LU decomposition and Cholesky factoization: LU decomposition: given A, find lowe/uppe tiangula matices L and U such that A = LU Cholesky factoization: given A = A T, find lowe tiangula matix L such that A = LL T symmetic peconditioning! vaiants with explicit diagonal matix D: A = LD 1 U o A = LD 1 L T, whee L = D + L and R = D + R with stict lowe/uppe tiangula L, R but: fo spase A, L and U may be non-spase Idea: disegad all fill-in duing factoization Conjugate Gadient Methods, Summe
40 Cholesky Factoization L 11 L 21 L 22 L 31 L 32 L 33 D 1 11 D 1 22 D 1 33 L T 11 L T 21 L T 31 A 11 A T 21 A T 31 L T 22 L T 32 = A 21 A 22 A T 32 L T 33 A 31 A 32 A 33 Deive the factoization algoithm:! assume that A 11 = L 11 D 1 11 LT 11 is aleady factoized let L 21 be a 1 k submatix, i.e., to compute next ow of L: L 21 D 1 11 LT! 11 = A 21 D 1 11 LT 11 uppe tiangula matix solve tiangula system fo L 21 by convention L 22 = D 22, which is computed fom: L 21 D 1 11 LT 21+L 22 D 1 22 L T 22! = A 22 L 22 = D 22 := A 22 L 21 D 1 11 LT 21 Conjugate Gadient Methods, Summe
41 Incomplete Cholesky Factoization Algoithm: A LD 1 L T ) initialize D := 0, L := 0 fo i = 1,..., n: 1. fo k = 1,..., i 1: if i, k) S then set L ik := A ik L ij D 1 2. set L ii = D ii := A ii L ij D 1 jj L ij j<i note: sums j<k and only conside non-zeo elements S j<i uses given patten S of non-zeo elements in the factoization fequent choice: use non-zeos of A fo S) Cholesky factoization computed in On) opeations fo spase matices with c n non-zeos) fequently used fo peconditioning j<k jj L kj Conjugate Gadient Methods, Summe
42 Liteatue/Refeences Conjugate Gadients: Shewchuk: An Intoduction to the Conjugate Gadient Method Without the Agonizing Pain. Hackbusch: Iteative Solution of Lage Spase Systems of Equations, Spinge M. Giebel: Multilevelmethoden als Iteationsvefahen ï 1 2 be Ezeugendensystemen, Teubne Skipten zu Numeik, 1994 M. Giebel: Multilevel algoithms consideed as iteative methods on semidefinite systems, SIAM Int. J. Sci. Stat. Comput. 153), Conjugate Gadient Methods, Summe
Conjugate Gradient Methods. Michael Bader. Summer term 2012
Gadient Methods Outlines Pat I: Quadatic Foms and Steepest Descent Pat II: Gadients Pat III: Summe tem 2012 Pat I: Quadatic Foms and Steepest Descent Outlines Pat I: Quadatic Foms and Steepest Descent
More informationConjugate Gradient (CG) Method
Optimization II Conugate Gaient CG Metho Anothe metho oes not equie explicit secon eivatives, an oes not even stoe appoximation to Hessian matix CG geneates sequence of conugate seach iections, implicitly
More informationMAC Module 12 Eigenvalues and Eigenvectors
MAC 23 Module 2 Eigenvalues and Eigenvectos Leaning Objectives Upon completing this module, you should be able to:. Solve the eigenvalue poblem by finding the eigenvalues and the coesponding eigenvectos
More informationMULTILAYER PERCEPTRONS
Last updated: Nov 26, 2012 MULTILAYER PERCEPTRONS Outline 2 Combining Linea Classifies Leaning Paametes Outline 3 Combining Linea Classifies Leaning Paametes Implementing Logical Relations 4 AND and OR
More informationEVALUATION OF PRECONDITIONED KRYLOV SUBSPACE METHODS FOR NAVIER-STOKES EQUATIONS
The Seventh Asia-Pacific Confeence on Wind Engineeing, Novembe 8-12, 2009, Taipei, Taiwan EVALUATION OF PRECONDITIONED KRYLOV SUBSPACE METHODS FOR NAVIER-STOKES EQUATIONS Aishe Zhang 1,2, Ming Gu 3, Cuilan
More information15.081J/6.251J Introduction to Mathematical Programming. Lecture 6: The Simplex Method II
15081J/6251J Intoduction to Mathematical Pogamming ectue 6: The Simplex Method II 1 Outline Revised Simplex method Slide 1 The full tableau implementation Anticycling 2 Revised Simplex Initial data: A,
More informationA Bijective Approach to the Permutational Power of a Priority Queue
A Bijective Appoach to the Pemutational Powe of a Pioity Queue Ia M. Gessel Kuang-Yeh Wang Depatment of Mathematics Bandeis Univesity Waltham, MA 02254-9110 Abstact A pioity queue tansfoms an input pemutation
More informationInformation Retrieval Advanced IR models. Luca Bondi
Advanced IR models Luca Bondi Advanced IR models 2 (LSI) Pobabilistic Latent Semantic Analysis (plsa) Vecto Space Model 3 Stating point: Vecto Space Model Documents and queies epesented as vectos in the
More informationMATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form
MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More informationMultifrontal sparse QR factorization on the GPU
Multifontal spase QR factoization on the GPU Tim Davis, Sanjay Ranka, Shaanyan Chetlu, Nui Yealan Univesity of Floida Feb 2012 GPU-based Multifontal QR factoization why spase QR? multifontal spase QR in
More informationTemporal-Difference Learning
.997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation
More information6 Matrix Concentration Bounds
6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationLight Time Delay and Apparent Position
Light Time Delay and ppaent Position nalytical Gaphics, Inc. www.agi.com info@agi.com 610.981.8000 800.220.4785 Contents Intoduction... 3 Computing Light Time Delay... 3 Tansmission fom to... 4 Reception
More information3. Electromagnetic Waves II
Lectue 3 - Electomagnetic Waves II 9 3. Electomagnetic Waves II Last time, we discussed the following. 1. The popagation of an EM wave though a macoscopic media: We discussed how the wave inteacts with
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationAppendix A. Appendices. A.1 ɛ ijk and cross products. Vector Operations: δ ij and ɛ ijk
Appendix A Appendices A1 ɛ and coss poducts A11 Vecto Opeations: δ ij and ɛ These ae some notes on the use of the antisymmetic symbol ɛ fo expessing coss poducts This is an extemely poweful tool fo manipulating
More informationNumerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 41 Digital Signal Pocessing Pof. Mak Fowle Note Set #31 Linea Phase FIR Design Optimum Equiipple (Paks-McClellan) Reading: Sect. 1.2.4 1.2.6 of Poakis & Manolakis 1/2 Motivation The window method and
More informationThe Strain Compatibility Equations in Polar Coordinates RAWB, Last Update 27/12/07
The Stain Compatibility Equations in Pola Coodinates RAWB Last Update 7//7 In D thee is just one compatibility equation. In D polas it is (Equ.) whee denotes the enineein shea (twice the tensoial shea)
More information18.06 Problem Set 4 Solution
8.6 Poblem Set 4 Solution Total: points Section 3.5. Poblem 2: (Recommended) Find the lagest possible numbe of independent vectos among ) ) ) v = v 4 = v 5 = v 6 = v 2 = v 3 =. Solution (4 points): Since
More informationCircuit Synthesizable Guaranteed Passive Modeling for Multiport Structures
Cicuit Synthesizable Guaanteed Passive Modeling fo Multipot Stuctues Zohaib Mahmood, Luca Daniel Massachusetts Institute of Technology BMAS Septembe-23, 2010 Outline Motivation fo Compact Dynamical Passive
More informationStanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012
Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,
More informationLocalization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matrix
Jounal of Sciences, Islamic Republic of Ian (): - () Univesity of Tehan, ISSN - http://sciencesutaci Localization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matix H Ahsani
More informationComplex Eigenvalues. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB
Pepaed by Vince Zaccone Fo ampus Leaning ssistance Sevices at USB omplex Numbes When solving fo the oots of a quadatic equation, eal solutions can not be found when the disciminant is negative. In these
More informationMinimization of the quadratic test function
Miimizatio of the quadatic test fuctio A quadatic fom is a scala quadatic fuctio of a vecto with the fom f ( ) A b c with b R A R whee A is assumed to be SPD ad c is a scala costat Note: A symmetic mati
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationThis is a very simple sampling mode, and this article propose an algorithm about how to recover x from y in this condition.
3d Intenational Confeence on Multimedia echnology(icm 03) A Simple Compessive Sampling Mode and the Recovey of Natue Images Based on Pixel Value Substitution Wenping Shao, Lin Ni Abstact: Compessive Sampling
More informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More information2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum
2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known
More informationEquilibrium Reconstruction with Grad-Shafranov Equation
Equilibium Reconstuction with Gad-Shafanov Equation Qian Peng (qp215) Apil 14, 211 1 Gad-Shafanov equation [1]Fo axial symmetic case, the tooidal B eld is B to = F φ. Whee 2π F (, θ) is the total poloidal
More informationFractional Tikhonov regularization for linear discrete ill-posed problems
BIT manuscipt No. (will be inseted by the edito) Factional Tikhonov egulaization fo linea discete ill-posed poblems Michiel E. Hochstenbach Lotha Reichel Received: date / Accepted: date Abstact Tikhonov
More informationA Crash Course in (2 2) Matrices
A Cash Couse in ( ) Matices Seveal weeks woth of matix algeba in an hou (Relax, we will only stuy the simplest case, that of matices) Review topics: What is a matix (pl matices)? A matix is a ectangula
More informationChapter 3 Optical Systems with Annular Pupils
Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The
More informationNumerical Solution of Boundary Value Problems for the Laplacian in R 3 in the Case of Complex Boundary Surface
Computational Applied Mathematics Jounal 5; (: 9-5 Published online Febuay 5 (http://www.aascit.og/jounal/camj Numeical Solution of Bounday Value Poblems fo the Laplacian in R in the Case of Complex Bounday
More information1 Spherical multipole moments
Jackson notes 9 Spheical multipole moments Suppose we have a chage distibution ρ (x) wheeallofthechageiscontained within a spheical egion of adius R, as shown in the diagam. Then thee is no chage in the
More informationAnalytical time-optimal trajectories for an omni-directional vehicle
Analytical time-optimal tajectoies fo an omni-diectional vehicle Weifu Wang and Devin J. Balkcom Abstact We pesent the fist analytical solution method fo finding a time-optimal tajectoy between any given
More informationPreconditioning of conjugate gradient iterations in interior point MPC method
MISUBISHI ELECRIC RESEARCH LABORAORIES http://wwwmelcom Peconditioning of conjugate gadient iteations in inteio point MPC method Malyshev, A; Quiynen, R; Knyazev, A R2018-120 August 25, 2018 Abstact hee
More informationTopic 4a Introduction to Root Finding & Bracketing Methods
/8/18 Couse Instucto D. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: cumpf@utep.edu Topic 4a Intoduction to Root Finding & Backeting Methods EE 4386/531 Computational Methods in EE Outline
More informationQuantum Fourier Transform
Chapte 5 Quantum Fouie Tansfom Many poblems in physics and mathematics ae solved by tansfoming a poblem into some othe poblem with a known solution. Some notable examples ae Laplace tansfom, Legende tansfom,
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationChapter 3: Theory of Modular Arithmetic 38
Chapte 3: Theoy of Modula Aithmetic 38 Section D Chinese Remainde Theoem By the end of this section you will be able to pove the Chinese Remainde Theoem apply this theoem to solve simultaneous linea conguences
More informationResearch Design - - Topic 17 Multiple Regression & Multiple Correlation: Two Predictors 2009 R.C. Gardner, Ph.D.
Reseach Design - - Topic 7 Multiple Regession & Multiple Coelation: Two Pedictos 009 R.C. Gadne, Ph.D. Geneal Rationale and Basic Aithmetic fo two pedictos Patial and semipatial coelation Regession coefficients
More informationQuestion 1: The dipole
Septembe, 08 Conell Univesity, Depatment of Physics PHYS 337, Advance E&M, HW #, due: 9/5/08, :5 AM Question : The dipole Conside a system as discussed in class and shown in Fig.. in Heald & Maion.. Wite
More informationMA557/MA578/CS557. Lecture 14. Prof. Tim Warburton. Spring
MA557/MA578/CS557 Lectue 4 Sping 3 Pof. Tim Wabuton timwa@math.unm.edu Matlab Notes To ceate a symbolic vaiable (say theta) use the command: theta = sym( theta ); Example manipulation: Matlab cont To evaluate
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationChapter 5 Linear Equations: Basic Theory and Practice
Chapte 5 inea Equations: Basic Theoy and actice In this chapte and the next, we ae inteested in the linea algebaic equation AX = b, (5-1) whee A is an m n matix, X is an n 1 vecto to be solved fo, and
More informationHomework # 3 Solution Key
PHYSICS 631: Geneal Relativity Homewok # 3 Solution Key 1. You e on you hono not to do this one by hand. I ealize you can use a compute o simply look it up. Please don t. In a flat space, the metic in
More informationMoment-free numerical approximation of highly oscillatory integrals with stationary points
Moment-fee numeical appoximation of highly oscillatoy integals with stationay points Sheehan Olve Abstact We pesent a method fo the numeical quadatue of highly oscillatoy integals with stationay points.
More informationLecture 8 - Gauss s Law
Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.
More informationValue Prediction with FA. Chapter 8: Generalization and Function Approximation. Adapt Supervised Learning Algorithms. Backups as Training Examples [ ]
Chapte 8: Genealization and Function Appoximation Objectives of this chapte:! Look at how expeience with a limited pat of the state set be used to poduce good behavio ove a much lage pat.! Oveview of function
More informationEQUATIONS OF MOTION LUCA GUIDO MOLINARI
EQUATIONS OF MOTION LUCA GUIDO MOLINARI 1. Equation of motion of destuction opeatos Conside a system of bosons o femions descibed by a Hamiltonian H = H 1 + H 2, whee H 1 and H 2 ae espectively the one
More informationPhysics 2B Chapter 22 Notes - Magnetic Field Spring 2018
Physics B Chapte Notes - Magnetic Field Sping 018 Magnetic Field fom a Long Staight Cuent-Caying Wie In Chapte 11 we looked at Isaac Newton s Law of Gavitation, which established that a gavitational field
More informationMODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ...
MODULE 5a and 5b (Stewat, Sections 2.2, 2.3) INTRO: In MATH 4 vectos wee witten eithe as ows (a, a2,..., an) o as columns a a 2... a n and the set of all such vectos of fixed length n was called the vecto
More informationPISCES II : 2.5-D RF Cavity Code
CAP'96 (Computational Acceleato Physics) Williamsbug, Viginia PISCES II :.5-D RF Cavity Code Yoshihisa Iwashita Acceleato Laboatoy, Nuclea Science Reseach Facility Institute fo Chemical Reseach, Kyoto
More informationA Power Method for Computing Square Roots of Complex Matrices
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 13, 39345 1997 ARTICLE NO. AY975517 A Powe Method fo Computing Squae Roots of Complex Matices Mohammed A. Hasan Depatment of Electical Engineeing, Coloado
More informationChapter 8: Generalization and Function Approximation
Chapte 8: Genealization and Function Appoximation Objectives of this chapte: Look at how expeience with a limited pat of the state set be used to poduce good behavio ove a much lage pat. Oveview of function
More informationSMT 2013 Team Test Solutions February 2, 2013
1 Let f 1 (n) be the numbe of divisos that n has, and define f k (n) = f 1 (f k 1 (n)) Compute the smallest intege k such that f k (013 013 ) = Answe: 4 Solution: We know that 013 013 = 3 013 11 013 61
More informationDetermining solar characteristics using planetary data
Detemining sola chaacteistics using planetay data Intoduction The Sun is a G-type main sequence sta at the cente of the Sola System aound which the planets, including ou Eath, obit. In this investigation
More informationReading Assignment. Problem Description for Homework #9. Read Chapters 29 and 30.
Reading Assignment Read Chaptes 29 and 30. Poblem Desciption fo Homewok #9 In this homewok, you will solve the inhomogeneous Laplace s equation to calculate the electic scala potential that exists between
More informationFE FORMULATIONS FOR PLASTICITY
G These slides ae designed based on the book: Finite Elements in Plasticity Theoy and Pactice, D.R.J. Owen and E. Hinton, 970, Pineidge Pess Ltd., Swansea, UK. Couse Content: A INTRODUCTION AND OVERVIEW
More informationPart V: Closed-form solutions to Loop Closure Equations
Pat V: Closed-fom solutions to Loop Closue Equations This section will eview the closed-fom solutions techniques fo loop closue equations. The following thee cases will be consideed. ) Two unknown angles
More informationDOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS
DOING PHYIC WITH MTLB COMPUTTIONL OPTIC FOUNDTION OF CLR DIFFRCTION THEORY Ian Coope chool of Physics, Univesity of ydney ian.coope@sydney.edu.au DOWNLOD DIRECTORY FOR MTLB CRIPT View document: Numeical
More informationQualifying Examination Electricity and Magnetism Solutions January 12, 2006
1 Qualifying Examination Electicity and Magnetism Solutions Januay 12, 2006 PROBLEM EA. a. Fist, we conside a unit length of cylinde to find the elationship between the total chage pe unit length λ and
More informationHua Xu 3 and Hiroaki Mukaidani 33. The University of Tsukuba, Otsuka. Hiroshima City University, 3-4-1, Ozuka-Higashi
he inea Quadatic Dynamic Game fo Discete-ime Descipto Systems Hua Xu 3 and Hioai Muaidani 33 3 Gaduate School of Systems Management he Univesity of suuba, 3-9- Otsua Bunyo-u, oyo -0, Japan xuhua@gssm.otsua.tsuuba.ac.jp
More informationComputational Methods of Solid Mechanics. Project report
Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate
More informationSecret Exponent Attacks on RSA-type Schemes with Moduli N = p r q
Secet Exponent Attacks on RSA-type Schemes with Moduli N = p q Alexande May Faculty of Compute Science, Electical Engineeing and Mathematics Univesity of Padebon 33102 Padebon, Gemany alexx@uni-padebon.de
More informationGeneral Solution of EM Wave Propagation in Anisotropic Media
Jounal of the Koean Physical Society, Vol. 57, No. 1, July 2010, pp. 55 60 Geneal Solution of EM Wave Popagation in Anisotopic Media Jinyoung Lee Electical and Electonic Engineeing Depatment, Koea Advanced
More informationMagnetic Dipoles Challenge Problem Solutions
Magnetic Dipoles Challenge Poblem Solutions Poblem 1: Cicle the coect answe. Conside a tiangula loop of wie with sides a and b. The loop caies a cuent I in the diection shown, and is placed in a unifom
More informationA STUDY OF HAMMING CODES AS ERROR CORRECTING CODES
AGU Intenational Jounal of Science and Technology A STUDY OF HAMMING CODES AS ERROR CORRECTING CODES Ritu Ahuja Depatment of Mathematics Khalsa College fo Women, Civil Lines, Ludhiana-141001, Punjab, (India)
More informationGeometry of the homogeneous and isotropic spaces
Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant
More informationGradient-based Neural Network for Online Solution of Lyapunov Matrix Equation with Li Activation Function
Intenational Confeence on Infomation echnology and Management Innovation (ICIMI 05) Gadient-based Neual Netwok fo Online Solution of Lyapunov Matix Equation with Li Activation unction Shiheng Wang, Shidong
More informationFractional Zero Forcing via Three-color Forcing Games
Factional Zeo Focing via Thee-colo Focing Games Leslie Hogben Kevin F. Palmowski David E. Robeson Michael Young May 13, 2015 Abstact An -fold analogue of the positive semidefinite zeo focing pocess that
More informationToday in Physics 218: radiation from moving charges
Today in Physics 218: adiation fom moving chages Poblems with moving chages Motion, snapshots and lengths The Liénad-Wiechet potentials Fields fom moving chages Radio galaxy Cygnus A, obseved by Rick Peley
More informationFresnel Diffraction. monchromatic light source
Fesnel Diffaction Equipment Helium-Neon lase (632.8 nm) on 2 axis tanslation stage, Concave lens (focal length 3.80 cm) mounted on slide holde, iis mounted on slide holde, m optical bench, micoscope slide
More informationPHYS 705: Classical Mechanics. Small Oscillations
PHYS 705: Classical Mechanics Small Oscillations Fomulation of the Poblem Assumptions: V q - A consevative system with depending on position only - The tansfomation equation defining does not dep on time
More informationCOMP Parallel Computing SMM (3) OpenMP Case Study: The Barnes-Hut N-body Algorithm
COMP 633 - Paallel Computing Lectue 8 Septembe 14, 2017 SMM (3) OpenMP Case Study: The Banes-Hut N-body Algoithm Topics Case study: the Banes-Hut algoithm Study an impotant algoithm in scientific computing»
More informationRelating Branching Program Size and. Formula Size over the Full Binary Basis. FB Informatik, LS II, Univ. Dortmund, Dortmund, Germany
Relating Banching Pogam Size and omula Size ove the ull Binay Basis Matin Saueho y Ingo Wegene y Ralph Wechne z y B Infomatik, LS II, Univ. Dotmund, 44 Dotmund, Gemany z ankfut, Gemany sauehof/wegene@ls.cs.uni-dotmund.de
More informationEE 5337 Computational Electromagnetics
Instucto D. Ramond Rumpf (95) 747 6958 cumpf@utep.edu EE 5337 Computational Electomagnetics Lectue # WEM Etas Lectue These notes ma contain copighted mateial obtained unde fai use ules. Distibution of
More informationGENLOG Multinomial Loglinear and Logit Models
GENLOG Multinomial Loglinea and Logit Models Notation This chapte descibes the algoithms used to calculate maximum-likelihood estimates fo the multinomial loglinea model and the multinomial logit model.
More informationMethod for Approximating Irrational Numbers
Method fo Appoximating Iational Numbes Eic Reichwein Depatment of Physics Univesity of Califonia, Santa Cuz June 6, 0 Abstact I will put foth an algoithm fo poducing inceasingly accuate ational appoximations
More informationClassical Mechanics Homework set 7, due Nov 8th: Solutions
Classical Mechanics Homewok set 7, due Nov 8th: Solutions 1. Do deivation 8.. It has been asked what effect does a total deivative as a function of q i, t have on the Hamiltonian. Thus, lets us begin with
More informationLINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE
LINEAR AND NONLINEAR ANALYSES O A WIND-TUNNEL INTRODUCTION BALANCE R. Kakehabadi and R. D. Rhew NASA LaRC, Hampton, VA The NASA Langley Reseach Cente (LaRC) has been designing stain-gauge balances fo utilization
More informationEXAM NMR (8N090) November , am
EXA NR (8N9) Novembe 5 9, 9. 1. am Remaks: 1. The exam consists of 8 questions, each with 3 pats.. Each question yields the same amount of points. 3. You ae allowed to use the fomula sheet which has been
More informationGoodness-of-fit for composite hypotheses.
Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test
More informationA Relativistic Electron in a Coulomb Potential
A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify
More informationMultipole Radiation. February 29, The electromagnetic field of an isolated, oscillating source
Multipole Radiation Febuay 29, 26 The electomagnetic field of an isolated, oscillating souce Conside a localized, oscillating souce, located in othewise empty space. We know that the solution fo the vecto
More informationNew Restricted Isometry Results for Noisy Low-rank Recovery
New Resticted Isomety Results fo Noisy Low-ank Recovey Kathik Mohan and Mayam Fazel Abstact he poblem of ecoveing a low-ank matix consistent with noisy linea measuements is a fundamental poblem with applications
More informationInternet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks
Intenet Appendix fo A Bayesian Appoach to Real Options: The Case of Distinguishing Between Tempoay and Pemanent Shocks Steven R. Genadie Gaduate School of Business, Stanfod Univesity Andey Malenko Gaduate
More informationIntroduction to Orbital-Free Density-Functional Theory. Ralf Gehrke FHI Berlin, February 8th 2005
Intoduction to Obital-Fee Density-Functional heoy Ralf Gehke FHI Belin, Febuay 8th 005 Outline Basics of functional deivatives I Pinciples of Obital-fee Density-Functional heoy basics of Density-Functional
More informationA FURTHER SUBSPACE METHOD OPTIMIZATION FOR TRACKING REAL VALUED SINUSOIDS IN NOISE
Électonique et tansmission de l infomation A FURHER SUBSPACE MEHOD OPIMIZAION FOR RACKING REAL VALUED SINUSOIDS IN NOISE ŞEFAN SLAVNICU, SILVIU CIOCHINĂ Key wods: Subspace tacking, Fequency estimation,
More informationCERFACS 42 av. Gaspard Coriolis, Toulouse, Cedex 1, France. Available at Date: April 2, 2008.
ON THE BLOCK TRIANGULAR FORM OF SYMMETRIC MATRICES IAIN S. DUFF and BORA UÇAR Technical Repot: No: TR/PA/08/26 CERFACS 42 av. Gaspad Coiolis, 31057 Toulouse, Cedex 1, Fance. Available at http://www.cefacs.f/algo/epots/
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationPhysics 235 Chapter 5. Chapter 5 Gravitation
Chapte 5 Gavitation In this Chapte we will eview the popeties of the gavitational foce. The gavitational foce has been discussed in geat detail in you intoductoy physics couses, and we will pimaily focus
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More information15 Solving the Laplace equation by Fourier method
5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More information