|
|
- Theodore Bridges
- 5 years ago
- Views:
Transcription
1 NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING Fernando Casas Universitat Jaume I, Castellón, Spain (on sabbatical leave at DAMTP, University of Cambridge) Edinburgh, June 2004
2 ...approach suggested by Arieh Iserles NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 2
3 Outline of the talk 1. Some (well known) Lie group methods for linear problems (Fer and Magnus expansions). 2. Schemes based on triangular matrices (splitting + solvability). 3. Some methods and practical issues in their construction NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 3
4 1 Lie group methods (linear problems) Let us consider a linear matrix ODE evolving in a Lie group G Y = A(t)Y, Y (t 0 ) = Y 0 G (0) with A : [t 0, [ G g smooth enough. g: Lie algebra associated with G Examples of G: SL(n), O(n), SU(n), Sp(n), SO(n),... Y (t) Lie group G if A(t) Lie algebra g There are several schemes preserving this feature (Magnus, Fer, Cayley,...) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 4
5 1.1 Magnus expansion For the equation Y = A(t)Y, Y (t 0 ) = I, Magnus (1954) proposed Y (t) = e Ω(t), Ω(t) = Ω k (t) (1) k=1 NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 5
6 1.1 Magnus expansion For the equation Y = A(t)Y, Y (t 0 ) = I, Magnus (1954) proposed Y (t) = e Ω(t), Ω(t) = Ω k (t) (1) k=1 with log(y (t)) satisfying Ω = d exp 1 Ω A(t) = k=0 B k k! adk Ω A(t), Ω(t 0) = 0, (2) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 6
7 1.1 Magnus expansion (II) Here and B k are Bernoulli numbers. ad 0 Ω A = A ad k Ω A = [Ω,adk 1 Ω A] [Ω,A] ΩA AΩ NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 7
8 1.1 Magnus expansion (III) First terms in the expansion (A i A(t i )): Ω 1 (t) = t Ω 2 (t) = 1 2 Ω 3 (t) = 1 6 t 0 A(t 1 )dt 1 t t 0 dt 1 t t1 t 0 dt 2 [A 1,A 2 ] t1 t2 dt 1 dt 2 dt 3 ([A 1,[A 2,A 3 ]] + [A 3,[A 2,A 1 ]]) t 0 t 0 t 0 e Ω(t) G even if the series Ω is truncated NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 8
9 1.1 Magnus expansion (III) First terms in the expansion (A i A(t i )): Ω 1 (t) = t Ω 2 (t) = 1 2 Ω 3 (t) = 1 6 t 0 A(t 1 )dt 1 t t 0 dt 1 t t1 t 0 dt 2 [A 1,A 2 ] t1 t2 dt 1 dt 2 dt 3 ([A 1,[A 2,A 3 ]] + [A 3,[A 2,A 1 ]]) t 0 t 0 t 0 e Ω(t) G even if the series Ω is truncated * Expansion widely used in Quantum Mechanics, NMR spectroscopy, infrared divergences in QED, control theory,... NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 9
10 1.1 Magnus expansion (IV) Magnus as a numerical integration method (Iserles & Nørsett, 1997) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 10
11 1.1 Magnus expansion (IV) Magnus as a numerical integration method (Iserles & Nørsett, 1997) Two critical factors in the computational cost of the resulting algorithms: NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 11
12 1.1 Magnus expansion (IV) Magnus as a numerical integration method (Iserles & Nørsett, 1997) Two critical factors in the computational cost of the resulting algorithms: (1) Evaluation of exp(ω) (Moler & Van Loan, Celledoni & Iserles,...) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 12
13 1.1 Magnus expansion (IV) Magnus as a numerical integration method (Iserles & Nørsett, 1997) Two critical factors in the computational cost of the resulting algorithms: (1) Evaluation of exp(ω) (Moler & Van Loan, Celledoni & Iserles,...) (2) Number of commutators involved in the expansion NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 13
14 1.1 Magnus expansion (IV) Magnus as a numerical integration method (Iserles & Nørsett, 1997) Two critical factors in the computational cost of the resulting algorithms: (1) Evaluation of exp(ω) (Moler & Van Loan, Celledoni & Iserles,...) (2) Number of commutators involved in the expansion To reduce this number is particularly useful the concept of graded free Lie algebra (Munthe-Kaas, Owren 1999) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 14
15 1.1 Magnus expansion (V) As a result, NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 15
16 1.1 Magnus expansion (V) As a result, * Numerical schemes based on Magnus up to order 8 have been constructed involving the minimum number of commutators in terms of quadratures and/or univariate integrals. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 16
17 1.1 Magnus expansion (V) As a result, * Numerical schemes based on Magnus up to order 8 have been constructed involving the minimum number of commutators in terms of quadratures and/or univariate integrals. * Efficient in applications NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 17
18 1.2 Other schemes Fer expansion. Built by Francis Fer (1958). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 18
19 1.2 Other schemes Fer expansion. Built by Francis Fer (1958). Referred erroneously in the (mathematical physics) literature (e.g., Wilcox 1967), but... NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 19
20 1.2 Other schemes Fer expansion. Built by Francis Fer (1958). Referred erroneously in the (mathematical physics) literature (e.g., Wilcox 1967), but... proposed (as an exercise!) by R. Bellman, Introduction to Matrix Analysis, 1960, page 204: NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 20
21 1.2 Other schemes Fer expansion. Built by Francis Fer (1958). Referred erroneously in the (mathematical physics) literature (e.g., Wilcox 1967), but... proposed (as an exercise!) by R. Bellman, Introduction to Matrix Analysis, 1960, page 204: The solution of dx/dt = Q(t)X, X(0) = I, can be put in the form e P e P 1 e Pn, where P = t 0 Q(s)ds, and P n = t 0 Q nds, with Q n = e P n 1 Q n 1 e P n The infinite product converges it t is sufficiently small. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 21 0 e sp n 1 Q n 1 e sp n 1 ds (See also Mathematical Reviews , review done by R. Bellman)
22 1.2 Other schemes (II) used as an (analytic) procedure in perturbation theory in Quantum Mechanics by Klarsfeld & Oteo (1989), but... NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 22
23 1.2 Other schemes (II) used as an (analytic) procedure in perturbation theory in Quantum Mechanics by Klarsfeld & Oteo (1989), but... numerical integration method built by Iserles (1984). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 23
24 1.2 Other schemes (II) used as an (analytic) procedure in perturbation theory in Quantum Mechanics by Klarsfeld & Oteo (1989), but... numerical integration method built by Iserles (1984). This class of methods can actually be built from Magnus. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 24
25 1.2 Other schemes (II) used as an (analytic) procedure in perturbation theory in Quantum Mechanics by Klarsfeld & Oteo (1989), but... numerical integration method built by Iserles (1984). This class of methods can actually be built from Magnus. They require the computation of several matrix exponentials. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 25
26 1.3 Methods based on the Cayley transform Let us suppose that Y = A(t)Y is defined in a J-orthogonal Lie group, J: constant matrix O J (n) = {A GL n (R) : A T JA = J}, NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 26
27 1.3 Methods based on the Cayley transform Let us suppose that Y = A(t)Y is defined in a J-orthogonal Lie group, O J (n) = {A GL n (R) : A T JA = J}, J: constant matrix Examples: orthogonal group (J = I), symplectic group, Lorentz group (J = diag(1, 1, 1, 1)). Solution: Y (t) = (I 12 C(t) ) 1 (I + 12 C(t) ) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 27
28 1.3 Methods based on the Cayley transform (II) with C(t) o J (n) satisfying (Iserles 2001) dc dt = A 1 2 [C,A] 1 4 CAC, t t 0, C(t 0 ) = 0. efficient methods without matrix exponentials! NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 28
29 1.3 Methods based on the Cayley transform (II) with C(t) o J (n) satisfying (Iserles 2001) dc dt = A 1 2 [C,A] 1 4 CAC, t t 0, C(t 0 ) = 0. efficient methods without matrix exponentials! In fact, one can also combine Magnus with Padé to avoid the use of matrix exponentials in J-orthogonal groups! NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 29
30 1.3 Methods based on the Cayley transform (II) with C(t) o J (n) satisfying (Iserles 2001) dc dt = A 1 2 [C,A] 1 4 CAC, t t 0, C(t 0 ) = 0. efficient methods without matrix exponentials! In fact, one can also combine Magnus with Padé to avoid the use of matrix exponentials in J-orthogonal groups! * It is possible to construct methods which are more efficient than those based on the Cayley transform (Blanes, C., Ros 2002). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 30
31 1.4 Summary These methods require the evaluation of one or several matrix exponentials NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 31
32 1.4 Summary These methods require the evaluation of one or several matrix exponentials They are expensive when n is very large NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 32
33 1.4 Summary These methods require the evaluation of one or several matrix exponentials They are expensive when n is very large In some cases, if the exponential is approximated by rational functions the method does not preserve the Lie group structure, in particular, when G = SL(n) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 33
34 1.4 Summary These methods require the evaluation of one or several matrix exponentials They are expensive when n is very large In some cases, if the exponential is approximated by rational functions the method does not preserve the Lie group structure, in particular, when G = SL(n) = Another class of methods is required. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 34
35 2 Solvability + splitting The procedure For the linear system Y = A(t)Y, Y (0) = I, we denote Y 0 Y, A 0 A and suppose that A 0 (t) = A 0+ (t) + A 0 (t), where A 0+ n is strictly upper-triangular A 0 n is weakly lower-triangular. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 35
36 2 Solvability + splitting (II) The idea is to write the solution as a product of upper and lower triangular matrices. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 36
37 2 Solvability + splitting (II) The idea is to write the solution as a product of upper and lower triangular matrices. More specifically, we propose the following factorization: Y 0 (t) = L 0 (t)u 0 (t)y 1 (t) such that L 0 = A 0 (t)l 0, L 0 (0) = I NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 37
38 2 Solvability + splitting (II) The idea is to write the solution as a product of upper and lower triangular matrices. More specifically, we propose the following factorization: Y 0 (t) = L 0 (t)u 0 (t)y 1 (t) such that L 0 = A 0 (t)l 0, L 0 (0) = I Observe then that L 0 (t) can be obtained by quadratures and L 0 (t) n. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 38
39 2 Solvability + splitting (III) Now we form the matrix C 0 = L 1 0 A 0 + L 0 NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 39
40 2 Solvability + splitting (III) Now we form the matrix C 0 = L 1 0 A 0 + L 0 which can also be split as where C 0 (t) = C 0+ (t) +C 0 (t), C 0+ n is weakly upper-triangular C 0 n is strictly lower-triangular. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 40
41 2 Solvability + splitting (IV) Next we choose U 0 as the solution of U 0 = C 0+ (t)u 0, U 0 (0) = I so that U 0 (t) can also be obtained by quadratures. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 41
42 2 Solvability + splitting (IV) Next we choose U 0 as the solution of U 0 = C 0+ (t)u 0, U 0 (0) = I so that U 0 (t) can also be obtained by quadratures. It is easy to show that Y 1 satisfies Y 1 = A 1 (t)y 1, Y 1 (0) = I, with A 1 = U 1 0 C 0 U 0. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 42
43 2 Solvability + splitting (V) This gives a single step of the solvable cycle, which we repeat with A 1. A 1 = A 1+ + A 1, A 1+ n, A 1 n Y 1 = L 1 U 1 Y 2 L 1 = A 1 L 1, L 1 (0) = I etc. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 43
44 2 Solvability + splitting (VI) In this way one has the following algorithm: with (k = 0,1,2,...) Y Y 0 = L 0 U 0 L 1 U 1 L k U k Y k+1 A k = A k+ + A k, A k+ n, A k n L k = A k L k, L k (0) = I C k Lk 1 A k+ L k = C k+ +C k C k+ n, C k n U k = C k+ U k, U k (0) = I NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 44
45 2 Solvability + splitting (VII) and finally A k+1 U 1 k C k U k, Y k+1 = A k+1 Y k+1 NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 45
46 2 Solvability + splitting (VII) and finally A k+1 U 1 k C k U k, Y k+1 = A k+1 Y k+1 Usually the factorization is truncated by taking Y k+1 = I. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 46
47 2 Solvability + splitting (VII) and finally A k+1 U 1 k C k U k, Y k+1 = A k+1 Y k+1 Usually the factorization is truncated by taking Y k+1 = I. In what follows we will analyse the main features of this procedure as a numerical integrator. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 47
48 2.1 Order of the method Suppose that A(t) = ε(a 0 + a 1 t + a 2 t 2 + ) for some parameter ε > 0. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 48
49 2.1 Order of the method Suppose that A(t) = ε(a 0 + a 1 t + a 2 t 2 + ) for some parameter ε > 0. Then A j = t n j ε n j ( εα 1 +t(εα 2 + ε 2 α 3 ) + O(t 2 ) ) A j+ = t m j ε m j ( εβ 1 +t(εβ 2 + ε 2 β 3 ) + O(t 2 ) ) for j = 1,2,..., so that L j (t) = I + 1 n j + 1 (tε)n j+1 α n j + 2 tn j+2 ε n j (εα 2 + ε 2 α 3 ) + U j (t) = I + 1 m j + 1 (tε)m j+1 β m j + 2 tm j+2 ε m j (εβ 2 + ε 2 β 3 ) + NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 49
50 2.1 Order of the method (II) Furthermore, n j+1 = n j + m j + 1 m j+1 = n j + 2m j + 2 j = 1,2,... NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 50
51 2.1 Order of the method (II) Furthermore, n j+1 = n j + m j + 1 m j+1 = n j + 2m j + 2 j = 1,2,... j n j m j NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 51
52 2.1 Order of the method (III) In consequence, NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 52
53 2.1 Order of the method (III) In consequence, (1) This algorithm could be useful for problems of the form Y = (B 0 + εb 1 )Y if the system Y = B 0 Y can be solved exactly. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 53
54 2.1 Order of the method (III) In consequence, (1) This algorithm could be useful for problems of the form Y = (B 0 + εb 1 )Y if the system Y = B 0 Y can be solved exactly. (2) The order of approximation is... NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 54
55 2.1 Order of the method (IV) Y 0 L 0 U 0 is order 1 Y 0 L 0 U 0 L 1 2 Y 0 L 0 U 0 L 1 U 1 4 Y 0 L 0 U 0 L 1 U 1 L 2 7 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 12 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 L 3 20 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 L 3 U 3 33 NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 55
56 2.1 Order of the method (IV) Y 0 L 0 U 0 is order 1 Y 0 L 0 U 0 L 1 2 Y 0 L 0 U 0 L 1 U 1 4 Y 0 L 0 U 0 L 1 U 1 L 2 7 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 12 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 L 3 20 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 L 3 U 3 33 With only 4 solvable cycles we get order 33! NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 56
57 2.1 Order of the method (IV) Y 0 L 0 U 0 is order 1 Y 0 L 0 U 0 L 1 2 Y 0 L 0 U 0 L 1 U 1 4 Y 0 L 0 U 0 L 1 U 1 L 2 7 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 12 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 L 3 20 Y 0 L 0 U 0 L 1 U 1 L 2 U 2 L 3 U 3 33 With only 4 solvable cycles we get order 33!...if we can compute L k and U k up to this order... NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 57
58 2.2 Questions Several problems involved NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 58
59 2.2 Questions Several problems involved (1) Does the approximate solution evolve in the Lie group if A is in the Lie algebra, i.e., is it a Lie group method? NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 59
60 2.2 Questions Several problems involved (1) Does the approximate solution evolve in the Lie group if A is in the Lie algebra, i.e., is it a Lie group method? (2) Solve explicitly the systems L k = A k L k and U k = C k + U k NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 60
61 2.2 Questions Several problems involved (1) Does the approximate solution evolve in the Lie group if A is in the Lie algebra, i.e., is it a Lie group method? (2) Solve explicitly the systems L k = A k L k and U k = C k + U k (3) Approximate efficiently the (multiple) integrals involved. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 61
62 3 Practical issues (1) Preservation of the Lie-group structure If A(t) sl(n), the algorithm provides by construction approximations to Y (t) in SL(n). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 62
63 3 Practical issues (1) Preservation of the Lie-group structure If A(t) sl(n), the algorithm provides by construction approximations to Y (t) in SL(n). Proof. A k = A k+ + A k, with A k+ n, A k n. In fact A k belongs to a solvable subalgebra of sl(n). Therefore the solution of L k = A k L k, L k (0) = I L k (t) SL(n) (in fact, a solvable subgroup of). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 63
64 3 Practical issues (1) Preservation of the Lie-group structure If A(t) sl(n), the algorithm provides by construction approximations to Y (t) in SL(n). Proof. A k = A k+ + A k, with A k+ n, A k n. In fact A k belongs to a solvable subalgebra of sl(n). Therefore the solution of L k = A k L k, L k (0) = I L k (t) SL(n) (in fact, a solvable subgroup of). Tr(A k+ ) = 0, and the trace is invariant under similarity, so that Tr(C k ) = Tr(Lk 1 A k+ L k ) = Tr(A k+ ) = 0 C k sl(n) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 64
65 3 Practical issues (II) Next, C k = C k+ +C k, with C k+ n, C k n and U k, solution of U k = C k + U k, U k (0) = I belongs to SL(n). Finally A k+1 U 1 k C k U k sl(n) and the process is repeated. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 65
66 3 Practical issues (II) Next, C k = C k+ +C k, with C k+ n, C k n and U k, solution of U k = C k + U k, U k (0) = I belongs to SL(n). Finally A k+1 U 1 k C k U k sl(n) and the process is repeated. Other properties (i.e., orthogonality) are preserved only up to the order of the method. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 66
67 3 Practical issues (III) (2a) Explicit solution of L k = A k L k NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 67
68 3 Practical issues (III) (2a) Explicit solution of L k = A k L k Consider k = 0 and denote A 0 (t) = (a i j ), i, j = 1,...,n, L 0 (t) = (L i j ), j i A ii (t) t 0 a ii (t 1 )dt 1. Then the solution of L 0 = A 0 (t)l 0, L 0 (0) = I is L ii (t) = e Aii(t), i = 1,...,n (3) ( ) t i 1 L i j (t) = e A ii(t) e A ii(t 1 ) a ik (t 1 )L k j (t 1 ) dt 1 0 k= j i = 2,...,n, j = 1,...,i 1. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 68
69 3 Practical issues (IV) (2b) Explicit solution of U k = C k + U k NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 69
70 3 Practical issues (IV) (2b) Explicit solution of U k = C k + U k Consider k = 0 and denote C 0 (t) = (c i j ), i, j = 1,...,n, U 0 (t) = (U i j ), j i C ii (t) t 0 c ii (t 1 )dt 1. Then the solution of U 0 = C 0 + (t)u 0, U 0 (0) = I is U ii (t) = e Cii(t), i = 1,...,n (4) ) U i j (t) = e C ii(t) t 0 ( j e C ii(t 1 ) c ik (t 1 )U k j (t 1 ) k=i+1 dt 1 i = 1,...,n 1, j = i + 1,...,n. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 70
71 3 Practical issues (V) Explicit expressions for the elements of L k and U k in terms of multivariate integrals. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 71
72 3 Practical issues (V) Explicit expressions for the elements of L k and U k in terms of multivariate integrals. They can be evaluated in sequence as follows: L ii i = 1,...,n U ii i = 1,...,n L i,i 1 i = 2,...,n U i,i+1 i = 1,...,n 1 L i,i 2 i = 3,...,n U i,i+2 i = 1,...,n 2.. L n1 U 1n NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 72
73 3 Practical issues (VI) In principle, the integrals appearing in L k and U k can be approximated by quadrature rules. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 73
74 3 Practical issues (VI) In principle, the integrals appearing in L k and U k can be approximated by quadrature rules. To minimise the computational cost this has to be done by using the minimum number of A evaluations in each integration step. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 74
75 3 Practical issues (VI) In principle, the integrals appearing in L k and U k can be approximated by quadrature rules. To minimise the computational cost this has to be done by using the minimum number of A evaluations in each integration step. Question: Is it possible to approximate all the nested integrals with the evaluations required to compute A ii = t 0 a ii (t 1 )dt 1, i.e., à la Magnus? NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 75
76 3 Practical issues (VI) In principle, the integrals appearing in L k and U k can be approximated by quadrature rules. To minimise the computational cost this has to be done by using the minimum number of A evaluations in each integration step. Question: Is it possible to approximate all the nested integrals with the evaluations required to compute A ii = t 0 a ii (t 1 )dt 1, i.e., à la Magnus? YES! NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 76
77 3.1 Example Illustration: method of order 4 with 2 A evaluations Step t = 0 t = h. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 77
78 3.1 Example Illustration: method of order 4 with 2 A evaluations Step t = 0 t = h. 1- Approximate A ii (h), i = 1,...,n up to order 4 A ii (h) = h 0 a ii (t)dt = h 3( aii (0) + 4a ii (h/2) + a ii (h) ) + O(h 5 ) Ã ii (h) + O(h 5 ) and A ii (h/2), i = 1,...,n 1, up to order 3 (necessary to approximate L i j ): A ii (h/2) = h 24 (5a ii(0) + 8a ii (h/2) a ii (h)) + O(h 4 ) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 78
79 3.1 Example (II) 2- L ii (h) = exp(ã ii (h)) + O(h 5 ) (i = 1,...,n) and L ii (h/2) = exp(ã ii (h/2)) + O(h 4 ) (i = 1,...,n 1). 3- Obtain an approximation to L i j (h), j < i, of order 4 and L i j (h/2) of order 3 with L i j (h) = e A ii(h) h 0 F i j (t)dt F i j (t) e A ii(t) i 1 a ik (t)l k j (t) k= j NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 79
80 3.1 Example (III) Then L i j (h) = h eãii(h) 3 (F i j(0) + 4F i j (h/2) + F i j (h)) + O(h 5 ) where F i j (0) = a i j (0) and F i j (h/2) and F i j (h) have to be approximated up to order h 3. The sequence of computation is (i = 2,...,n): (a) F i,i 1 (h/2) = e à ii (h/2) a i,i 1 (h/2)l i 1,i 1 (h/2) + O(h 4 ) (b) F i,i 1 (h) = e à ii (h/2) a i,i 1 (h)l i 1,i 1 (h) + O(h 5 ) (c) L i,i 1 (h), i = 2,...,n up to order 4 NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 80
81 3.1 Example (IV) (d) L i,i 1 (h/2) = eãii(h/2) h 24 (5a i,i 1(0) + 8F i,i 1 (h/2) F i,i 1 (h)) + O(h 4 ) (e) L i,i 2 (h), i = 3,...,n, up to order 4 and L i,i 2 (h/2) up to order 3...and so on. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 81
82 3.1 Example (IV) (d) L i,i 1 (h/2) = eãii(h/2) h 24 (5a i,i 1(0) + 8F i,i 1 (h/2) F i,i 1 (h)) + O(h 4 ) (e) L i,i 2 (h), i = 3,...,n, up to order 4 and L i,i 2 (h/2) up to order 3...and so on. In this way we have L 0 (h) computed up to order O(h 5 ) and also L 0 (h/2) up to order O(h 4 ) with 2 evaluations of A(t). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 82
83 3.1 Example (V) 3- Next we compute C 0 : C 0 (0) = A 0+ (0) error O(h 5 ) C 0 (h/2) = L 1 0 (h/2)a 0 + (h/2)l 0 (h/2) error O(h 4 ) C 0 (h) = L 1 0 (h)a 0 + (h)l 0 (h) error O(h 5 ) 4- C ii (h) = h 3( cii (0) + 4c ii (h/2) + c ii (h) ) + O(h 5 ) C ii (h/2) = h 24 (5c ii(0) + 8c ii (h/2) c ii (h)) + O(h 4 ) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 83
84 3.1 Example (VI) 5- U i,i+1 (h), i = 1,...,n 1, up to order O(h 5 ); U i,i+1 (h/2), i = 1,...,n 1, up to order O(h 4 ); U i,i+2 (h), i = 1,...,n 2, up to order O(h 5 ); U i,i+2 (h), i = 1,...,n 2, up to order O(h 4 );... and so on. Thus we compute U 0 (h) with error O(h 5 ) and also U 0 (h/2) with error O(h 4 ). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 84
85 3.1 Example (VII) 6- A 1 : A 1 (0) = C 0 (0) error O(h 5 ) A 1 (h/2) = U 1 0 (h/2)c 0 (h/2)u 0 (h/2) error O(h 4 ) A 1 (h) = U 1 0 (h)c 0 (h)u 0 (h) error O(h 5 )... and the process is repeated again for the second cycle NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 85
86 3.1 Example (VII) 6- A 1 : A 1 (0) = C 0 (0) error O(h 5 ) A 1 (h/2) = U 1 0 (h/2)c 0 (h/2)u 0 (h/2) error O(h 4 ) A 1 (h) = U 1 0 (h)c 0 (h)u 0 (h) error O(h 5 )... and the process is repeated again for the second cycle it is possible to construct a method of order 4 with only 2 A(t) evaluations (3 for the first step). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 86
87 3.2 Other possibilities One could use other quadrature rules instead, for instance Gauss Legendre, but... NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 87
88 3.2 Other possibilities One could use other quadrature rules instead, for instance Gauss Legendre, but... in that case there are not enough nodes to approximate all the multivariate integrals up to the required order. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 88
89 3.2 Other possibilities One could use other quadrature rules instead, for instance Gauss Legendre, but... in that case there are not enough nodes to approximate all the multivariate integrals up to the required order. Remark. This is not the case with Newton Cotes rules, although the error introduced may be important for higher order (negative coefficients) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 89
90 3.2 Other possibilities One could use other quadrature rules instead, for instance Gauss Legendre, but... in that case there are not enough nodes to approximate all the multivariate integrals up to the required order. Remark. This is not the case with Newton Cotes rules, although the error introduced may be important for higher order (negative coefficients) Solution: use G L with matrix evaluations in the previous/next step. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 90
91 3.2 Other possibilities One could use other quadrature rules instead, for instance Gauss Legendre, but... in that case there are not enough nodes to approximate all the multivariate integrals up to the required order. Remark. This is not the case with Newton Cotes rules, although the error introduced may be important for higher order (negative coefficients) Solution: use G L with matrix evaluations in the previous/next step. method of order 4 with 2 evaluations (and 1 in the next step) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 91
92 3.3 Some methods Order 4 Y L 0 U 0 L 1 U 1 Quadratures NC / GL, 2 matrix evaluations per step Order 6 Y L 0 U 0 L 1 U 1 L 2 order 6 with a 5 points NC quadrature (4 evaluations per step) order 7 with a 7 points NC (6 evaluations) Order 12 Y L 0 U 0 L 1 U 1 L 2 U 2 with a 11 points NC (or GL involving several steps). NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 92
93 3.4 Variable step size Local extrapolation technique is trivial to implement in this setting. NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 93
94 3.4 Variable step size Local extrapolation technique is trivial to implement in this setting. For instance, Y 1 L 0 U 0 L 1 Ŷ 1 L 0 U 0 L 1 U 1 = Y 1 U 1 Then Ŷ 1 Y 1 = Y 1 U 1 Y 1 = Y 1 (U 1 I) and Ŷ 1 Y 1 can be used as a measure of the error NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 94
95 3.5 Future work Analyse the convergence of the procedure NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 95
96 3.5 Future work Analyse the convergence of the procedure Consider numerical examples in SL(n) with (very) large n NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 96
97 3.5 Future work Analyse the convergence of the procedure Consider numerical examples in SL(n) with (very) large n Highly oscillatory problems (with special quadratures) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 97
98 3.5 Future work Analyse the convergence of the procedure Consider numerical examples in SL(n) with (very) large n Highly oscillatory problems (with special quadratures) Analyse in practice the preservation of other structures (Blanes & Moan) NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 98
99 3.5 Future work Analyse the convergence of the procedure Consider numerical examples in SL(n) with (very) large n Highly oscillatory problems (with special quadratures) Analyse in practice the preservation of other structures (Blanes & Moan) Try to generalize to nonlinear problems NEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING 99
100 The End Copyright c 2004 Fernando Casas Universitat Jaume I, Castellón, Spain Fernando.Casas@uji.es (on sabbatical leave at DAMTP, University of Cambridge) Made with the ujislides document class c Sergio Barrachina (barrachi@icc.uji.es)
IMPROVED HIGH ORDER INTEGRATORS BASED ON MAGNUS EPANSION S. Blanes ;, F. Casas, and J. Ros 3 Abstract We build high order ecient numerical integration
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Improved High Order Integrators Based on Magnus Expansion S. Blanes, F. Casas and J. Ros DAMTP 999/NA8 August, 999 Department of Applied Mathematics and
More informationCOMPOSITION MAGNUS INTEGRATORS. Sergio Blanes. Lie Group Methods and Control Theory Edinburgh 28 June 1 July Joint work with P.C.
COMPOSITION MAGNUS INTEGRATORS Sergio Blanes Departament de Matemàtiques, Universitat Jaume I Castellón) sblanes@mat.uji.es Lie Group Methods and Control Theory Edinburgh 28 June 1 July 24 Joint work with
More informationS. Sánchez, F. Casas & A. Fernández
New analytic approximations based on the Magnus expansion S. Sánchez, F. Casas & A. Fernández Journal of Mathematical Chemistry ISSN 259-9791 Volume 49 Number 8 J Math Chem (211) 49:1741-1758 DOI 1.17/s191-11-9855-
More informationTHE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS
THE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS MARIANNA KHANAMIRYAN Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road,
More informationHigh-order actions and their applications to honor our friend and
High-order actions and their applications to honor our friend and collaborator Siu A. Chin Universitat Politecnica de Catalunya, Barcelona, Spain Lecture: Magnus Expansion and Suzuki s Method Outline of
More informationComparative numerical solutions of stiff ordinary differential equations using Magnus series expansion method
NTMSCI, No. 1, 5-45 (215) 5 New Trends in Mathematical Sciences http://www.ntmsci.com Comparative numerical solutions of stiff ordinary differential equations using Magnus series expansion method Mehmet
More informationPhysics Reports 470 (2009) Contents lists available at ScienceDirect. Physics Reports. journal homepage:
Physics Reports 47 9 151 38 Contents lists available at ScienceDirect Physics Reports journal homepage: www.elsevier.com/locate/physrep The Magnus expansion and some of its applications S. Blanes a, F.
More informationExponential integrators
Exponential integrators Marlis Hochbruck Heinrich-Heine University Düsseldorf, Germany Alexander Ostermann University Innsbruck, Austria Helsinki, May 25 p.1 Outline Time dependent Schrödinger equations
More informationEfficient methods for time-dependence in semiclassical Schrödinger equations
Efficient methods for time-dependence in semiclassical Schrödinger equations Philipp Bader, Arieh Iserles, Karolina Kropielnicka & Pranav Singh November 25, 25 Abstract We build an efficient and unitary
More informationA pedagogical approach to Magnus expansion
A pedagogical approach to Magnus expansion S Blanes F Casas J A Oteo and J. Ros Instituto de Matemática Multidisciplinar, Universidad Politécnica de Valencia, E-46022 Valencia, Spain Departament de Matemàtiques,
More informationCOMPUTING THE MATRIX EXPONENTIAL WITH AN OPTIMIZED TAYLOR POLYNOMIAL APPROXIMATION
COMPUTING THE MATRIX EXPONENTIAL WITH AN OPTIMIZED TAYLOR POLYNOMIAL APPROXIMATION PHILIPP BADER, SERGIO BLANES, FERNANDO CASAS Abstract. We present a new way to compute the Taylor polynomial of the matrix
More informationLie algebraic quantum control mechanism analysis
Lie algebraic quantum control mechanism analysis June 4, 2 Existing methods for quantum control mechanism identification do not exploit analytic features of control systems to delineate mechanism (or,
More informationApplications of Lie group integrators and exponential schemes
Applications of Lie group integrators and exponential schemes Workshop on Lie group methods and control theory, Edinburgh, 28.6-1.7, 2004, Brynjulf Owren Dept of Math Sci, NTNU Applications of Lie group
More informationMethods in Computer Vision: Introduction to Matrix Lie Groups
Methods in Computer Vision: Introduction to Matrix Lie Groups Oren Freifeld Computer Science, Ben-Gurion University June 14, 2017 June 14, 2017 1 / 46 Definition and Basic Properties Definition (Matrix
More informationLie-group methods in geometric numerical integration
UNIVERSITY OF CAMBRIDGE CENTRE FOR MATHEMATICAL SCIENCES DEPARTMENT OF APPLIED MATHEMATICS & THEORETICAL PHYSICS Lie-group methods in geometric numerical integration Arieh Iserles 1 OUTLINE 1. Why geometric
More informationLECTURE 25-26: CARTAN S THEOREM OF MAXIMAL TORI. 1. Maximal Tori
LECTURE 25-26: CARTAN S THEOREM OF MAXIMAL TORI 1. Maximal Tori By a torus we mean a compact connected abelian Lie group, so a torus is a Lie group that is isomorphic to T n = R n /Z n. Definition 1.1.
More information1 Smooth manifolds and Lie groups
An undergraduate approach to Lie theory Slide 1 Andrew Baker, Glasgow Glasgow, 12/11/1999 1 Smooth manifolds and Lie groups A continuous g : V 1 V 2 with V k R m k open is called smooth if it is infinitely
More informationFavourable time integration methods for. Non-autonomous evolution equations
Favourable time integration methods for non-autonomous evolution equations International Conference on Scientific Computation and Differential Equations (SciCADE) Session on Numerical integration of non-autonomous
More informationA numerical approximation to some specific nonlinear differential equations using magnus series expansion method
NTMSCI 4, No. 1, 125-129 (216) 125 New Trends in Mathematical Sciences http://dx.doi.org/1.2852/ntmsci.21611566 A numerical approximation to some specific nonlinear differential equations using magnus
More informationThree approaches for the design of adaptive time-splitting methods
Three approaches for the design of adaptive time-splitting methods Mechthild Thalhammer Leopold Franzens Universität Innsbruck, Austria Workshop on Geometric Integration and Computational Mechanics Organisers:
More informationSymplectic propagators for the Kepler problem with time-dependent mass
Celestial Mechanics and Dynamical Astronomy manuscript No. (will be inserted by the editor) Symplectic propagators for the Kepler problem with time-dependent mass Philipp Bader Sergio Blanes Fernando Casas
More informationUnitary transformations depending on a small parameter
Unitary transformations depending on a small parameter Fernando Casas Fernando.Casas@mat.uji.es www.gicas.uji.es Departament de Matemàtiques Institut de Matemàtiques i Aplicacions de Castelló (IMAC) Universitat
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 informationSplitting methods (with processing) for near-integrable problems
Splitting methods (with processing) for near-integrable problems Fernando Casas casas@mat.uji.es www.gicas.uji.es Departament de Matemàtiques Institut de Matemàtiques i Aplicacions de Castelló (IMAC) Universitat
More informationEnergy-Preserving Runge-Kutta methods
Energy-Preserving Runge-Kutta methods Fasma Diele, Brigida Pace Istituto per le Applicazioni del Calcolo M. Picone, CNR, Via Amendola 122, 70126 Bari, Italy f.diele@ba.iac.cnr.it b.pace@ba.iac.cnr.it SDS2010,
More informationMagnus series expansion method for solving nonhomogeneous stiff systems of ordinary differential equations
Kuwait J. Sci. 43 (1) pp. 25-38, 2016 Magnus series expansion method for solving nonhomogeneous stiff systems of ordinary differential equations Mehmet T. Atay 1, Aytekin Eryılmaz 2, Sure Köme 2,* 1 Department
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS (Semester I: AY 2014-15) Time Allowed: 2 Hours INSTRUCTIONS TO CANDIDATES 1. Please write your student number only. 2. This examination
More informationSymmetric Spaces Toolkit
Symmetric Spaces Toolkit SFB/TR12 Langeoog, Nov. 1st 7th 2007 H. Sebert, S. Mandt Contents 1 Lie Groups and Lie Algebras 2 1.1 Matrix Lie Groups........................ 2 1.2 Lie Group Homomorphisms...................
More informationBasic Concepts of Group Theory
Chapter 1 Basic Concepts of Group Theory The theory of groups and vector spaces has many important applications in a number of branches of modern theoretical physics. These include the formal theory of
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationSymplectic time-average propagators for the Schödinger equation with a time-dependent Hamiltonian
Symplectic time-average propagators for the Schödinger equation Symplectic time-average propagators for the Schödinger equation with a time-dependent Hamiltonian Sergio Blanes,, a Fernando Casas, 2, b
More informationTHEORY OF GROUP REPRESENTATIONS AND APPLICATIONS
THEORY OF GROUP REPRESENTATIONS AND APPLICATIONS ASIM 0. BARUT Institute for Theoretical Physics, University of Colorado, Boulder, Colo., U.S.A. RYSZARD RATJZKA Institute for Nuclear Research, Warszawa,
More informationLECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups
LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS 1. Lie groups A Lie group is a special smooth manifold on which there is a group structure, and moreover, the two structures are compatible. Lie groups are
More informationSplittings, commutators and the linear Schrödinger equation
UNIVERSITY OF CAMBRIDGE CENTRE FOR MATHEMATICAL SCIENCES DEPARTMENT OF APPLIED MATHEMATICS & THEORETICAL PHYSICS Splittings, commutators and the linear Schrödinger equation Philipp Bader, Arieh Iserles,
More informationMATRIX LIE GROUPS. Claudiu C Remsing. Dept of Mathematics (Pure and Applied) Rhodes University Grahamstown September Maths Seminar 2007
MATRIX LIE GROUPS Claudiu C Remsing Dept of Mathematics (Pure and Applied) Rhodes University Grahamstown 6140 26 September 2007 Rhodes Univ CCR 0 TALK OUTLINE 1. What is a matrix Lie group? 2. Matrices
More informationAttempts at relativistic QM
Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and
More informationAccurate approximation of stochastic differential equations
Accurate approximation of stochastic differential equations Simon J.A. Malham and Anke Wiese (Heriot Watt University, Edinburgh) Birmingham: 6th February 29 Stochastic differential equations dy t = V (y
More informationQuantum Field Theory III
Quantum Field Theory III Prof. Erick Weinberg January 19, 2011 1 Lecture 1 1.1 Structure We will start with a bit of group theory, and we will talk about spontaneous symmetry broken. Then we will talk
More informationThese notes are incomplete they will be updated regularly.
These notes are incomplete they will be updated regularly. LIE GROUPS, LIE ALGEBRAS, AND REPRESENTATIONS SPRING SEMESTER 2008 RICHARD A. WENTWORTH Contents 1. Lie groups and Lie algebras 2 1.1. Definition
More informationExplicit Wei-Norman formulæ for matrix Lie groups via Putzer s method
Explicit Wei-Norman formulæ for matrix Lie groups via Putzer s method Claudio Altafini SISSA-ISAS International School for Advanced Studies via Beirut 2-4, 34014 Trieste, Italy; Abstract The Wei-Norman
More informationLecture 2: Computing functions of dense matrices
Lecture 2: Computing functions of dense matrices Paola Boito and Federico Poloni Università di Pisa Pisa - Hokkaido - Roma2 Summer School Pisa, August 27 - September 8, 2018 Introduction In this lecture
More informationQUATERNIONS AND ROTATIONS
QUATERNIONS AND ROTATIONS SVANTE JANSON 1. Introduction The purpose of this note is to show some well-known relations between quaternions and the Lie groups SO(3) and SO(4) (rotations in R 3 and R 4 )
More informationBlack Hole thermodynamics
Black Hole thermodynamics I Black holes evaporates I Black holes have a partition function For a Schwarzschild black hole, the famous Bekenstein-Hawking results are: T = 1 8 M S = A 4G = 4 r + 4G Note
More informationResearch Article The Approximate Solution of Fredholm Integral Equations with Oscillatory Trigonometric Kernels
Applied Mathematics, Article ID 72327, 7 pages http://ddoiorg/055/204/72327 Research Article The Approimate Solution of Fredholm Integral Equations with Oscillatory Trigonometric Kernels Qinghua Wu School
More informationLinear ODE s with periodic coefficients
Linear ODE s with periodic coefficients 1 Examples y = sin(t)y, solutions Ce cos t. Periodic, go to 0 as t +. y = 2 sin 2 (t)y, solutions Ce t sin(2t)/2. Not periodic, go to to 0 as t +. y = (1 + sin(t))y,
More informationCompact schemes for laser-matter interaction in Schrödinger equation
Compact schemes for laser-matter interaction in Schrödinger equation Arieh Iserles, Karolina Kropielnicka & Pranav Singh April 21, 218 AMS Mathematics Subject Classification: Primary 65M7, Secondary 35Q41,
More informationInterpolation in Special Orthogonal Groups
Interpolation in Special Orthogonal Groups Tatiana Shingel DAMTP, Centre for Mathematical Sciences University of Cambridge August 7, 7 Abstract The problem of constructing smooth interpolating curves in
More informationMulti-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems
Multi-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems Jürgen Geiser and Siu A. Chin geiser@mathematik.hu-berlin.de Abstract. In this paper we discuss the
More informationSince a matrix is invertible if and only if its determinant is nonzero, we can also write. GL(n) = {A R n n : deta 0}. (2)
MATHEMATICS 7302 (Analytical Dynamics) YEAR 2016 2017, TERM 2 HANDOUT #13: INTRODUCTION TO THE ROTATION GROUP 1 Lie groups As mentioned in Handout #11, the general framework for studying symmetry in mathematics
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationNumerical Methods I: Numerical Integration/Quadrature
1/20 Numerical Methods I: Numerical Integration/Quadrature Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 30, 2017 umerical integration 2/20 We want to approximate the definite integral
More informationMatrix Lie groups. and their Lie algebras. Mahmood Alaghmandan. A project in fulfillment of the requirement for the Lie algebra course
Matrix Lie groups and their Lie algebras Mahmood Alaghmandan A project in fulfillment of the requirement for the Lie algebra course Department of Mathematics and Statistics University of Saskatchewan March
More informationAutomorphic Equivalence Within Gapped Phases
1 Harvard University May 18, 2011 Automorphic Equivalence Within Gapped Phases Robert Sims University of Arizona based on joint work with Sven Bachmann, Spyridon Michalakis, and Bruno Nachtergaele 2 Outline:
More informationSplitting and composition methods for the time dependent Schrödinger equation
Splitting and composition methods for the time dependent Schrödinger equation S. Blanes, joint work with F. Casas and A. Murua Instituto de Matemática Multidisciplinar Universidad Politécnica de Valencia,
More information(Kac Moody) Chevalley groups and Lie algebras with built in structure constants Lecture 1. Lisa Carbone Rutgers University
(Kac Moody) Chevalley groups and Lie algebras with built in structure constants Lecture 1 Lisa Carbone Rutgers University Slides will be posted at: http://sites.math.rutgers.edu/ carbonel/ Video will be
More informationSplitting Methods for Non-autonomous Hamiltonian Equations
Journal of Computational Physics 170, 205 230 2001 doi:10.1006/jcph.2001.6733, available online at http://www.idealibrary.com on Splitting Methods for Non-autonomous Hamiltonian Equations S. Blanes and
More informationLinearized methods for ordinary di erential equations
Applied Mathematics and Computation 104 (1999) 109±19 www.elsevier.nl/locate/amc Linearized methods for ordinary di erential equations J.I. Ramos 1 Departamento de Lenguajes y Ciencias de la Computacion,
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationSchrödinger equation
Splitting methods for the time dependent Schrödinger equation S. Blanes, Joint work with F. Casas and A. Murua Instituto de Matemática Multidisciplinar Universidad Politécnica de Valencia, SPAIN Workshop:
More informationAn homomorphism to a Lie algebra of matrices is called a represetation. A representation is faithful if it is an isomorphism.
Lecture 3 1. LIE ALGEBRAS 1.1. A Lie algebra is a vector space along with a map [.,.] : L L L such that, [αa + βb,c] = α[a,c] + β[b,c] bi linear [a,b] = [b,a] Anti symmetry [[a,b],c] + [[b,c],a][[c,a],b]
More informationThe Toda Lattice. Chris Elliott. April 9 th, 2014
The Toda Lattice Chris Elliott April 9 th, 2014 In this talk I ll introduce classical integrable systems, and explain how they can arise from the data of solutions to the classical Yang-Baxter equation.
More informationTopics in Representation Theory: Lie Groups, Lie Algebras and the Exponential Map
Topics in Representation Theory: Lie Groups, Lie Algebras and the Exponential Map Most of the groups we will be considering this semester will be matrix groups, i.e. subgroups of G = Aut(V ), the group
More informationSANGRADO PAGINA 17CMX24CM. PhD Thesis. Splitting methods for autonomous and non-autonomous perturbed equations LOMO A AJUSTAR (AHORA 4CM)
SANGRADO A3 PAGINA 17CMX24CM LOMO A AJUSTAR (AHORA 4CM) We have considered the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients.
More informationAngular matrix integrals
Montreal, 25 August 2008 1 Angular matrix integrals Montreal, 25 August 2008 J.-B. uber A. Prats Ferrer, B. Eynard, P. Di Francesco, J.-B... Correlation Functions of Harish-Chandra Integrals over the Orthogonal
More information1: Lie groups Matix groups, Lie algebras
Lie Groups and Bundles 2014/15 BGSMath 1: Lie groups Matix groups, Lie algebras 1. Prove that O(n) is Lie group and that its tangent space at I O(n) is isomorphic to the space so(n) of skew-symmetric matrices
More informationCourse 216: Ordinary Differential Equations
Course 16: Ordinary Differential Equations Notes by Chris Blair These notes cover the ODEs course given in 7-8 by Dr. John Stalker. Contents I Solving Linear ODEs 1 Reduction of Order Computing Matrix
More information2.1 Definition and general properties
Chapter 2 Gaussian states Gaussian states are at the heart of quantum information processing with continuous variables. The basic reason is that the vacuum state of quantum electrodynamics is itself a
More informationParallel Methods for ODEs
Parallel Methods for ODEs Levels of parallelism There are a number of levels of parallelism that are possible within a program to numerically solve ODEs. An obvious place to start is with manual code restructuring
More informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
More informationAnton ARNOLD. with N. Ben Abdallah (Toulouse), J. Geier (Vienna), C. Negulescu (Marseille) TU Vienna Institute for Analysis and Scientific Computing
TECHNISCHE UNIVERSITÄT WIEN Asymptotically correct finite difference schemes for highly oscillatory ODEs Anton ARNOLD with N. Ben Abdallah (Toulouse, J. Geier (Vienna, C. Negulescu (Marseille TU Vienna
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationValidated Explicit and Implicit Runge-Kutta Methods
Validated Explicit and Implicit Runge-Kutta Methods Alexandre Chapoutot joint work with Julien Alexandre dit Sandretto and Olivier Mullier U2IS, ENSTA ParisTech 8th Small Workshop on Interval Methods,
More informationThe Spinor Representation
The Spinor Representation Math G4344, Spring 2012 As we have seen, the groups Spin(n) have a representation on R n given by identifying v R n as an element of the Clifford algebra C(n) and having g Spin(n)
More informationGROUP THEORY PRIMER. New terms: so(2n), so(2n+1), symplectic algebra sp(2n)
GROUP THEORY PRIMER New terms: so(2n), so(2n+1), symplectic algebra sp(2n) 1. Some examples of semi-simple Lie algebras In the previous chapter, we developed the idea of understanding semi-simple Lie algebras
More informationHigh-order time-splitting methods for Schrödinger equations
High-order time-splitting methods for Schrödinger equations and M. Thalhammer Department of Mathematics, University of Innsbruck Conference on Scientific Computing 2009 Geneva, Switzerland Nonlinear Schrödinger
More informationMANIFOLD STRUCTURES IN ALGEBRA
MANIFOLD STRUCTURES IN ALGEBRA MATTHEW GARCIA 1. Definitions Our aim is to describe the manifold structure on classical linear groups and from there deduce a number of results. Before we begin we make
More informationCommutator-free Magnus Lanczos methods for the linear Schrödinger equation
Commutator-free Magnus Lanczos methods for the linear Schrödinger equation Arieh Iserles, Karolina Kropielnicka & Pranav Singh October 3, 27 AMS Mathematics Subject Classification: Primary 65M7, Secondary
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationThe Hamiltonian operator and states
The Hamiltonian operator and states March 30, 06 Calculation of the Hamiltonian operator This is our first typical quantum field theory calculation. They re a bit to keep track of, but not really that
More informationImplicit-explicit exponential integrators
Implicit-explicit exponential integrators Bawfeh Kingsley Kometa joint work with Elena Celledoni MaGIC 2011 Finse, March 1-4 1 Introduction Motivation 2 semi-lagrangian Runge-Kutta exponential integrators
More informationReal, Complex, and Quarternionic Representations
Real, Complex, and Quarternionic Representations JWR 10 March 2006 1 Group Representations 1. Throughout R denotes the real numbers, C denotes the complex numbers, H denotes the quaternions, and G denotes
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationLyapunov-based control of quantum systems
Lyapunov-based control of quantum systems Symeon Grivopoulos Bassam Bamieh Department of Mechanical and Environmental Engineering University of California, Santa Barbara, CA 936-57 symeon,bamieh@engineering.ucsb.edu
More informationAccuracy, Precision and Efficiency in Sparse Grids
John, Information Technology Department, Virginia Tech.... http://people.sc.fsu.edu/ jburkardt/presentations/ sandia 2009.pdf... Computer Science Research Institute, Sandia National Laboratory, 23 July
More informationOn the Harish-Chandra Embedding
On the Harish-Chandra Embedding The purpose of this note is to link the Cartan and the root decompositions. In addition, it explains how we can view a Hermitian symmetric domain as G C /P where P is a
More informationCHERN-SIMONS THEORY AND WEYL QUANTIZATION Răzvan Gelca
CHERN-SIMONS THEORY AND WEYL QUANTIZATION Răzvan Gelca CHERN-SIMONS THEORY AND WEYL QUANTIZATION Răzvan Gelca this talk is based on joint work and discussions with Alejandro Uribe, Alastair Hamilton, Charles
More informationMASLOV INDEX. Contents. Drawing from [1]. 1. Outline
MASLOV INDEX J-HOLOMORPHIC CURVE SEMINAR MARCH 2, 2015 MORGAN WEILER 1. Outline 1 Motivation 1 Definition of the Maslov Index 2 Questions Demanding Answers 2 2. The Path of an Orbit 3 Obtaining A 3 Independence
More informationFermionic coherent states in infinite dimensions
Fermionic coherent states in infinite dimensions Robert Oeckl Centro de Ciencias Matemáticas Universidad Nacional Autónoma de México Morelia, Mexico Coherent States and their Applications CIRM, Marseille,
More informationFourier series on triangles and tetrahedra
Fourier series on triangles and tetrahedra Daan Huybrechs K.U.Leuven Cambridge, Sep 13, 2010 Outline Introduction The unit interval Triangles, tetrahedra and beyond Outline Introduction The topic of this
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationExercises on chapter 1
Exercises on chapter 1 1. Let G be a group and H and K be subgroups. Let HK = {hk h H, k K}. (i) Prove that HK is a subgroup of G if and only if HK = KH. (ii) If either H or K is a normal subgroup of G
More informationRepresentations of Matrix Lie Algebras
Representations of Matrix Lie Algebras Alex Turzillo REU Apprentice Program, University of Chicago aturzillo@uchicago.edu August 00 Abstract Building upon the concepts of the matrix Lie group and the matrix
More informationLIE ALGEBRAS MICHAEL PENKAVA
LIE ALGEBRAS MICHAEL PENKAVA 1. Introduction In this text, we will be studying Lie algebra structures on a vector space V over a field K. Most of the time, the field will be C, and the vector space will
More informationarxiv: v1 [math.na] 12 Oct 2018
arxiv:181.5412v1 [math.na] 12 Oct 218 Sixth-order schemes for laser matter interaction in Schrödinger equation Pranav Singh October 15, 218 AMS Mathematics Subject Classification: Primary 65M7, Secondary
More informationNumerical methods for matrix functions
Numerical methods for matrix functions SF2524 - Matrix Computations for Large-scale Systems Lecture 13 Numerical methods for matrix functions 1 / 26 Reading material Lecture notes online Numerical methods
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationTheoretical Analysis of Nonlinear Differential Equations
Stephen F. Austin State University SFA ScholarWorks Electronic Theses and Dissertations 2-2018 Theoretical Analysis of Nonlinear Differential Equations Emily Jean Weymier Stephen F Austin State University,
More informationarxiv:math-ph/ v1 31 May 2000
An Elementary Introduction to Groups and Representations arxiv:math-ph/0005032v1 31 May 2000 Author address: Brian C. Hall University of Notre Dame, Department of Mathematics, Notre Dame IN 46556 USA E-mail
More informationNotes on nilpotent orbits Computational Theory of Real Reductive Groups Workshop. Eric Sommers
Notes on nilpotent orbits Computational Theory of Real Reductive Groups Workshop Eric Sommers 17 July 2009 2 Contents 1 Background 5 1.1 Linear algebra......................................... 5 1.1.1
More information