From the invention of the Schwarz method to the Best Current Methods for Oscillatory Problems: Part 2
|
|
- Osborn Bryant
- 5 years ago
- Views:
Transcription
1 and From the invention of the Schwarz method to the Best Current Methods for Oscillatory Problems: Part 2 martin.gander@unige.ch University of Geneva Woudschoten, October 2014
2 Are There Two Equations? ( +k 2 )u = f or ( η)u = f, η > 0? and
3 Are There Two Equations? ( +k 2 )u = f or ( η)u = f, η > 0? Textbook by Zauderer (1989) and
4 Are There Two Equations? ( +k 2 )u = f or ( η)u = f, η > 0? Textbook by Zauderer (1989) Leslie, Bryant, McAveney (1973) Rosmond and Faulkner (1976) and
5 Are There Two Equations? ( +k 2 )u = f or ( η)u = f, η > 0? Textbook by Zauderer (1989) (1859), Hertz... Leslie, Bryant, McAveney (1973) Rosmond and Faulkner (1976) and
6 Are There Two Equations? ( +k 2 )u = f or ( η)u = f, η > 0? Textbook by Zauderer (1989) (1859), Hertz... Leslie, Bryant, McAveney (1973) Rosmond and Faulkner (1976) and
7 Are There Two Equations? ( +k 2 )u = f or ( η)u = f, η > 0? Textbook by Zauderer (1989) (1859), Hertz... Leslie, Bryant, McAveney (1973) Rosmond and Faulkner (1976) and
8 Are There Two Equations? ( +k 2 )u = f or ( η)u = f, η > 0? Textbook by Zauderer (1989) (1859), Hertz... Leslie, Bryant, McAveney (1973) Rosmond and Faulkner (1976) and
9 Are There Two Equations? ( +k 2 )u = f or ( η)u = f, η > 0? Textbook by Zauderer (1989) (1859), Hertz... Leslie, Bryant, McAveney (1973) Rosmond and Faulkner (1976) x Solution of the equation 0.5 y x y x Solution of Laplaces equation x and
10 with ILU, Multigrid or Schwarz Example: Open cavity with point source at the center QMR ILU( 0 ) ILU(1e-2) k it Mflops it Mflops it Mflops > 2000 > > 2000 > and
11 with ILU, Multigrid or Schwarz Example: Open cavity with point source at the center QMR ILU( 0 ) ILU(1e-2) k it Mflops it Mflops it Mflops > 2000 > > 2000 > k Smoothing steps 2.5π 5π 10π 20π Iterative ν = 2 7 div div div Preconditioner ν = Iterative ν = 5 7 stag div div Preconditioner ν = and
12 with ILU, Multigrid or Schwarz Example: Open cavity with point source at the center QMR ILU( 0 ) ILU(1e-2) k it Mflops it Mflops it Mflops > 2000 > > 2000 > k Smoothing steps 2.5π 5π 10π 20π Iterative ν = 2 7 div div div Preconditioner ν = Iterative ν = 5 7 stag div div Preconditioner ν = k Overlap 10π 20π 40π 80π 160π Iterative h div div div div div Preconditioner h Iterative fixed div div div div div Preconditioner fixed and
13 with ILU, Multigrid or Schwarz Example: Open cavity with point source at the center QMR ILU( 0 ) ILU(1e-2) k it Mflops it Mflops it Mflops > 2000 > > 2000 > k Smoothing steps 2.5π 5π 10π 20π Iterative ν = 2 7 div div div Preconditioner ν = Iterative ν = 5 7 stag div div Preconditioner ν = k Overlap 10π 20π 40π 80π 160π Iterative h div div div div div Preconditioner h Iterative fixed div div div div div Preconditioner fixed and Why it is Difficult to Solve Problems with Classical Iterative Methods, Ernst and MJG 2012
14 New Method 1: Preconditioner Enquist, Ying 2010: Preconditioner for the Equation The paper introduces the sweeping preconditioner, which is highly efficient for iterative solutions of the variable coefficient equation including very high frequency problems. The first central idea of this novel approach is to construct an approximate factorization of the discretized equation by sweeping the domain layer by layer, starting from an absorbing layer or boundary condition. and
15 Numerical Experiment (Enquist, Ying 2010) and
16 New Method 2: DDM Chen, Xiang 2012: A Domain Decomposition Method for Equations in Unbounded Domain The method is based on the decomposition of the domain into non-overlapping layers and the idea of source transfer which transfers the sources equivalently layer by layer so that the solution in the final layer can be solved using a PML method defined locally outside the last two layers. and
17 Numerical Experiment (Chen, Xiang 2012) and
18 New Method 3: Stolk 2013: A rapidly converging domain decomposition method for the equation A new domain decomposition method is introduced for the heterogeneous 2-D and 3-D equations. Transmission conditions based on the perfectly matched layer (PML) are derived that avoid artificial reflections and match incoming and outgoing waves at the subdomain interfaces. Our most remarkable finding concerns the situation where the domain is split into many thin layers along one of the axes, say J subdomains numbered from 1 to J. Following [3] we will also call these quasi 2-D subdomains. Generally, an increase in the number of subdomains leads to an increase in the number of iterations required for convergence. Here we propose and study a method where the number of iterations is essentially independent of the number of subdomains. and
19 Numerical Experiments (Stolk 2013) and
20 Numerical Experiments (Stolk 2013) and
21 Underlying Fundamental Algorithms 1 An optimal discrete algorithm based on matrix factorization and Au = y ( +k 2 )u = f u 1m u 11 D 1 L 1,2 L 2,1 D Ln 1,n L n,n 1 in Ω = (0,1) (0,π) D n u 1 u 2. u n u nm u n1 x = f 1 f 2. f n
22 Block Factorization The block LU factorization of the block tridiagional matrix is D 1 L 1,2 L 2,1 D 2 L 2,3 A = = LU L n,n 1 D n T 1 L 2,1 T 2 = L n,n 1 T n I T 1 1 L 1,2 I T 1 2 L 2, where the diagonal blocks T i satisfy the recurrence { D1 i = 1, T i = D i L i,i 1 T 1 i 1 L i 1,i 1 < i n. I and
23 Forward and Backward Substitution Au = LUu = f Lv = f, Uu = v Forward substitution: T 1 L 2,1 T L n,n 1 T n v 1 =T 1 1 f 1 v 1 v 2. v n f 1 = f 2. f n v 2 =T 1 2 (f 2 L 2,1 v 1 )=T 1 2 (f 2 L 2,1 T 1 1 f 1 )=T 1 2 f 2 v 3 =T 1 3 (f 3 L 3,2 v 2 )=T 1 3 (f 3 L 3,2 T 1 2 f 2 )=T 1 3 f 3... and
24 Forward and Backward Substitution Au = LUu = f Lv = f, Uu = v Forward substitution: T 1 L 2,1 T L n,n 1 T n v 1 =T 1 1 f 1 v 1 v 2. v n f 1 = f 2. f n v 2 =T 1 2 (f 2 L 2,1 v 1 )=T 1 2 (f 2 L 2,1 T 1 1 f 1 )=T 1 2 f 2 v 3 =T 1 3 (f 3 L 3,2 v 2 )=T 1 3 (f 3 L 3,2 T 1 2 f 2 )=T 1 3 f 3.. = Source transfer: f j = f j L j,j 1 T 1 j 1 f j 1 and v n = u n. and
25 Forward and Backward Substitution Au = LUu = f Lv = f, Uu = v Forward substitution: T 1 L 2,1 T L n,n 1 T n v 1 =T 1 1 f 1 v 1 v 2. v n f 1 = f 2. f n v 2 =T 1 2 (f 2 L 2,1 v 1 )=T 1 2 (f 2 L 2,1 T 1 1 f 1 )=T 1 2 f 2 v 3 =T 1 3 (f 3 L 3,2 v 2 )=T 1 3 (f 3 L 3,2 T 1 2 f 2 )=T 1 3 f 3.. = Source transfer: f j = f j L j,j 1 T 1 j 1 f j 1 and v n = u n = preconditioner: forward and backward sweep. and
26 Underlying Fundamental Algorithms 2 π y Ω 1 Γ Ω Nataf et al (1994): algorithm ( +k 2 )u1 n = f in Ω 1 x u1 n +DtN 2(u1 n) = xu2 n 1 +DtN 2 (u2 n 1 ) on Γ ( +k 2 )u2 n = f in Ω 2 x u2 n DtN 1(u2 n) = xu1 n 1 DtN 1 (u1 n 1 ) on Γ This algorithm converges in two iterations, x and
27 Underlying Fundamental Algorithms 2 π y Ω 1 Γ 21 Γ 12 Ω Nataf et al (1994): algorithm ( +k 2 )u1 n = f in Ω 1 x u1 n +DtN 2(u1 n) = xu2 n 1 +DtN 2 (u2 n 1 ) on Γ 12 ( +k 2 )u2 n = f in Ω 2 x u2 n DtN 1(u2 n) = xu1 n 1 DtN 1 (u1 n 1 ) on Γ 21 This algorithm converges in two iterations, independantly of the overlap! x and
28 Underlying Fundamental Algorithms 2 π y Γ 12 Γ 23 Γ 34 Γ 45 Ω 1 Ω 2 Ω 3 Ω 4 Ω Nataf et al (1994): algorithm ( +k 2 )ui n = f in Ω i x ui n +DtN i,i+1 (ui n) = xu n 1 i+1 +DtN i,i+1(u n 1 i+1 ) on Γ i,i+1 x ui n DtN i,i 1 (ui n) = xu n 1 i 1 DtN i,i 1(u n 1 i 1 ) on Γ i 1,i With N subdomains, it converges in N iterations, x and
29 Underlying Fundamental Algorithms 2 π y Γ 12 Γ 23 Γ 34 Γ 45 Ω 1 Ω 2 Ω 3 Ω 4 Ω Nataf et al (1994): algorithm ( +k 2 )ui n = f in Ω i x ui n +DtN i,i+1 (ui n) = xu n 1 i+1 +DtN i,i+1(u n 1 i+1 ) on Γ i,i+1 x ui n DtN i,i 1 (ui n) = xu n 1 i 1 DtN i,i 1(u n 1 i 1 ) on Γ i 1,i With N subdomains, it converges in N iterations, or in two when sweeping back and forth once, independently of N = Stolk DDM. x and
30 Underlying Fundamental Algorithms 2 π y Γ 12 Γ 23 Γ 34 Γ 45 Ω 1 Ω 2 Ω 3 Ω 4 Ω Nataf et al (1994): algorithm ( +k 2 )ui n = f in Ω i x ui n +DtN i,i+1 (ui n) = xu n 1 i+1 +DtN i,i+1(u n 1 i+1 ) on Γ i,i+1 x ui n DtN i,i 1 (ui n) = xu n 1 i 1 DtN i,i 1(u n 1 i 1 ) on Γ i 1,i With N subdomains, it converges in N iterations, or in two when sweeping back and forth once, independently of N = Stolk DDM. x and It is the continuous formulation of block LU!
31 1 DtN operator is non local. Fourier symbol is DtN(u) = k 2 ξ 2 Use approximations for practical algorithms: absorbing boundary conditions (ABCs) k 2 ξ 2 p +qξ 2. and perfectly matched layers (PMLs) k 2 ξ 2 P(ξ) Q(ξ). and
32 1 DtN operator is non local. Fourier symbol is DtN(u) = k 2 ξ 2 Use approximations for practical algorithms: absorbing boundary conditions (ABCs) k 2 ξ 2 p +qξ 2. and perfectly matched layers (PMLs) k 2 ξ 2 P(ξ) Q(ξ). This leads to Optimized (OSM): based on subdomain decomposition can run in parallel or sweeping and
33 1 DtN operator is non local. Fourier symbol is DtN(u) = k 2 ξ 2 Use approximations for practical algorithms: absorbing boundary conditions (ABCs) k 2 ξ 2 p +qξ 2. and perfectly matched layers (PMLs) k 2 ξ 2 P(ξ) Q(ξ). This leads to Optimized (OSM): based on subdomain decomposition can run in parallel or sweeping Examples: OO0 and OO2 (Japhet, Nataf 1999), OSM (many in , MJG, Kwok 2010), (Chen, Xiang 2012), DDM (Stolk 2013) and
34 2 Block factorization: need to approximate dense matrices T i { D1 i = 1, T i = D i L i,i 1 T 1 i 1 L i 1,i 1 < i n. = Use algebraic (compute and approximate the limit by a sparse matrix) or continuous techniques (as in OSM) and
35 2 Block factorization: need to approximate dense matrices T i { D1 i = 1, T i = D i L i,i 1 T 1 i 1 L i 1,i 1 < i n. = Use algebraic (compute and approximate the limit by a sparse matrix) or continuous techniques (as in OSM) This leads to Approximate Factorization Preconditioners: based on grid layers instead of subdomains sweep forward and backward Examples: Frequency Filtering (Wittum 1992), AILU (G, Nataf 2000), preconditioner (Enquist, Ying 2010) and
36 2 Block factorization: need to approximate dense matrices T i { D1 i = 1, T i = D i L i,i 1 T 1 i 1 L i 1,i 1 < i n. = Use algebraic (compute and approximate the limit by a sparse matrix) or continuous techniques (as in OSM) This leads to Approximate Factorization Preconditioners: based on grid layers instead of subdomains sweep forward and backward Examples: Frequency Filtering (Wittum 1992), AILU (G, Nataf 2000), preconditioner (Enquist, Ying 2010) Disadvantage: can not use adaptive mesh refinement! and MJG, Zhang 2014: Iterative Solvers for the Equation: Factorizations, Preconditioners and
37 s Equations James C., in Faradays Lines of Force (1856): All the mathematical sciences are founded on relations between physical laws and laws of numbers, so that the aim of exact science is to reduce the problems of nature to the determination of quantities by operations with numbers. The original s equations (20 equations in 20 unknowns) are nowadays written in the form ε E t curl H+σE = J, µ H Related form using only the electric field: t +curl E = 0 ε 2 E t 2 +σ E t +curl 1 curl E = J µ t and
38 Time Harmonic and Time Discretized Variants Time harmonic s equations: E(x,t) = Re(E(x)exp(iωt)) H(x,t) = Re(H(x)exp(iωt)) = (σ +iωε)e curl H = J, iωµh+curl E = 0. Time discretized s equations: ( ε En+1 E n t curl H n+1 ) ( +H n 2 +σ E n+1 ) +E n 2 = J, ( µ Hn+1 H n t +curl E n+1 ) +E n 2 = 0. = (σ +ε η)e curl H = J, µ ηh+curl E = g, where (E,H) := (E n+1,h n+1 ), η := 2 t, J := J ηεe n +2σE n curl H n, g = ηµh n curl E n. and
39 Are There Two s Equations? Eliminating H, we obtain in the time harmonic case ( ω 2 ε+iωσ)e+curl 1 curl E = iωj. µ and in the time discretized case (ε η +σ)e+curl 1 µ curl E = J curl η Many solvers for the time discretized case: 1 µ η g. Domain decomposition: Toselli 1997, Toselli, Widlund, Wohlmuth 1998, Zou 2011,... Multigrid: Hiptmair 1999, Hiptmair, Xu 2007,... Very few solvers for the time harmonic case with large wave number! and
40 for s Equation x z Γ 21 Γ 12 Ω 2 Ω 1 iωεe 1,n +curl H 1,n σe 1,n =J in Ω 1 iωµh 1,n +curl E 1,n =0 in Ω 1 (Bn 1 +S 1 Bn 2 )(E 1,n,H 1,n )=(Bn 1 +S 1 Bn 2 )(E 2,n 1,H 2,n 1 )on Γ 12 iωεe 2,n +curl H 2,n σe 2,n =J in Ω 2 iωµh 2,n +curl E 2,n =0 in Ω 2 (Bn 2 +S 2 Bn 1 )(E 2,n,H 2,n )=(Bn 2 +S 2 Bn 1 )(E 1,n 1,H 1,n 1 )on Γ 21 where Bn is the impedance operator Ω Bn(E,H) := n E +n (H n) = s. Z y and
41 Result for One New Variant Theorem (Dolean, G, Gerardo-Giorda (2009)) For σ = 0, if S l have the Fourier symbol [ k 2 σ l := F(S l ) = γ y kz 2 2k y k z l 2k y k z kz 2 ky 2 ], γ l C(k z,k y ), then the convergence factor of the new algorithm is ρ= ( k 2 ω 2 i ω) 2 1 γ 1 ( k ( k 2 ω 2 +i ω) 2 1 γ 2 ( 2 ω 2 +i ω) 2 1 γ 1 ( k 2 ω 2 i ω) 2 1 γ 2 ( Theorem ( Method (2009)) 1 k 2 ω 2 +i ω) 2 k 2 ω 2 2 L k 2 ω 2 i ω) 2e 2 The choice γ l = 1/( k 2 ω 2 +i ω) 2 is optimal, since then ρ 0, for all Fourier modes k; the method converges in two iterations. Similar results for two other variants. and
42 Optimized and 1. taking γ 1 = γ 2 = 0: classical algorithm 2. taking γ 1 = γ 2 = 1 s iω εµ k 2 s+iω εµ ( ρ 2 (ω,ε,µ,l, k,s) = with s C: k 2 ω 2 εµ s k 2 ω 2 εµ+s 1 3. taking γ 1 = γ 2 = k 2 2ω 2 εµ+2iω εµs ρ 3 (ω,ε,µ,l, k,s)= 4. taking γ l = 1 s l iω εµ k 2 k 2 ω 2 εµ iω εµ k 2 ω 2 εµ+iω εµ ) 2 e 2 s l +iω εµ, l = 1,2, s 1 s 2. with s C: k 2 ω 2 εµl 1 2 ρ 2(ω,ε,µ,L, k,s) 5. taking γ l = 1 k 2 2ω 2 εµ+2iω εµs l, l = 1,2, s 1 s 2.
43 Optimal Choice Based on Asymptotic Analysis Based on equioscillation, we find for the choice s l = p l (1+i) the optimized results: with overlap, L = h without overlap, L = 0 Case ρ parameters ρ parameters 1 1 k + ω 2 h none 1 none 2 1 2C 1 6 ωh 1 3 p = C 1 3ω 2 h (k+ 2 ω 2 ) h 3 p = (k2 + ω2 ) h p 5 C 10 ω h 1 1 = C 5ω h p 2 = C 5ω h 5 3 p (k 2 + ω 2 ) = (k2 + ω2 ) 2 5 h h 1 5 p 2= (k2 + ω2 ) h C 1 4ω C h p = CC 1 4ω 2 h 2(k ω 2 ) 4 1 C(k 2 C h p = + ω 2 ) h 3 p 1 = C 8ω C C 8ω h C (k2 + ω2 ) 1 8 where ω = ω εµ and C ω = min ( k 2 + ω2, ω 2 k 2 ). C 1 4 h h p 2 = C 8ω C 3 4 h 4 3 p 1= (k2 + ω2 ) 8 3 C h 1 4 p 2= (k2 + ω2 ) 8 1 C 3 4 h 3 4 and
44 Transverse Electric Waves Joint work with Dolean and Gerardo-Giorda (2009) We consider the transverse electric waves problem (TE) in Ω = (0,1) 2 : Subdomains Ω 1 = (0,β) (0,1) and Ω 2 = (α,1) (0,1) Overlap L = β α pulsation ω = 2π relative residual reduction of 10 6 with overlap, L = h without overlap, L = 0 h 1/32 1/64 1/128 1/256 1/32 1/64 1/128 1/256 Case 1 27(21) 47(27) 72(33) 118(45) -(73) -(100) -(138) -(181) Case 2 19(14) 19(15) 22(17) 26(19) 41(26) 57(34) 79(40) 111(47) Case 3 13(13) 14(14) 17(17) 21(18) 41(23) 56(25) 80(28) 115(35) Case 4 19(14) 21(16) 24(18) 27(19) 31(24) 35(28) 41(30) 47(33) Case 5 13(13) 15(15) 17(18) 19(19) 56(26) 64(30) 76(32) 76(35) and
45 Comparison of the Asymptotics and iterations 10 2 case 1 case 2 case 3 case 4 case5 O(h 1 ) O(h 1/3 ) O(h 1/3 ) O(h 1/5 ) O(h 1/5 ) iterations 10 2 case 2 case 3 case 4 case5 O(h 1/2 ) O(h 1/2 ) O(h 1/4 ) O(h 1/4 ) h overlapping case 10 2 non-overlapping case h
46 The Chicken Problem Joint work with Dolean (2009) and Heating a chicken in our Whirlpool Talent Combi 4 microwave oven
47 The Chicken Problem The electric field intensity in the chicken and
48 Scattering of a plane wave by a dielectric cylinder Joint Work with Dolean, El Bouajaji and Lanteri (2010) Method L 2 err. E z L 2 err. E z N s # iter # iter Classical Optimized Classical Optimized DGTH-P ( 6.1) ( 4.7) DGTH-P (10.7) ( 6.7) DGTH-P (15.0) ( 8.2) DGTH-P (19.5) (10.3) DGTH-P pk ( 7.2) ( 5.1) EZ and
49 Exposure of head tissue to plane wave 1800 MHz Joint with Dolean, El Bouajaji, Lanteri, Perrussel (2010) Z Y and REX X # d.o.f Ns RAM LU time LU iter Elapsed time 26,854, /3.1 GB 496 sec sec /1.2 GB 132 sec sec 44,491, /3.2 GB 528 sec sec /1.3 GB 142 sec sec DGTH-P1, single precision MUMPS LU, contour lines of Ex Skin (εr = and σ = 1.23 S/m), skull (εr = and σ = 0.43 S/m), CSF (Cerebro Spinal Fluid) (εr = and σ = 2.92 S/m), brain (εr = and σ = 1.15 S/m)
50 COBRA Cavity Joint work with Dolean, Lanteri, Lee and Peng (2014): 33 subdomains, h = λ 0 /4, DOFs plane wave incident opon the cavity aperture f = 17.6GHz Orthogonal incidence: 18 iterations to reach 1e-3 (29 without optimization) Oblique incidence: 23 iterations to reach 1e-3 (51 without optimization) and
51 COBRA Cavity Results 1 and
52 COBRA Cavity Results 2 and
53 s New preconditioners based on one fundamental idea: 1. Block factorizations 2. Schwarz methods for strip decompositions Both are direct solvers, if DtN or Schur complements are used: = methods Approximations lead to practical algorithms: = methods Examples: OO0, OO2, Generalized Schwarz, Robin Schwarz,, Frequency filtering, AILU, Preconditioner, and many more Preprints are available at gander and MJG, Zhang 2014: Iterative Solvers for the Equation: Factorizations, Preconditioners and
Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping
Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping M. El Bouajaji 1, V. Dolean 2, M. J. Gander 3 and S. Lanteri 1 1 Introduction Transmission conditions
More informationOptimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients
Optimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients Victorita Dolean, Martin Gander, Erwin Veneros To cite this version: Victorita Dolean, Martin Gander, Erwin Veneros. Optimized
More informationDG discretization of optimized Schwarz methods for Maxwell s equations
DG discretization of optimized Schwarz methods for Maxwell s equations Mohamed El Bouajaji, Victorita Dolean,3, Martin J. Gander 3, Stéphane Lanteri, and Ronan Perrussel 4 Introduction In the last decades,
More informationIntroduction to Domain Decomposition Methods in the numerical approximation of PDEs. Luca Gerardo-Giorda
Introduction to Domain Decomposition Methods in the numerical approximation of PDEs Luca Gerardo-Giorda Plan 0 Motivation Non-overlapping domain decomposition methods Multi-domain formulation Variational
More informationHybrid iterative-direct domain decomposition based solvers for the time-harmonic Maxwell equations
Hybrid iterative-direct domain decomposition based solvers for the time-harmonic Maxwell equations 1st workshop of the joint INRIA-UIUC laboratory for Petascale Computing Victorita Dolean 1, Mohamed El
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J SCI COMPUT Vol 34, No 4, pp A48 A7 c Society for Industrial and Applied Mathematics OPTIMIZED SCHWARZ METHODS FOR THE TIME-HARMONIC MAXWELL EQUATIONS WITH DAMPING M EL BOUAJAJI, V DOLEAN, M J GANDER,
More informationAdditive Schwarz method for scattering problems using the PML method at interfaces
Additive Schwarz method for scattering problems using the PML method at interfaces Achim Schädle 1 and Lin Zschiedrich 2 1 Zuse Institute, Takustr. 7, 14195 Berlin, Germany schaedle@zib.de 2 Zuse Institute,
More informationOptimized Schwarz methods for the time-harmonic Maxwell equations with damping
Optimized Schwarz methods for the time-harmonic Maxwell equations with damping Mohamed El Bouajaji, Victorita Dolean, Martin Gander, Stephane Lanteri To cite this version: Mohamed El Bouajaji, Victorita
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 informationOptimized Schwarz Methods for Maxwell equations
Optimized Schwarz Methods for Maxwell equations Victorita Dolean, Martin Gander, Luca Gerardo-Giorda To cite this version: Victorita Dolean, Martin Gander, Luca Gerardo-Giorda. Optimized Schwarz Methods
More informationTwo-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations
Two-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations Marcella Bonazzoli 2, Victorita Dolean 1,4, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier
More informationEfficient domain decomposition methods for the time-harmonic Maxwell equations
Efficient domain decomposition methods for the time-harmonic Maxwell equations Marcella Bonazzoli 1, Victorita Dolean 2, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier 4 1 Inria Saclay (Defi
More informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationOn domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption
On domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Emerson Roxar,
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationOptimized Schwarz methods with overlap for the Helmholtz equation
Optimized Schwarz methods with overlap for the Helmholtz equation Martin J. Gander and Hui Zhang Introduction For the Helmholtz equation, simple absorbing conditions of the form n iω were proposed as transmission
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationMulti-Domain Approaches for the Solution of High-Frequency Time-Harmonic Propagation Problems
Académie universitaire Wallonie Europe Université de Liège Faculté des Sciences Appliquées Collège de doctorat en Électricité, électronique et informatique Multi-Domain Approaches for the Solution of High-Frequency
More informationShifted Laplace and related preconditioning for the Helmholtz equation
Shifted Laplace and related preconditioning for the Helmholtz equation Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Schlumberger Gould Research), Martin Gander (Geneva) Douglas
More informationLinear and Non-Linear Preconditioning
and Non- martin.gander@unige.ch University of Geneva June 2015 Invents an Iterative Method in a Letter (1823), in a letter to Gerling: in order to compute a least squares solution based on angle measurements
More informationDeflated Krylov Iterations in Domain Decomposition Methods
305 Deflated Krylov Iterations in Domain Decomposition Methods Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 1 Introduction The goal of this research is an investigation of some advanced versions of
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationA PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY
More informationA Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation
A Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation Y. Boubendir, X. Antoine, C. Geuzaine June 28, 2011 Abstract This paper presents a new non-overlapping domain decomposition
More informationA Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation
A Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation Y. Boubendir, X. Antoine, C. Geuzaine May 8, 2012 Abstract This paper presents a new non-overlapping domain decomposition
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationWhy it is Difficult to Solve Helmholtz Problems with Classical Iterative Methods
Why it is Difficult to Solve Helmholtz Problems with Classical Iterative Methods Oliver G. Ernst and Martin J. Gander Abstract In contrast to the positive definite Helmholtz equation, the deceivingly similar
More informationDirichlet-Neumann and Neumann-Neumann Methods
Dirichlet-Neumann and Neumann-Neumann Methods Felix Kwok Hong Kong Baptist University Introductory Domain Decomposition Short Course DD25, Memorial University of Newfoundland July 22, 2018 Outline Methods
More informationOn the performance of the Algebraic Optimized Schwarz Methods with applications
On the performance of the Algebraic Optimized Schwarz Methods with applications Report 12-1-18 Department of Mathematics Temple University October 212 This report is available at http://www.math.temple.edu/~szyld
More informationScalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems
Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Pierre Jolivet, F. Hecht, F. Nataf, C. Prud homme Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann INRIA Rocquencourt
More informationLinear and Non-Linear Preconditioning
and Non- martin.gander@unige.ch University of Geneva July 2015 Joint work with Victorita Dolean, Walid Kheriji, Felix Kwok and Roland Masson (1845): First Idea of After preconditioning, it takes only three
More informationComputational Electromagnetics Definitions, applications and research
Computational Electromagnetics Definitions, applications and research Luis E. Tobón Pontificia Universidad Javeriana Seminario de investigación Departamento de Electrónica y Ciencias de la Computación
More informationarxiv: v1 [math.na] 11 Jul 2018
An Additive Overlapping Domain Decomposition Method for the Helmholtz Equation Wei Leng Lili Ju arxiv:1807.04180v1 [math.na] 11 Jul 2018 July 12, 2018 Abstract In this paper, we propose and analyze an
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 informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationFrom Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D
From Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D Luc Giraud N7-IRIT, Toulouse MUMPS Day October 24, 2006, ENS-INRIA, Lyon, France Outline 1 General Framework 2 The direct
More informationA High-Performance Parallel Hybrid Method for Large Sparse Linear Systems
Outline A High-Performance Parallel Hybrid Method for Large Sparse Linear Systems Azzam Haidar CERFACS, Toulouse joint work with Luc Giraud (N7-IRIT, France) and Layne Watson (Virginia Polytechnic Institute,
More informationOpen-source finite element solver for domain decomposition problems
1/29 Open-source finite element solver for domain decomposition problems C. Geuzaine 1, X. Antoine 2,3, D. Colignon 1, M. El Bouajaji 3,2 and B. Thierry 4 1 - University of Liège, Belgium 2 - University
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 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 informationEigenvalue Problems in Surface Acoustic Wave Filter Simulations
Eigenvalue Problems in Surface Acoustic Wave Filter Simulations S. Zaglmayr, J. Schöberl, U. Langer FWF Start Project Y-192 3D hp-finite Elements: Fast Solvers and Adaptivity JKU + RICAM, Linz in collaboration
More informationSome examples in domain decomposition
Some examples in domain decomposition Frédéric Nataf nataf@ann.jussieu.fr, www.ann.jussieu.fr/ nataf Laboratoire J.L. Lions, CNRS UMR7598. Université Pierre et Marie Curie, France Joint work with V. Dolean
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 informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationHybrid (DG) Methods for the Helmholtz Equation
Hybrid (DG) Methods for the Helmholtz Equation Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology Contributions by
More informationOptimized Schwarz Waveform Relaxation: Roots, Blossoms and Fruits
Optimized Schwarz Waveform Relaxation: Roots, Blossoms and Fruits Laurence Halpern 1 LAGA, Université Paris XIII, 99 Avenue J-B Clément, 93430 Villetaneuse, France, halpern@math.univ-paris13.fr 1 Introduction:
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationDomain Decomposition Method for Electromagnetic Scattering Problems: Application to EUV Lithography
Domain Decomposition Method for Electromagnetic Scattering Problems: Application to EUV Lithography Lin Zschiedrich, Sven Burger, Achim Schädle, Frank Schmidt Zuse Institute Berlin, JCMwave GmbH NUSOD,
More informationOPTIMIZED SCHWARZ METHODS
SIAM J. NUMER. ANAL. Vol. 44, No., pp. 699 31 c 006 Society for Industrial and Applied Mathematics OPTIMIZED SCHWARZ METHODS MARTIN J. GANDER Abstract. Optimized Schwarz methods are a new class of Schwarz
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
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 informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
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 informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationUne méthode parallèle hybride à deux niveaux interfacée dans un logiciel d éléments finis
Une méthode parallèle hybride à deux niveaux interfacée dans un logiciel d éléments finis Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI Victorita Dolean (Univ. Nice
More information27. A Non-Overlapping Optimized Schwarz Method which Converges with Arbitrarily Weak Dependence on h
Fourteent International Conference on Domain Decomposition Metods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 003 DDM.org 7. A Non-Overlapping Optimized Scwarz Metod wic Converges
More informationFast evaluation of boundary integral operators arising from an eddy current problem. MPI MIS Leipzig
1 Fast evaluation of boundary integral operators arising from an eddy current problem Steffen Börm MPI MIS Leipzig Jörg Ostrowski ETH Zürich Induction heating 2 Induction heating: Time harmonic eddy current
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationITERATIVE METHODS FOR SPARSE LINEAR SYSTEMS
ITERATIVE METHODS FOR SPARSE LINEAR SYSTEMS YOUSEF SAAD University of Minnesota PWS PUBLISHING COMPANY I(T)P An International Thomson Publishing Company BOSTON ALBANY BONN CINCINNATI DETROIT LONDON MADRID
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationNewton-Krylov-Schwarz Method for a Spherical Shallow Water Model
Newton-Krylov-Schwarz Method for a Spherical Shallow Water Model Chao Yang 1 and Xiao-Chuan Cai 2 1 Institute of Software, Chinese Academy of Sciences, Beijing 100190, P. R. China, yang@mail.rdcps.ac.cn
More informationOn complex shifted Laplace preconditioners for the vector Helmholtz equation
On complex shifted Laplace preconditioners for the vector Helmholtz equation C. Vuik, Y.A. Erlangga, M.B. van Gijzen, C.W. Oosterlee, D. van der Heul Delft Institute of Applied Mathematics c.vuik@tudelft.nl
More informationMultispace and Multilevel BDDC. Jan Mandel University of Colorado at Denver and Health Sciences Center
Multispace and Multilevel BDDC Jan Mandel University of Colorado at Denver and Health Sciences Center Based on joint work with Bedřich Sousedík, UCDHSC and Czech Technical University, and Clark R. Dohrmann,
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationComparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses
Journal of the Korean Physical Society, Vol. 58, No. 4, April 2011, pp. 1014 1020 Comparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses Jian-Bo Chen, Yan Shen,
More informationIndefinite and physics-based preconditioning
Indefinite and physics-based preconditioning Jed Brown VAW, ETH Zürich 2009-01-29 Newton iteration Standard form of a nonlinear system F (u) 0 Iteration Solve: Update: J(ũ)u F (ũ) ũ + ũ + u Example (p-bratu)
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 informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationarxiv: v1 [math.na] 28 Feb 2008
BDDC by a frontal solver and the stress computation in a hip joint replacement arxiv:0802.4295v1 [math.na] 28 Feb 2008 Jakub Šístek a, Jaroslav Novotný b, Jan Mandel c, Marta Čertíková a, Pavel Burda a
More informationParallelizing large scale time domain electromagnetic inverse problem
Parallelizing large scale time domain electromagnetic inverse problems Eldad Haber with: D. Oldenburg & R. Shekhtman + Emory University, Atlanta, GA + The University of British Columbia, Vancouver, BC,
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 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 informationFEM and Sparse Linear System Solving
FEM & sparse system solving, Lecture 7, Nov 3, 2017 1/46 Lecture 7, Nov 3, 2015: Introduction to Iterative Solvers: Stationary Methods http://people.inf.ethz.ch/arbenz/fem16 Peter Arbenz Computer Science
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
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 informationRapid evaluation of electrostatic interactions in multi-phase media
Rapid evaluation of electrostatic interactions in multi-phase media P.G. Martinsson, The University of Colorado at Boulder Acknowledgements: Some of the work presented is joint work with Mark Tygert and
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationTHE OPTIMISED SCHWARZ METHOD AND THE TWO-LAGRANGE MULTIPLIER METHOD FOR HETEROGENEOUS PROBLEMS. by Neil Greer
THE OPTIMISED SCHWARZ METHOD AND THE TWO-LAGRANGE MULTIPLIER METHOD FOR HETEROGENEOUS PROBLEMS by Neil Greer Submitted for the degree of Doctor of Philosophy Department of Mathematics School of Mathematical
More informationBlock-tridiagonal matrices
Block-tridiagonal matrices. p.1/31 Block-tridiagonal matrices - where do these arise? - as a result of a particular mesh-point ordering - as a part of a factorization procedure, for example when we compute
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 informationPseudospectral and High-Order Time-Domain Forward Solvers
Pseudospectral and High-Order Time-Domain Forward Solvers Qing H. Liu G. Zhao, T. Xiao, Y. Zeng Department of Electrical and Computer Engineering Duke University DARPA/ARO MURI Review, August 15, 2003
More informationAlgebraic Coarse Spaces for Overlapping Schwarz Preconditioners
Algebraic Coarse Spaces for Overlapping Schwarz Preconditioners 17 th International Conference on Domain Decomposition Methods St. Wolfgang/Strobl, Austria July 3-7, 2006 Clark R. Dohrmann Sandia National
More informationMéthodes de raffinement espace-temps
Méthodes de raffinement espace-temps Caroline Japhet Université Paris 13 & University of Maryland Laurence Halpern Université Paris 13 Jérémie Szeftel Ecole Normale Supérieure, Paris Pascal Omnes CEA,
More informationSolving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation
Introduction Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation Plane waves as basis functions Peter Monk 1 Tomi Huttunen 2 1 Department of Mathematical Sciences University
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 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 informationPerformance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
March 4-5, 2015 Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain M. Bonnasse-Gahot 1,2, H. Calandra 3, J. Diaz 1 and S. Lanteri 2
More informationTechnology USA; 3 Dept. of Mathematics, University of California, Irvine. SUMMARY
Adrien Scheuer, Leonardo Zepeda-Núñez,3*), Russell J. Hewett 2 and Laurent Demanet Dept. of Mathematics and Earth Resources Lab, Massachusetts Institute of Technology; 2 Total E&P Research & Technology
More informationCLASSICAL ITERATIVE METHODS
CLASSICAL ITERATIVE METHODS LONG CHEN In this notes we discuss classic iterative methods on solving the linear operator equation (1) Au = f, posed on a finite dimensional Hilbert space V = R N equipped
More informationAN INTRODUCTION TO DOMAIN DECOMPOSITION METHODS. Gérard MEURANT CEA
Marrakech Jan 2003 AN INTRODUCTION TO DOMAIN DECOMPOSITION METHODS Gérard MEURANT CEA Introduction Domain decomposition is a divide and conquer technique Natural framework to introduce parallelism in the
More informationCoupling Stokes and Darcy equations: modeling and numerical methods
Coupling Stokes and Darcy equations: modeling and numerical methods Marco Discacciati NM2PorousMedia 24 Dubrovnik, September 29, 24 Acknowledgment: European Union Seventh Framework Programme (FP7/27-23),
More informationSpace-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure
Space-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure Elyes Ahmed, Caroline Japhet, Michel Kern INRIA Paris Maison de la Simulation ENPC University Paris
More informationA decade of fast and robust Helmholtz solvers
A decade of fast and robust Helmholtz solvers Werkgemeenschap Scientific Computing Spring meeting Kees Vuik May 11th, 212 1 Delft University of Technology Contents Introduction Preconditioning (22-28)
More information