Total Overlapping Schwarz Preconditioners for Elliptic PDEs ( a ) ( b )
|
|
- Erica Jones
- 5 years ago
- Views:
Transcription
1 1 Seminaire, PARIS XIII, 05/12/08 Total Overlapping Schwarz Preconditioners for Elliptic PDEs ( a ) ( b ) F. Ben Belgacem (LMAC-UTCompiègne) N. Gmati (LAMSIN-ENITunis) F. Jelassi (LMAC-UTCompiègne) a Main ideas originated in the Ph.D Thesis of F. Jelassi, b The Helmholtz equation is considered as elliptic.
2 2 Schwarz Method (1870) : A Brief Description Find u such that u = 0 in ] 1, 1[, u( 1) = u(+1) = 0, 1 α β + 1 Find u 2m+1 such that (u 2m+1 ) = 0, ] 1, β[, (u 2m+1 )( 1) = 0, (u 2m+1 )(β) = (u 2m )(β) Find u 2m+2 such that (u 2m+1 ) = 0, ]α, +1[, (u 2m+2 )(+1) = 0, (u 2m+2 )(α) = (u 2m+1 )(α)
3 3 Convergence This is a Mathematical Proof) CL Dirichlet x = 1 CL Dirichlet x = CL Dirichlet x = +1 CL Dirichlet x = CL Dirichlet x = 1 CL Dirichlet x = Linear convergence
4 4 Motivations Construction of Sub-structuring Preconditioners MSM, ASM, FETI,... (A Huge Literature). Coupling Methods, Models, or Programs. Matching Meshes,... (Great Number of papers) Parallelism (A lot of work)
5 5 Bibliography H. A. Schwarz, 1870 maximum principle K. Miller, 1965 efficient computational tool A. M. Matsokin, S. V. Nepomnyaschikh, 1985 P. E. Bjørstad, O. B. Widlund, 1986 P.-L. Lions, 1988 variationnal approach!!! M. Dryja, O. B. Widlund, 1989, 1990 (Additive Schwarz ) B. Smith, P. Bjørstad, W. Gropp, 1996 C. Farhat, F. X. Roux, 1991 (FETI) A. Quarteroni, Valli, 1999 A. Toselli, O. Widlund, 2005 BenDali, Boubendir, Fares, Collino, 2003
6 The Total Overlapping Schwarz Method (TOSM) 6
7 7 Bibliography (2) : post 2000 Liu, Jin, Adaptive absorbing BC/FEM (2001), Preconditioners for FEM/Boundary Integral (2002) Maury, Ismaïl, The Fat boundary Method (2001) Pironeau, Lions, Brezzi, Apoung-Kamga, Hecht Chimera Method, Numerical Zooming (2005) Glowinski, He, Lozinsky, Rappaz, Wagner, Multi-Domain (Multi-Scale) Method, FEM with patches (2004)
8 8 A Cracked Domain (I) γ + ω γ Σ Ω γ γ Find ϕ such that ϕ = f in Ω γ, ϕ = 0 on Γ, n ϕ = 0 in γ + γ
9 9 A Cracked Domain (II) Variational Principle (= ) Find ϕ such that ϕ ψ dx = fψ dx, Ω γ Ω γ ψ. Finite Element Approximation (= ) Find ϕ h such that A h ϕ h = f h. The stiffness matrix A h sees the cracks (ϕ h suffers from jumps) (= ) The condition number degrades rapidly with the number and of the geometrical form of the cracks. (= ) Necessity of Domain Decomposition devices in the computations.
10 10 The Mile Stone of (TOSM) : Decompose the Problem One sub-problem (= ) Find η in Ω γ such that η = f in Ω γ, [η] = [χ] on γ, [ n η] = 0 on γ, η = 0 on Γ. The other (= ) Find χ in ω γ such that χ = f in ω γ, χ = η on Σ n χ = 0 on γ + γ. We have that η = ϕ and χ = ϕ ωγ Look at (η χ) to be convinced.
11 11 Uncouple the Transmission (= ) (TOSM) Assume (χ m, η m ) is known, η m+1 is solution of η m+1 = f in Ω γ, [η m+1 ] = [χ m ] on γ, [ n η m+1 ] = 0 on γ, η m+1 = 0 on Γ, and χ m+1 satisfies χ m+1 = f in ω γ, χ m+1 = η m+1 on Σ n χ m+1 = 0 on γ + γ.
12 12 If η m+1 wants not to see the crack at all γ = any regular curve embracing the crack γ and located inside ω γ. f = is the trivial extension of f inside γ. Replace η m+1 by η m+1 solution of and χ m+1 satisfies η m+1 = f in Ω \ γ, [ η m+1 ] = χ m on γ, [ n η m+1 ] = n χ m on γ, η m+1 = 0 on Γ. χ m+1 = f in ω γ, χ m+1 = η m+1 on Σ n χ m+1 = 0 on γ + γ.
13 13 Advantages Mesh for η that do not see the crack (= ) The Stiffness Matrix is related to the safe domain Ω (= ) Better conditioned and Easier to solve. Deconnection between multiple cracks (= ) Solve local subproblems on χ at the vicinity of the individual cracks (small subproblems). Uncoupled two-scale computations : we may use Standard FEM for η and XFEM( a ) for χ (Blind to each other). a extended FEM
14 14 Higher Advantages for Exterior Problems (I) Ω γ = R d γ (= ) η m+1 the solution of η m+1 = f in R d γ, [η m+1 ] = [χ m ] on γ, [ n η m+1 ] = 0 on γ, η m+1 = 0 at, can be obtained thanks to an Integral Representation (G : Green function). η m+1 (x) = fg(x, ) dy + ([χ m ]) n G(x, ) dγ = T(χ m ). R d γ
15 15 Higher Advantages for exterior problems (II) (= ) Solve iteratively χ m+1 = f in ω γ, χ m+1 = (η m+1 ) = T γ (χ m ) on Σ n χ m+1 = 0 on γ + γ. (= ) Uncoupling FEM/Boundary Integral. T γ is a Regular Kernel Operator Allows Programs dedicated to Bounded Domains to handle Exterior Domains without considering Singular Kernels.
16 16 In what is it different of Schwarz Method (for cracks)? Assume (χ m, η m ) is known, η m+1 is solution of and χ m+1 satisfies η m+1 = f in Ω γ, η m+1 = χ m on γ + γ, η m+1 = 0 on Γ, χ m+1 = f in ω γ, χ m+1 = η m+1 on Σ n χ m+1 = 0 on γ + γ. Bounded domains : Both approaches, when written on γ, are economically equivalent if they have the same convergence rate (not the case). Unbounded domains : The Total Overlapping approach is pretty more economical.
17 17 Perforated Domains γ Ω γ Σ Ω γ ω γ Find ϕ such that ϕ = f in Ω γ, ϕ = 0 on Γ, n ϕ = 0 in γ
18 18 Decompose the Problem Extend f by zero in Ω γ = f One sub-problem (= ) Find η Ω γ Ω γ such that η = f in Ω γ Ω γ, [η] = χ on γ, [ n η] = 0 on γ, η = 0 on Γ. The other (= ) Find χ in ω γ such that We have that χ = f in ω γ, χ = η on Σ, n χ = 0 on γ. η Ωγ = ϕ Ωγ, η Ω γ = 0 and χ = ϕ ωγ
19 19 Toal Overlapping Schwarz Method Find η m+1 Ω γ Ω γ such that Find χ m+1 in ω γ such that η m+1 = f in Ω γ Ω γ, [η m+1 ] = χ m on γ, [ n η m+1 ] = 0 on γ, η m+1 = 0 on Γ, χ m+1 = f in ω γ, χ m+1 = η m+1 on Σ, n χ m+1 = 0 on γ. The Fat Boundary Method (Bertrand Maury 2001) (=) A Penalized (TO)-Schwarz Method ( =) Wrong Transmission Cdts.
20 20 Differences with the Classical Schwarz Method Find η m+1 Ω γ such that Find χ m+1 in ω γ such that η m+1 = f in Ω γ, η m+1 = χ m on γ, η m+1 = 0 on Γ, χ m+1 = f in ω γ, χ m+1 = η m+1 on Σ, n χ m+1 = 0 on γ. Classical Schwarz method (= ) Compute η m+1 in the perforated domain. (TO)-Schwarz method (= ) Compute η m+1 in the filled domain ( ).
21 21 Convergence? Classical Schwarz γ,σ, Γ : Circles with radius 1,, and 1 <. The Schwarz iterates are given by (χ m, η m )(r, θ) = (χ m p, ηp m )(r) cos(pθ), p=0 η m p = a m r p, r χ m p = b m r p + d m r p, 1 r Linear Convergence (= ) µ p = 2 2p 1 + 2p
22 22 Convergence? (TOSM) (χ m, η m )(r, θ) = (χ m p, ηp m )(r) cos(pθ), p=0 ηp m = a i rp, r 1 a e r p, r 1 χ m p = b m r p + d m r p Linear Convergence (= ) µ p = 2p (= ) Twice faster than the Cl. Sch p 2,5 2,5 2, ,5 1 1,5 1 1,5 1 0,5 0,5 0, ,5 1 1,5 2 2,5 3 3, ,5 1 1,5 2 2,5 3 3,5-0,
23 23 Convergence of (TOSM) Denote by e m = (ϕ χ m ), τ m = (ϕ η m ). Théo. 1 (fj, 2005) (TOSM) converges with a linear rate e m L2 (ω γ ) + τ m L2 (Ω γ ) be am m. a > 0 and increases when the distance between Σ and γ grows.
24 24 Numerical Examples: Exterior Problems MELINA fortran 90 (D. Martin, IRMAR) Implementation of specific procedures (F. Jelassi, LMAC)
25 25 Rankine s Oval (1) Potential flow : Source-sink dipole + Uniform flow f(z) = 5π 4 z + log ( z 0.2 z ) Rankine oval: Streamlines in ω γ
26 26 Rankine Oval (2) Ecart entre deux iteres successifs (l,w) = (0.6,0.2) (l,w) = (1,0.5) (l,w) = (2,1) (l,w) = (3,1.5) Ecart entre deux iteres successifs (l,w) = (0.6,0.2) (l,w) = (1,0.5) (l,w) = (2,1) (l,w) = (3,1.5) Nombre d iterations Nombre d iterations Consecutive iterates gap : Potential (left) and Stream function (right).
27 27 A Wing Profile (1)
28 28 A Wing Profile (2) Wing profile : NACA0012. Angle of Attack = π 12. Erreurs relatives 10-1 Carre Ellipse Spline Erreurs relatives 10-1 Carre Ellipse Spline Nombre d iteration Nombre d iteration Consecutive iterates gap: l 2 (left) et l (right).
29 29 Electrical Dipole : Non-connected Σ (1) The local problems are dis-connecte. They are solved individually
30 30 Electrical Dipole : Non-connected Σ (2) Circular Conductors radius = 0.5. Distance = 4. Electric Charge = πσ of opposite signs. Non-Connected Calculus. Erreur relative R = 0.7 R = 1 R = 1.5 Erreur relative R = 0.7 R = 1 R = Iterations Iterations Gaps (u sc u co ) : l 2 (left) et l (right).
31 31 Electric Quadrupole : Non-Connected Σ Quadrupole. Symmetric et Dis-symmetric Cases.
32 32 Electric Quadrupole in 3D : Non-Connected Electrostatic Conductors (Spheres of radius 1 U). Uniformly charged on their skins and Distant of 4 U from each other Erreurs relatives Cas symetrique Cas symetrique Cas non-symetrique Cas non-symetrique Nombre d iteration Gap between Consecutive Iterates with l 2 -norm (empty symbols) and l (full symbols).
33 33 Flat Quadrupole (1) (Cracked Domain) Equipotentials of the non-symmetric flat quadrupole. Charges concentration are higher along the sides where equipotentials are more dense
34 34 Flat Quadrupole (2) 10 0 Residual Schwarz GMRES Iterations (TO)-Schwarz method and Preconditioned GMRES residual curves for the cracked domain: the potential for the non symmetric quadrupole.
35 35 Electric Dipole within Heterogeneous Media GMRES Computations (left) and (TO)-Schwarz Method (right).
36 36 To overcome this Weakness (= ) Krylov Sub-Spaces Solvers Take back the coupled problem (= ) Find η λ Ω γ Ω γ such that η λ = f in Ω γ Ω γ, [η λ ] = χ λ on γ, [ n η λ ] = 0 on γ, η λ = 0 on Γ. Find χ λ in ω γ such that χ λ = f in ω γ, χ λ = λ( = η λ ) on Σ, n χ λ = 0 on γ. The condensed problem to solve λ (η λ ) Σ = 0 on Σ.
37 37 Static Condensation (I) The Constant Part of the Condensed Equation Define χ such that χ = f in ω γ, χ = 0 on Σ n χ = 0 on γ + γ. Define η that satisfies η = f in Ω γ, [η] = [χ] on γ, [ n η] = 0 on γ, η = 0 on Γ. Set G = (η) Σ. No real coupling between both sub-problems.
38 38 Static Condensation (II) The linear Part Define χ λ such that χ λ = 0 in ω γ, χ λ = λ on Σ n χ λ = 0 on γ + γ, Define η λ that satisfies η λ = 0 in Ω γ, [η λ ] = [χ λ ] on γ, [ n η λ ] = 0 on γ, η λ = 0 on Γ. Set B Σ (λ) = (η λ ) Σ
39 39 (TOS) Preconditioner The suitable λ(= χ) solves the equation χ λ = η λ on Σ, (I + B Σ )(λ) = G ( ) λ (η λ ) Σ = (η) Σ. Richardson method (= Jacobi Algorithm) λ m+1 + B Σ (λ m ) = G. ( ) (TOSM) with χ m = χ λ m + χ. B Σ is a small perturbation of I because all the singular values of it are smaller than one ( =) (TOSM) converges. B Σ is a highly Regular Kernel operator (= ) the singular values of B Σ decrease exponentially fast toward zero. (= ) Krylov Subspaces Method will dramatically accelerate the convergence.
40 40 Convergence of Krylov Methods Denote by r m = G (I + B Σ )(λ m ) = (χ λ η λ ) m. Théo. 2 (fbb-ng-fj, 2008) Krylov Methods converge with a quadratic rate in two dimensions that is r m L2 (Σ) be am2 m. The convergence rate in three dimensions is three halves since we have that r m L2 (Σ) be am m m. In both estimates, a > 0 and increases when the distance between Σ and γ grows.
41 41 Circular Quadrupole (1) (= ) Perforated Domain) The non-symmetric circular quadrupole..
42 42 Circular Quadrupole (2) Equipotentials of the symmetric circular quadrupole (in ω γ ).
43 43 Circular Quadrupole (3) Residual Schwarz GMRES Residual Schwarz GMRES Iterations Iterations Convergence curves for the quadrupole with circular conductors. The symmetric is in the left diagram and non-symmetric one in the right diagram. Full symbols are for connected ω γ while empty symbols are for non-connected ω γ. Circular symbols correspond to the TOSM and squared symbols are for the precondit. GMRES.
44 44 Back to the Flat Quadrupole (2b) Residual Schwarz GMRES Iterations (TO)-Schwarz method and Preconditioned GMRES residual curves for the cracked domain: the potential for the non symmetric quadrupole.
45 Wave Diffraction 45
46 46 Helmholtz Problem Ω : (Gray) Bounded Obstacle = Ω : medium surrounding Ω. Incident wave Scattered wave Find u such that u + κ 2 u = 0 in Ω, n u = g on γ, ( r iκ)u = e iκr O( 1 ) at, r2 Sommerfeld Condition at = Existence and Uniqueness
47 47 Exact Non-Reflecting Boundary Condition (FEM/BEM) Find u such that u + κ 2 u = 0 in ω γ, n u = g on γ, ( n + β)u = T β K (u) on Σ. The operator T β K is defined by (K β = ( n + β)g, Regular Green Function ) T β K (u)(x) = ( (u) n K β (x, y) + (g)k β (x, y) ) dγ y, on Σ γ Variational Formulation = Find u (A β + T β )u = G in ω γ. Exterior and Bounded problems are equivalent for β C
48 48 (TOSM) Find v m+1 such that Find u m+1 such that I(β) 0 is compulsory. v m+1 + κ 2 v m+1 = 0 in Ω Ω, [v m+1 ] = [u m ] on γ, [ n v m+1 ] = g on γ, ( r iκ)v m+1 = e iκr O( 1 ) at, r2 u m+1 + κ 2 u m+1 = 0 in ω γ, n u m+1 = g on Γ, ( n + β)u m+1 = ( n + β)v m+1 on Σ.
49 49 (TOSM) = (keep u m and eliminate v m ) Find u such that u m+1 + κ 2 u m+1 = 0 in ω γ, n u m+1 = g on Γ, ( n + β)u m+1 = T β K (um ) on Σ. Uncoupled variational Formulation (= ) u m known, find u m+1 such that A β u m+1 + T β u m = G. (1) I(β) 0 = A β is invertible. Exact non-reflection at = IneExact Non-Reflection at Σ + Iterations. (TOSM) = Liu and Jin Method (IEEE Trans. Antennas Propagat. 2001) = Preconditioned Richardson Method applied to (1).
50 50 Analytical Investigation Γ, Σ : Circles with radius 1, et R 1 R. Sound Hard Obstacle (a Neumann Condition) g(θ) = p Z g p exp(ipθ). Iterates of (TOSM) are provided by u m (r, θ) = p Z ϕ m (r)exp(ipθ), Attenuation Cœfficients of (TOSM) 1 τ p (β) = 1 H(1) p (κ) J p (κ) (κj p βj p )(κr) (κh p (1) βh p (1) )(κr).
51 51 Attenuation Cœfficients τ p (β) (1) Choice of β = iκ, Wave Length : λ = 2π κ, Thickness of the Computational Domain = eλ, 0,4 0,5 τ 0,3 κ = 30 κ = 10 κ = 1 τ 0,4 κ = 30 κ = 10 κ = 1 0,3 0,2 0,2 0,1 0, p p Attenuation Cœfficients τ p for e = 1 and e = 0.25.
52 52 Attenuation Cœfficients τ p (β) (2) τ 1 0,8 κ = 30 κ = 10 κ = 1 partie imaginaire 1 0,5 0 κ = 1 κ = 10 κ = 30 0,6-0,5 0,4 0, p -1,5-1 -0,5 0 0,5 1 1,5 partie reelle Attenuation Cœfficients for e = 0.01 (left). Eigenvalue of the Operator of Iteration (right). Relaxation of Richardson with a Small Parameter = It works.
53 53 (TOS) Preconditioner for Krylov Methods (1) Find u such that (I + (A β ) 1 T β ) u = (A β ) 1 G. by GMRES or Bi-Conjugate Gradient. Théo. 3 Assume that the fictitious boundary Σ is piecewise analytic. Then, for any real number a > 0, there exists c = c(a) > 0, such that a(m log m) r m L2 (Σ) ce m. Rem. 1 There is no reason why we have Less sharp convergence than for the Poisson Problem (= ) Hope to improve our Mathematical proofs.
54 Numerical Examples (Melina) 54
55 55 Non-Convex Guide (1) Slit Length = 0.4. Wave length = 1 or 0.1. Thickness of the computational domain = The truncating Σ and the boundary of the scatterer Γ have similar shapes. To the right is a zoom on the mesh we used.
56 56 Non-Convex Guide (2) 8 8 Im (eigenv) 6 Im (eigenv) 6 4 λ = 1 4 λ = 1 λ = λ = Re (eigenv) Re (eigenv) The distribution of the eigenvalues of non preconditioned system (to the left). The effect of the (TOS) preconditioner (to the right).
57 57 Non-Convex Guide (3) Residuals λ = 1 λ = 0.1 Residuals λ = 1 λ = Iterations Iterations Convergence curves of the GMRES with no preconditioner (to the left) With the Schwarz preconditioner (to the right). After capturing the first non-converging modes the Precond. GMRES exhibits a quadratic convergence.
58 58 Non-Convex 3D Guide Residuals λ = 1 λ = 0.5 λ = Iterations The shape of the scatterer (left panel) and (TOS) precondtioned GMRES convergence curves for an incident planar wave with wavelengths λ = 1, 0.5 and 0.2. The thickness of the computational domain e = 0.05.
59 59 Conclusion (TOSM) improves the Schwarz. Perspectives (1) More Computations (3D). (2) Many possible Extensions. (3) More Mathematics.
Capacitance 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 informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
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 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 informationLecture 18 Classical Iterative Methods
Lecture 18 Classical Iterative Methods MIT 18.335J / 6.337J Introduction to Numerical Methods Per-Olof Persson November 14, 2006 1 Iterative Methods for Linear Systems Direct methods for solving Ax = b,
More informationProblems of Corner Singularities
Stuttgart 98 Problems of Corner Singularities Monique DAUGE Institut de Recherche MAthématique de Rennes Problems of Corner Singularities 1 Vertex and edge singularities For a polyhedral domain Ω and a
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 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 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 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 informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
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 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 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 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 informationCIMPA Summer School on Current Research in Finite Element Methods
CIMPA Summer School on Current Research in Finite Element Methods Local zooming techniques for the finite element method Alexei Lozinski Laboratoire de mathématiques de Besançon Université de Franche-Comté,
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 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 informationDomain Decomposition Methods for Mortar Finite Elements
Domain Decomposition Methods for Mortar Finite Elements Dan Stefanica Courant Institute of Mathematical Sciences New York University September 1999 A dissertation in the Department of Mathematics Submitted
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 information1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl)
DDM Preprint Editors: editor1, editor2, editor3, editor4 c DDM.org 1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Bernhard Hientzsch 1 1. Introduction.
More informationMaster Thesis Literature Study Presentation
Master Thesis Literature Study Presentation Delft University of Technology The Faculty of Electrical Engineering, Mathematics and Computer Science January 29, 2010 Plaxis Introduction Plaxis Finite Element
More informationFast iterative BEM for high-frequency scattering
Fast iterative BEM for high-frequency scattering elastodynamic problems In 3D. Marion Darbas LAMFA UMR CNRS 7352, Amiens, France S. Chaillat (POEMS, CNRS-INRIA-ENSTA) and F. Le Louër (LMAC, UTC) Conference
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 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 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 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 Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems
Contemporary Mathematics Volume 8, 998 B 0-88-0988--03030- A Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems Daoqi Yang. Introduction
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationAdditive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems
Additive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems Atle Loneland 1, Leszek Marcinkowski 2, and Talal Rahman 3 1 Introduction In this paper
More informationSharp condition number estimates for the symmetric 2-Lagrange multiplier method
Sharp condition number estimates for the symmetric -Lagrange multiplier method Stephen W. Drury and Sébastien Loisel Abstract Domain decomposition methods are used to find the numerical solution of large
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 informationNatural Boundary Integral Method and Its Applications
Natural Boundary Integral Method and Its Applications By De-hao Yu State Key Laboratory of Scientific and Engineering Computing Institute of Computational Mathematics and Scientific/Engineering Computing
More informationAn Introduction to Domain Decomposition Methods: algorithms, theory and parallel implementation
An Introduction to Domain Decomposition Methods: algorithms, theory and parallel implementation Victorita Dolean, Pierre Jolivet, Frédéric ataf To cite this version: Victorita Dolean, Pierre Jolivet, Frédéric
More informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
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 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 information9.1 Preconditioned Krylov Subspace Methods
Chapter 9 PRECONDITIONING 9.1 Preconditioned Krylov Subspace Methods 9.2 Preconditioned Conjugate Gradient 9.3 Preconditioned Generalized Minimal Residual 9.4 Relaxation Method Preconditioners 9.5 Incomplete
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 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 informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationTwo-level multiplicative domain decomposition algorithm for recovering the Lamé coefficient in biological tissues
Two-level multiplicative domain decomposition algorithm for recovering the Lamé coefficient in biological tissues Si Liu and Xiao-Chuan Cai Department of Applied Mathematics, University of Colorado at
More informationAn eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More informationSpectral analysis of complex shifted-laplace preconditioners for the Helmholtz equation
Spectral analysis of complex shifted-laplace preconditioners for the Helmholtz equation C. Vuik, Y.A. Erlangga, M.B. van Gijzen, and C.W. Oosterlee Delft Institute of Applied Mathematics c.vuik@tudelft.nl
More information33 RASHO: A Restricted Additive Schwarz Preconditioner with Harmonic Overlap
Thirteenth International Conference on Domain Decomposition ethods Editors: N. Debit,.Garbey, R. Hoppe, J. Périaux, D. Keyes, Y. Kuznetsov c 001 DD.org 33 RASHO: A Restricted Additive Schwarz Preconditioner
More informationThe 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 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 informationAn additive average Schwarz method for the plate bending problem
J. Numer. Math., Vol. 10, No. 2, pp. 109 125 (2002) c VSP 2002 Prepared using jnm.sty [Version: 02.02.2002 v1.2] An additive average Schwarz method for the plate bending problem X. Feng and T. Rahman Abstract
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
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 11: CMSC 878R/AMSC698R. Iterative Methods An introduction. Outline. Inverse, LU decomposition, Cholesky, SVD, etc.
Lecture 11: CMSC 878R/AMSC698R Iterative Methods An introduction Outline Direct Solution of Linear Systems Inverse, LU decomposition, Cholesky, SVD, etc. Iterative methods for linear systems Why? Matrix
More informationSolving Ill-Posed Cauchy Problems in Three Space Dimensions using Krylov Methods
Solving Ill-Posed Cauchy Problems in Three Space Dimensions using Krylov Methods Lars Eldén Department of Mathematics Linköping University, Sweden Joint work with Valeria Simoncini February 21 Lars Eldén
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 informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
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 informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
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 informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
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 informationXIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods.
Contemporary Mathematics Volume 00, 0000 Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA Abstract. In this paper, we study several overlapping
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Iteration basics Notes for 2016-11-07 An iterative solver for Ax = b is produces a sequence of approximations x (k) x. We always stop after finitely many steps, based on some convergence criterion, e.g.
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 informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
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 informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 11, pp. 1-24, 2000. Copyright 2000,. ISSN 1068-9613. ETNA NEUMANN NEUMANN METHODS FOR VECTOR FIELD PROBLEMS ANDREA TOSELLI Abstract. In this paper,
More informationA THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS
A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor
More informationError analysis and fast solvers for high-frequency scattering problems
Error analysis and fast solvers for high-frequency scattering problems I.G. Graham (University of Bath) Woudschoten October 2014 High freq. problem for the Helmholtz equation Given an object Ω R d, with
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 informationNested Krylov methods for shifted linear systems
Nested Krylov methods for shifted linear systems M. Baumann, and M. B. van Gizen Email: M.M.Baumann@tudelft.nl Delft Institute of Applied Mathematics Delft University of Technology Delft, The Netherlands
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 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 informationFactorization method in inverse
Title: Name: Affil./Addr.: Factorization method in inverse scattering Armin Lechleiter University of Bremen Zentrum für Technomathematik Bibliothekstr. 1 28359 Bremen Germany Phone: +49 (421) 218-63891
More informationImproving the Robustness of a Surface Integral Formulation for Wideband Impendance Extraction of 3D Structures
Improving the Robustness of a Surface Integral Formulation for Wideband Impendance Extraction of 3D Structures Zhenhai Zhu, Jingfang Huang, Ben Song, Jacob White Department of Electrical Engineering and
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 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 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 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 informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
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 informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
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 informationDomain decomposition for the Jacobi-Davidson method: practical strategies
Chapter 4 Domain decomposition for the Jacobi-Davidson method: practical strategies Abstract The Jacobi-Davidson method is an iterative method for the computation of solutions of large eigenvalue problems.
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 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 informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationA HIGH-ORDER ACCURATE ACCELERATED DIRECT SOLVER FOR ACOUSTIC SCATTERING FROM SURFACES
A HIGH-ORDER ACCURATE ACCELERATED DIRECT SOLVER FOR ACOUSTIC SCATTERING FROM SURFACES JAMES BREMER,, ADRIANNA GILLMAN, AND PER-GUNNAR MARTINSSON Abstract. We describe an accelerated direct solver for the
More informationNewton additive and multiplicative Schwarz iterative methods
IMA Journal of Numerical Analysis (2008) 28, 143 161 doi:10.1093/imanum/drm015 Advance Access publication on July 16, 2007 Newton additive and multiplicative Schwarz iterative methods JOSEP ARNAL, VIOLETA
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationThe Deflation Accelerated Schwarz Method for CFD
The Deflation Accelerated Schwarz Method for CFD J. Verkaik 1, C. Vuik 2,, B.D. Paarhuis 1, and A. Twerda 1 1 TNO Science and Industry, Stieltjesweg 1, P.O. Box 155, 2600 AD Delft, The Netherlands 2 Delft
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 informationCoupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations
Coupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations P. Bochev, J. Cheung, M. D Elia, A. Frishknecht, K. Kim, M. Parks, M. Perego CM4 summer school, Stanford,
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
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 informationAn Accelerated Block-Parallel Newton Method via Overlapped Partitioning
An Accelerated Block-Parallel Newton Method via Overlapped Partitioning Yurong Chen Lab. of Parallel Computing, Institute of Software, CAS (http://www.rdcps.ac.cn/~ychen/english.htm) Summary. This paper
More informationInverse Obstacle Scattering
, Göttingen AIP 2011, Pre-Conference Workshop Texas A&M University, May 2011 Scattering theory Scattering theory is concerned with the effects that obstacles and inhomogenities have on the propagation
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 informationMultigrid and Domain Decomposition Methods for Electrostatics Problems
Multigrid and Domain Decomposition Methods for Electrostatics Problems Michael Holst and Faisal Saied Abstract. We consider multigrid and domain decomposition methods for the numerical solution of electrostatics
More information