Parallel KS Block-Step Method. Sverre Aarseth. Institute of Astronomy, Cambridge

Size: px
Start display at page:

Download "Parallel KS Block-Step Method. Sverre Aarseth. Institute of Astronomy, Cambridge"

Transcription

1 Parallel KS Block-Step Method Sverre Aarseth Institute of Astronomy, Cambridge Code Overview Hermite KS Prediction & Correction Iteration Time-Steps Decision-Making Binary Project

2 Code Overview Directories Ncode, Docs, Chain, Nchain, ARchain, ARint, Block Standard make nbody6 in dir Ncode nbody6 Block-steps make nbody6 in dir Block nbody6b NBODY7 make -f Makefile7 nbody6 in Ncode nbody7 NBODY7b make -f Makefile7 nbody6 in Block nbody7b GPU make gpu in dir GPU2 produces nbody6.gpu GPUb make -f Makefile6 gpu nbody6b.gpu ARC make -f Makefile7 gpu nbody7.gpu ARCb make -f Makefile8 gpu nbody7b.gpu

3 Program Flow KS predict Predict u, u to order u (5) Perturbers Predict r, ṙ and obtain F, Ḟ Corrector Temporary correct u,u & h Iteration One or more, same F, Ḟ Final stage Form u (4), u (5), copy to common Time Derivatives up to t (5) or t (6) Time-step Set τ, convert to t & truncate Invert Newton Raphson gives τ bs

4 Hermite KS Standard KS New notation u h t = 2 hu + 2 R LT F kl = 2u L T F kl = u u P = F kl F u = u Q = L T P, Basic equations F u = 2 hu + 2 RQ h t = 2u Q = u u Hermite F, F formulation F u = 2 hu + 2 RQ F u = 2 (h u + hu + R Q + RQ ) = 2u Q h = 2F u Q + 2u Q t = u u h

5 Hermite KS with block-steps Predict u, u, h to highest order u pred = u + τu + τ2 2 u + τ3 6 u + τ4 24 u 4 + τ5 20 u u pred = u + τu + τ2 2 u + τ3 6 u 4 + τ4 24 u 5, h pred = h + τh + τ 2 h + τ 6 h + τ 24 h. Force evaluation P = L T uf, P = L T u F + rl T u F. (2) Derivatives for corrector u = (hu + rp), 2 u = 2 (h u + hu + r P + rp ), h = 2(u P), h = 2(u P + u P ). (3)

6 Fourth-order corrector (Hut, Makino & McMillan 995) u corr = u old + τ 2 (u corr + u old ) τ2 2 (u new u old ), u corr = u old + τ 2 (u new + u old ) τ2 2 (u new u old ), h corr = h old + τ 2 (h new + h old ) τ2 2 (h new h old ). (4) Velocity derivatives for prediction u 4 mid = u new u τ u 5 mid = 2 τ 2 old, ( u new + u old 2 u new u ) old. (5) τ Shifted from midpoint u 4 new = u 4 mid + τ 2 u 5 mid, u 5 new = u 5 mid, h new = h mid + τ h new = h 2 h mid, mid. (6)

7 Taylor series for physical time t = u u, t = 2(u u ), t = 2(u u + u u ), t 4 = 2(u u ) + 6(u u ), t 5 = 2(u u 4 + 8(u u ) + 6(u u ), t 6 = 2(u u 5 + 0(u u (u u ). (7) Next look-up time t next = t + 6 n= t n τ n. (8) n! Regularized time-step by inversion t + 6 n= t n τ n n! = 0. (9) Newton-Raphson iteration to solve for τ x i+ = x i f(x i) f (x i ).

8 Initial guess τ = t bs t τ orig, with quantized time-step t bs (= 2 n t orig ). Commensurability condition mod (t sys, t bs ) = 0. Solve equation (9) to obtain τ bs.

9 0.000 ECC = 0.96 e-05 DT e-06 e e Time

10 0.03 ECC = DTU e Time

11 subint.f!!!!!! # %&# ()! # () +(, (./0!!!!!!! :;58<!=7>!6? 6Α6Β396Χ!!!!!!!!!!!!!() Ε#2!Φ479978ΓΧΗΦ!!!!!! &ΙΙ&)ϑΚ ΛΜ%ϑ!!( Μ +ΚΙΜΝ0!!!!!!! 33!ΟΗ3ΒΗ3>!Β7!:ΠΘ:843!Κ!67Ρ?Β5786!:Β!83Ο! Ρ74; 6Β3ΣΧ!!!!!!(Τ!+)ΛΜ(% ΧΥ Χ/0! ς )!!!!!!!Ι:;3!:!Ρ56Β!7=!:ΡΡ!939 3>6!Π?3!?Σ!Β7!38Π!7=! Ρ74; 6Β3ΣΧ!!!!!!!!!!(.!Ω!.!!!!!!!!!!)Κ!Ω!/!!!!!!!!!!2&!Ξ!(Λ!Ω!. )ΛΜ(%!!!!!!!!!!!!(.!Ω!(.!Ψ!Ξ!!!!!!!!!!!!(Τ!+ /+(.0!Ψ! Λ+(.0ΧΕ Χ Ε& Κ0! ς )!!!!!!!!!!!!!!)Κ!Ω!)Κ!Ψ!.!!!!!!!!!!!!!!Ε( Κ +)Κ0!Ω!(.!!!!!!!!!!!! )2!(Τ!!!!Ξ!!!!! &) ()#!!!!!!!23Β3>9583!Ζ?:8Β5[3Π!Θ:Ρ?3!7=! Ι()Χ!!!!!!!!! Ι()!Ω! Ε& Κ!!!!!!!!!!2&!Γ!Ε!Ω!. )Κ!!!!!!!!!!!!(.!Ω!Ε( Κ +Ε0!!!!!!!!!!!! Ι()!Ω!Ι()+ /+(.0!Ψ! Λ+(.0 Ι()0!!!!Γ!!!!! &) ()#!!!!!!!#ΣΠ:Β3!Β593!:8Π!=7>9!Ρ56Β!7=! 58:>536!Π?3!:Β! Ι()Χ!!!!!!!!!! (Ι!Ω! Ι()!!!!!!!!!!Ε )Υ ς!ω!/!!!!!!!!!!2&!]!ε!ω!. )Κ!!!!!!!!!!!!(.!Ω!Ε( Κ +Ε0!!!!!!!!!!!!(Τ!+ /+(.0!Ψ! Λ+(.0Χ,Χ Ι()0! ς )!!!!!!!!!!!!!!Ε )Υ ς!ω!ε )Υ ς!ψ!.!!!!!!!!!!!!!!) Ν Ε+Ε )Υ ς0!ω!(.!!!!!!!!!!!! )2!(Τ!!!!]!!!!! &) ()#!!!!!!! Η34;!Σ7665 Ρ3!4Η:58!58Β3<>:Β578!78![3>7!Ε )Υ ςχ!!!!!!!!!!(τ!+ε )Υ ςχ,χ/0! ς )!!!!!!!!!!!!Υ&! &!Γ/!!!!!!!!!!! )2!(Τ!!!!!!! Ρ3:8!4?>>38Β!Ρ38<ΒΗ!7=!( Μ Χ!!!!!!!!!!2&!!Ε!Ω!. Ε )Υ ς!!!!!!!!!!!!( Μ +Ε0!Ω!/!!!!!!!!! &) ()#!!!!!!! 3<58!Κ!58Β3<>:Β578!+63Ζ?38Β5:Ρ!7>!Σ:>:ΡΡ3Ρ0Χ!!!!!!!!!!(Τ!+Ε )Υ ςχε Χ)ΕΙΜΝ0! ς )!!!!!!!!!!!!2&!./!Ε(!Ω!. Ε )Υ ς!!!!!!!!!!!!!!(.!ω!) Ν Ε+Ε(0!!!!!!!!!!!!!! ΜΕΕ!Κ () +(. Ε(0!!!./!!!!!!! &) ()#!!!!!!! 3:>4Η!878 [3>7!=Ρ:<6!+( Μ!_!/!93:86!47ΡΡ565780Χ!!!!!!!!!!!!(Τ!+(,Χ,Χ/0! ς )!!!!!!!!!!!!!!2&!.Ξ!Ε(!Ω!. Ε )Υ ς!!!!!!!!!!!!!!!!(τ!+( Μ +Ε(0Χ) Χ/0! ς )!!!!!!!!!!!!!!!!!!(,!Ω!( Μ +Ε(0!!!!!!!!!!!!!!!!!!(./!Ω!) Ν Ε+Ε(0!!!!!!!!!!!!!!!!!!Κ ΛΜ(%!Ω!Κ +(./0!!!!!!!)7Β3!(ΛςΜ!9?6Β! 3!Π3=583Π!=7>!4Η:58Χ!!!!!!!!!!!!!!!!!!(ΛςΜ!Ω!(,!!!!!!!!!!!!!!!!!!Υ&! &!.!!!!!!!!!!!!!!!! )2!(Τ!!!!.Ξ!!!!!!!!! &) ()#!!!!!!!!!!!! )2!(Τ!!!!!!! >3:Β!47ΡΡ56578!3αΣΡ545ΒΡΑ! 3=7>3!Ζ?5ΒΒ58<Χ!!!.!!!!!!!(Τ!+(,ΧΕ Χ/0! ς )

12 subint.f 2!!!!!!!!!!!!!!(Λ!Ω!Κ +(./0!!!!!!!!!!!!!!Κ ΛΜ(%!Ω!(Λ!!!!!!!!!!!!!!,Λ %(!Ω!%+(Λ0!!!!!!!!!!!!!!(, &ΕΕ!Ω! Ξ!!!!!!!!!!!!!! ΜΕΕ! Ι &2β+,Λ %( Ξ0!!!!!!! 86?>3! Ρ74; 6Β3Σ!Β593!=7>!87>9:Ρ!478Β58?:Β578Χ!!!!!!!!!!!!!! (Ι!Ω! Ε& Κ!!!!!!!)7Β3!>36Β!7=!ΒΗ3! Ρ74;!Ο5ΡΡ! 3!Π783!83αΒ!Β593Χ!!!!!!!!!!!!!!Υ&! &!Γ/!!!!!!!!!!!! )2!(Τ!!!!!!! 79ΣΡ3Β3!:ΡΡ! Ρ74; 6Β3Σ6! 3=7>3!6Σ345:Ρ!Σ>743Π?>3Χ!!!!!!!!!!!!(Τ!+ Ι()ΧΕ Χ Ε& Κ0!Υ&! &!!!!!!!!!!!!! (Ι!Ω! Ε& Κ!!!!!!!!!! Ε!!!!!!!Λ3>=7>9!ΒΗ3!Σ:>:ΡΡ3Ρ!3α34?Β578!Ρ77ΣΧ χδ79σ!σ:>:ρρ3ρ!π7!σ>5θ:β3+ε( (.0!!!!!!!!!!!!2&!Ξ/!Ε(!Ω!. Ε )Υ ς!!!!!!!!!!!!!!(.!ω!) Ν Ε+Ε(0!!!!!!!!!!!!!! ΜΕΕ!Κ () Λ+(. Ε(0!!!Ξ/!!!!!!! &) ()# χδ79σ!38π!σ:>:ρρ3ρ!π7!!!!!!!!!!!!) Λ#!Ω!) Λ#!Ψ!Ε )Υ ς!!!!!!! 3:>4Η!878 [3>7!=Ρ:<6!+( Μ!_!/!=7>!47ΡΡ565780Χ!!!!!!!!!!!!(Τ!+(,Χ,Χ/0! ς )!!!!!!!!!!!!!!2&!Ξ!Ε(!Ω!. Ε )Υ ς!!!!!!!!!!!!!!!!(τ!+( Μ +Ε(0Χ) Χ/0! ς )!!!!!!!!!!!!!!!!!!(,!Ω!( Μ +Ε(0!!!!!!!!!!!!!!!!!!(./!Ω!) Ν Ε+Ε(0!!!!!!!!!!!!!!!!!!(ΛςΜ!Ω!(,!!!!!!!!!!!!!!!!!!Κ ΛΜ(%!Ω!Κ +(./0!!χ!9:Α! 3!833Π3ΠΧ!!!!!!!!!!!!!!!!!!Υ&! &!Ξ]!!!!!!!!!!!!!!!! )2!(Τ!!!!Ξ!!!!!!!!! &) ()#!!!!!!!!!!!! )2!(Τ!!!!!!! >3:Β!47ΡΡ56578!3αΣΡ545ΒΡΑ! 3=7>3!Ζ?5ΒΒ58<Χ!!!Ξ]!!!!!!!(Τ!+(,ΧΕ Χ/0! ς )!!!!!!!!!!!!!!(Λ!Ω!Κ +(./0!!!!!!!!!!!!!!Κ ΛΜ(%!Ω!(Λ!!!!!!!!!!!!!!,Λ %(!Ω!%+(Λ0!!!!!!!!!!!!!!(, &ΕΕ!Ω! Ξ!!!!!!!!!!!!!! ΜΕΕ! Ι &2β+,Λ %( Ξ0!!!!!!!)7Β3!>36Β!7=!ΒΗ3! Ρ74;!Ο5ΡΡ! 3!Π783!83αΒ!Β593Χ!!!!!!!!!!!!!!Υ&! &!Γ/!!!!!!!!!!!! )2!(Τ!!!!!!! Η34;!97>3! Ρ74; 6Β3Σ!Ρ3Θ3Ρ6!?8Ρ366!47ΡΡ56578Χ!!!!!!!!!!!!(Τ!+ Ι()ΧΕ Χ Ε& Κ0!Υ&! &!!!!!!!!!!! )2!(Τ!!!!!!! 7ΣΑ!7>5<58:Ρ! Ρ74;!Β593!:Β!38Π!7=!Κ!Β>3:Β938ΒΧ!!!!!!!!!! (Ι!Ω! Ε& Κ!!!!!! )2!(Τ!!!Γ/!% #%)!!!!!! )2

13 Massive KS Binary Initial condition e, a, m, m 2, c GR coalescence R coal = 6(m + m 2 ) c 2 PN perturbation PN & 2.5PN Isolated binary e = 0.99, t/t K = 0 5 Time-step τ = η u 2 h /2, N = 2π η u η u = 0., it = 3, E/E = η u = 0.3, it = 3, E/E = 0 8 η u = 0.3, it = 4, E/E = 3 0

14 NBODY7 Binary Project Initial condition e, a, m, m 2, c Perturbation PN & 2.5PN in pnpert.f Energy loss Perturbed KS or Peters 964 Accumulated Hermite, previous Ė GR, Ë GR Integration E GR = 2 (Ė 0 + Ė) t + 2 (Ë 0 Ë) t 2 Identity check h < 4 h

15 PN PN+2.5PN Peters SEMI e-05 e-06 e Time

16 Conclusions Performance Significant speed-up Accuracy One or two extra iterations Memory No new variables needed Time-steps Similar to current method Usage Special procedures require attention

Post-Newtonian N-body Codes. Sverre Aarseth. Institute of Astronomy, Cambridge

Post-Newtonian N-body Codes. Sverre Aarseth. Institute of Astronomy, Cambridge Post-Newtonian N-body Codes Sverre Aarseth Institute of Astronomy, Cambridge Introduction N-body tools Three-body formulation PN implementations Numerical examples Discussion Hermite Integration Taylor

More information

. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Friday April 1 ± ǁ

. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Friday April 1 ± ǁ . α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω Friday April 1 ± ǁ 1 Chapter 5. Photons: Covariant Theory 5.1. The classical fields 5.2. Covariant

More information

and in each case give the range of values of x for which the expansion is valid.

and in each case give the range of values of x for which the expansion is valid. α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Mathematics is indeed dangerous in that it absorbs students to such a degree that it dulls their senses to everything else P Kraft Further Maths A (MFPD)

More information

ASSIGNMENT COVER SHEET omicron

ASSIGNMENT COVER SHEET omicron ASSIGNMENT COVER SHEET omicron Name Question Done Backpack Ready for test Drill A differentiation Drill B sketches Drill C Partial fractions Drill D integration Drill E differentiation Section A integration

More information

. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Wednesday March 30 ± ǁ

. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Wednesday March 30 ± ǁ . α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω Wednesday March 30 ± ǁ 1 Chapter 5. Photons: Covariant Theory 5.1. The classical fields 5.2. Covariant

More information

K e sub x e sub n s i sub o o f K.. w ich i sub s.. u ra to the power of m i sub fi ed.. a sub t to the power of a

K e sub x e sub n s i sub o o f K.. w ich i sub s.. u ra to the power of m i sub fi ed.. a sub t to the power of a - ; ; ˆ ; q x ; j [ ; ; ˆ ˆ [ ˆ ˆ ˆ - x - - ; x j - - - - - ˆ x j ˆ ˆ ; x ; j κ ˆ - - - ; - - - ; ˆ σ x j ; ˆ [ ; ] q x σ; x - ˆ - ; J -- F - - ; x - -x - - x - - ; ; 9 S j P R S 3 q 47 q F x j x ; [ ]

More information

4sec 2xtan 2x 1ii C3 Differentiation trig

4sec 2xtan 2x 1ii C3 Differentiation trig A Assignment beta Cover Sheet Name: Question Done Backpack Topic Comment Drill Consolidation i C3 Differentiation trig 4sec xtan x ii C3 Differentiation trig 6cot 3xcosec 3x iii C3 Differentiation trig

More information

A2 Assignment zeta Cover Sheet. C3 Differentiation all methods. C3 Sketch and find range. C3 Integration by inspection. C3 Rcos(x-a) max and min

A2 Assignment zeta Cover Sheet. C3 Differentiation all methods. C3 Sketch and find range. C3 Integration by inspection. C3 Rcos(x-a) max and min A Assignment zeta Cover Sheet Name: Question Done Backpack Ready? Topic Comment Drill Consolidation M1 Prac Ch all Aa Ab Ac Ad Ae Af Ag Ah Ba C3 Modulus function Bb C3 Modulus function Bc C3 Modulus function

More information

Beauty Contests and Iterated Expectations in Asset Markets. Franklin Allen Stephen Morris Hyun Song Shin

Beauty Contests and Iterated Expectations in Asset Markets. Franklin Allen Stephen Morris Hyun Song Shin Beauty Contests and Iterated Expectations in Asset Markets Franklin Allen Stephen Morris Hyun Song Shin 1 professional investment may be likened to those newspaper competitions in which the competitors

More information

Numerical Approximation of Phase Field Models

Numerical Approximation of Phase Field Models Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School

More information

Further Maths A2 (M2FP2D1) Assignment ψ (psi) A Due w/b 19 th March 18

Further Maths A2 (M2FP2D1) Assignment ψ (psi) A Due w/b 19 th March 18 α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω The mathematician s patterns, like the painter s or the poet s, must be beautiful: the ideas, like the colours or the words, must fit together in a harmonious

More information

b e i ga set s oane f dast heco mm on n va ns ing lo c u soft w section

b e i ga set s oane f dast heco mm on n va ns ing lo c u soft w section 66 M M Eq: 66 - -I M - - - -- -- - - - - -- - S I T - I S W q - I S T q ] q T G S W q I ] T G ˆ Gα ˆ ˆ ] H Z ˆ T α 6H ; Z - S G W [6 S q W F G S W F S W S W T - I ] T ˆ T κ G Gα ±G κ α G ˆ + G > H O T

More information

CSE 1400 Applied Discrete Mathematics Definitions

CSE 1400 Applied Discrete Mathematics Definitions CSE 1400 Applied Discrete Mathematics Definitions Department of Computer Sciences College of Engineering Florida Tech Fall 2011 Arithmetic 1 Alphabets, Strings, Languages, & Words 2 Number Systems 3 Machine

More information

A2 Assignment lambda Cover Sheet. Ready. Done BP. Question. Aa C4 Integration 1 1. C4 Integration 3

A2 Assignment lambda Cover Sheet. Ready. Done BP. Question. Aa C4 Integration 1 1. C4 Integration 3 A Assignment lambda Cover Sheet Name: Question Done BP Ready Topic Comment Drill Mock Exam Aa C4 Integration sin x+ x+ c 4 Ab C4 Integration e x + c Ac C4 Integration ln x 5 + c Ba C Show root change of

More information

License: Creative Commons Attribution Non-commercial No Derivatives

License: Creative Commons Attribution Non-commercial No Derivatives Title Type URL Design in the Time of Policy Problems Article Date 2016 Citation Creators http://ualresearchonline.arts.ac.uk/9542/ Kimbell, Lucy (2016) Design in the Time of Policy Problems. Proceedings

More information

Differential Equations

Differential Equations Differential Equations Definitions Finite Differences Taylor Series based Methods: Euler Method Runge-Kutta Methods Improved Euler, Midpoint methods Runge Kutta (2nd, 4th order) methods Predictor-Corrector

More information

Rømer Science Mission Plan

Rømer Science Mission Plan Institute of Physics and Astronomy, University of Aarhus Rømer Science Mission Plan Danish Small Satellite Programme Document No.(issue): MONS/IFA/MAN/PLN/0001(1) Date: 2001-05-29 Prepared by: Jørgen Christensen-Dalsgaard,

More information

Post-Newtonian N-body Dynamics

Post-Newtonian N-body Dynamics Minor Master Project Post-Newtonian N-body Dynamics Author: Emiel Por por@strw.leidenuniv.nl Supervisors: Simon F. Portegies Zwart Adrian S. Hamers October 2014 1 Post-Newtonian N-body Dynamics Emiel Por

More information

Numerical Methods in Physics and Astrophysics

Numerical Methods in Physics and Astrophysics Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation

More information

Dynamic Asset Allocation - Identifying Regime Shifts in Financial Time Series to Build Robust Portfolios

Dynamic Asset Allocation - Identifying Regime Shifts in Financial Time Series to Build Robust Portfolios Downloaded from orbit.dtu.dk on: Jan 22, 2019 Dynamic Asset Allocation - Identifying Regime Shifts in Financial Time Series to Build Robust Portfolios Nystrup, Peter Publication date: 2018 Document Version

More information

Y1 Double Maths Assignment λ (lambda) Exam Paper to do Core 1 Solomon C on the VLE. Drill

Y1 Double Maths Assignment λ (lambda) Exam Paper to do Core 1 Solomon C on the VLE. Drill α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Nature is an infinite sphere of which the centre is everywhere and the circumference nowhere Blaise Pascal Y Double Maths Assignment λ (lambda) Tracking

More information

Mathematics Review Exercises. (answers at end)

Mathematics Review Exercises. (answers at end) Brock University Physics 1P21/1P91 Mathematics Review Exercises (answers at end) Work each exercise without using a calculator. 1. Express each number in scientific notation. (a) 437.1 (b) 563, 000 (c)

More information

Numerical Methods in Physics and Astrophysics

Numerical Methods in Physics and Astrophysics Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation

More information

A Method for Solving and Estimating Heterogenous Agent Macro Models (by Thomas Winberry)

A Method for Solving and Estimating Heterogenous Agent Macro Models (by Thomas Winberry) A Method for Solving and Estimating Heterogenous Agent Macro Models (by Thomas Winberry) Miguel Angel Cabello Macro Reading Group October, 2018 Miguel Cabello Perturbation and HA Models 1 / 95 Outline

More information

PRELIMINARIES FOR GENERAL TOPOLOGY. Contents

PRELIMINARIES FOR GENERAL TOPOLOGY. Contents PRELIMINARIES FOR GENERAL TOPOLOGY DAVID G.L. WANG Contents 1. Sets 2 2. Operations on sets 3 3. Maps 5 4. Countability of sets 7 5. Others a mathematician knows 8 6. Remarks 9 Date: April 26, 2018. 2

More information

Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs

Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Barbara Kaltenbacher, Tuomo Valkonen,

More information

Lattice Bhatnagar Gross Krook model for the Lorenz attractor

Lattice Bhatnagar Gross Krook model for the Lorenz attractor Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,

More information

STOP, a i+ 1 is the desired root. )f(a i) > 0. Else If f(a i+ 1. Set a i+1 = a i+ 1 and b i+1 = b Else Set a i+1 = a i and b i+1 = a i+ 1

STOP, a i+ 1 is the desired root. )f(a i) > 0. Else If f(a i+ 1. Set a i+1 = a i+ 1 and b i+1 = b Else Set a i+1 = a i and b i+1 = a i+ 1 53 17. Lecture 17 Nonlinear Equations Essentially, the only way that one can solve nonlinear equations is by iteration. The quadratic formula enables one to compute the roots of p(x) = 0 when p P. Formulas

More information

ECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations

ECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent

More information

A Dynamic Model for Investment Strategy

A Dynamic Model for Investment Strategy A Dynamic Model for Investment Strategy Richard Grinold Stanford Conference on Quantitative Finance August 18-19 2006 Preview Strategic view of risk, return and cost Not intended as a portfolio management

More information

Symbols and dingbats. A 41 Α a 61 α À K cb ➋ à esc. Á g e7 á esc. Â e e5 â. Ã L cc ➌ ã esc ~ Ä esc : ä esc : Å esc * å esc *

Symbols and dingbats. A 41 Α a 61 α À K cb ➋ à esc. Á g e7 á esc. Â e e5 â. Ã L cc ➌ ã esc ~ Ä esc : ä esc : Å esc * å esc * Note: Although every effort ws tken to get complete nd ccurte tble, the uhtor cn not be held responsible for ny errors. Vrious sources hd to be consulted nd MIF hd to be exmined to get s much informtion

More information

Computational Astrophysics

Computational Astrophysics Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2: The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear

More information

Nonlinear wave-wave interactions involving gravitational waves

Nonlinear wave-wave interactions involving gravitational waves Nonlinear wave-wave interactions involving gravitational waves ANDREAS KÄLLBERG Department of Physics, Umeå University, Umeå, Sweden Thessaloniki, 30/8-5/9 2004 p. 1/38 Outline Orthonormal frames. Thessaloniki,

More information

Review Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester

Review Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE

More information

Parameterized Expectations Algorithm

Parameterized Expectations Algorithm Parameterized Expectations Algorithm Wouter J. Den Haan London School of Economics c by Wouter J. Den Haan Overview Two PEA algorithms Explaining stochastic simulations PEA Advantages and disadvantages

More information

AN OVERVIEW. Numerical Methods for ODE Initial Value Problems. 1. One-step methods (Taylor series, Runge-Kutta)

AN OVERVIEW. Numerical Methods for ODE Initial Value Problems. 1. One-step methods (Taylor series, Runge-Kutta) AN OVERVIEW Numerical Methods for ODE Initial Value Problems 1. One-step methods (Taylor series, Runge-Kutta) 2. Multistep methods (Predictor-Corrector, Adams methods) Both of these types of methods are

More information

Data Analysis Question Sheet

Data Analysis Question Sheet MSE 11: Mathematics: Sam Cooper and David Dye Data Analysis Question Sheet Comments and corrections to david.dye@imperial.ac.uk Exercises Logs and curve sketching 1. Simplify the following expressions:

More information

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction

More information

Asynchronous Training in Wireless Sensor Networks

Asynchronous Training in Wireless Sensor Networks u t... t. tt. tt. u.. tt tt -t t t - t, t u u t t. t tut t t t t t tt t u t ut. t u, t tt t u t t t t, t tt t t t, t t t t t. t t tt u t t t., t- t ut t t, tt t t tt. 1 tut t t tu ut- tt - t t t tu tt-t

More information

1 Integration of Rational Functions Using Partial Fractions

1 Integration of Rational Functions Using Partial Fractions MTH Fall 008 Essex County College Division of Mathematics Handout Version 4 September 8, 008 Integration of Rational Functions Using Partial Fractions In the past it was far more usual to simplify or combine

More information

Outline. Logic. Definition. Theorem (Gödel s Completeness Theorem) Summary of Previous Week. Undecidability. Unification

Outline. Logic. Definition. Theorem (Gödel s Completeness Theorem) Summary of Previous Week. Undecidability. Unification Logic Aart Middeldorp Vincent van Oostrom Franziska Rapp Christian Sternagel Department of Computer Science University of Innsbruck WS 2017/2018 AM (DCS @ UIBK) week 11 2/38 Definitions elimination x φ

More information

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction

More information

Chapter 7: Natural Convection

Chapter 7: Natural Convection 7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection

More information

5.1 Classical Harmonic Oscillator

5.1 Classical Harmonic Oscillator Chapter 5 Harmonic Oscillator 5.1 Classical Harmonic Oscillator m l o l Hooke s Law give the force exerting on the mass as: f = k(l l o ) where l o is the equilibrium length of the spring and k is the

More information

T m / A. Table C2 Submicroscopic Masses [2] Symbol Meaning Best Value Approximate Value

T m / A. Table C2 Submicroscopic Masses [2] Symbol Meaning Best Value Approximate Value APPENDIX C USEFUL INFORMATION 1247 C USEFUL INFORMATION This appendix is broken into several tables. Table C1, Important Constants Table C2, Submicroscopic Masses Table C3, Solar System Data Table C4,

More information

Nonlinear equations. Norms for R n. Convergence orders for iterative methods

Nonlinear equations. Norms for R n. Convergence orders for iterative methods Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector

More information

Milnor s Exotic 7-Spheres Jhan-Cyuan Syu ( June, 2017 Introduction

Milnor s Exotic 7-Spheres Jhan-Cyuan Syu ( June, 2017 Introduction 許展銓 1 2 3 R 4 R 4 R n n 4 R 4 ξ : E π M 2n J : E E R J J(v) = v v E ξ ξ ξ R ξ : E π M M ξ : Ē π M ξ ξ R = ξ R i : E Ē i(cv) = ci(v) c C v E ξ : E π M n M c(ξ) = c i (ξ) H 2n (M, Z) i 0 c 0 (ξ) = 1 c i

More information

From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction

From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy

More information

Jim Lambers MAT 772 Fall Semester Lecture 21 Notes

Jim Lambers MAT 772 Fall Semester Lecture 21 Notes Jim Lambers MAT 772 Fall Semester 21-11 Lecture 21 Notes These notes correspond to Sections 12.6, 12.7 and 12.8 in the text. Multistep Methods All of the numerical methods that we have developed for solving

More information

MEMORY EFFECT IN HOMOGENIZATION OF VISCOELASTIC KELVIN-VOIGT MODEL WITH TIME DEPENDENT COEFFICIENTS

MEMORY EFFECT IN HOMOGENIZATION OF VISCOELASTIC KELVIN-VOIGT MODEL WITH TIME DEPENDENT COEFFICIENTS Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company MEMORY EFFECT IN HOMOGENIZATION OF VISCOELASTIC KELVIN-VOIGT MODEL WITH TIME DEPENDENT COEFFICIENTS ZOUHAIR ABDESSAMAD,

More information

(1) u (t) = f(t, u(t)), 0 t a.

(1) u (t) = f(t, u(t)), 0 t a. I. Introduction 1. Ordinary Differential Equations. In most introductions to ordinary differential equations one learns a variety of methods for certain classes of equations, but the issues of existence

More information

Persistence of Plummer-Distributed Small Globular Clusters as a Function of Primordial- Binary Population Size

Persistence of Plummer-Distributed Small Globular Clusters as a Function of Primordial- Binary Population Size Persistence of Plummer-Distributed Small Globular Clusters as a Function of Primordial- Binary Population Size Jack K. Horner P.O. Box 266 Los Alamos NM 87544 USA CSC 2013 Abstract Globular stellar clusters

More information

Wavelets For Computer Graphics

Wavelets For Computer Graphics {f g} := f(x) g(x) dx A collection of linearly independent functions Ψ j spanning W j are called wavelets. i J(x) := 6 x +2 x + x + x Ψ j (x) := Ψ j (2 j x i) i =,..., 2 j Res. Avge. Detail Coef 4 [9 7

More information

A class of non-convex polytopes that admit no orthonormal basis of exponentials

A class of non-convex polytopes that admit no orthonormal basis of exponentials A class of non-convex polytopes that admit no orthonormal basis of exponentials Mihail N. Kolountzakis and Michael Papadimitrakis 1 Abstract A conjecture of Fuglede states that a bounded measurable set

More information

Hoover Institution, Stanford University and NBER. University of Alicante and Hoover Institution, Stanford University

Hoover Institution, Stanford University and NBER. University of Alicante and Hoover Institution, Stanford University Supplementary Material Supplement to Numerically stable and accurate stochastic simulation approaches for solving dynamic economic models : Appendices (Quantitative Economics, Vol. 2, No. 2, July 2011,

More information

SHANGHAI JIAO TONG UNIVERSITY LECTURE

SHANGHAI JIAO TONG UNIVERSITY LECTURE Lecture 7 SHANGHAI JIAO TONG UNIVERSITY LECTURE 7 017 Anthony J. Leggett Department of Physics University of Illinois at Urbana-Champaign, USA and Director, Center for Complex Physics Shanghai Jiao Tong

More information

2 Post-Keplerian Timing Parameters for General Relativity

2 Post-Keplerian Timing Parameters for General Relativity 1 Introduction General Relativity has been one of the pilars of modern physics for over 100 years now. Testing the theory and its consequences is therefore very important to solidifying our understand

More information

Application of the perturbation iteration method to boundary layer type problems

Application of the perturbation iteration method to boundary layer type problems DOI 10.1186/s40064-016-1859-4 RESEARCH Open Access Application of the perturbation iteration method to boundary layer type problems Mehmet Pakdemirli * *Correspondence: mpak@cbu.edu.tr Applied Mathematics

More information

Stellar Dynamics and Structure of Galaxies

Stellar Dynamics and Structure of Galaxies Stellar Dynamics and Structure of Galaxies in a given potential Vasily Belokurov vasily@ast.cam.ac.uk Institute of Astronomy Lent Term 2016 1 / 59 1 Collisions Model requirements 2 in spherical 3 4 Orbital

More information

Do not turn over until you are told to do so by the Invigilator.

Do not turn over until you are told to do so by the Invigilator. UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.

More information

3.4 Numerical orbit integration

3.4 Numerical orbit integration 196 Chapter 3: The Orbits of Stars about the center of the galaxy. For small-amplitude librations, the libration frequency is p, consistent with our assumption that the oscillation frequency is of order

More information

QFT Perturbation Theory

QFT Perturbation Theory QFT Perturbation Theory Ling-Fong Li (Institute) Slide_04 1 / 43 Interaction Theory As an illustration, take electromagnetic interaction. Lagrangian density is The combination L = ψ (x ) γ µ ( i µ ea µ

More information

Analysis Methods in Atmospheric and Oceanic Science

Analysis Methods in Atmospheric and Oceanic Science Analysis Methods in Atmospheric and Oceanic Science AOSC 652 Week 7, Day 1 13 Oct 2014 1 Student projects: 20% of the final grade: you will receive a numerical score for the project and final grade will

More information

A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit

A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit David Vener Department of Mathematics, MIT May 5, 3 Introduction In 8.366, we discussed the relationship

More information

VARIATIONAL PRINCIPLE FOR THE ENTROPY

VARIATIONAL PRINCIPLE FOR THE ENTROPY VARIATIONAL PRINCIPLE FOR THE ENTROPY LUCIAN RADU. Metric entropy Let (X, B, µ a measure space and I a countable family of indices. Definition. We say that ξ = {C i : i I} B is a measurable partition if:

More information

R k. t + 1. n E t+1 = ( 1 χ E) W E t+1. c E t+1 = χ E Wt+1 E. Γ E t+1. ) R E t+1q t K t. W E t+1 = ( 1 Γ E t+1. Π t+1 = P t+1 /P t

R k. t + 1. n E t+1 = ( 1 χ E) W E t+1. c E t+1 = χ E Wt+1 E. Γ E t+1. ) R E t+1q t K t. W E t+1 = ( 1 Γ E t+1. Π t+1 = P t+1 /P t R k E 1 χ E Wt E n E t+1 t t + 1 n E t+1 = ( 1 χ E) W E t+1 c E t+1 = χ E Wt+1 E t + 1 q t K t Rt+1 E 1 Γ E t+1 Π t+1 = P t+1 /P t W E t+1 = ( 1 Γ E t+1 ) R E t+1q t K t Π t+1 Γ E t+1 K t q t q t K t j

More information

Initial value problems for ordinary differential equations

Initial value problems for ordinary differential equations Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP

More information

ASSIGNMENT BOOKLET. Numerical Analysis (MTE-10) (Valid from 1 st July, 2011 to 31 st March, 2012)

ASSIGNMENT BOOKLET. Numerical Analysis (MTE-10) (Valid from 1 st July, 2011 to 31 st March, 2012) ASSIGNMENT BOOKLET MTE-0 Numerical Analysis (MTE-0) (Valid from st July, 0 to st March, 0) It is compulsory to submit the assignment before filling in the exam form. School of Sciences Indira Gandhi National

More information

Numerical Methods - Initial Value Problems for ODEs

Numerical Methods - Initial Value Problems for ODEs Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value

More information

Time Evolution of Matter States

Time Evolution of Matter States Tie Evolution of Matter States W. M. Hetherington February 15, 1 The Tie-Evolution Operat The tie-evolution of a wavefunction is deterined by the effect of a tie evolution operat through the relation Ψ

More information

CHAPTER 6 : LITERATURE REVIEW

CHAPTER 6 : LITERATURE REVIEW CHAPTER 6 : LITERATURE REVIEW Chapter : LITERATURE REVIEW 77 M E A S U R I N G T H E E F F I C I E N C Y O F D E C I S I O N M A K I N G U N I T S A B S T R A C T A n o n l i n e a r ( n o n c o n v e

More information

P E R E N C O - C H R I S T M A S P A R T Y

P E R E N C O - C H R I S T M A S P A R T Y L E T T I C E L E T T I C E I S A F A M I L Y R U N C O M P A N Y S P A N N I N G T W O G E N E R A T I O N S A N D T H R E E D E C A D E S. B A S E D I N L O N D O N, W E H A V E T H E P E R F E C T R

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 14: Formulation of the Stability Problem

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 14: Formulation of the Stability Problem .65, MHD Theory of Fusion Systems Prof. Freidberg Lecture 4: Formulation of the Stability Problem Hierarchy of Formulations of the MHD Stability Problem for Arbitrary 3-D Systems. Linearized equations

More information

Input to state Stability

Input to state Stability Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part III: Lyapunov functions and quantitative aspects ISS Consider with

More information

DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science

DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI097 TU) Thursday January 014, 18:0-1:0

More information

(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by

(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by 1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How

More information

Simulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients

Simulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients Journal of mathematics and computer Science 8 (2014) 28-32 Simulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients Behrouz Fathi Vajargah Department of statistics,

More information

Variations on Backpropagation

Variations on Backpropagation 2 Variations on Backpropagation 2 Variations Heuristic Modifications Moentu Variable Learning Rate Standard Nuerical Optiization Conjugate Gradient Newton s Method (Levenberg-Marquardt) 2 2 Perforance

More information

K(ζ) = 4ζ 2 x = 20 Θ = {θ i } Θ i=1 M = {m i} M i=1 A = {a i } A i=1 M A π = (π i ) n i=1 (Θ) n := Θ Θ (a, θ) u(a, θ) E γq [ E π m [u(a, θ)] ] C(π, Q) Q γ Q π m m Q m π m a Π = (π) U M A Θ

More information

Research Article On Existence, Uniform Decay Rates, and Blow-Up for Solutions of a Nonlinear Wave Equation with Dissipative and Source

Research Article On Existence, Uniform Decay Rates, and Blow-Up for Solutions of a Nonlinear Wave Equation with Dissipative and Source Abstract and Applied Analysis Volume, Article ID 65345, 7 pages doi:.55//65345 Research Article On Existence, Uniform Decay Rates, and Blow-Up for Solutions of a Nonlinear Wave Equation with Dissipative

More information

Periodic orbits high above the ecliptic plane in the solar sail 3-body problem

Periodic orbits high above the ecliptic plane in the solar sail 3-body problem Periodic orbits high above the ecliptic plane in the solar sail 3-body problem Thomas Waters Department of Mechanical Engineering, University of Strathclyde, Glasgow. 4th October 2006 In conjunction with

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference

More information

MAS114: Exercises. October 26, 2018

MAS114: Exercises. October 26, 2018 MAS114: Exercises October 26, 2018 Note that the challenge problems are intended to be difficult! Doing any of them is an achievement. Please hand them in on a separate piece of paper if you attempt them.

More information

Lecture 3 Optimization methods for econometrics models

Lecture 3 Optimization methods for econometrics models Lecture 3 Optimization methods for econometrics models Cinzia Cirillo University of Maryland Department of Civil and Environmental Engineering 06/29/2016 Summer Seminar June 27-30, 2016 Zinal, CH 1 Overview

More information

High-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation

High-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y

More information

ON ESTIMATION OF TEMPERATURE UNCERTAINTY USING THE SECOND ORDER ADJOINT PROBLEM

ON ESTIMATION OF TEMPERATURE UNCERTAINTY USING THE SECOND ORDER ADJOINT PROBLEM ON ESTIMATION OF TEMPERATURE UNCERTAINTY USING THE SECOND ORDER ADJOINT PROBLEM Aleksey K. Alekseev a, Michael I. Navon b,* a Department of Aerodynamics and Heat Transfer, RSC, ENERGIA, Korolev (Kaliningrad),

More information

Theoretical Statistical Correlation for Biometric Identification Performance

Theoretical Statistical Correlation for Biometric Identification Performance Theoretical Statistical Correlation for Biometric Identification Performance Michael E. Schuckers schuckers@stlawu.edu Theoretical Statistical Correlation for Biometric Identification Performance p. 1/20

More information

Talk 3: Examples: Star exponential

Talk 3: Examples: Star exponential Akira Yoshioka 2 7 February 2013 Shinshu Akira Yoshioka 3 2 Dept. Math. Tokyo University of Science 1. Star exponential 1 We recall the definition of star exponentials. 2 We show some construction of star

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme

More information

Candidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.

Candidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request. UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book

More information

Massive star clusters

Massive star clusters Massive star clusters as a host of compact binaries Michiko Fujii ( 藤井通子 ) The University of Tokyo Outline Massive star clusters and compact binaries Dynamical evolution of star clusters Distribution of

More information

Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US

Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US Gerdie Everaert 1, Lorenzo Pozzi 2, and Ruben Schoonackers 3 1 Ghent University & SHERPPA 2 Erasmus

More information

A sharp diffuse interface tracking method for approximating evolving interfaces

A sharp diffuse interface tracking method for approximating evolving interfaces A sharp diffuse interface tracking method for approximating evolving interfaces Vanessa Styles and Charlie Elliott University of Sussex Overview Introduction Phase field models Double well and double obstacle

More information

Weak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij

Weak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij Weak convergence and Brownian Motion (telegram style notes) P.J.C. Spreij this version: December 8, 2006 1 The space C[0, ) In this section we summarize some facts concerning the space C[0, ) of real

More information

Electromagnetically Induced Flows in Water

Electromagnetically Induced Flows in Water Electromagnetically Induced Flows in Water Michiel de Reus 8 maart 213 () Electromagnetically Induced Flows 1 / 56 Outline 1 Introduction 2 Maxwell equations Complex Maxwell equations 3 Gaussian sources

More information

Planet formation (aero)dynamics project

Planet formation (aero)dynamics project Planet formation (aero)dynamics project Chris Ormel, Carsten Dominik, Jacob Arcangeli February 2, 2016 1 Goals This project consist of two steps: 1. Program a small N-body code and test several algorithms

More information

Waveform inversion and time-reversal imaging in attenuative TI media

Waveform inversion and time-reversal imaging in attenuative TI media Waveform inversion and time-reversal imaging in attenuative TI media Tong Bai 1, Tieyuan Zhu 2, Ilya Tsvankin 1, Xinming Wu 3 1. Colorado School of Mines 2. Penn State University 3. University of Texas

More information

SECOND GENERAL REPORT - APRIL 1954 35 During the same period, France and the Saar increased their hard coal purchases from other countries by 689,000 metric tons, or 14.4 %, with 5,486,000 tons in 1953

More information

Mark Gales October y (x) x 1. x 2 y (x) Inputs. Outputs. x d. y (x) Second Output layer layer. layer.

Mark Gales October y (x) x 1. x 2 y (x) Inputs. Outputs. x d. y (x) Second Output layer layer. layer. University of Cambridge Engineering Part IIB & EIST Part II Paper I0: Advanced Pattern Processing Handouts 4 & 5: Multi-Layer Perceptron: Introduction and Training x y (x) Inputs x 2 y (x) 2 Outputs x

More information