Inexact Data-Sparse BETI Methods by Ulrich Langer. (joint talk with G. Of, O. Steinbach and W. Zulehner)
|
|
- Winfred Johnston
- 5 years ago
- Views:
Transcription
1 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 and Institute for Computational Mathematics Johannes Kepler University Linz, Austria This paper was supported by the Austrian Sciene Fund (FWF) and by the German Research Foundation (DFG)
2 Contents 1. Introduction 2. Data-Sparse Boundary Element Domain Decompostion Methods 3. Inexact BETI Methods 4. Preconditioning and Solution of Two-fold Saddle Point Problems 5. Data-Sparse Boundary Element Preconditioners 6. Numerical Results 7. Summary and Work in Progress
3 1. Introduction A lot of FETI papers have been published since Farhat and Roux proposed this DD technique in 1991: new applications: structural mechanics, contact problems, Helmholtz,... new versions: FETI-2, FETI-DP, Algebraic versions etc. FETI Analysis History: Mandel and Tezaur (1996): Dirichlet FETI preconditioner (PC) C: Klawonn and Widlund (2000,2001): New FETI PCs; Inexact versions: Brenner (2002, 2003): Schwarz analysis! Mandel and Tezaur (2001), Klawonn, Widlund, Dryja (2002),...: FETI-DP κ(c 1 F ) O((1 + log(h/h)) 2 ) and robust!
4 BETI = Boundary Element counterpart of FETI: Langer and Steinbach (2002/2003): BETI Langer and Steinbach (2003 / DD15): BETI-FETI-Coupling Idea: Start from the skeleton formulation of elliptic BVP: p p (S i u)(x)v(x)ds x = [(N i f)(x) (S i g)(x)]v(x)ds x i=1 Γ i=1 i Γ i (1) q p q p Ri Si,hR FEM i u + Ri Si,hR BEM i u = Ri f FEM R i i Si,hR BEM i g (2) i=1 i=q+1 i=1 i=q+1
5 In this talk: = Inexact (no elimination of any unkowns!) = Data-sparse approximations of all BE matrices (V, K, K T, D): Ṽ V, K K, D D: X X = O(discretization error) and property preserving! = Two-fold saddle point problem for determining t, u, λ = Subspace iteration for n-fold saddle point problems (n = 2) = Data-sparse preconditioners for Ṽi, S i, F = B T i S + i B i Ultimate Goal for Inexact Data-Sparse BETI = asymptotically almost optimal complexity + robustness! Complexity = O(N S ln q (N S )) = O(h (d 1) ln q (h 1 ))!
6 2. Data-Sparse Boundary Element DDM Model boundary value problem (BVP): d = 2, 3 div(α û) = 0 in Ω R d and û = g on Γ = Ω = Γ D (3) Non-overlapping quasi-regular domain decomposition (DD): Ω Ω Ω 7 Ω 8 9 Ω α Ω Ω Ω Ω n 1 α Ω = Ω 1 Ω q Ω q+1 Ω p Ω i Ω j = for i j, Γ ij = Γ i Γ j, Γ i = Ω i, Γ S = p i=1 Γ i Γ C = Γ S \ Γ D = Γ S \ Γ diam(ω i ) = O(H), α(x) = α i > 0, x Ω i, i = 1,..., p (4)
7 Therefore, the weak form of our model BVP (3): Find û H 1 (Ω) with given û Γ = g H 1/2 (Γ) such that: p α i i=1 Ω i û(x) v(x)dx = 0 v H 1 0 (Ω), (5) can be rewritten in the form p α i û(x)v(x)ds x = 0 for all v H n 0 1 (Ω) (6) i i=1 Γ i by the use of Green s first formula and the fact that û(x) = 0 for x Ω i. (7)
8 The solutions of the local subproblems (7) are given by the representation formula that holds for x Ω i : û(x) = E(x, y) û(y)ds y E(x, y)û(y)ds y (8) n i n i Γ i where E(x, y) = (1/(2π)) log x y resp. E(x, y) = (1/(4π)) x y 1 is the fundamental solution of in 2D resp. 3D. Taking the trace on Γ i and the normal derivative gives Caldaron s representation of the Cauchy date on Γ i : 1 2 u i(x) := 1 1 2û Γ i (x) =... 2 t i(x) := 1 u i (x) =..., x Γ 2 n i (9) i Γ i
9 The complete Cauchy data (u i, t i = u i / n i ) on Γ i satisfy the Calderon equations ( ) ( 1 ui = 2 I K ) ( ) i V i ui 1 t i D i 2 I + (10) K i t i where (V i t i )(x) := E(x, y)t i (y)ds y, x Γ i, Γ i (K i u i )(x) := E(x, y)u i (y)ds y, x Γ i, n y Γ i (K it i )(x) := E(x, y)t i (y)ds y, x Γ i, and n x Γ i (D i u i )(x) := E(x, y)u i (y)ds y, x Γ i, n x n y Γ i
10 Inserting the second boundary integral equation of (10) into (6) we have to find û H 1/2 (Γ S ) with û Γ = g such that p [ α i D i û Γi, v Γi Γi + ( 1 ] 2 I + K i )t i, v Γi Γi = 0 (11) i=1 for all test functions v H 1/2 0 (Γ S ) = {v H 1/2 (Γ S ) : v Γ = 0}, where t i H 1/2 (Γ i ) are the unique solutions of the local variational problems [ α i V i t i, τ i Γi ( 1 ] 2 I + K i)û Γi, τ i Γi = 0 (12) for all test functions τ i H 1/2 (Γ i ). REF.: Hsiao and Wendland (1991), Costabel (1987)
11 After homogenization of the Dirichlet boundary condition via the ansatz û = ĝ + u with ĝ Γ = g and u Γ = 0, we can rewrite (11) and (12) as mixed variational problem to find t = (t 1, t 2,..., t p ) T = T 1 T 2... T p = H 1/2 (Γ 1 ) H 1/2 (Γ 2 )... H 1/2 (Γ p ) and u U = H 1/2 0 (Γ S ) such that [ α i τ i, V i t i Γi τ i, ( 1 ] 2 I + K i)u Γi Γi = τ i, g i for all τ i T i (13) for i = 1, 2,..., p and p [ α i ( 1 ] 2 I + K i)t i, v Γi Γi D i u Γi, v Γi Γi i=1 where the right-hand sides are given by the relations = F, v v U, (14) τ i, g i = α i τ i, ( 1 2 I + K i)ĝ Γi Γi i = 1, 2,..., p, (15) F, v = p α i f i, v Γi Γi = i=1 p α i D i ĝ Γi, v Γi Γi. (16) i=1
12 The Galerkin discretization of the symmetric mixed variational formulation (13) and (14) now gives immediately the mixed boundary element equations to find t h = (t 1,h, t 2,h,..., t p,h ) T h = T 1,h T 2,h... T p,h (e.g. piecewise constant) and u h U h (e.g. continuous piecewise linear) such that [ α i τ i,h, V i t i,h Γi τ i,h, ( 1 ] 2 I + K i)u h Γi Γi = τ i,h, g i (17) for all τ i,h T i,h, i = 1, 2,..., p, and p α i [ ( 1 ] 2 I + K i)t i,h, v h Γi Γi D i u h Γi, v h Γi Γi + i=1 = F, v h (18) for all v h U h. { } u u h 2 U + t t h 2 T c inf u v h 2 U + inf t τ h 2 T. (19) v h U h τ h T h
13 The symmetric mixed boundary element equations (17)-(18) now immediately gives the symmetric, but indefinite Galerkin system with dense block matrix entries: α 1 V 1,h α 1 K 1,h R 1,h t 1. g α p V p,h α p K p,h R p,h t pu =. g p, (20) α 1 R1,h K 1,h... α prp,h K p,h D h f where K i,h = [ ψl, i ( 1 2 I + K i)ϕ i m L2 (Γ i )], (21) p D h = α i Ri,hD i,h R i,h, (22) f = i=1 p Ri,hf i. (23) i=1
14 Data-Sparse Approximation via symmetric FMM: α 1 Ṽ 1,h α 1 K 1,h R 1,h.... α p Ṽ p,h α 1 R1,h K 1,h... α prp,h K p,h Results: Data-sparse and property preserving! α p Kp,h R p,h D h t 1. t pũ = g 1. g p f. (24) Quasi-optimal error estimate: { } u ũ h 2 U + t t h 2 T c inf u v h 2 U + inf t τ h 2 T. (25) v h U h τ h T h
15 3. Inexact BETI In order to avoid assambled matrices and vectors, we tear the global potential (displacements) vector u on the subdomain boundaries Γ i by introducing the individual local unknows u i = R i,h u. (26) Conversely, the global continuity of the potentials is enforced by the constraints p B i u i = 0 (27) i=1 again interconnecting the local potential vectors across the subdomain boundaries as usual in FETI.
16 By introducing Lagrange multipliers λ R L, we easily see that the one-fold saddle point problem (24) is equivalent to the twofold saddle point problem: Find t = (t 1,..., t p ) R M, u = (u 1,..., u p ) R N, and λ = (λ 1,..., λ L ) R L such that V K 0 K D B 0 B 0 t u λ = g f 0, (28) where V = blockdiag(α i Ṽ i,h ), D = blockdiag(α i Di,h ), K = blockdiag( α i Ki,h ), B = (B 1,..., B p ), g = ( g 1,..., g p ) R M and f = ( f 1,..., f p ) R N are given..
17 For the floating subdomains Ω j, j = 1,..., q, we have D j,h e j = 0, Kj,h e j = 0, Sj,h e j = 0, (29) where S j,h e j = α j ( D j,h + K j,hṽ j,h 1 K j,h ) are the FMM BE SC. If Dj,h s psd is replaced by D j,h + β j e j e j spd then there exists a unique solution v j of the regularized equations α j ( D j,h + β j e j e j )v j = f j + B j λ α j K j,h t j (30) If the right-hand side of (30) fulfils the solvability conditions for the original equations without the regularization term e j ( f j + B j λ α j K i,h t j ) = e j ( f j + B j λ) = 0, (31) then the solution v j of (30) is orthogonal to ker D j,h, i.e. e j v j = 0.
18 Now the solution u = (u 1,..., u p ) of the original system (28) can be recovered by the formulas u i = v i + γ i e i, (32) with γ i = 0 for i = q + 1,..., p and appropriatly choosen γ i for i = 1,..., q. Taking into account (29) - (32), we easily observe that the unique solution t R M, v R N, λ R L and γ R q of the three-fold saddle point problem V K 0 0 K D B 0 0 B 0 G 0 0 G 0 t v λ γ = g f 0 e (33) immediately yields the solution (t, u, λ) of the two-fold saddle point problem (28), where now D is SPD, since we replaced the singular blocks D j,h by the regularized ones for the floating subdomains. The L q matrix G and the vector e R q are defined by the relations G = (B 1 e 1,..., B q e 1 ) and e = (e 1 f 1,..., e q f q ). (34)
19 Let us now introduce the orthogonal projection P = I G(G G) 1 G : Λ := R M Λ 0 = ker G = (range G) Λ (35) Since P G = 0, the application of P to the third equation of (33) gives P Bv = 0 that excludes γ from the first three equations of (33). Let us now represent λ as λ = T 0 λ 0 + λ e (36) with known λ e = G(G G) 1 e (ker G ) = range G, fulfilling the constraints G λ e = e, and unknown T 0 λ 0 ker G, i.e. G T 0 λ 0 = 0, where λ 0 R L 0 and L 0 = dim Λ 0. The columns of T 0 span a basis of Λ 0. Now we can define t, v and λ 0 from the two-fold saddle point problem V K 0 K D B P T 0 0 T0 P B 0 t v λ 0 = g d 0, (37) where here d = f B λ e. Once t, v, λ 0 are defined from (37), we get λ from (36), γ = (G G) 1 G Bv from the third equation of (33), and, finally, u from (32).
20 4. Preconditioning and Solution of Two-fold SPP Let us rewrite the two-fold saddle point problem (37) in the form A 1 B1 0 B 1 A 2 B2 =, (38) 0 B 2 0 x 1 x 2 x 3 where A 1 = V, A 2 = D, B 1 = K and B 2 = T 0 P B. The following matrices are SPD: S 1 = A 1 = V = blockdiag(α i Ṽ i,h ) A 2 = D = blockdiag(α i Di,h + β j e j e j ) S 2 = A 2 + B 1 A 1 1 B 1 = D + K V 1 K S 3 = B 2 S 1 2 B 2 = T 0 P B(D + K V 1 K) 1 B P T 0 = F b 1 b 2 b 3
21 It is easy to see that the system matrix A 1 B1 0 S B 1 A 2 B2 1 I S 1 = B 1 S 2 I 0 B 2 0 B 2 S 3 1 BT 1 S 1 2 BT 2 I of (38) possesses a block LU-decomposition K 3 = L 3 U 3 (39) that motivates the preconditioner (denoted by hats ˆX) ( ) τ τ 2 τ 1 Ŝ 1 ˆL 3 = 2 ˆL2 = τ B 2 E 2 Ŝ 2 B 1 τ 2 Ŝ 2 (40) 3 B 2 Ŝ 3
22 Theorem 4.1 If the parameter τ 1 is chosen such that 0 < τ 1 < λ min (Ŝ 1 1 A 1) then the matrix ˆL 1 2 K 2 is symmetric and positive definite in the scalar product (u, v) 2 = (D 2 u, v) l2 with ( ) ( τ ˆL 2 = 1 Ŝ 1 A1 B, K 2 = T ) ) 1 (A, D B 1 Ŝ2 B 1 A 2 = 1 τ 1 Ŝ 1. 2 Ŝ2 ˆL 1 If the parameter τ 2 is chosen such that 0 < τ 2 < λ min ( 2 K 2) then the preconditioned matrix ˆL 1 3 K 3 is SPD in the scalar product (u, v) 3 = (D 3 u, v) l2 with ( ) D D 3 = 2 ( ˆL 1 2 K 2 τ 2 I) 0. 0 Ŝ 3
23 Theorem 4.2 If σ 1 Ŝ 1 S 1 σ 1 Ŝ 1, σ 2 Ŝ 2 S 2 σ 2 Ŝ 2, σ 3 Ŝ 3 S 3 σ 3 Ŝ 3 and the parameters τ 1 and τ 2 are chosen such that τ 1 < λ 1, τ 2 < λ 2, then the matrices ˆL 1 2 K 2 and ˆL 1 3 K 3 are symmetric and positive definite with rerspect to the scalar products (.,.) 2 and (.,.) 3 and λ min ( ˆL 1 3 K 3) λ 3, λ max ( ˆL 1 3 K 3) λ 3, κ( with explicit formulas for λ 1, λ 2, λ 3, λ 3 = λ(σ)! ˆL 1 3 K 3) λ 3 /λ 3 = Use PCG-Solver with PC ˆL 3 and SP (D 3 u, v) l2!
24 5. Data-Sparse BE Preconditioners ˆX X Ŝ1 = ˆV V, i.e. ˆVi,h Ṽi,h (local!) = Artificial Multilevel Single Layer Preconditioner proposed by Steinbach (2003): σ 1 /σ 1 = O((1 + log(h/h)) 2 ) Ŝ2 S 2 = D + K V 1 K = S BE, i.e. Ŝ BE,i S BE,i (local!) = Opposite order preconditioner Ŝ 1 BE,i = M 1 i,hṽ i,hm 1 i,h proposed by Steinbach & Wendland (1998) (dense): σ 2 /σ 2 = O(1) Ŝ3 = ˆF S 3 = F = T 0 P B(D + K V 1 K) 1 B P T 0 = Data-sparse scaled hypersingular BETI preconditioner proposed by Langer and Steinbach (2002) in our BETI paper: ˆF 1 = Q DQ: σ 3 /σ 3 = O((1 + log(h/h)) 2 )
25 Theorem 5.1 (Final Result) If σ 1 /σ 1 = O(1) (e.g. Data-sparse GMG- or AMG-PC), σ 2 /σ 2 = O(1) (e.g. opposite order preconditioner), σ 3 /σ 3 = O((1 + log(h/h)) 2 ) (e.g. scaled data-sparse hypersingular BETI preconditioner), then the PCG acceleration gives the complexity estimates: = Iter(ε) = O((1 + log H h ) log 1 ε ) = Ops parallel (ε) = O(( H h )d 1 (1 + log H h )3 log 1 ε ) Robust with respect to coefficient jumps!
26 6. Numerical Results We solve the Dirichlet potential problem div(α û) = 0 in Ω with given Dirichlet data g on Γ = Ω:
27 L N M M c N i t 1 t 2 It. D-Error e e e e e e e-6 Table 6.1 BETI two-fold saddle point formulation for 8 cubes.
28 BETI Schur BETI spp BETI 2spp L t 2 it t 2 it t 2 it 0 0 3(2) (10) (15) (17) (19) (21) Table 6.2 Comparison of different BETI versions for 8 cubes.
29 L N M M c N i t 1 t 2 It. D-error e e e e e e e-12 Table 6.3 BETI two-fold SP formulation for 8 cubes with α = 10 5.
30 7. Summary and Work in Progress We propose Inexact Data-Sparse BETI methods the complexity (w.r.t. arithmetical costs and memory) of which is up to some polylogarithmic factor proportional to the nubmber of skeleton unknowns!!! Parallel Implementation: Ops parallel (ε) = O(( H h )d 1 (1 + log H h )3 log 1 ε ) Inexact Data-Sparse BETI-FETI, especially, for non-linear problems in the FEM subdomains Application to Elasticity, Helmholtz, Maxwell etc. equations
31 I am hoping very much to see you at the DD17 held in St. Wolfgang (Austria), July 3-7, More information can be found at the preliminary DD17 homepage Homepages of our university institutes: Institute of Computational Mathematics, JKU Linz Institute of Mathematics D, TU Graz IANS, University of Stuttgart Thank you very much for your attention!
The All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationInexact Data Sparse Boundary Element Tearing and Interconnecting Methods
Inexact Data Sparse Boundary Element Tearing and Interconnecting Methods U. Langer 1,2, G. Of 3, O. Steinbach 4, W. Zulehner 2 1 Radon Institute for Computational and Applied Mathematics, Austrian Academy
More informationTechnische Universität Graz
Technische Universität Graz Inexact Data Sparse Boundary Element Tearing and Interconnecting Methods U. Langer, G. Of, O. Steinbach, W. Zulehner Berichte aus dem Institut für Mathematik D Numerik und Partielle
More informationApplication of Preconditioned Coupled FETI/BETI Solvers to 2D Magnetic Field Problems
Application of Preconditioned Coupled FETI/BETI Solvers to 2D Magnetic Field Problems U. Langer A. Pohoaţǎ O. Steinbach 27th September 2004 Institute of Computational Mathematics Johannes Kepler University
More informationCoupled FETI/BETI for Nonlinear Potential Problems
Coupled FETI/BETI for Nonlinear Potential Problems U. Langer 1 C. Pechstein 1 A. Pohoaţǎ 1 1 Institute of Computational Mathematics Johannes Kepler University Linz {ulanger,pechstein,pohoata}@numa.uni-linz.ac.at
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 informationAll-Floating Coupled Data-Sparse Boundary and Interface-Concentrated Finite Element Tearing and Interconnecting Methods
All-Floating Coupled Data-Sparse Boundary and Interface-Concentrated Finite Element Tearing and Interconnecting Methods Ulrich Langer 1,2 and Clemens Pechstein 2 1 Institute of Computational Mathematics,
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 informationTechnische Universität Graz
Technische Universität Graz Robust boundary element domain decomposition solvers in acoustics O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2009/9 Technische Universität
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 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 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 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 informationAdaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers
Adaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers Jan Mandel University of Colorado at Denver Bedřich Sousedík Czech Technical University
More informationTechnische Universität Graz
Technische Universität Graz The All-Floating Boundary Element Tearing and Interconnecting Method G. Of, O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/3 Technische Universität
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 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 informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationScalable BETI for Variational Inequalities
Scalable BETI for Variational Inequalities Jiří Bouchala, Zdeněk Dostál and Marie Sadowská Department of Applied Mathematics, Faculty of Electrical Engineering and Computer Science, VŠB-Technical University
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 informationFETI Methods for the Simulation of Biological Tissues
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz FETI Methods for the Simulation of Biological Tissues Ch. Augustin O. Steinbach
More informationNumerical Approximation Methods for Elliptic Boundary Value Problems
Numerical Approximation Methods for Elliptic Boundary Value Problems Olaf Steinbach Numerical Approximation Methods for Elliptic Boundary Value Problems Finite and Boundary Elements Olaf Steinbach Institute
More informationExtending the theory for domain decomposition algorithms to less regular subdomains
Extending the theory for domain decomposition algorithms to less regular subdomains Olof Widlund Courant Institute of Mathematical Sciences New York University http://www.cs.nyu.edu/cs/faculty/widlund/
More informationSubstructuring for multiscale problems
Substructuring for multiscale problems Clemens Pechstein Johannes Kepler University Linz (A) jointly with Rob Scheichl Marcus Sarkis Clark Dohrmann DD 21, Rennes, June 2012 Outline 1 Introduction 2 Weighted
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
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 informationNonoverlapping Domain Decomposition Methods with Simplified Coarse Spaces for Solving Three-dimensional Elliptic Problems
Nonoverlapping Domain Decomposition Methods with Simplified Coarse Spaces for Solving Three-dimensional Elliptic Problems Qiya Hu 1, Shi Shu 2 and Junxian Wang 3 Abstract In this paper we propose a substructuring
More informationSelecting Constraints in Dual-Primal FETI Methods for Elasticity in Three Dimensions
Selecting Constraints in Dual-Primal FETI Methods for Elasticity in Three Dimensions Axel Klawonn 1 and Olof B. Widlund 2 1 Universität Duisburg-Essen, Campus Essen, Fachbereich Mathematik, (http://www.uni-essen.de/ingmath/axel.klawonn/)
More informationCONVERGENCE ANALYSIS OF A BALANCING DOMAIN DECOMPOSITION METHOD FOR SOLVING A CLASS OF INDEFINITE LINEAR SYSTEMS
CONVERGENCE ANALYSIS OF A BALANCING DOMAIN DECOMPOSITION METHOD FOR SOLVING A CLASS OF INDEFINITE LINEAR SYSTEMS JING LI AND XUEMIN TU Abstract A variant of balancing domain decomposition method by constraints
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 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 informationConvergence analysis of a balancing domain decomposition method for solving a class of indefinite linear systems
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 6 [Version: 2002/09/18 v1.02] Convergence analysis of a balancing domain decomposition method for solving a class of indefinite
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationParallel scalability of a FETI DP mortar method for problems with discontinuous coefficients
Parallel scalability of a FETI DP mortar method for problems with discontinuous coefficients Nina Dokeva and Wlodek Proskurowski University of Southern California, Department of Mathematics Los Angeles,
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 information18. Balancing Neumann-Neumann for (In)Compressible Linear Elasticity and (Generalized) Stokes Parallel Implementation
Fourteenth nternational Conference on Domain Decomposition Methods Editors: smael Herrera, David E Keyes, Olof B Widlund, Robert Yates c 23 DDMorg 18 Balancing Neumann-Neumann for (n)compressible Linear
More informationGroupe de Travail Méthodes Numériques
Groupe de Travail Méthodes Numériques Laboratoire Jacques-Louis Lions, Paris 6 Monday, May 11th, 2009 Domain Decomposition with Lagrange Multipliers: A continuous framework for a FETI-DP+Mortar method
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
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 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 informationMultispace and Multilevel BDDC
Multispace and Multilevel BDDC Jan Mandel Bedřich Sousedík Clark R. Dohrmann February 11, 2018 arxiv:0712.3977v2 [math.na] 21 Jan 2008 Abstract BDDC method is the most advanced method from the Balancing
More informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationOn the Use of Inexact Subdomain Solvers for BDDC Algorithms
On the Use of Inexact Subdomain Solvers for BDDC Algorithms Jing Li a, and Olof B. Widlund b,1 a Department of Mathematical Sciences, Kent State University, Kent, OH, 44242-0001 b Courant Institute of
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
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 informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More informationSome Domain Decomposition Methods for Discontinuous Coefficients
Some Domain Decomposition Methods for Discontinuous Coefficients Marcus Sarkis WPI RICAM-Linz, October 31, 2011 Marcus Sarkis (WPI) Domain Decomposition Methods RICAM-2011 1 / 38 Outline Discretizations
More informationON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000
ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Abstract In the Dual-Primal FETI method, introduced by Farhat et al [5], the domain is decomposed into
More informationGRUPO DE GEOFÍSICA MATEMÁTICA Y COMPUTACIONAL MEMORIA Nº 8
GRUPO DE GEOFÍSICA MATEMÁTICA Y COMPUTACIONAL MEMORIA Nº 8 MÉXICO 2014 MEMORIAS GGMC INSTITUTO DE GEOFÍSICA UNAM MIEMBROS DEL GGMC Dr. Ismael Herrera Revilla iherrerarevilla@gmail.com Dr. Luis Miguel de
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 informationCapacitance Matrix Method
Capacitance Matrix Method Marcus Sarkis New England Numerical Analysis Day at WPI, 2019 Thanks to Maksymilian Dryja, University of Warsaw 1 Outline 1 Timeline 2 Capacitance Matrix Method-CMM 3 Applications
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 informationAn Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03024-7 An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors Michel Lesoinne
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationIncomplete Cholesky preconditioners that exploit the low-rank property
anapov@ulb.ac.be ; http://homepages.ulb.ac.be/ anapov/ 1 / 35 Incomplete Cholesky preconditioners that exploit the low-rank property (theory and practice) Artem Napov Service de Métrologie Nucléaire, Université
More informationA DECOMPOSITION RESULT FOR BIHARMONIC PROBLEMS AND THE HELLAN-HERRMANN-JOHNSON METHOD
Electronic ransactions on Numerical Analysis. Volume 45, pp. 257 282, 2016. Copyright c 2016,. ISSN 1068 9613. ENA A DECOMPOSIION RESUL FOR BIHARMONIC PROBLEMS AND HE HELLAN-HERRMANN-JOHNSON MEHOD WOLFGANG
More informationFETI-DPH: A DUAL-PRIMAL DOMAIN DECOMPOSITION METHOD FOR ACOUSTIC SCATTERING
Journal of Computational Acoustics, c IMACS FETI-DPH: A DUAL-PRIMAL DOMAIN DECOMPOSITION METHOD FOR ACOUSTIC SCATTERING Charbel Farhat, Philip Avery and Radek Tezaur Department of Mechanical Engineering
More informationDomain Decomposition solvers (FETI)
Domain Decomposition solvers (FETI) a random walk in history and some current trends Daniel J. Rixen Technische Universität München Institute of Applied Mechanics www.amm.mw.tum.de rixen@tum.de 8-10 October
More informationA smooth variant of the fictitious domain approach
A smooth variant of the fictitious domain approach J. Haslinger, T. Kozubek, R. Kučera Department of umerical Mathematics, Charles University, Prague Department of Applied Mathematics, VŠB-TU, Ostrava
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 informationPreconditioning of Saddle Point Systems by Substructuring and a Penalty Approach
Preconditioning of Saddle Point Systems by Substructuring and a Penalty Approach Clark R. Dohrmann 1 Sandia National Laboratories, crdohrm@sandia.gov. Sandia is a multiprogram laboratory operated by Sandia
More informationSOME PRACTICAL ASPECTS OF PARALLEL ADAPTIVE BDDC METHOD
Conference Applications of Mathematics 2012 in honor of the 60th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 2012 SOME PRACTICAL ASPECTS OF PARALLEL ADAPTIVE BDDC METHOD Jakub Šístek1,2,
More informationfor three dimensional problems are often more complicated than the quite simple constructions that work well for problems in the plane; see [23] for a
DUAL-PRIMAL FETI METHODS FOR THREE-DIMENSIONAL ELLIPTIC PROBLEMS WITH HETEROGENEOUS COEFFICIENTS AEL KLAWONN Λ, OLOF B. WIDLUND y, AND MAKSYMILIAN DRYJA z Abstract. In this paper, certain iterative substructuring
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 informationOVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY TR
OVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY MINGCHAO CAI, LUCA F. PAVARINO, AND OLOF B. WIDLUND TR2014-969 Abstract. Low order finite element discretizations of the linear
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationTwo new enriched multiscale coarse spaces for the Additive Average Schwarz method
346 Two new enriched multiscale coarse spaces for the Additive Average Schwarz method Leszek Marcinkowski 1 and Talal Rahman 2 1 Introduction We propose additive Schwarz methods with spectrally enriched
More informationTechnische Universität Graz
Technische Universität Graz Schur complement preconditioners for the biharmonic Dirichlet boundary value problem L. John, O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2013/4
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
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 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 informationProjected Schur Complement Method for Solving Non-Symmetric Saddle-Point Systems (Arising from Fictitious Domain Approaches)
Projected Schur Complement Method for Solving Non-Symmetric Saddle-Point Systems (Arising from Fictitious Domain Approaches) Jaroslav Haslinger, Charles University, Prague Tomáš Kozubek, VŠB TU Ostrava
More informationUncertainty analysis of large-scale systems using domain decomposition
Center for Turbulence Research Annual Research Briefs 2007 143 Uncertainty analysis of large-scale systems using domain decomposition By D. Ghosh, C. Farhat AND P. Avery 1. Motivation and objectives A
More informationParallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients
Parallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients Nina Dokeva and Wlodek Proskurowski Department of Mathematics, University of Southern California, Los Angeles,
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 informationMultipréconditionnement adaptatif pour les méthodes de décomposition de domaine. Nicole Spillane (CNRS, CMAP, École Polytechnique)
Multipréconditionnement adaptatif pour les méthodes de décomposition de domaine Nicole Spillane (CNRS, CMAP, École Polytechnique) C. Bovet (ONERA), P. Gosselet (ENS Cachan), A. Parret Fréaud (SafranTech),
More informationConferences Fictitious domain methods for numerical realization of unilateral problems. J.H. (jointly with T. Kozubek and R. Kučera) p.
p. /2 Conferences 28 Fictitious domain methods for numerical realization of unilateral problems J.H. (jointly with T. Kozubek and R. Kučera) Lyon 23-24.6. 28 p. 2/2 Outline Fictitious domain methods -
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 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 informationParallel Scalable Iterative Substructuring: Robust Exact and Inexact FETI-DP Methods with Applications to Elasticity
Parallel Scalable Iterative Substructuring: Robust Exact and Inexact FETI-DP Methods with Applications to Elasticity Oliver Rheinbach geboren in Hilden Fachbereich Mathematik Universität Duisburg-Essen
More informationBoundary Value Problems - Solving 3-D Finite-Difference problems Jacob White
Introduction to Simulation - Lecture 2 Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Reminder about
More informationBath Institute For Complex Systems
BICS Bath Institute for Complex Systems Energy Minimizing Coarse Spaces for Two-level Schwarz Methods for Multiscale PDEs J. Van lent, R. Scheichl and I.G. Graham Bath Institute For Complex Systems Preprint
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationEFFICIENT NUMERICAL TREATMENT OF HIGH-CONTRAST DIFFUSION PROBLEMS
EFFICIENT NUMERICAL TREATMENT OF HIGH-CONTRAST DIFFUSION PROBLEMS A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial Fulfillment of the Requirements
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 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 informationPANM 17. Marta Čertíková; Jakub Šístek; Pavel Burda Different approaches to interface weights in the BDDC method in 3D
PANM 17 Marta Čertíková; Jakub Šístek; Pavel Burda Different approaches to interface weights in the BDDC method in 3D In: Jan Chleboun and Petr Přikryl and Karel Segeth and Jakub Šístek and Tomáš Vejchodský
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 informationSOME NONOVERLAPPING DOMAIN DECOMPOSITION METHODS
SIAM REV. c 1998 Society for Industrial and Applied Mathematics Vol. 40, No. 4, pp. 857 914, December 1998 004 SOME NONOVERLAPPING DOMAIN DECOMPOSITION METHODS JINCHAO XU AND JUN ZOU Abstract. The purpose
More informationGRUPO DE GEOFÍSICA MATEMÁTICA Y COMPUTACIONAL MEMORIA Nº 3
GRUPO DE GEOFÍSICA MATEMÁTICA Y COMPUTACIONAL MEMORIA Nº 3 MÉXICO 2014 MEMORIAS GGMC INSTITUTO DE GEOFÍSICA UNAM MIEMBROS DEL GGMC Dr. Ismael Herrera Revilla iherrerarevilla@gmail.com Dr. Luis Miguel de
More informationIntroduction to Domain Decomposition Methods. Alfio Quarteroni
. Aachen, EU Regional School Series February 18, 2013 Introduction to Domain Decomposition Methods Alfio Quarteroni Chair of Modeling and Scientific Computing Mathematics Institute of Computational Science
More informationAn Iterative Domain Decomposition Method for the Solution of a Class of Indefinite Problems in Computational Structural Dynamics
An Iterative Domain Decomposition Method for the Solution of a Class of Indefinite Problems in Computational Structural Dynamics Charbel Farhat a, and Jing Li a a Department of Aerospace Engineering 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 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 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 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 informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
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 information