Accuracy of Symmetric Partitioned Runge-Kutta Methods for Differential Equations on Lie-Groups
|
|
- Bethany Hardy
- 5 years ago
- Views:
Transcription
1 Accuracy of Symmetric Partitioned Runge-Kutta Method for Differential Equation on Lie-Group WUB / BUW-IMACM 11-19, Michael Günther, Franceco Knechtli and Michèle Wandelt Bergiche Univerität Wuppertal, Faculty of Mathematic and Natural Science, Germany {triebel, guenther, wandelt}@math.uni-wuppertal.de, knechtli@phyik.uni-wuppertal.de Computer imulation in QCD are baed on the dicretization of the theory on a Euclidean lattice. To compute the mean value of an obervable, uually the Hybrid Monte Carlo method i applied. Here equation of motion, derived from an Hamiltonian, have to be olved numerically. Commonly the Leapfrog (Stoermer-Verlet) method or plitting method with multiple timecale à la Sexton-Weingarten are ued to integrate the dynamical ytem, defined on a Lie group. Here we formulate time-reverible higher order integrator baed on implicit partitioned Runge- Kutta cheme and how that they allow for larger tep-ize than the Leapfrog method. Since thee method are baed on an infinite erie of exponential function, we concentrate on the truncation of thi erie with repect to the global error and accuracy. Finally, we ee that the global error of a SPRK cheme i alway even uch that a convergence order of one i gained for method with odd convergence order. The XXIX International Sympoium on Lattice Field Theory - Lattice 211 July 1-16, 211 Squaw Valley, Lake Tahoe, California Speaker. Thi work wa upported by the Deutche Forchunggemeinchaft through the Collaborative Reearch Centre SFB-TR 55 Hadron phyic from Lattice QCD. Furthermore, acknowledge the Marie Curie Initial Training Network STRONGnet Strong Interaction Supercomputing Training Network for travel upport. c Copyright owned by the author() under the term of the Creative Common Attribution-NonCommercial-ShareAlike Licence.
2 1. Introduction and Motivation In the molecular dynamic tep of the Hybrid Monte Carlo method [1], Hamiltonian equation of motion have to be olved. Thee equation form coupled ytem of matrix differential equation of the form H ([U],[A]) U ν = = A ν U ν, A ν (1.1a) H ([U],[A]) Ȧ ν = = g(u ν ), U ν for ν = 1,...,n. (1.1b) In thi notation, U ν i an element of a matrix Lie group G and A ν an element of it aociated Lie algebra g. Thu, [U] can be imagined a a vector of n matrix Lie group element U 1,U 2,...,U n, and [A] a a vector of n Lie algebra element A 1,A 2,...,A n. In pure lattice gauge theory, the element U ν can be een a the link U x,µ between the lattice ite x and x + a ˆµ. Thu, A µ i it aociated momentum P x,µ time the complex i. In thi context, the vector of matrice [U] and [A] are the whole configuration of the link and it momenta. The equation of motion have to be olved in a Lie group, repectively in a Lie algebra with a timereverible and area-preerving cheme. In a recent paper [2], we have invetigated the potentialitie of higher order partitioned Runge-Kutta cheme for olving the equation of motion uch that the deired propertie are met. We found out that ymmetric partitioned Runge-Kutta method baed on the Magnu and Munthe-Kaa approach can be time-reverible. So far, area-preervation i not fulfilled and mut be corrected in the acceptance tep. Furthermore, the global error of thi cheme i alway even and invetigated in detail in thi paper. In doing o, we tart with a hort derivation of ymmetric partitioned Runge-Kutta cheme baed on the idea of Magnu and Munthe-Kaa. Afterward, we focu on the global error and accuracy of the method and how ome numerical reult. 2. Numerical Integration The differential equation (1.1) become an initial value problem (IVP) by precribed initial value: U ν () := U ν, and A ν () := A ν, for ν = 1,...,n. Thereby, the initial value U ν () have to be in the Lie group and the element A ν () in the Lie algebra. Conidering the tructure of equation (1.1), (1.1a) i a differential equation in a Lie group uch that it ha to be olved with a numerical cheme that guarantee a olution in the Lie group a decribed in paragraph 2.1 For the econd equation (1.1b), no pecial treatment ha to be applied. It i an equation in the Lie algebra g, which i a linear pace. Thu, thi equation can be olved with any time-reverible and areapreerving numerical cheme. For convenience, we leave out the index ν from now on. Thi mean we invetigate jut one coupled differential equation for a pecial but arbitrary index ν: U = A U and Ȧ = g(u). The reult can then be extended traightforward to the whole vector [U] and [A]. 2.1 Differential Equation in Lie Group Concerning equation (1.1a), we follow the idea of Magnu and Munthe-Kaa. Magnu [3] tated that the differential equation (1.1a) in the Lie group can be replaced through a differential equation in the Lie algebra. Thi new differential equation can be olved directly due to the linearity 2
3 of the Lie algebra. The trategy i a follow: Identify U(t) with exp(ω(t)) uch that the variable change from U to Ω. Ω i the olution of the differential equation Ω = d exp 1 Ω (A), (2.1) with Ω(t) g and initial value Ω() :=. The derivative of the invere exponential map (2.1) i given by an infinite erie a d expω 1 (A) = k k! adk Ω (A). In thi erie, the variable are the k-th Bernoulli number and the adjoint operator adω k i a mapping in the Lie algebra g given by ad Ω (A) := [Ω,A] = ΩA AΩ. It follow the convention ad Ω (A) = A and adk Ω (A) = [Ω,adk 1(A)]. Thi mean, Ω i the olution of the differential equation Ω Ω = k k! adk Ω (A). Knowing Ω, the olution U of (1.1a) can be attained via U = exp(ω)u. In total, we record that the initial value problem (1.1) i equivalent to Ω = k= k! adk Ω (A), Ȧ = g(u) with U = exp(ω)u (2.2) and U() := U G, A() := A g and Ω() := g. Thi tranformed problem can now be olved directly by a Runge-Kutta method without detroying the Lie group tructure: A the Lie algebra g i a vector pace, the analytic olution (Ω(t),A(t)) a well a it approximation (Ω 1,A 1 ) attained by a numerical integration cheme both are element of the Lie algebra g. Furthermore, a for any a g the matrix exponential exp(a) i in the aociated matrix Lie group G, alo U i in G. 2.2 Symmetric Partitioned Runge-Kutta cheme The problem in olving (2.2) i that Ω i given a infinite erie which ha to be uitably truncated after q + 1 term. Thi mean, the truncation index q of Ω ha to be choen properly uch that a numerical integration cheme meet a precribed convergence order p. Thereby, the convergence order of a numerical integration method i p if the deviation between the exact olution and it numerical approximation after one tep i of order p + 1 in a uitable norm. Here, the idea of Munthe-Kaa come into play. He tate in [4] that the truncation index q of Ω ha to be choen a a value larger than the deired convergence order p minu one. Conequently, for a Runge-Kutta cheme of convergence order p, Ω i et a a function depending on the truncation q p 2 of the aforementioned infinite erie, i. e. Ω = q k= k! adk Ω (A) =: f q(ω,a). (2.3) All in all, the exact olution of (2.2) i approximated through an integration cheme of order p of the truncated model Ω = q=p 2 k= k! adk Ω(Â), Â = g(û) with Û = exp( Ω)Û, (2.4) 3
4 Û() := U G, Â() := A g and Ω() := g. Thereby, thi model can be olved with higher order time-reverible ymmetric partitioned Runge-Kutta (SPRK) cheme derived in [2] a follow: Compute the approximation Ω 1 = h i=1 b i K i, A 1 = A + h i=1 bi L i, (2.5) with increment K i = f q ( Ω i,ā i ) and L i = g(ū i ) for i = 1,...,. In the coure of thi, the internal tage are defined a Ω i = h j=1 α i j K j, Ā i = A + h j=1 α i j L j, Ū i = exp( X i )exp ( 1 2 Ω 1) U, X i = h j=1 γ i j K j. At the end, the olution U 1 i attained via U 1 = exp(ω 1 )U. In thi cheme, the coefficient b i, b i,α i j, α i j and γ i j for i, j = 1,..., can be determined to guarantee time-reveribility (and ymmetry). Their value for convergence order p = 3 can be found in [2]. 3. Global Error and Accuracy of the SPRK Method For the local error, the olution of the integration method after one tep ha to be compared with the exact olution U(t + h),a(t + h) of the differential equation (1.1). The SPRK method (2.5) i of convergence order p if U(t + h) U 1 = O(h p+1 ) and A(t + h) A 1 = O(h p+1 ) (3.1) hold. Since the approximation to U 1 i computed from evaluating the matrix exponential (we aume here that we can evaluate thi exactly), which i Lipchitz on every cloed interval, it uffice to demand Ω(t + h) Ω 1 = O(h p+1 ) and A(t + h) A 1 = O(h p+1 ) with exact olution Ω(t + h),a(t + h). According to Munthe-Kaa, the approximation Ω 1 and A 1 of the exact olution of the uitably truncated problem (2.4) can alo be interpreted a approximation to the original problem (2.2). With the ame argument, we can even formulate a tronger tatement on the local accuracy (3.1). A the method i ymmetric, theorem 3.2 in [6, II.3] applie, which tate that the maximal convergence order p of a ymmetric method i even, which mean that the local error i alway odd. Hence, the SPRK method developed a a method of an odd convergence order p i of order p + 1. The global error of a numerical integrated cheme i computed a the um of the local error. Thi mean, for the computation of a trajectory with length τ, an integration method with fixed tep ize h i applied N = τ/h time. Hence, the global error i of the order local error minu one. Thu, the SPRK method of convergence order p ha at leat a global error of order p. Again, in cae of an odd p, the global error i of order p + 1. The accuracy depend on the truncation k = q = p 2 of the erie given in (2.2) with p being the convergence order of the Runge-Kutta method that i ued to olve the problem numerically. We recall the baic tep of the proof given in [6] for a deeper undertanding of the choice of the truncation parameter q. For thi purpoe, we retrict to an uncoupled Lie algebra problem Ω = F (Ω) := k= k! adk Ω (A), with Ω() = g, (3.2) 4
5 which arie from the Magnu approach mentioned in paragraph 2.1. The truncation of the erie in (3.2) at k = q yield the truncated problem Ω = F q ( Ω), with Ω() =, (3.3) uch that F (Ω(t)) F q (Ω(t)) = k=q+1 k! adk Ω(t) (A). For ufficiently mooth A we recognize ad k Ω(t) (A) = O(tk+1 ), which i due to the neted tructure of the ad-operator [6]. Hence, we have F q ( Ω(t)) = F ( Ω(t)) +C(t) with C(t) = c 1 t q+2 + c 2 t q+3 + and contant value c 1,c 2,... g. For fixed h > and t [,h] the exact olution Ω(t) of (3.2) and Ω(t) of (3.3) atify Ω(t) Ω(t) = Ω() + F (Ω(τ))dτ ( = F (Ω(τ))dτ F (Ω(τ)) F ( Ω(τ)) dτ + ( Ω() + ) F q ( Ω(τ))dτ F ( Ω(τ)) +C(τ)dτ) C(τ) dτ. The function F i Lipchitz continuou on every cloed interval for ufficiently mooth A. We aume that for an interval where both Ω(t) and Ω(t) reide in for t [,h], the Lipchitz contant i L R, i. e., Furthermore, for t [,h] we ee that Hence from (3.4) it follow that F (Ω(t)) F ( Ω(t)) L Ω(t) Ω(t). (3.4) C(t) c 1 h q+2 + c 2 h q := c(h) R. (3.5) Ω(t) Ω(t) c(h) t + L Ω(τ) Ω(τ) dτ, uch that the requirement of the "Gronwall lemma" [5] are atified by which Ω(t) Ω(t) c(h) ( ) e L h 1 L = 1 ( ( (1 c1 h q+2 + c 2 h q+3 + ) + L h + 1 L 2! (L h) 2 + ) 1) = ( ( c 1 h q+2 + c 2 h q+3 + ) h L h ) 3! L2 h
6 Thu, the difference between the exact olution of the full problem (3.2) and the truncated problem (3.3) i Ω(h) Ω(h) = O(h q+3 ) (3.6) after one time tep h. Applying a one tep method of convergence order p on the truncated problem (3.3) mean to calculate an approximation Ω 1 to the exact value Ω(h) uch that Ω(h) Ω 1 = O(h p+1 ). (3.7) Finally, we interpret Ω 1 a an approximation to the exact olution of the original problem (3.2). The quality of thi approximation i determined by the deviation (3.6) introduced by the modeling and the dicretization error (3.7): Ω(h) Ω 1 Ω(h) Ω(h) + Ω(h) Ω 1 = O(h q+3 ) + O(h p+1 ). Thi clearly indicate that Ω 1 i a numerical approximation to Ω(h) of convergence order p, i. e., 4. Numerical Tet Ω(h) Ω 1 = O(h p+1 ) if q + 3 p + 1, i. e., q p 2. We conider a pure lattice gauge theory in SU(2,C) with Wilon action and compare the SPRK method decribed in (2.5) with the Leapfrog method. For thi purpoe, we invetigate a ymmetric partitioned Runge-Kutta cheme of convergence order p = 4 which contain the truncated function Ω = f q (Ω,A) given in (2.3). Becaue of the ymmetry, the method ha an even convergence order uch that the choice p = 3,q = 1 already lead to a local error of order 5. Thi mean, we ue the equation Ω = f 1 (Ω,A) = A 1 [Ω,A] (4.1) 2 according to (2.3) and perform imulation on a 2-dimenional lattice with lattice ize L = T = 32. There are 2 reult hown in figure 1: On the left ide, the convergence order of the different method can be een. For thi purpoe, we conider the energy change H of two ucceive configuration after a whole trajectory of length 1 and take the mean of 5 configuration. The Leapfrog SPRK h H 1 5 h 4 Determinant Leapfrog SPRK.995 1/32 1/16 1/8 1/4 1/3 Step ize 1/32 1/16 1/8 1/4 1/3 Step ize Figure 1: Left: Convergence order. Right: Area-preervation. 6
7 tatitical error are o mall that they are not viible in the plot. Since the energy change H deviate from zero jut becaue of the numerical error of the integration cheme, the violation of the energy preervation give the global error. We ee that for a given H the SPRK allow for larger tep ize. On the right ide of figure 1, we ee the violation of the area-preervation in dependence of the tep ize choen in the numerical method. Area-preervation (up to roundoff error) i given if the determinant of (Ω 1,A 1 )/ (Ω,A ) ha exactly the value 1. Here, the determinant i numerically approximated by firt order difference quotient. 5. Concluion We invetigated the accuracy of the time-reverible ymmetric partitioned Runge-Kutta cheme (2.5). The order of accuracy conit of two component: On the one hand, the convergence order depend of coure on the order p of the method itelf. A the method i ymmetric, the local error i alway odd, i. e., the cheme ha a local error of order p+1 for an even convergence order p. On the other hand, the SPRK cheme contain one truncated erie (2.3). The truncation index q ha to be larger than or equal to p 2 to meet the precribed convergence order. All in all, chooing an SPRK method with an even local error hould be preferred ince a convergence order of one i gained by the ymmetry. We performed imulation for an SPRK cheme of convergence order 4 and ee that the global error given in the numerical reult ha order 4 a theoretically expected. In the development of the SPRK method, area-preervation ha not been conidered. Thu, it i not urpriing, that area-preervation i not met applying thi cheme. Thi property ha to be invetigated in future work. Reference [1] S. Duane, A. D. Kennedy, B. J. Pendleton and D. Roweth, Hybrid Monte Carlo, Phyic Letter B, , 195, [2] M. Wandelt, M. Günther, F. Knechtli, M. Striebel, Symmetric Partitioned Runge-Kutta Method, 211, arxiv: [hep-lat]. [3] Wilhelm Magnu, On the exponential olution of differential equation for a linear operator, Communication on Pure and Applied Mathematic, 4, , doi:1.12/cpa [4] Han Munthe-Kaa, High order Runge-Kutta method on manifold, Appl. Numer. Math, 29, , [5] Hairer, E and Nørett, S. P. and Wanner, G., Solving Ordinary Differential Equation I Nontiff Problem, Springer, econd revied, 2. [6] E.Hairer, C.Lubich and G.Wanner, Geometric Numerical Integration Structure-Preerving Algorithm for Ordinary Differential Equation, Springer Ser. Comput. Math. 31, 2nd ed., Springer, 26. [7] W. Kamleh and M. Peardon, Polynomial Filtered HMC an algorithm for lattice QCD with dynamical quark, arxiv: [hep-lat]. [8] M. A. Clark, B. Joó, A. D. Kennedy, P. J. Silva, Improving dynamical lattice QCD imulation through integrator tuning uing Poion bracket and a force-gradient integrator, arxiv: [hep-lat]. 7
Bogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationChapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem
Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,
More informationMATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:
MATEMATIK Datum: 20-08-25 Tid: eftermiddag GU, Chalmer Hjälpmedel: inga A.Heintz Telefonvakt: Ander Martinon Tel.: 073-07926. Löningar till tenta i ODE och matematik modellering, MMG5, MVE6. Define what
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationAvoiding Forbidden Submatrices by Row Deletions
Avoiding Forbidden Submatrice by Row Deletion Sebatian Wernicke, Jochen Alber, Jen Gramm, Jiong Guo, and Rolf Niedermeier Wilhelm-Schickard-Intitut für Informatik, niverität Tübingen, Sand 13, D-72076
More information696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of
Vol 12 No 7, July 2003 cfl 2003 Chin. Phy. Soc. 1009-1963/2003/12(07)/0695-05 Chinee Phyic and IOP Publihing Ltd Lie ymmetrie and conerved quantitie of controllable nonholonomic dynamical ytem Fu Jing-Li(ΛΠ±)
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationAdvanced methods for ODEs and DAEs
Lecture : Implicit Runge Kutta method Bojana Roić, 9. April 7 What you need to know before thi lecture numerical integration: Lecture from ODE iterative olver: Lecture 5-8 from ODE 9. April 7 Bojana Roić
More informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014
Phyic 7 Graduate Quantum Mechanic Solution to inal Eam all 0 Each quetion i worth 5 point with point for each part marked eparately Some poibly ueful formula appear at the end of the tet In four dimenion
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationApproximate Analytical Solution for Quadratic Riccati Differential Equation
Iranian J. of Numerical Analyi and Optimization Vol 3, No. 2, 2013), pp 21-31 Approximate Analytical Solution for Quadratic Riccati Differential Equation H. Aminikhah Abtract In thi paper, we introduce
More informationModeling the scalar wave equation with Nyström methods
GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 200 ; P. T151 T158, 7 FIGS. 10.1190/1.2335505 Modeling the calar wave equation with Nytröm method Jing-Bo Chen 1 ABSTRACT High-accuracy numerical cheme for
More informationTRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL
GLASNIK MATEMATIČKI Vol. 38583, 73 84 TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL p-laplacian Haihen Lü, Donal O Regan and Ravi P. Agarwal Academy of Mathematic and Sytem Science, Beijing, China, National
More informationDesign of Digital Filters
Deign of Digital Filter Paley-Wiener Theorem [ ] ( ) If h n i a caual energy ignal, then ln H e dω< B where B i a finite upper bound. One implication of the Paley-Wiener theorem i that a tranfer function
More informationFeedback Control Systems (FCS)
Feedback Control Sytem (FCS) Lecture19-20 Routh-Herwitz Stability Criterion Dr. Imtiaz Huain email: imtiaz.huain@faculty.muet.edu.pk URL :http://imtiazhuainkalwar.weebly.com/ Stability of Higher Order
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationNumerical Time Integration in Hybrid Monte Carlo Schemes for Lattice Quantum Chromo Dynamics
Numerical Time Integration in Hybrid Monte Carlo Schemes for Lattice Quantum Chromo Dynamics M. Günther and M. Wandelt Bergische Universität Wuppertal IMACM Dortmund, March 12, 2018 Contents Lattice QCD
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationThe Hassenpflug Matrix Tensor Notation
The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationRobustness analysis for the boundary control of the string equation
Routne analyi for the oundary control of the tring equation Martin GUGAT Mario SIGALOTTI and Mariu TUCSNAK I INTRODUCTION AND MAIN RESULTS In thi paper we conider the infinite dimenional ytem determined
More informationNew bounds for Morse clusters
New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationarxiv: v1 [math.ac] 30 Nov 2012
ON MODULAR INVARIANTS OF A VECTOR AND A COVECTOR YIN CHEN arxiv:73v [mathac 3 Nov Abtract Let S L (F q be the pecial linear group over a finite field F q, V be the -dimenional natural repreentation of
More informationCompact finite-difference approximations for anisotropic image smoothing and painting
CWP-593 Compact finite-difference approximation for aniotropic image moothing and painting Dave Hale Center for Wave Phenomena, Colorado School of Mine, Golden CO 80401, USA ABSTRACT Finite-difference
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationLecture 8: Period Finding: Simon s Problem over Z N
Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing
More information5.5 Application of Frequency Response: Signal Filters
44 Dynamic Sytem Second order lowpa filter having tranfer function H()=H ()H () u H () H () y Firt order lowpa filter Figure 5.5: Contruction of a econd order low-pa filter by combining two firt order
More informationMarch 18, 2014 Academic Year 2013/14
POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of
More informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
More informationAssignment for Mathematics for Economists Fall 2016
Due date: Mon. Nov. 1. Reading: CSZ, Ch. 5, Ch. 8.1 Aignment for Mathematic for Economit Fall 016 We now turn to finihing our coverage of concavity/convexity. There are two part: Jenen inequality for concave/convex
More informationMoment of Inertia of an Equilateral Triangle with Pivot at one Vertex
oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.
More informationReformulation of Block Implicit Linear Multistep Method into Runge Kutta Type Method for Initial Value Problem
International Journal of Science and Technology Volume 4 No. 4, April, 05 Reformulation of Block Implicit Linear Multitep Method into Runge Kutta Type Method for Initial Value Problem Muhammad R., Y. A
More informationResearch Article Existence for Nonoscillatory Solutions of Higher-Order Nonlinear Differential Equations
International Scholarly Reearch Network ISRN Mathematical Analyi Volume 20, Article ID 85203, 9 page doi:0.502/20/85203 Reearch Article Exitence for Nonocillatory Solution of Higher-Order Nonlinear Differential
More informationUnbounded solutions of second order discrete BVPs on infinite intervals
Available online at www.tjna.com J. Nonlinear Sci. Appl. 9 206), 357 369 Reearch Article Unbounded olution of econd order dicrete BVP on infinite interval Hairong Lian a,, Jingwu Li a, Ravi P Agarwal b
More informationarxiv: v1 [math.mg] 25 Aug 2011
ABSORBING ANGLES, STEINER MINIMAL TREES, AND ANTIPODALITY HORST MARTINI, KONRAD J. SWANEPOEL, AND P. OLOFF DE WET arxiv:08.5046v [math.mg] 25 Aug 20 Abtract. We give a new proof that a tar {op i : i =,...,
More informationRELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS
www.arpapre.com/volume/vol29iue1/ijrras_29_1_01.pdf RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS Sevcan Demir Atalay 1,* & Özge Elmataş Gültekin
More informationLie series An old concept going back to Sophus Lie, but already used by Newton and made rigorous by Cauchy. Widely exploited, e.g.
ÄÁ Ë ÊÁ Ë Æ ÄÁ ÌÊ ÆË ÇÊÅË Lie erie An old concept going back to Sophu Lie, but already ued by Newton and made rigorou by Cauchy Widely exploited, eg, in differential geometry Ued a a method for numerical
More informationarxiv: v2 [math.nt] 30 Apr 2015
A THEOREM FOR DISTINCT ZEROS OF L-FUNCTIONS École Normale Supérieure arxiv:54.6556v [math.nt] 3 Apr 5 943 Cachan November 9, 7 Abtract In thi paper, we etablih a imple criterion for two L-function L and
More informationA SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho
A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES Sanghyun Cho Abtract. We prove a implified verion of the Nah-Moer implicit function theorem in weighted Banach pace. We relax the
More informationHELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES
15 TH INTERNATIONAL CONFERENCE ON GEOMETRY AND GRAPHICS 0 ISGG 1-5 AUGUST, 0, MONTREAL, CANADA HELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES Peter MAYRHOFER and Dominic WALTER The Univerity of Innbruck,
More informationLinear Motion, Speed & Velocity
Add Important Linear Motion, Speed & Velocity Page: 136 Linear Motion, Speed & Velocity NGSS Standard: N/A MA Curriculum Framework (006): 1.1, 1. AP Phyic 1 Learning Objective: 3.A.1.1, 3.A.1.3 Knowledge/Undertanding
More informationDimensional Analysis A Tool for Guiding Mathematical Calculations
Dimenional Analyi A Tool for Guiding Mathematical Calculation Dougla A. Kerr Iue 1 February 6, 2010 ABSTRACT AND INTRODUCTION In converting quantitie from one unit to another, we may know the applicable
More informationFinite Element Truss Problem
6. rue Uing FEA Finite Element ru Problem We tarted thi erie of lecture looking at tru problem. We limited the dicuion to tatically determinate tructure and olved for the force in element and reaction
More informationp. (The electron is a point particle with radius r = 0.)
- pin ½ Recall that in the H-atom olution, we howed that the fact that the wavefunction Ψ(r) i ingle-valued require that the angular momentum quantum nbr be integer: l = 0,,.. However, operator algebra
More informationPHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS B - HW #6 Spring 4, Solution by David Pace Any referenced equation are from Griffith Problem tatement are paraphraed. Problem. from Griffith Show that the following, A µo ɛ o A V + A ρ ɛ o Eq..4 A
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationSTOCHASTIC DIFFERENTIAL GAMES:THE LINEAR QUADRATIC ZERO SUM CASE
Sankhyā : he Indian Journal of Statitic 1995, Volume 57, Serie A, Pt. 1, pp.161 165 SOCHASIC DIFFERENIAL GAMES:HE LINEAR QUADRAIC ZERO SUM CASE By R. ARDANUY Univeridad de Salamanca SUMMARY. hi paper conider
More informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
More informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationNULL HELIX AND k-type NULL SLANT HELICES IN E 4 1
REVISTA DE LA UNIÓN MATEMÁTICA ARGENTINA Vol. 57, No. 1, 2016, Page 71 83 Publihed online: March 3, 2016 NULL HELIX AND k-type NULL SLANT HELICES IN E 4 1 JINHUA QIAN AND YOUNG HO KIM Abtract. We tudy
More informationQuark-Gluon Plasma in Proton-Proton Scattering at the LHC?
Quark-Gluon Plama in Proton-Proton Scattering at the LHC? K. Werner (a), Iu. Karpenko (b), T. Pierog (c) (a) SUBATECH, Univerity of Nante INP/CNRS EMN, Nante, France (b) Bogolyubov Intitute for Theoretical
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationTHE DIVERGENCE-FREE JACOBIAN CONJECTURE IN DIMENSION TWO
THE DIVERGENCE-FREE JACOBIAN CONJECTURE IN DIMENSION TWO J. W. NEUBERGER Abtract. A pecial cae, called the divergence-free cae, of the Jacobian Conjecture in dimenion two i proved. Thi note outline an
More informationAn Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem
An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent
More informationTuning of High-Power Antenna Resonances by Appropriately Reactive Sources
Senor and Simulation Note Note 50 Augut 005 Tuning of High-Power Antenna Reonance by Appropriately Reactive Source Carl E. Baum Univerity of New Mexico Department of Electrical and Computer Engineering
More informationChapter 7: The Laplace Transform Part 1
Chapter 7: The Laplace Tranform Part 1 王奕翔 Department of Electrical Engineering National Taiwan Univerity ihwang@ntu.edu.tw November 26, 213 1 / 34 王奕翔 DE Lecture 1 Solving an initial value problem aociated
More informationUSPAS Course on Recirculated and Energy Recovered Linear Accelerators
USPAS Coure on Recirculated and Energy Recovered Linear Accelerator G. A. Krafft and L. Merminga Jefferon Lab I. Bazarov Cornell Lecture 6 7 March 005 Lecture Outline. Invariant Ellipe Generated by a Unimodular
More informationLinear System Fundamentals
Linear Sytem Fundamental MEM 355 Performance Enhancement of Dynamical Sytem Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Content Sytem Repreentation Stability Concept
More informationON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS. Volker Ziegler Technische Universität Graz, Austria
GLASNIK MATEMATIČKI Vol. 1(61)(006), 9 30 ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS Volker Ziegler Techniche Univerität Graz, Autria Abtract. We conider the parameterized Thue
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More informationNotes on the geometry of curves, Math 210 John Wood
Baic definition Note on the geometry of curve, Math 0 John Wood Let f(t be a vector-valued function of a calar We indicate thi by writing f : R R 3 and think of f(t a the poition in pace of a particle
More information(3) A bilinear map B : S(R n ) S(R m ) B is continuous (for the product topology) if and only if there exist C, N and M such that
The material here can be found in Hörmander Volume 1, Chapter VII but he ha already done almot all of ditribution theory by thi point(!) Johi and Friedlander Chapter 8. Recall that S( ) i a complete metric
More informationFast Convolutional Sparse Coding (FCSC)
Fat Convolutional Spare Coding (FCSC) Bailey ong Department of Computer Science Univerity of California, Irvine bhkong@ic.uci.edu Charle C. Fowlke Department of Computer Science Univerity of California,
More informationLINEAR ALGEBRA METHOD IN COMBINATORICS. Theorem 1.1 (Oddtown theorem). In a town of n citizens, no more than n clubs can be formed under the rules
LINEAR ALGEBRA METHOD IN COMBINATORICS 1 Warming-up example Theorem 11 (Oddtown theorem) In a town of n citizen, no more tha club can be formed under the rule each club have an odd number of member each
More informationLaplace Transformation
Univerity of Technology Electromechanical Department Energy Branch Advance Mathematic Laplace Tranformation nd Cla Lecture 6 Page of 7 Laplace Tranformation Definition Suppoe that f(t) i a piecewie continuou
More informationWeber Schafheitlin-type integrals with exponent 1
Integral Tranform and Special Function Vol., No., February 9, 47 53 Weber Schafheitlin-type integral with exponent Johanne Kellendonk* and Serge Richard Univerité de Lyon, Univerité Lyon, Intitut Camille
More informationonline learning Unit Workbook 4 RLC Transients
online learning Pearon BTC Higher National in lectrical and lectronic ngineering (QCF) Unit 5: lectrical & lectronic Principle Unit Workbook 4 in a erie of 4 for thi unit Learning Outcome: RLC Tranient
More informationLecture 9: Shor s Algorithm
Quantum Computation (CMU 8-859BB, Fall 05) Lecture 9: Shor Algorithm October 7, 05 Lecturer: Ryan O Donnell Scribe: Sidhanth Mohanty Overview Let u recall the period finding problem that wa et up a a function
More information1 Routh Array: 15 points
EE C28 / ME34 Problem Set 3 Solution Fall 2 Routh Array: 5 point Conider the ytem below, with D() k(+), w(t), G() +2, and H y() 2 ++2 2(+). Find the cloed loop tranfer function Y () R(), and range of k
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationCOHOMOLOGY AS A LOCAL-TO-GLOBAL BRIDGE
COHOMOLOGY AS A LOCAL-TO-GLOBAL BRIDGE LIVIU I. NICOLAESCU ABSTRACT. I dicu low dimenional incarnation of cohomology and illutrate how baic cohomological principle lead to a proof of Sperner lemma. CONTENTS.
More informationHyperbolic Partial Differential Equations
Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationGiven the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is
EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -
More informationLecture 7: Testing Distributions
CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting
More informationThe Laplace Transform , Haynes Miller and Jeremy Orloff
The Laplace Tranform 8.3, Hayne Miller and Jeremy Orloff Laplace tranform baic: introduction An operator take a function a input and output another function. A tranform doe the ame thing with the added
More information1. Preliminaries. In [8] the following odd looking integral evaluation is obtained.
June, 5. Revied Augut 8th, 5. VA DER POL EXPASIOS OF L-SERIES David Borwein* and Jonathan Borwein Abtract. We provide concie erie repreentation for variou L-erie integral. Different technique are needed
More informationControl of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach
Proceeding of the 7th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING (ICOSSSE '8) Control of Delayed Integrating Procee Uing Two Feedback Controller R MS Approach LIBOR
More informationDesign By Emulation (Indirect Method)
Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal
More informationTHE THERMOELASTIC SQUARE
HE HERMOELASIC SQUARE A mnemonic for remembering thermodynamic identitie he tate of a material i the collection of variable uch a tre, train, temperature, entropy. A variable i a tate variable if it integral
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationSOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU. I will collect my solutions to some of the exercises in this book in this document.
SOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU CİHAN BAHRAN I will collect my olution to ome of the exercie in thi book in thi document. Section 2.1 1. Let A = k[[t ]] be the ring of
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationExercises for lectures 19 Polynomial methods
Exercie for lecture 19 Polynomial method Michael Šebek Automatic control 016 15-4-17 Diviion of polynomial with and without remainder Polynomial form a circle, but not a body. (Circle alo form integer,
More informationSECTION x2 x > 0, t > 0, (8.19a)
SECTION 8.5 433 8.5 Application of aplace Tranform to Partial Differential Equation In Section 8.2 and 8.3 we illutrated the effective ue of aplace tranform in olving ordinary differential equation. The
More informationLaplace Adomian Decomposition Method for Solving the Nonlinear Volterra Integral Equation with Weakly Kernels
Studie in Nonlinear Science (4): 9-4, ISSN -9 IDOSI Publication, Laplace Adomian Decompoition Method for Solving the Nonlinear Volterra Integral Equation with Weakly Kernel F.A. Hendi Department of Mathematic
More informationMulticolor Sunflowers
Multicolor Sunflower Dhruv Mubayi Lujia Wang October 19, 2017 Abtract A unflower i a collection of ditinct et uch that the interection of any two of them i the ame a the common interection C of all of
More informationDigital Control System
Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital
More informationSymmetry Lecture 9. 1 Gellmann-Nishijima relation
Symmetry Lecture 9 1 Gellmann-Nihijima relation In the lat lecture we found that the Gell-mann and Nihijima relation related Baryon number, charge, and the third component of iopin. Q = [(1/2)B + T 3 ]
More information