Parallel Multipatch Space-Time IGA Solvers for Parabloic Initial-Boundary Value Problems
|
|
- Malcolm Byrd
- 5 years ago
- Views:
Transcription
1 Parallel Multipatch Space-Time IGA Solvers for Parabloic Initial-Boundary Value Problems Ulrich Langer Institute of Computational Mathematics (NuMa) Johannes Kepler University Linz (RICAM) Austrian Academy of Sciences (ÖAW) Linz, Austria AANMPDE12, Särkisaari, Finland, August 6 10, 2018
2 Joint work with my collaborators Christoph Hofer (JKU, DK) Martin Neumüller (JKU, NuMa) Ioannis Toulopoulos (RICAM, CM4PDE) Main results have just been published in [1] C. Hofer, U. Langer, M. Neumüller, I. Toulopoulos. Time-multipatch discontinuous Galerkin space-time isogeometric analysis of parabolic evolution problems. ETNA, 2018, v. 49, [2] C. Hofer, U. Langer, M. Neumüller. Robust preconditioning for space-time isogeometric analysis of parabolic evolution problems. [math.na], 2018, arxiv: , arxiv.org. Doctoral Program Computational Mathematics Numerical Analysis and Symbolic Computation
3 Outline 1 Introduction 2 Time Multi-Patch Space-Time IgA 3 Space-Time Solvers 4 Conclusions & Outlooks
4 Parabolic Initial-Boundary Value Model Problem Let us consider the IBVP problem: Find u : Q R such that t u u = f in Q := Ω (0, T ), u = u D := 0 on Σ := Ω (0, T ), u = u 0 on Σ 0 := Ω {0}, (1) as the typical model problem for a linear parabolic evolution equation posed in the space-time cylinder Q = Ω [0, T ]. Our space-time technology can be generalized to more general parabolic equations like div x (A(x, t) x u) + b(x, t) x u + c(x, t) t u + a(x, t)u = f, eddy-current problems, non-linear problems etc.
5 Parabolic Initial-Boundary Value Model Problem Let us consider the IBVP problem: Find u : Q R such that t u u = f in Q := Ω (0, T ), u = u D := 0 on Σ := Ω (0, T ), u = u 0 on Σ 0 := Ω {0}, (1) as the typical model problem for a linear parabolic evolution equation posed in the space-time cylinder Q = Ω [0, T ]. Our space-time technology can be generalized to more general parabolic equations like div x (A(x, t) x u) + b(x, t) x u + c(x, t) t u + a(x, t)u = f, eddy-current problems, non-linear problems etc.
6 Standard Weak Space-Time Variational Formulation Find u H 1,0 0 (Q) such that a(u, v) = l(v) v H 1,1 (Q), (2) 0,0 with the bilinear form a(u, v) = u(x, t) t v(x, t)dxdt + Q Q x u(x, t) x v(x, t)dxdt and the linear form l(v) = f (x, t)v(x, t)dxdt + u 0 (x)v(x, 0)dx, Q Ω For solvability and regularity results, we refer to original papers by the Leningrad School of Mathematical Physics in the 50s, summerized in the monographs by Ladyzhenskaya, Solonnikov & Uralceva (1967) and Ladyzhenskaya (1973)!
7 Parabolic Solvability and Regularity Results If f L 2,1 (Q T ) := {v : T 0 v(, t) L 2 (Ω) dt < } and u 0 L 2 (Ω), then there exists a unique generalized (weak) solution u H 1,0 1,0 0 (Q) of (2) that even belongs to V2,0 (Q T ). If f L 2 (Q T ) and u 0 H0 1 (Ω), then the generalized solution u of (2) belongs to space H,1 0 (Q T ) = W,1 2,0 (Q T ), and u continuously depends on t in the norm of the space H0 1(Ω), where W,1 2,0 (Q T ) = {v H0 1(Q T ) : x v L 2 (Q T )}. Maximal parabolic regularity: t u div x (A(x, t) x u) = f t u X + div x (A(x, t) x u) X C f X, where X = L p ((0, T ), L q (Ω)) = L p,q (Q T ), 1 < p, q <, u 0 = 0, and a = 0, b = 0, c = 1.
8 Parabolic Solvability and Regularity Results If f L 2,1 (Q T ) := {v : T 0 v(, t) L 2 (Ω) dt < } and u 0 L 2 (Ω), then there exists a unique generalized (weak) solution u H 1,0 1,0 0 (Q) of (2) that even belongs to V2,0 (Q T ). If f L 2 (Q T ) and u 0 H0 1 (Ω), then the generalized solution u of (2) belongs to space H,1 0 (Q T ) = W,1 2,0 (Q T ), and u continuously depends on t in the norm of the space H0 1(Ω), where W,1 2,0 (Q T ) = {v H0 1(Q T ) : x v L 2 (Q T )}. Maximal parabolic regularity: t u div x (A(x, t) x u) = f t u X + div x (A(x, t) x u) X C f X, where X = L p ((0, T ), L q (Ω)) = L p,q (Q T ), 1 < p, q <, u 0 = 0, and a = 0, b = 0, c = 1.
9 Some References to Time-parallel methods: Gander (2015): Nice historical overview on 50 years time-parallel method Parareal introduced by Lions, Maday, Turinici (2001) Time-parallel multigrid: Hackbusch (1984),...,Vandewalle (1993),..., Gander & Neumüller (2014): smart time-parallel multigrid,..., Neumüller & Smears (2018),... Space-time methods for parabolic evolution problems: Babuška & Janik (1989,1990), Behr (2008), Schwab & Stevenson (2009), Neumüller & Steinbach (2011), Neumüller (2013), Andreev (2013), Bank & Metti (2013), Mollet (2014), Urban & Patera (2014), Schwab & Stevenson (2014), Bank & Vassilevski (2014), Karabelas & Neumüller (2015), Langer & Moore & Neumüller (2016), Bank & Vassilevski & Zikatanov (2016), Larsson & Molteni (2017), Steinbach & Yang (2017), Langer & Neumüller & Schafelner (2018),...
10 Some References to Time-parallel methods: Gander (2015): Nice historical overview on 50 years time-parallel method Parareal introduced by Lions, Maday, Turinici (2001) Time-parallel multigrid: Hackbusch (1984),...,Vandewalle (1993),..., Gander & Neumüller (2014): smart time-parallel multigrid,..., Neumüller & Smears (2018),... Space-time methods for parabolic evolution problems: Babuška & Janik (1989,1990), Behr (2008), Schwab & Stevenson (2009), Neumüller & Steinbach (2011), Neumüller (2013), Andreev (2013), Bank & Metti (2013), Mollet (2014), Urban & Patera (2014), Schwab & Stevenson (2014), Bank & Vassilevski (2014), Karabelas & Neumüller (2015), Langer & Moore & Neumüller (2016), Bank & Vassilevski & Zikatanov (2016), Larsson & Molteni (2017), Steinbach & Yang (2017), Langer & Neumüller & Schafelner (2018),...
11 References Time-parallel methods: Gander (2015): Nice historical overview on 50 years time-parallel methods Space-time methods for parabolic evolution problems: Steinbach & Yang (2018): Nice overview on space-time methods in the space-time book that will be published in RSCAM by de Gruyter
12 Single-Patch Space-Time IgA: LMN 16 [3] U. Langer, S. Moore, M. Neumüller. Space-time isogeometric analysis of parabolic evolution equations. Comput. Methods Appl. Mech. Engrg., v. 306, pp , t T 0 0 a Q K Ω b x Φ 1 Φ ˆt 1 K Q 0 Ω 1 ˆx t T Q ( ) n x n t Φ 1 1 ˆx2 x2 Ω Φ 0 0 Ω 0 x1 0 1 ˆt ˆQ ˆx1 Space-Time IgA paraphernalia: Q R d+1 ; d = 1 (l) and d = 2 (r). In this talk: Generalization to the multipatch case + Solvers!
13 Outline 1 Introduction 2 Time Multi-Patch Space-Time IgA 3 Space-Time Solvers 4 Conclusions & Outlooks
14 Time Multi-Patch Decomposition of Q We decompose the space-time cylinder Q = Ω (0, T ) into N non-overlapping space-time subcylinder b n = 1, 2,..., N, such that Qn = Ω (tn 1, tn ) = Φn (Q), S Q= N n=1 Q n with the time faces Σn = Q n+1 Q n = Ω {tn }, ΣN = ΣT : T Introduction Time Multi-Patch Space-Time IgA Space-Time Solvers Conclusions & Outlooks
15 Time Multipatch dg IgA spaces V 0h We look for an approximate solution u h to the IBVP (2) in the globally discontinuous, but patch-wise smooth IgA (B-spline, NURBS) spaces V 0h = {v h H 1,0 0 (Q) : v h n := v h Qn B Ξ d+1(q n ), n = 1,..., N} = {v h L 2 (Q) : v n h V n 0h, n = 1,..., N} = span{ϕ i} i I, V0h n = {v h n B Ξ d+1 (Q n n ) : vh n = 0 on Σ} = span{ϕ n,i} i In where B Ξ d+1(q n n ) is the smooth (depending of the polynomial degrees and multiciplicity of the knots) IgA space corresponding to the knot vector Ξ d+1 n = Ξ d+1 n (n n 1,.., n n d+1 ; pn 1,.., p n d+1 ) =... n
16 Stable Time Multipatch dg IgA Scheme Multiplying the PDE t u x u = f by v h + θ n h n t v h, integrating over Q n, integrating by parts, suming over n, and using that the jumps [ u ] across Σ n are 0 at the solution u H,1 0 (Q), we get the consistency identity where a h (u, v h ) = l h (v h ) v h V 0h, N a h (u, v h ) = ( tuv h + θ nh n tu tv h + x u x v h + θ nh n x u x tv h ) dxdt n=1 Q n N + u n 1 v n 1 h,+ dx, n=1 Σ n 1 N l h (v h ) = f [v h + θ nh n tv h ] dxdt + u 0 vh,+ 0 dx. n=1 Q n Σ 0
17 Stable Time Multipatch dg IgA Scheme Multiplying the PDE t u x u = f by v h + θ n h n t v h, integrating over Q n, integrating by parts, suming over n, and using that the jumps [ u ] across Σ n are 0 at the solution u H,1 0 (Q), we get the multi-patch space-time scheme: Find u h V 0h such that where a h (u h, v h ) = l h (v h ) v h V 0h, N a h (u h, v h ) = ( tu h v h + θ nh n tu h tv h + x u h x v h + θ nh n x u h x tv h ) dxdt n=1 Q n N + u n 1 h v n 1 h,+ dx, n=1 Σ n 1 N l h (v h ) = f [v h + θ nh n tv h ] dxdt + u 0 vh,+ 0 dx. n=1 Q n Σ 0
18 Space-Time IgA Scheme and System of IgA Eqns Hence, we look for the solution u h V 0h of the IgA scheme in the form of a h (u h, v h ) = l h (v h ) v h V 0h (3) u h (x, t) = u h (x 1,..., x d, x d+1 ) = i I u iϕ i (x, t) where u h := [u i ] i I R N h= I is the unknown solution vector of control points defined by the solution of the linear system L h u h = f h (4) with huge, non-symmetric, but positive definite system matrix L h.
19 Road Map of the Numerical Analysis v h 2 h = N n=1 ( x v h 2 L 2 (Q + n) θn hn tv h 2 L 2 (Q + 1 ) n) 2 v h n 1 2 L 2 (Σ n 1 ) N v 2 h, = v 2 h + (θ nh n) 1 v 2 L 2 (Q + N n) n=1 n=2 v n 1 2 L 2 (Σ n 1 ) v h 2 L 2 (Σ N ) Coercivity: a h (v h, v h ) µ c v h 2 h, v h V 0h,! u h u h : (4) Boundedness: a h (u, v h ) µ b u h, v h h, u V 0h + H,1 0 (Q), v h V 0h, Consistency: a h (u, v h ) = l h (v h ) v h V 0h, u H 1,0 0 (Q) H,1 (Q): (2) GO: a h (u u h, v h ) = 0 v h V 0h, Cea-like: u u h h (1 + µ b /µ c) inf vh V 0h u v h h,... Convergence rates: u u h h ch p u H p+1 (Q)
20 V 0h Coercivity of the bilinear form a h (, ) We now introduce the mesh-dependent norm v h 2 h = N n=1 ( x v h 2 L 2 (Q + n) θn hn tv h 2 L 2 (Q + 1 ) n) 2 v h n 1 2 L 2 (Σ n 1 ) Lemma (Coercivity / Ellipticity on V 0h ) v h 2 L 2 (Σ N ). The bilinear form a h (, ) : V 0h V 0h R is V 0h coercive wrt the norm h, i.e., there exists a constant µ c = 1/2 such that a h (v h, v h ) µ c v h 2 h, v h V 0h. (5) provided that θ n c 2 inv,0, where v h 2 L 2 ( E) c inv,0hn 1 v h 2 L 2 (E) This lemma immediately yields uniqueness and existence of the solution u h V 0h and u h R N h of (9) and (4), respectively.
21 Uniform Boundedness of a h (, ) on V 0h, V 0h Let us introduce the space V 0h, = V 0h + H,1 0 (Q) equipped with the norm v h, = ( v N 2 h + (θ nh n) 1 v 2 L 2 (Q + N n) v n L 2 (Σ n 1 )). (6) n=1 n=2 Lemma (Boundedness) The bilinear form a h (, ) is uniformly bounded on V 0h, V 0h : a h (u, v h ) µ b u h, v h h, u V 0h,, v h V 0h, (7) with µ b = max(c inv,1 θ max, 2), where θ max = max n {θ n } c 2 inv,0. and c inv,k = c inv,k (p) are constants in the inverse inequalities t xi v h 2 L 2 (E) c inv,1h 2 n x i v h 2 L 2 (E) and v h 2 L 2 ( E) c inv,0h 1 n v h 2 L 2 (E).
22 Consistency and Galerkin orthogonality Lemma (Consitency) If the solution u H 1,0 0 (Q) of the variational problem (2) belongs to H,1 (Q), then it satisfies the consistency identity a h (u, v h ) = l h (v h ) v h V 0h. (8) Lemma (Galerkin orthogonality) If the solution u H 1,0 0 (Q) of the variational problem (2) belongs to H,1 (Q), then the Galerkin orthogonality a h (u u h, v h ) = 0 v h V 0h. (9) holds, where u h V 0h is the space-time dg IgA solution.
23 Cea-like Discretization Error Estimate Theorem (Cea-like Estimate) Let the exact solution u of (2) belong to H 1,0 0 (Q) H,1 (Q), and let u h be the solution of the space-time IgA scheme (9). Then we get ( where v h, = ( N v h = n=1 u u h h (1 + µ b µ c ) inf v h V 0h u v h h,, (10) v 2 h + N n=1 (θnhn) 1 v 2 L 2 (Q + N n) n=2 ( 1 2 x v 2 L 2 (Qn) + θn hn tv 2 L 2 (Qn) v 2 L 2 (Σn 1)) v 2 L 2 (Σ N )). Proof: u u h h u v h h + v h u h h µ c v h u h 2 h a h(v h u h, v h u h ) = a h (v h u, v h u h ) µ b u v h h, v h u h h v n 1 2 L 2 (Σ n 1 )) 1 2 and
24 Approximation Error Estimate Theorem (Approximation Theorem) Let p n + 1 l n 2 and p n + 1 m n 1 be integers, and let u L 2 (Q) such that the restriction u n := u Qn belongs to H ln,mn (Q n) for n = 1,..., N. Then there exists a quasi-interpolant Π h u V 0h such that ( N u Π h u 2 h, = ( x (u Π n h u) 2 L 2 (Q + n) θn hn t(u Πn h u) 2 L 2 (Q n) n=1 + 1 ) 2 (u Π hu) n 1 2 L 2 (Σ n 1 ) + 1 ) 2 u ΠN h u 2 L 2 (Σ N ) N 1 + u Π n h θ n=1 nh u 2 L 2 (Q + N n) n n=2 N ( (C n hn 2(ln 1) + θ nhn 2ln 1 + hn 2mn 1 n=1 N ( ( C n hn 2(ln 1) + h 2(mn 2 1 ) ) n u 2 Hln,mn (Qn) n=1 (u Π n 1 h u) n 1 2 L 2 (Σ n 1 ) ) + θ nhn 2mn 1 u 2 Hln,mn (Qn)
25 A priori Discretization Error Estimate Theorem (A priori Disscretization Error Estimate) u u h h (1 + µ b µ c ) N n=1 ( ( ) C n hn 2(ln 1) + h 2(mn 1 2 ) n u 2 H ln,mn (Q n) We remark that for the case of highly smooth solutions, i.e., p + 1 min(l n, m n), the above estimate takes the form u u h h C N hn p u H p+1,p+1 (Q n) n=1 C h p u H p+1,p+1 (Q) where the last estimate holds if u H p+1,p+1 (Q) and h = max{h n} is assumed.
26 (sp cg, mp dg) IgA for d = 1 and p = 2, 3, error in dg-norm dg meshsize h p=2 p=3 p=4 O(h 2 ) O(h 3 ) O(h 4 )
27 (sp cg, mp dg) IgA for d = 1 and p = 2, 3, 4 Error in the dg-norm and convergence rate for the exact solution u(x, t) = sin(πx) sin( π 2 (t + 1)) and for B-Spline degrees 2, 3 and 4 p = 2 p = 3 p = 4 refinement error eoc error eoc error eoc E E E E E E E E E E E E E E E E E E E E E E E E
28 (mp cg, mp dg) IgA for d = 1 and p = 2, 3 ref Introduction p=2 error e e-05 Time Multi-Patch Space-Time IgA eoc p=3 error e e e e-08 Space-Time Solvers eoc Conclusions & Outlooks
29 Outline 1 Introduction 2 Time Multi-Patch Space-Time IgA 3 Space-Time Solvers 4 Conclusions & Outlooks
30 One huge system of IgA equations Once the basis is chosen, the IgA scheme (9) can be rewritten as a huge system of algebraic equations of the form L h u h = f h (11) for determining the vector u h = ((u 1,i ) i I1,..., (u N,i ) i IN ) R N h of the control points of the IgA solution u h (x, t) = i I n u n,i ϕ n,i (x, t), (x, t) Q n, n = 1,..., N, solving the IgA scheme (9). The system matrix L h is the usual Galerkin (stiffness) matrix, and f h is the rhs (load) vector. The matrix L h is non-symmetric, but positive definite!
31 System Matrix L h The Galerkin matrix L h can be rewritten in the block form A 1 B 2 A 2 L h = B 3 A , B N A N with the matrices A n := M n,x K n,t + K n,x M n,t for n = 1,..., N, B n := M n,x N n,t for n = 2,..., N.
32 Parallel Space-Time Multigrid Solvers Solve A 1 B 2 A 2 L h u h = f h, with L h = B N A N by the time-parallel MGM proposed by Gander & Neumüller ( 16): Ingredients: Time-Restriction and Prolongation Smoother
33 Time-Restriction Two time-slabs are restricted to a single time-slab. Assumption: Same basis on each slab
34 Smoother Inexact damped block Jacobi smoother: u h (k+1) = u h (k) + ω t D 1 h [ f h L h u h (k) ] for k = 1, 2,... with ω t = 1 2 and ω t = 1 2 D h := diag{a n } n=1,...,n Parallel w.r.t. time Replace D 1 1 h by multigrid D h w.r.t space parallel in space
35 Parallel Solver Studies for d = 3 and p = 1 Figure: Computational spatial domain Ω decomposed into 4096 elements (left) and distributed over 32 processors (right). The IgA solution u h (x, t) u(x, t) = sin(πx 1 ) sin(πx 2 ) sin(πx 3 ) sin(πt) is plotted at t = 0.5.
36 Parallel Solver Studies for d = 3 and p = 1 Convergence results in the h - norm for the regular solution u(x, t) = sin(πx 1 ) sin(πx 2 ) sin(πx 3 ) sin(πt) as well as iteration numbers and solving times for the parallel space-time multigrid preconditioned GMRES method on Vulcan N overall dof u u h h eoc c x c t cores iter time [s] E E E E E E E
37 Parallel Solver Studies for d = 3 and p = 1 Convergence results in the norm h for the low regularity solution u(x, t) = cos(βx 1 ) cos(βx 2 ) cos(βx 3 )(1 t) α H s,α+ 1 2 ε (Q), with α = 0.75 and β = 0.3, for an arbitrary s 2 and for an arbitrary small ε > 0, as well as iteration numbers and solving times for the parallel space-time multigrid preconditioned GMRES method on Vulcan BlueGene/Q at LLNL N overall dof u u h h eoc c x c t cores iter time [s] E E E E E E E
38 New Space-Time Multigrid Solver / Preconditioner Utilizing the tensor product structure A = A n = K t M x + M t K x, to construct a cheap approximation  1 to A 1, i.e., D 1 to D 1 : K t K T t, M t M T t K t and M t are small ( 1d - matrices) compared to K x and M x ( 1d, 2d, 3d - SPD-matrices). Idea: Perform decomposition of M 1 t K t 1 Fast Diagonalization: M 1 t K t = XDX 1, X, D C nt nt 2 Complex-Schur: M 1 t K t = QTQ, Q, T C nt nt 3 Real-Schur: M 1 t K t = QTQ T, Q, T R nt nt 4 Ref.: Sangalli & Tani (2016), Tani (2017) for elliptic BVP
39 New Space-Time Multigrid Solver / Preconditioner Utilizing the tensor product structure A = A n = K t M x + M t K x, to construct a cheap approximation  1 to A 1, i.e., D 1 to D 1 : K t K T t, M t M T t K t and M t are small ( 1d - matrices) compared to K x and M x ( 1d, 2d, 3d - SPD-matrices). Idea: Perform decomposition of M 1 t K t 1 Fast Diagonalization: M 1 t K t = XDX 1, X, D C nt nt 2 Complex-Schur: M 1 t K t = QTQ, Q, T C nt nt 3 Real-Schur: M 1 t K t = QTQ T, Q, T R nt nt 4 Ref.: Sangalli & Tani (2016), Tani (2017) for elliptic BVP
40 Fast Diagonalization Eigenvalue decomposition M 1 t K t = XDX 1, i.e. K t X = M t XD D = diag(λ i ), λ i C (Eigenvalues) X C nt nt (Eigenvectors), X 1 X! Defining Y := (M t X) 1 gives M t = Y 1 X 1 K t = Y 1 DX 1 Hence, we obtain A 1 = (K t M x + M t K x ) 1 = ( (Y 1 I) (D M x + I K x ) (X 1 I) ) 1 = (X I) (D M x + I K x ) 1 (Y I) }{{} diag((k x +λ i M x ) 1 )
41 Fast Diagonalization Eigenvalue decomposition M 1 t K t = XDX 1, i.e. K t X = M t XD D = diag(λ i ), λ i C (Eigenvalues) X C nt nt (Eigenvectors), X 1 X! Defining Y := (M t X) 1 gives M t = Y 1 X 1 K t = Y 1 DX 1 Hence, we obtain A 1 = (K t M x + M t K x ) 1 = ( (Y 1 I) (D M x + I K x ) (X 1 I) ) 1 = (X I) (D M x + I K x ) 1 (Y I) }{{} diag((k x +λ i M x ) 1 )
42 Solver for K x + λ i M x Case 1: λ i R + : K x + λ i M x... SPD matrix solvers available, e.g., DD, MG Case 2: λ i = α i + β i i C, α i > 0: (K x + λ i M x ) H K x + λ i M x Rewrite as symmetric, indefinite problem: ( ) Kx + α A i := i M x βm x βm x (K x + α i M x ) Robust preconditioner 1 cond 2 (P 1 i A i ) 2: ( ) Kx + (α P i := i + β i )M x 0 0 K x + (α i + β i )M x K x + (α i + β i )M x is a SPD matrix DD, MG 1 W. Zulehner, SIAM J. Matrix Anal. & Appl., 32(2), , 2011
43 Solver for K x + λ i M x Case 1: λ i R + : K x + λ i M x... SPD matrix solvers available, e.g., DD, MG Case 2: λ i = α i + β i i C, α i > 0: (K x + λ i M x ) H K x + λ i M x Rewrite as symmetric, indefinite problem: ( ) Kx + α A i := i M x βm x βm x (K x + α i M x ) Robust preconditioner 1 cond 2 (P 1 i A i ) 2: ( ) Kx + (α P i := i + β i )M x 0 0 K x + (α i + β i )M x K x + (α i + β i )M x is a SPD matrix DD, MG 1 W. Zulehner, SIAM J. Matrix Anal. & Appl., 32(2), , 2011
44 Solver for K x + λ i M x Case 1: λ i R + : K x + λ i M x... SPD matrix solvers available, e.g., DD, MG Case 2: λ i = α i + β i i C, α i > 0: (K x + λ i M x ) H K x + λ i M x Rewrite as symmetric, indefinite problem: ( ) Kx + α A i := i M x βm x βm x (K x + α i M x ) Robust preconditioner 1 cond 2 (P 1 i A i ) 2: ( ) Kx + (α P i := i + β i )M x 0 0 K x + (α i + β i )M x K x + (α i + β i )M x is a SPD matrix DD, MG 1 W. Zulehner, SIAM J. Matrix Anal. & Appl., 32(2), , 2011
45 Numerical Tests Parallel application for each i = 1,..., n t Unfortunately, cond 2 (X) >> 1 (Eigenvectors) Reliable approximation  only for small n t. n t p p = 2 p = 3 p = 4 p = 5 p = 6 p = 7 p = e e+06 6e+06 3e+07 1e e+06 1e+07 6e+07 4e+08 7e+09 1e e+07 4e+07 3e+08 3e+09 6e+10 3e+11 1e e+08 1e+09 1e+10 3e+11 2e+13 5e+13 4e+14 Table: condition number of X : θ = 0.01 and t i+1 t i = 0.1.
46 Complex Schur decomposition The complex Schur decomposition gives M 1 t K t = QTQ T = uppertriag(t ij ) C, T ii = λ i (Eigenvalues) Q C nt nt with Q 1 = Q cond 2 (Q) = 1. Again, by defining Y := (M t Q ) 1, we obtain M t = Y 1 Q K t = Y 1 TQ. Hence, we have A 1 = (M x K t + K x M t ) 1 = (Q I) (T M x + I K x ) 1 (Y I)
47 Complex Schur decomposition The complex Schur decomposition gives M 1 t K t = QTQ T = uppertriag(t ij ) C, T ii = λ i (Eigenvalues) Q C nt nt with Q 1 = Q cond 2 (Q) = 1. Again, by defining Y := (M t Q ) 1, we obtain M t = Y 1 Q K t = Y 1 TQ. Hence, we have A 1 = (M x K t + K x M t ) 1 = (Q I) (T M x + I K x ) 1 (Y I)
48 Complex Schur decomposition (cont.) T M x + I K x has the following structure K x + λ 1 M x T 12 M x... 0 K x + λ 2 M x T 23 M x Tnt 1n t M x K x + λ nt M x Application is done staggered K x + λ i M x has the same structure as in the diagonal case. Real Schur decomposition works similar only real arithmetic.
49 YETI-footprint: Space mp cg & time mp dg IgA
50 Numerical Tests: Spatial Solver: K x + λ i M x λ i C MinRes method with tolerance ε = 10 8 (K x + (α i + β i )M x ) 1 Direct solver Maximum number of iterations for i = 1,..., n t. ref. x and t p = 2 p = 3 p = 4 p = 5 p =
51 Numerical Tests: Sparse Direct Solver Degree p x, p t = (3, 3); θ = 0.01 and t i t i 1 = 0.1 A i MinRes with ε = 10 4 (It. 10). (K x + (α i + β i )M x ) 1 Direct solver Direct Solver: PARDISO Tolerance Multigrid ε = 10 8 #dofs ref #slabs MG-It Direct Diag x t Setup Solving Setup Solving #dofs ref #slabs MG-It C-Schur R-Schur
52 Numerical Tests: Parallelization in time - Real Schur Degree p x, p t = (3, 3), refinement r x, r t = (2, 3) θ = 0.01 and t i t i 1 = 0.1; (K x + (α i + β i )M x ) 1 Direct solver Direct Solver: PARDISO Tolerance ε = 10 8 #slabs = #Processors. Weak Scaling #dofs #slabs MG-It. Setup Solving 4.5 Setup 4 Solving Procs Time (s)
53 Parallelization in Space & Time Complex Schur Spatial domain YETI-Footprint (84 patches) Degree p x, p t = (3, 3), θ = 0.01 and t [0, 4] Parallel MG as preconditioner/solver for spatial problems #dofs refs #sl. c total c x c t it Setup Solving
54 Parallelization in Space & Time Real Schur Spatial domain YETI-Footprint (84 patches) Degree p x, p t = (3, 3), θ = 0.01 and t [0, 4] Parallel MG as preconditioner/solver for spatial problems #dofs refs #sl. c total c x c t it Setup Solving
55 Parallelization in Space & Time Real Schur Spatial domain YETI-Footprint (84 patches) Degree p x, p t = (3, 3), θ = 0.01 and t [0, 4] Parallel MG as preconditioner/solver for spatial problems #dofs refs #sl. c total c x c t it Setup Solving serial c t =
56 Numerical Tests: Parallelization in Space & Time Alternative Version non-symmetric IETI-DP for A n Spatial domain Ω: 4 4 square grid Degree p x, p t = (3, 3), θ = 0.01 and t i t i 1 = 0.1. GMRes-IETI-DP in the smoother (3 Iterations). Coarsening only in time! #dofs ref x #sl. c total c x c t It. Setup Solving
57 Outline 1 Introduction 2 Time Multi-Patch Space-Time IgA 3 Space-Time Solvers 4 Conclusions & Outlooks
58 Conclusions & Ongoing Work & Outlook Space-time IgA: space singlepatch cg and time-multipatch dg Space-time IgA: space multipatch cg and time-multipatch dg Space-time IgA: Fast generation via tensor techniques Space-time IgA: Fast parallel solvers Functional a posteriori estimates and THB-Spline adaptivity = joint work with S. Matculevich and S. Repin = talk by Svetlana Matculevich on Thursday! space-time multipatch dg IgA + adaptivity + fast generation + efficient parallel solvers??? Adaptive Space-Time FEM = talk by Andreas Schafelner on Tuesday!
59 Conclusions & Ongoing Work & Outlook Space-time IgA: space singlepatch cg and time-multipatch dg Space-time IgA: space multipatch cg and time-multipatch dg Space-time IgA: Fast generation via tensor techniques Space-time IgA: Fast parallel solvers Functional a posteriori estimates and THB-Spline adaptivity = joint work with S. Matculevich and S. Repin = talk by Svetlana Matculevich on Thursday! space-time multipatch dg IgA + adaptivity + fast generation + efficient parallel solvers??? Adaptive Space-Time FEM = talk by Andreas Schafelner on Tuesday!
60 Conclusions & Ongoing Work & Outlook Space-time IgA: space singlepatch cg and time-multipatch dg Space-time IgA: space multipatch cg and time-multipatch dg Space-time IgA: Fast generation via tensor techniques Space-time IgA: Fast parallel solvers Functional a posteriori estimates and THB-Spline adaptivity = joint work with S. Matculevich and S. Repin = talk by Svetlana Matculevich on Thursday! space-time multipatch dg IgA + adaptivity + fast generation + efficient parallel solvers??? Adaptive Space-Time FEM = talk by Andreas Schafelner on Tuesday!
61 Numerical Results: Parallel Solution for p = 1 2d fixed spatial domain Ω = (0, 1) 2 yielding Q = Ω (0, T ) = (0, 1) 3 Exact solution: u(x, t) = sin(πx 1 ) sin(πx 2 ) sin(πt) Figure shows space-time decomposition with 64 subdomains Table shows parallel performance of the parallel AMG hypre preconditioned GMRES stopped after the relative residual error reduction by 10 10! Dofs u u h L2 (Q) Rate iter time [s] cores e e e e e e e e e e e This example was computed on the supercomputer Vulcan BlueGene/Q in Livermore by M. Neumüller.
62 Supercomputing Results on Vulcan: (3+1)D, p = 1 3d fixed spatial domain Ω = (0, 1) 3 yielding Q = Ω (0, T ) = (0, 1) 4 Exact solution: u(x, t) = sin(πx 1 ) sin(πx 2 ) sin(πx 3 ) sin(πt) Table shows parallel performance of the parallel AMG hypre preconditioned GMRES stopped after the relative residual error reduction by 10 6! dofs u u h L2 (Q) rate iter time [s] cores e e e e e e e e e This example was computed on the supercomputer Vulcan BlueGene/Q in Livermore by M. Neumüller.
63 THANK YOU VERY MUCH! Solution 3.32e Introduction Time Multi-Patch Space-Time IgA Space-Time Solvers Conclusions & Outlooks
Multipatch Space-Time Isogeometric Analysis of Parabolic Diffusion Problems. Ulrich Langer, Martin Neumüller, Ioannis Toulopoulos. G+S Report No.
Multipatch Space-Time Isogeometric Analysis of Parabolic Diffusion Problems Ulrich Langer, Martin Neumüller, Ioannis Toulopoulos G+S Report No. 56 May 017 Multipatch Space-Time Isogeometric Analysis of
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationSpace-time isogeometric analysis of parabolic evolution equations
www.oeaw.ac.at Space-time isogeometric analysis of parabolic evolution equations U. Langer, S.E. Moore, M. Neumüller RICAM-Report 2015-19 www.ricam.oeaw.ac.at SPACE-TIME ISOGEOMETRIC ANALYSIS OF PARABOLIC
More informationJOHANNES KEPLER UNIVERSITY LINZ. Institute of Computational Mathematics
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Multipatch Space-Time Isogeometric Analysis of Parabolic Diffusion Problems Ulrich Langer Institute of Computational Mathematics,
More informationDual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations
Dual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations Christoph Hofer and Ulrich Langer Doctoral Program Computational Mathematics Numerical
More informationSpace time finite element methods in biomedical applications
Space time finite element methods in biomedical applications Olaf Steinbach Institut für Angewandte Mathematik, TU Graz http://www.applied.math.tugraz.at SFB Mathematical Optimization and Applications
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationIsogEometric Tearing and Interconnecting
IsogEometric Tearing and Interconnecting Christoph Hofer and Ludwig Mitter Johannes Kepler University, Linz 26.01.2017 Doctoral Program Computational Mathematics Numerical Analysis and Symbolic Computation
More informationSpace time finite and boundary element methods
Space time finite and boundary element methods Olaf Steinbach Institut für Numerische Mathematik, TU Graz http://www.numerik.math.tu-graz.ac.at based on joint work with M. Neumüller, H. Yang, M. Fleischhacker,
More informationInexact Data-Sparse BETI Methods by Ulrich Langer. (joint talk with G. Of, O. Steinbach and W. Zulehner)
Inexact Data-Sparse BETI Methods by Ulrich Langer (joint talk with G. Of, O. Steinbach and W. Zulehner) Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences http://www.ricam.oeaw.ac.at
More informationParallel Discontinuous Galerkin Method
Parallel Discontinuous Galerkin Method Yin Ki, NG The Chinese University of Hong Kong Aug 5, 2015 Mentors: Dr. Ohannes Karakashian, Dr. Kwai Wong Overview Project Goal Implement parallelization on Discontinuous
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationElectronic Transactions on Numerical Analysis Volume 49, 2018
Electronic Transactions on Numerical Analysis Volume 49, 2018 Contents 1 Adaptive FETI-DP and BDDC methods with a generalized transformation of basis for heterogeneous problems. Axel Klawonn, Martin Kühn,
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
More informationMesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain Stephen Edward Moore Johann Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences,
More informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
More informationA new approach for Kirchhoff-Love plates and shells
A new approach for Kirchhoff-Love plates and shells Walter Zulehner Institute of Computational Mathematics JKU Linz, Austria AANMPDE 10 October 2-6, 2017, Paleochora, Crete, Greece Walter Zulehner (JKU
More informationFunctional type error control for stabilised space-time IgA approximations to parabolic problems
www.oeaw.ac.at Functional type error control for stabilised space-time IgA approximations to parabolic problems U. Langer, S. Matculevich, S. Repin RICAM-Report 2017-14 www.ricam.oeaw.ac.at Functional
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationNew Multigrid Solver Advances in TOPS
New Multigrid Solver Advances in TOPS R D Falgout 1, J Brannick 2, M Brezina 2, T Manteuffel 2 and S McCormick 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O.
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationJOHANNES KEPLER UNIVERSITY LINZ. A Robust FEM-BEM Solver for Time-Harmonic Eddy Current Problems
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics A Robust FEM-BEM Solver for Time-Harmonic Eddy Current Problems Michael Kolmbauer Institute of Computational Mathematics, Johannes
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationAuxiliary space multigrid method for elliptic problems with highly varying coefficients
Auxiliary space multigrid method for elliptic problems with highly varying coefficients Johannes Kraus 1 and Maria Lymbery 2 1 Introduction The robust preconditioning of linear systems of algebraic equations
More informationCondition number estimates for matrices arising in the isogeometric discretizations
www.oeaw.ac.at Condition number estimates for matrices arising in the isogeometric discretizations K. Gahalaut, S. Tomar RICAM-Report -3 www.ricam.oeaw.ac.at Condition number estimates for matrices arising
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationA Chebyshev-based two-stage iterative method as an alternative to the direct solution of linear systems
A Chebyshev-based two-stage iterative method as an alternative to the direct solution of linear systems Mario Arioli m.arioli@rl.ac.uk CCLRC-Rutherford Appleton Laboratory with Daniel Ruiz (E.N.S.E.E.I.H.T)
More informationA Jacobi Davidson Method with a Multigrid Solver for the Hermitian Wilson-Dirac Operator
A Jacobi Davidson Method with a Multigrid Solver for the Hermitian Wilson-Dirac Operator Artur Strebel Bergische Universität Wuppertal August 3, 2016 Joint Work This project is joint work with: Gunnar
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationJOHANNES KEPLER UNIVERSITY LINZ. Multiharmonic Finite Element Analysis of a Time-Periodic Parabolic Optimal Control Problem
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Multiharmonic Finite Element Analysis of a Time-Periodic Parabolic Optimal Control Problem Ulrich Langer Institute of Computational
More informationINTRODUCTION TO MULTIGRID METHODS
INTRODUCTION TO MULTIGRID METHODS LONG CHEN 1. ALGEBRAIC EQUATION OF TWO POINT BOUNDARY VALUE PROBLEM We consider the discretization of Poisson equation in one dimension: (1) u = f, x (0, 1) u(0) = u(1)
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationRobust solution of Poisson-like problems with aggregation-based AMG
Robust solution of Poisson-like problems with aggregation-based AMG Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Paris, January 26, 215 Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationAggregation-based algebraic multigrid
Aggregation-based algebraic multigrid from theory to fast solvers Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire CEMRACS, Marseille, July 18, 2012 Supported by the Belgian FNRS
More informationToward black-box adaptive domain decomposition methods
Toward black-box adaptive domain decomposition methods Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI joint work with Victorita Dolean (Univ. Nice Sophia-Antipolis)
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationJOHANNES KEPLER UNIVERSITY LINZ. A Robust Preconditioned MinRes-Solver for Time-Periodic Eddy Current Problems
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics A Robust Preconditioned MinRes-Solver for Time-Periodic Eddy Current Problems Michael Kolmbauer Institute of Computational Mathematics,
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationJOHANNES KEPLER UNIVERSITY LINZ. Institute of Computational Mathematics
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Schur Complement Preconditioners for Multiple Saddle Point Problems of Block Tridiagonal Form with Application to Optimization Problems
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationOptimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms
Optimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms Marcus Sarkis Worcester Polytechnic Inst., Mass. and IMPA, Rio de Janeiro and Daniel
More informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More informationMultilevel and Adaptive Iterative Substructuring Methods. Jan Mandel University of Colorado Denver
Multilevel and Adaptive Iterative Substructuring Methods Jan Mandel University of Colorado Denver The multilevel BDDC method is joint work with Bedřich Sousedík, Czech Technical University, and Clark Dohrmann,
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More information1. Fast Iterative Solvers of SLE
1. Fast Iterative Solvers of crucial drawback of solvers discussed so far: they become slower if we discretize more accurate! now: look for possible remedies relaxation: explicit application of the multigrid
More informationA BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationA Comparison of Adaptive Algebraic Multigrid and Lüscher s Inexact Deflation
A Comparison of Adaptive Algebraic Multigrid and Lüscher s Inexact Deflation Andreas Frommer, Karsten Kahl, Stefan Krieg, Björn Leder and Matthias Rottmann Bergische Universität Wuppertal April 11, 2013
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationJOHANNES KEPLER UNIVERSITY LINZ. A Robust Preconditioner for Distributed Optimal Control for Stokes Flow with Control Constraints
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics A Robust Preconditioner for Distributed Optimal Control for Stoes Flow with Control Constraints Marus Kollmann Doctoral Program "Computational
More informationEFFICIENT MULTIGRID BASED SOLVERS FOR ISOGEOMETRIC ANALYSIS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK EFFICIENT MULTIGRID BASED SOLVERS FOR ISOGEOMETRIC
More informationJOHANNES KEPLER UNIVERSITY LINZ. Institute of Computational Mathematics
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Robust Multigrid for Isogeometric Analysis Based on Stable Splittings of Spline Spaces Clemens Hofreither Institute of Computational
More informationMultigrid Method ZHONG-CI SHI. Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China. Joint work: Xuejun Xu
Multigrid Method ZHONG-CI SHI Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China Joint work: Xuejun Xu Outline Introduction Standard Cascadic Multigrid method(cmg) Economical
More informationPreconditioners for reduced saddle point systems arising in elliptic PDE-constrained optimization problems
Zeng et al. Journal of Inequalities and Applications 205 205:355 DOI 0.86/s3660-05-0879-x RESEARCH Open Access Preconditioners for reduced saddle point systems arising in elliptic PDE-constrained optimization
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationA STOKES INTERFACE PROBLEM: STABILITY, FINITE ELEMENT ANALYSIS AND A ROBUST SOLVER
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) O. Nevanlinna and R. Rannacher (assoc. eds.) Jyväskylä,
More informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationSome Geometric and Algebraic Aspects of Domain Decomposition Methods
Some Geometric and Algebraic Aspects of Domain Decomposition Methods D.S.Butyugin 1, Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 Abstract Some geometric and algebraic aspects of various domain decomposition
More informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More informationAlgebraic Multigrid Methods for the Oseen Problem
Algebraic Multigrid Methods for the Oseen Problem Markus Wabro Joint work with: Walter Zulehner, Linz www.numa.uni-linz.ac.at This work has been supported by the Austrian Science Foundation Fonds zur Förderung
More informationOn the interplay between discretization and preconditioning of Krylov subspace methods
On the interplay between discretization and preconditioning of Krylov subspace methods Josef Málek and Zdeněk Strakoš Nečas Center for Mathematical Modeling Charles University in Prague and Czech Academy
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationContents. Preface... xi. Introduction...
Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More information