Adaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers
|
|
- Benjamin Cummings
- 5 years ago
- Views:
Transcription
1 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 and University of Colorado at Denver Parts of this presentation are based on joint work with Clark R. Dohrmann, Sandia National Laboratories Radek Tezaur, Stanford University DD16, January 2005 Supported by Sandia National Laboratories and National Science Foundation Grant CNS
2 Synopsis 1. BDDC and FETI-DP formulation and condition number bound following Jan Mandel, Clark Dohrmann, and Radek Tezaur, An Algebraic Theory for Primal and Dual Substructuring Methods by Constraints, Applied Numerical Mathematics, in print (6th IMACS Conference on Iterative Methods, Denver, 2003) 2. Numerical computation of the condition number bound from a global eigenvalue problem: in experiments, reasonably sharp 3. Estimate of the condition number bound from local (face based) eigenvalue problems: in experiments, reasonably accurate 4. Optimal selection of face averages from local eigenvectors to decrease the condition number estimate = practical and efficient adaptive algorithm
3 The FETI-DP and BDDC Methods Both the BDDC and FETI-DP methods are build from same components. The methods use different algebraic algorithms. FETI-DP (Farhat, Lesoinne, Pierson 2000): substructuring method of the FETI family (Farhat, Roux, 1992), based on Lagrange Multipliers. Convergence proof Mandel and Tezaur (2001). BDDC (Dohrmann 2002): method from the Balancing Domain Decomposition (BDD) family (Mandel 1993), based on Additive Schwarz framework of Neumann Neumann type (Dryja, Widlund 1995); convergence analysis Mandel and Dohrmann For corner constraints only, Cros the eigenvalues of the preconditioned operators of FETI-DP and BDDC are identical (Mandel, Dohrmann, Tezaur 2004)
4 Development of Adaptive FETI-DP and BDDC Algorithms Continue algebraic analysis as long as possible before switching to calculus (FEM, functional analysis) arguments, abstract from a specific FEM framework to widen the scope of the methods Algebraic bounds on the condition number allow to select components of the methods adaptively Select method components adaptively to force condition number estimate under a specified value: control number of iterations but each iteration gets more expensive. Future: minimize estimated total computational work.
5 Substructuring and Reduction to Interfaces K i is the stiffness matrix for substructure i, symmetric positive problem in decomposed form 1 2 vt Kv v T f min, v = v 1. v N K = K continuity of dofs between substructures partition the dofs in each subdomain i into internal and interface (boundary) and eliminate interior dofs K i = [ K ii i Ki ib T Ki ib K bb i ], v i = [ ] v i i = decomposed problem reduced to interfaces v b i, f i = [ f i i 1 2 wt Sw w T g min, S = diag(s i ), S i = Ki bb Ki ib +continuity of dofs between substructures f b i ]. K N T K ii 1 i K ib i
6 Enforcing Intersubdomain Continuity primal methods (BDD,N-N): U is the space of global dofs R i : U W i restriction to substructure i, R = R 1. R N continuity of dofs between substructures: w = Ru for some u U dual methods (FETI): continuity of dofs between [ substructures: ] Bw = B = [B 1,..., B N ] = : W Λ.... by construction: R i R T i = I, range R = null B
7 Constraints and Coarse Degrees of Freedom Choose matrix Q T P that selects coarse dofs: u c = Q T Pu (e.g. values of global dofs u at corners, averages on sides) Define R ci : all constraint values values that can be nonzero on substructure i, define C i = R ci Q T P RT i : u c = 0 C i w i = 0 i Assume the generalized coarse dofs define interpolation: w U u c i : C i R i w = R ci u c a coarse dof can only involve nodes adjacent to the same set of substructures Define subspace of vectors with coarse dofs continuous across substructures: W = {w W : u c i : C i w i = R ci u c } Dual approach (FETI): construct Q D so that w W Q T D Bw = 0
8 Intersubdomain Averaging and Weight matrices primary diagonal weight matrices D Pi : W i W i, D P = diag(d Pi ) decomposition of unity: R T D P R = I purpose: average between subdomains to get global dofs: u = R T D P w dual weight matrices D Di : Λ Λ, defined by d α ij = dα j - d α j : diagonal entry of D Pj associated with the same global dof α - d α ij : diagonal entry of D Di for Ω i and Ω j and the same global dof α define B D = [D D1 B 1,... D DN B N ], then B T D is a generalized inverse of B: BBT D B = B associated projections are complementary: B T D B + RRT D P = I Rixen, Farhat, Tezaur, Mandel 1998, Rixen, Farhat 1999, Klawonn, Widlund, Dryja 2001, 2002, Fragakis, Papadrakakis 2003 same identities hold also for other versions of the operators
9 FETI-DP Approach enforce equality of coarse dofs directly : Q T D Bw = 0 enforce the rest of the constraints Bw = 0 (or all) by multipliers λ } = saddle point problem: min w W max λ L(w, λ) = max λ min w W L(w, λ) L(w, λ) = 1 2 wt Sw w T f + w T B T λ W = { w : Q T Bw = 0 } functions with coarse dofs same between substructures = dual problem: F(λ) λ = 0, F(λ) = min w W preconditioner M = B D SBD T L(w, λ) Q T D Bw = 0 enforces continuity - of values across crosspoints (Farhat, Lesoinne, Le Tallec, Pierson, Rixen 2001) - also of averages across edges and faces (for 3D) (Farhat, Lesoinne, Pierson 2000, Pierson thesis 2000, Klawonn, Widlund, Dryja 2002)
10 BDDC: Balancing Domain Decomposition Based on Constraints The system operator of the BDDC method is the assembled Schur complement The preconditioner P is defined by A = R T SR. where u c is the solution of SPD coarse problems and z is the solution of Pr = R T D P (Ψu c + z) (1) Ψ T SΨu c = Ψ T DP T Rr (2) Sz + C T µ = D T P Rr Cz = 0 (independent substructure problems because of block structure of S and C) Here, the coarse basis functions Ψ are defined by energy minimization [ ] [ ] [ ] S C T Ψ 0 =. (4) C 0 Λ R c (3)
11 Algebraic Bound on Condition Number Asssume that S is positive definite on W : there is enough constraints to eliminate mechanisms between substructures in the whole model. Theorem. κ BDDC = κ FETI DP ω = sup w W B T D Bw 2 S w 2 S Theorem. Eigenvalues or the preconditioned operator are the same except for zero eigenvalue in FETI-DP and different multiplicity of eigenvalue one. Zero eigenvalues of FETI-DP are from redundant constraints in B. All other eigenvalues are 1.
12 Related Work: Theoretical Asymptotic Condition Number Bounds The algebraic bound + standard substructuring arguments (Klawonn, Widlund, Dryja 2002) = ω C ( 1 + log 2 H ) h However, the arguments are quite subtle, especially for elasticity. A similar asymptotic bound (with h, H independent constants) in terms of BD T B was given by Klawonn and Widlund 2001 for some other versions of FETI and BDD. We carry on our algebra a bit longer so our algebraic bound here has no undetermined constants. It also helps that the analysis is simpler in the case of FETI-DP and BDDC. This introduction was all taken from Mandel, Dohrmann, Tezaur The following material is new.
13 Computation of Condition Number Bound Lemma. If Π is an orthogonal projection, a 0, and ω 0, then ΠTΠu = ωπsπu u rangeπ ΠTΠu = ω (ΠSΠ + a (I Π)) u Theorem. a > 0, κ BDDC = κ FETI DP ω = sup w W B T D Bw 2 S w 2 S = sup w w T ΠB T B D SB T D BΠw w T (ΠSΠ + a (I Π)) w where w T (ΠSΠ + a (I Π)) w >0 and Π is the orthogonal projection onto W: Π = I [ C T ] [ CC 0 T ] 1 [ ] R c C Rc T 0 0 Approximate max eigenvalue and eigenvector can be computed easily e.g. by LOBPCG (Knyazev, 2001). Or just run the iterations and compute the extreme eigenvalues from the Lanczos sequence generated by PCG.
14 Heuristic Local Condition Number Estimates We estimate the convergence bound by the maximum of the same bounds computed from two adjacent substructures i, j at a time: w ij = [ ] wi w j, S ij = [ Si 0 0 S j ], C ij = [ Ci 0 0 C j etc. Define B ij as the submatrix of [ B i B j ] consisting of all rows that have exactly one +1 and one 1, define W ij accordingly. We propose where ω ij = sup w ij Wij ω ω = max ij ω ij BDij T B 2 ijw ij S ij = sup 2 w ij w S ij Wij ij w T ij BT ij B DijS ij B T Dij B ijw ij ] w T ij S ijw ij We assume that there is already enough constraints between the two substructures to eliminate mechanisms (i.e., relative motion).
15 Computation of Local Condition Number Estimates We have restricted all computations on two substructures i and j with a common face and assumed that there is already enough constraints to eliminate relative motion between the two substructures. Theorem. wij T ω ij = sup Π ijbij T B DijS ij BDij T B ijπ ij w ij [ w ij Πij S ij Π ij + a ( )] I Π ij wij w T ij Problem: both matrices are still typically singular because of rigid body modes that move substructures i and j as a whole. This makes computing the eigenvalues harder. Need to add to the right hand side matrix a projection on the complement of the nullspace of [ Π ij S ij Π ij + a ( I Π ij )].
16 Computation of Local Condition Number Estimates (cont.) Find matrices Z i, Z j that generate a superspace of rigid body modes of the two substructures: null S i range Z i, null S j range Z j, Z ij = [ ] Zi (the coarse basis functions may be used, Z i = Ψ i because span of coarse basis functions contains rigid body modes, but there are too many of them; RBM often available directly) Z j. construct generators of nullspace by solving a smaller eigenvalue problem range Q = null ( Π ij S ij Π ij + a ( I Π ij )) = null ( Z T ij [ Πij S ij Π ij + a ( I Π ij )] Zij )
17 Computation of Local Condition Number Estimates (cont.) Theorem. where ω ij = sup w ij w T ij X ijw ij w T ij Y ijw ij X ij = Π ij B T ij B DijS ij B T Dij B ijπ ij, Y ij = ( ( Π ij Πij S ij Π ij + a ( )) I Π ij Πij + a ( )) I Π ij > 0 ( and Π ij = I QQ T ) is the orthogonal projection onto range ( Π ij S ij Π ij + a ( I Π ij )). This is a symmetric generalized matrix eigenvalue problem with one of the matrices positive definite, easy to solve.
18 Optimal Constraints to Decrease the Rayleigh Quotient The condition number estimate is the maximum of Rayleigh quotient over space W ij. We can add constraints to make W ij smaller. Lemma (well known) Let X, Y be symmetric positive semidefinite operators in n-dimensional space V, Y positive definite, and Xu k = λ k Y u k, λ 1 λ 2... λ n all eigenvalues and eigenvectors. Then for any V k V such that dim V k = k, with equality if and only if max {J(u) : v, Xu = 0 v V k } λ k+1 V k = span {u 1,..., u k } k constraints from eigenvectors eliminates the k largest eigenvalues This will allow us to add constraints that decrease the estimates ω ij optimally.
19 Correspondence Between Dual and Primal Definitions of Coarse DOFs Dual view: W ij = { w ij : Q T Dij B ijw ij = 0 } is the space defined by averages of continuity constraints: w ij orthogonal to something Primal view: w ij W ij the values of the coarse dofs q P R T i w i are same for both substructures Entries of B ij are +1 for i and 1 for j = the value of the coarse dof is the common value on both sides: B ij w ij = I ij w i I ji w j (I ij selects from w i the DOFs on interface between substructures i and j) w ij W ij Q T Dij I ijw i = Q T Dij I jiw j Each column of q D of Q Dij defines a coarse degree of freedom which can be written using the corresponding column q P of Q P as q T P RT i w i = q T D I ijw i R i and I ij are 0-1 q P is a scattering of the entries of the q D.
20 Construction of Optimal Additional Face Constraints Recall that ω ij is the maximal eigenvalue λ ij,k of B T ij B DijS ij B T Dij B iju ij,k = λ ij,k S ij u ij,k Given target τ add columns q D of Q D so that u T ij,k BT ij B DijS ij BDij T B ij w ij = 0 }{{} qd T λ ij,k > τ = guaranteed that the heuristic condition estimate ω ij τ The additional (primal) coarse degrees of freedom are then defined by q T P RT i w i = q T D I ijw i
21 Adaptive Coarse DOFs for Plane Elasticity λ = 1 µ = 2 H/h Ndof τ Nc ω κ it < < < < < < < < Ndof: number of degrees of freedom τ condition number target: coarse dofs added for all λ ij > τ; corners only if not specified ω = max ω ij : heuristic condition estimate κ approximate condition number estimate from PCG it: number of iterations for stopping tolerance 10 8
22 Adaptive Constraints for H/h = 16 Eigenvalues ω ij on interfaces sub i sub j λ ij,1 λ ij,2 λ ij,3 λ ij,4 λ ij, Number of added constraints τ = Numbering of substructures:
23 Adaptive Coarse DOFs for Almost Incompressible Elasticity λ = 1000 µ = 2 H/h Ndof τ Nc ω κ it < < < < < < < < NA < < < NA < < < Ndof: number of degrees of freedom τ condition number target: coarse dofs added for all λ ij > τ; corners only if not specified ω = max ω ij : heuristic condition estimate κ approximate condition number estimate from PCG it: number of iterations for stopping tolerance 10 8
24 Related Work: Selection of Constraints for Elasticity Lesoinne 2001: selection of corner constraints to prevent global mechanisms (= guarantee that S is positive definite on W). Examples that S positive definite on {u : Cu = 0} (each substructure sufficiently constrained by coarse dofs to prevent rigid body modes) is not enough. Klawonn, Widlund, 2004: how to select a minimal number of corner constraints to prevent global mechanisms and to get rigorous C ( 1 + log 2 H h ) condition number bounds. Here we assume that S is positive definite on W and even more: there is already enough constraints to prevent mechanism between any two substructures with a common face. Dohrmann 2004: additional constraint (flux) for almost incompressible elasticity; Jing Li 2002: edge average constraints on all velocity components for Stokes problem. Here we find the additional constraints automatically.
25 Conclusion Demonstrated that condition number is estimated accurately and reliably Developed efficient control of condition number by additional constraints: a robust method Shown that additional constraints concentrated on the problem spot: low computational complexity Works also for almost incompressible elasticity. No other information needed besides the matrix in substructure form and what are the corners. Applicable to both FETI-DP and BDDC. Future: tests and development on real problems; adaptive algorithms to minimize total work
Multilevel 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 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 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 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 informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 31, pp. 384-402, 2008. Copyright 2008,. ISSN 1068-9613. ETNA ON THE EQUIVALENCE OF PRIMAL AND DUAL SUBSTRUCTURING PRECONDITIONERS BEDŘICH SOUSEDÍK
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationParallel Sums and Adaptive BDDC Deluxe
249 Parallel Sums and Adaptive BDDC Deluxe Olof B. Widlund 1 and Juan G. Calvo 2 1 Introduction There has recently been a considerable activity in developing adaptive methods for the selection of primal
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 informationINSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC. Adaptive-Multilevel BDDC and its parallel implementation
my Cze of Scie ch R n epu ces blic INSTITUTE of MATHEMATICS Aca de ACADEMY of SCIENCES of the CZECH REPUBLIC INSTITU TE THEMATICS MA of Adaptive-Multilevel BDDC and its parallel implementation Bedřich
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 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 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 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 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 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 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 informationTR THREE-LEVEL BDDC IN THREE DIMENSIONS
R2005-862 REE-LEVEL BDDC IN REE DIMENSIONS XUEMIN U Abstract BDDC methods are nonoverlapping iterative substructuring domain decomposition methods for the solution of large sparse linear algebraic systems
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 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 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 informationON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000 Revised April 2000
ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Revised April 2000 Abstract In te Dual-Primal FETI metod, introduced by Farat et al [5], te domain is
More informationAN ADAPTIVE CHOICE OF PRIMAL CONSTRAINTS FOR BDDC DOMAIN DECOMPOSITION ALGORITHMS
lectronic Transactions on Numerical Analysis. Volume 45, pp. 524 544, 2016. Copyright c 2016,. ISSN 1068 9613. TNA AN ADAPTIV CHOIC O PRIMAL CONSTRAINTS OR BDDC DOMAIN DCOMPOSITION ALGORITHMS JUAN G. CALVO
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 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 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 informationParallel Implementation of BDDC for Mixed-Hybrid Formulation of Flow in Porous Media
Parallel Implementation of BDDC for Mixed-Hybrid Formulation of Flow in Porous Media Jakub Šístek1 joint work with Jan Březina 2 and Bedřich Sousedík 3 1 Institute of Mathematics of the AS CR Nečas Center
More informationA Neumann-Dirichlet Preconditioner for FETI-DP 2 Method for Mortar Discretization of a Fourth Order 3 Problems in 2D 4 UNCORRECTED PROOF
1 A Neumann-Dirichlet Preconditioner for FETI-DP 2 Method for Mortar Discretization of a Fourth Order 3 Problems in 2D 4 Leszek Marcinkowski * 5 Faculty of Mathematics, University of Warsaw, Banacha 2,
More informationInexact Data-Sparse BETI Methods by Ulrich Langer. (joint talk with G. Of, O. Steinbach and W. Zulehner)
Inexact Data-Sparse BETI Methods by Ulrich Langer (joint talk with G. Of, O. Steinbach and W. Zulehner) Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences http://www.ricam.oeaw.ac.at
More 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 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 informationAN ANALYSIS OF A FETI DP ALGORITHM ON IRREGULAR SUBDOMAINS IN THE PLANE TR
AN ANALYSIS OF A FETI DP ALGORITHM ON IRREGULAR SUBDOMAINS IN THE PLANE TR27 889 AXEL KLAWONN, OLIVER RHEINBACH, AND OLOF B. WIDLUND Abstract. In the theory for domain decomposition algorithms of the iterative
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 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 information(a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? Solution: dim N(A) 1, since rank(a) 3. Ax =
. (5 points) (a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? dim N(A), since rank(a) 3. (b) If we also know that Ax = has no solution, what do we know about the rank of A? C(A)
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 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 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 informationcfl Jing Li All rights reserved, 22
NYU-CS-TR83 Dual-Primal FETI Methods for Stationary Stokes and Navier-Stokes Equations by Jing Li A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationA FETI-DP method for the parallel iterative solution of indefinite and complex-valued solid and shell vibration problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.1282 A FETI-DP method for the parallel iterative solution
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 informationFETI-DP for Elasticity with Almost Incompressible 2 Material Components 3 UNCORRECTED PROOF. Sabrina Gippert, Axel Klawonn, and Oliver Rheinbach 4
1 FETI-DP for Elasticity with Almost Incompressible 2 Material Components 3 Sabrina Gippert, Axel Klawonn, and Oliver Rheinbach 4 Lehrstuhl für Numerische Mathematik, Fakultät für Mathematik, Universität
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
More informationAdaptive Coarse Spaces and Multiple Search Directions: Tools for Robust Domain Decomposition Algorithms
Adaptive Coarse Spaces and Multiple Search Directions: Tools for Robust Domain Decomposition Algorithms Nicole Spillane Center for Mathematical Modelling at the Universidad de Chile in Santiago. July 9th,
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 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 informationDOMAIN DECOMPOSITION
ČESKÉ VYSOKÉ UČENÍ TECHNICKÉ V PRAZE Fakulta stavební Doktorský studijní program: STAVEBNÍ INŽENÝRSTVÍ Studijní obor: Matematika ve stavebním inženýrství Ing. Bedřich Sousedík SROVNÁNÍ NĚKTERÝCH METOD
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 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 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 informationOn the Robustness and Prospects of Adaptive BDDC Methods for Finite Element Discretizations of Elliptic PDEs with High-Contrast Coefficients
On the Robustness and Prospects of Adaptive BDDC Methods for Finite Element Discretizations of Elliptic PDEs with High-Contrast Coefficients Stefano Zampini Computer, Electrical and Mathematical Sciences
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationMultilevel Methods for Eigenspace Computations in Structural Dynamics
Multilevel Methods for Eigenspace Computations in Structural Dynamics Ulrich Hetmaniuk & Rich Lehoucq Sandia National Laboratories, Computational Math and Algorithms, Albuquerque, NM Joint work with Peter
More informationA hybrid domain decomposition method based on one-level FETI and BDDC algorithms
A hybrid domain decomposition method based on one-level FET and BDDC algorithms Jungho Lee Computer Science and Mathematics Division Oak Ridge National Laboratory July 26, 2 Abstract A three-level domain
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 informationPreconditioned Locally Minimal Residual Method for Computing Interior Eigenpairs of Symmetric Operators
Preconditioned Locally Minimal Residual Method for Computing Interior Eigenpairs of Symmetric Operators Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University
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 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 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 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 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 informationAccomodating Irregular Subdomains in Domain Decomposition Theory
Accomodating Irregular Subdomains in Domain Decomposition Theory Olof B. Widlund 1 Courant Institute, 251 Mercer Street, New York, NY 10012, USA widlund@cims.nyu.edu Summary. In the theory for domain decomposition
More informationA Unified Framework for Adaptive BDDC
www.oeaw.ac.at A Unified Framework for Adaptive BDDC C. Pechstein, C.R. Dohrmann RICAM-Report 2016-20 www.ricam.oeaw.ac.at A UNIFIED FRAMEWORK FOR ADAPTIVE BDDC CLEMENS PECHSTEIN 1 AND CLARK R. DOHRMANN
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 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 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 informationFETI domain decomposition method to solution of contact problems with large displacements
FETI domain decomposition method to solution of contact problems with large displacements Vít Vondrák 1, Zdeněk Dostál 1, Jiří Dobiáš 2, and Svatopluk Pták 2 1 Dept. of Appl. Math., Technical University
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 informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More informationLINEAR ALGEBRA 1, 2012-I PARTIAL EXAM 3 SOLUTIONS TO PRACTICE PROBLEMS
LINEAR ALGEBRA, -I PARTIAL EXAM SOLUTIONS TO PRACTICE PROBLEMS Problem (a) For each of the two matrices below, (i) determine whether it is diagonalizable, (ii) determine whether it is orthogonally diagonalizable,
More informationAN ADAPTIVE CHOICE OF PRIMAL CONSTRAINTS FOR BDDC DOMAIN DECOMPOSITION ALGORITHMS TR
AN ADAPTIVE CHOICE OF PRIMAL CONSTRAINTS FOR BDDC DOMAIN DECOMPOSITION ALGORITHMS JUAN G. CALVO AND OLOF B. WIDLUND TR2015-979 Abstract. An adaptive choice for primal spaces, based on parallel sums, is
More informationRESEARCH ARTICLE. A strategy of finding an initial active set for inequality constrained quadratic programming problems
Optimization Methods and Software Vol. 00, No. 00, July 200, 8 RESEARCH ARTICLE A strategy of finding an initial active set for inequality constrained quadratic programming problems Jungho Lee Computer
More informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
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 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 informationGRAPH PARTITIONING WITH MATRIX COEFFICIENTS FOR SYMMETRIC POSITIVE DEFINITE LINEAR SYSTEMS
GRAPH PARTITIONING WITH MATRIX COEFFICIENTS FOR SYMMETRIC POSITIVE DEFINITE LINEAR SYSTEMS EUGENE VECHARYNSKI, YOUSEF SAAD, AND MASHA SOSONKINA Abstract. Prior to the parallel solution of a large linear
More informationA BDDC ALGORITHM FOR RAVIART-THOMAS VECTOR FIELDS TR
A BDDC ALGORITHM FOR RAVIART-THOMAS VECTOR FIELDS DUK-SOON OH, OLOF B WIDLUND, AND CLARK R DOHRMANN TR013-951 Abstract A BDDC preconditioner is defined by a coarse component, expressed in terms of primal
More informationA Numerical Algorithm for Block-Diagonal Decomposition of Matrix -Algebras, Part II: General Algorithm
A Numerical Algorithm for Block-Diagonal Decomposition of Matrix -Algebras, Part II: General Algorithm Takanori Maehara and Kazuo Murota May 2008 / May 2009 Abstract An algorithm is proposed for finding
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 informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More information18.06SC Final Exam Solutions
18.06SC Final Exam Solutions 1 (4+7=11 pts.) Suppose A is 3 by 4, and Ax = 0 has exactly 2 special solutions: 1 2 x 1 = 1 and x 2 = 1 1 0 0 1 (a) Remembering that A is 3 by 4, find its row reduced echelon
More informationConceptual Questions for Review
Conceptual Questions for Review Chapter 1 1.1 Which vectors are linear combinations of v = (3, 1) and w = (4, 3)? 1.2 Compare the dot product of v = (3, 1) and w = (4, 3) to the product of their lengths.
More informationROBUST BDDC PRECONDITIONERS FOR REISSNER-MINDLIN PLATE BENDING PROBLEMS AND MITC ELEMENTS
ROBUST BDDC PRECONDITIONERS FOR REISSNER-MINDLIN PLATE BENDING PROBLEMS AND MITC ELEMENTS L. BEIRÃO DA VEIGA, C. CHINOSI, C. LOVADINA, AND L. F. PAVARINO Abstract. A Balancing Domain Decomposition Method
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 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 informationMa/CS 6b Class 23: Eigenvalues in Regular Graphs
Ma/CS 6b Class 3: Eigenvalues in Regular Graphs By Adam Sheffer Recall: The Spectrum of a Graph Consider a graph G = V, E and let A be the adjacency matrix of G. The eigenvalues of G are the eigenvalues
More information