A hybrid domain decomposition method based on one-level FETI and BDDC algorithms
|
|
- Sandra Blair
- 5 years ago
- Views:
Transcription
1 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 decomposition is considered. Bodies in contact with each other are divided into subdomains, which in turn are the union of elements. Using an approach based purely on FET (finite element tearing and interconnecting) algorithms with only Lagrange multipliers as unknowns, which has been developed by the engineering community, does not lead to a scalable algorithm with respect to the number of subdomains in each body. nstead, we consider a new method based on the saddle point formulation of the FET methods with both displacement vectors and Lagrange multipliers as unknowns. The resulting system is solved with a block-diagonal preconditioner which combines the one-level FET and the BDDC (balancing domain decomposition by constraints) methods. We show that this new method is scalable with respect to the number of subdomains. A model contact problem is solved by a nonlinear algorithm which combines the new method and a primal-dual active set method. Keywords domain decomposition, scalable algorithms, FET, BDDC, contact problems, primaldual active set method AMS Subject Classification 65F, 65K, 65N55 ntroduction We consider contact problems without friction. n 8, 2, we considered the FET-FET method, which is a domain decomposition method for a linearized contact problem. We showed that the FET-FET method, which has been used in the engineering community 2,, has a condition number estimate which grows linearly with the number of subdomains, or processors. This paper is a sequel to 8; we introduce a scalable alternative to the FET-FET method, which we call a hybrid method. This method combines the one-level FET and the BDDC methods. n this paper, we assume the use of an active set method to solve our contact problem. n each step of an active set method, the active set is updated, and a imization problem on the current active set is approximately solved, until a desired accuracy is achieved. Thus, an active set method requires (nonlinear) outer iterations in which the active set is updated and (linear) inner iterations in which a imization problem is solved. We use the primal-dual active set strategy viewed as a semismooth Newton method, proposed and analyzed in 3, 4,, 9, to detere the outer iterations. This paper is organized as follows. n Section 2, we review the one-level FET, FET-DP (dualprimal FET), and BDDC methods. We introduce a nonlinear model problem and briefly describe the This author s work was supported in part by the U.S. Department of Energy under contracts DE-FG2-6ER2578 and DE-FC2-ER25482 and in part by National Science Foundation grant DMS This submission was sponsored by a contractor of the United States Government under contract DE-AC5-OR22725 with the United States Department of Energy. The United States Government retains, and the publisher, by accepting this submission for publication, acknowledges that the United States Government retains, a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this submission, or allow others to do so, for United States Government purposes. An earlier version of the material contained in this paper already appeared in the Ph.D. thesis of the author, see 2.
2 FET-FET method in Section 3. We introduce the hybrid method and provide an eigenvalue analysis of its preconditioned operator in Sections 4 and 5. n Section 6 we present numerical results which confirm the scalability of the hybrid method. n Section 7, we solve a model problem using a combination of a primal-dual active set method and our hybrid method. 2 Building blocks of the Hybrid method n this section, we briefly review some popular domain decomposition methods with nonoverlapping subdomains. 2. A model problem and notation We consider a second-order scalar elliptic problem on a bounded domain Ω R n,n = 2,3. We denote the boundary of Ω by Ω, and assume that homogeneous Dirichlet boundary conditions are imposed on Ω D Ω, which is a subset of Ω with a positive measure. Let Ω N := Ω \ Ω D be its complement. The corresponding Sobolev space in which the solution will be found is H (Ω, Ω D ) := {v H (Ω) : u = on Ω D }. We find u H (Ω, Ω D ) such that where a(u,v) = f(v), v H (Ω, Ω D ), (2.) a(u,v) := Ω ρ(x) u v, f(v) = Note that (2.) is equivalent to the following imization problem: Ω fv. (2.2) a(u,u) f(u). (2.3) u H (Ω, Ω D 2 ) We decompose Ω into N nonoverlapping subdomains Ω i,i =,,N, each of which is the union of shape-regular elements with the finite element nodes on the boundaries of neighboring subdomains matching across the interface := ( N i= Ω i) \ Ω D. is the union of faces, edges and vertices in three dimensions: faces, regarded as open subsets of, are shared by two subdomains. Edges, regarded as open subsets of the boundaries of the faces, are shared by more than two subdomains. Vertices are endpoints of edges. edges and vertices in two dimensions: edges, regarded as open subsets of, are shared by two subdomains. Vertices, as in three dimensions, are endpoints of edges. We assume that ρ(x) = ρ i ρ >, x Ω i,i =,,N. We also introduce the corresponding set of interface nodes h := ( N i= Ω i,h) \ Ω h, where Ω h and Ω i,h are the sets of finite element nodes on Ω and Ω i, respectively. We also define local bilinear forms and linear functionals, a (u,v) := ρ(x) u v, Ω i f (v) = fv. Ω i (2.4) n the rest of this section, we discuss the choice of the space of finite element functions in onelevel FET, FET-DP, and BDDC methods. We denote a standard finite element space of continuous, piecewise linear functions on Ω i by W. We will always assume that these functions vanish on Ω D. Each W is decomposed into a subdomain interior part W W = W W. and a subdomain interface part W : We denote the associated product spaces by W := N i= W,W := N, and W := N. The functions in W and W are in general discontinuous across the interface, whereas the finite element solutions are continuous across the interface. Therefore we introduce continuous subspaces i= W i= W 2
3 of W and W by Ŵ and Ŵ, respectively. For the FET-DP and BDDC methods, we will also need a subspace W W, intermediate between W and Ŵ, which consists of finite element functions which satisfy certain continuity constraints. The corresponding interface space is denoted by W. We introduce the following decomposition of W : ( N ) W = W ŴΠ = ŴΠ, where ŴΠ is a subspace of continuous functions and W is a subspace of functions which are allowed to be discontinuous across the interface. More precisely, ŴΠ is spanned by subdomain vertex nodal basis functions, i.e., consists of functions which are nonzero only at subdomain vertices, in the twodimensional case. Accordingly, W W consist of functions which are zero at the vertices of the subdomain Ω i. n other words, W consists of functions that are continuous at subdomain vertices. n the three-dimensional case, such vertex constraints are not enough to ensure scalability and we need to enforce the edge averages to be continuous; for instance, see 6. See Figure for a depiction of W, W, and Ŵ in the two-dimensional case. For each subdomain Ω i,i =,,N, we assemble local stiffness matrices A and local load vectors f obtained by integrating appropriate expressions over individual subdomains. We also introduce scaling factors δ i (x) for each node x h Ω i,h,i =,,N: for γ /2, ), δ i (x) = ρ γ i j N x ρ γ j i= W, x Ω i,h h. Here, N x is the set of indices j of the subdomains such that x Ω j,h. Figure : W, W and Ŵ (a) W: One-Level (b) W: FET-DP (c) Ŵ: BDDC 2.2 The One-Level FET method n this subsection, we review the one-level FET method, following 26, Section 6.3. We use the finite element functions in the space W to discretize the imization problem (2.3). Since the functions in W are in general discontinuous across the interface, we need to enforce the continuity condition explicitly: u W 2 ut Au f T u, subject to Bu =, (2.5) 3
4 where A = A () A (N), f = f (). f (N) Bu = represents continuity constraints across the interface, where B = B (),B (2),,B (N) is a matrix consisting of elements,, such that Bu = if and only if all the values of u associated with more than one subdomain boundary coincide. The columns of B which correspond to the interior nodes of Ω i are zero. Thus, B = B when the interior degrees of freedom are ordered first. We call B a jump operator. ntroducing a vector of Lagrange multipliers λ to enforce the continuity constraint Bu =, we obtain the following Karush-Kuhn-Tucker (KKT) system: Find (u, λ) W range(b), such that. Au + B T λ = f Bu =. (2.6) λ is unique only up to an additive element of ker(b T ). The space of Lagrange multipliers, U, is therefore chosen as range(b). Eliating the interior unknowns in each subdomain, we obtain the following: Find (u,λ) W range(b ), such that where S = S () g = g (). g (N) S (N), Su + B T λ = g B u =, (2.7), S = A A g = f A A A A T,i =,,N, f,i =,,N and B = B (),B(2),,B(N) is obtained by removing the zero columns of B that correspond to the interior nodes of individual subdomains, resulting in Bu = B u where B = B and u T = u T u T. n all FET methods, we reduce the KKT system (2.7) to an equation of λ alone, by solving the first equation of (2.7) for u. The matrices A in (2.6) and S in (2.7), however, are generally singular, when there are subdomains with boundaries which do not intersect the Dirichlet boundary Ω D. We call such subdomains floating. n such a case the solution of the first equation of (2.7) exists if and only if g B T λ range(s); this requirement leads to the introduction of a projection P, which will be introduced shortly. First, we introduce a matrix R such that range(r) = ker(s): R () R =, where R consists of the null vectors of S,i =,,N. Subdomains with nonsingular stiffness matrices do not contribute to the matrix R, i.e., R is an empty matrix if the subdomain Ω i intersects the Dirichlet boundary Ω D. We now solve the first equation of (2.7) for u : R (N) u = S (g B T λ) + Rα if g B T λ range(s) = ker(s) = range(r), (2.8) 4
5 where S is a pseudoinverse of S and α has to be detered. Substituting (2.8) into the second equation of (2.7), we obtain B S B T λ = B S g + B Rα. (2.9) We introduce the notation F := B S B T,d := B S g,g := B R,e := R T g and P := G(G T G) G T. Note that P is a projection operator with its range orthogonal to G. We apply this P to (2.9) to eliate the term with α and rewrite the orthogonality condition in (2.8) to obtain the following: { PFλ = Pd G T (2.) λ = e. We define the space V := {µ U : B T µ range(s)} = ker(g T ), which we call the space of admissible increments, following Chen and Mandel 7. The one-level FET method is a preconditioned conjugate gradient method applied to PFλ = Pd, λ λ + V (2.) where λ is chosen such that G T λ = e. Here, we only consider the Dirichlet preconditioner M D := B D, SBD, T, where B D, = B () D, BN D, is a scaled jump operator. B D, is obtained as follows: each nonzero entry of B is related to the Lagrange multiplier enforcing the continuity at a node x Ω i Ω j and is multiplied by δ j (x) to produce the corresponding entry of B D,. With this choice of preconditioner, the preconditioned operator of the one-level FET method has the following condition number bound: K C( + log(h/h)) 2, (2.2) where K denotes the condition number of the preconditioned operator in the appropriate subspace. For a proof of (2.2), see 24 or 26, Section 6.3. Thus the convergence rate of the one-level FET method depends only polylogarithmically on the number of degrees of freedom of a subdomain. 2.3 The FET-DP method n this subsection, we closely follow the notation of 2. For more details on various FET-DP methods, see, e.g., 2, 6,, 2, 26 and the references therein. n the FET-DP method, we use finite element functions in W to discretize (2.3). We first note that the local stiffness matrices A and the local load vectors f can be written as follows: A = A A A Π A T A A Π A T Π AT Π A ΠΠ, f = f f f Π, (2.3) where,, and Π indicate the index sets corresponding to the interior nodes, dual nodes, i.e., those of W, and primal nodes, i.e., those of W Π, respectively. We introduce the matrix Ã, which can be thought of as the restriction of A, defined for the functions in W, to the subspace W: Ã = A () A ()T A () A () Ã () Π Ã ()T Π Ã ()T Π A (N) A (N) Ã () Π Ã (N) Π A (N)T A (N) Ã (N) Π Ã (N)T Π Ã (N)T Π Ã ΠΠ. (2.4) 5
6 Here, and Ã Π = RT Π A Π, à ΠΠ = Ã Π = RT Π A Π, N i= R T Π A ΠΠ R Π, i =,,N, where R Π : ŴΠ W Π,i =,,N, is a restriction operator which extracts the relevant subdomain component belonging to W Π from a vector in ŴΠ. As in the one-level FET method, we introduce a vector of Lagrange multipliers and obtain the following saddle point problem: Find (u,λ) W range( B), such that Ãu + B T λ = f Bu =. (2.5) Again, B is a jump operator such that Bu =,u W if and only if the values of u associated with more than one subdomain coincide. Eliating the interior unknowns of each subdomain from the system (2.5), we obtain: Find (u,λ) W range( B ), such that S u + B T λ = g B u =. (2.6) S can also be regarded as the restriction of S, defined on W, to the subspace W : S = R T S R, where R : W W is a direct sum of restriction operators that extract the subdomain part belonging to W from a vector in W. The matrix Ã, and therefore also S, are nonsingular, so we can solve the first equation of (2.6) for u and substitute the resulting expression into the second equation of (2.6): B S B T λ = B S g. (2.7) The Dirichlet preconditioner used in the FET-DP algorithms to solve the equation (2.7) is B D, S BT D,. () (N) B D, = B D,,, B D, is obtained in exactly the same manner as B D, in Section 2.2. With this choice of preconditioner, the preconditioned operator for the FET-DP method also has the condition number bound (2.2). For a proof of this convergence bound for the two-dimensional case, see, e.g., 25. For three-dimensional scalar elliptic problems and linear elasticity problems, see 7 and 6, respectively. 2.4 The BDDC method n this subsection, we review the BDDC method, following 27. The discretized problem on the entire domain Ω is: Find (u,u ) (W,Ŵ), such that ( ) ( ) A A T u = A A u The equation (2.8) can be rewritten as A () R ()T A () A ()T R(N) T A (N) A (N) N i= A (N)T R T R () R (N) A 6 R ( f f u () Ị ). (2.8). u (N) u f ()... f (N) N i= RT f, (2.9)
7 R where : Ŵ W are restriction operators which extract the subdomain parts. Eliating the interior unknowns of each subdomain, i.e., eliating the upper left block of (2.9), we obtain where and Ŝ = = N i= N i= Ŝ u = g, (2.2) R T (A A R T S R A A T ) R = R T S R, (2.2) g = N i= R T (f A A f ). R : Ŵ W and S are direct sums of R and S, respectively. From (2.2), we can see that Ŝ can be regarded as the restriction of S, defined on W = N i= W, to the continuous subspace Ŵ. We can also view Ŝ as the restriction of S to Ŵ: Ŝ = R T S R, where R : Ŵ W. We introduce a few more restriction operators that constitute R. R : Ŵ W extracts the part that belongs to W from a vector in Ŵ. R Π : Ŵ ŴΠ is defined similarly. Thus R = R ()T R(N)T RT Π. We also define R D,, which are scaled versions of R ; each row of R has exactly one nonzero entry corresponding to a node x on the subdomain interface. Multiplying each such entry with δ i (x) results in the scaled version R D,. n the BDDC method, we use M = R D, T S R D, as the preconditioner. With this choice of preconditioner, the BDDC and the FET-DP algorithms have the same set of eigenvalues; see 22, 23, 2. 3 A Model Problem 3. Motivation and Notations Our ultimate goal is to solve multi-body contact problems without friction. Contact problems are characterized by an active area of contact, which is unknown a priori, and inequality constraints such as non-penetration conditions; e.g., see 2,. We consider two slightly different approaches to derive the same saddle point formulation we are going to solve. n the first approach the method of attack is an active set method from the outset, in which we solve a sequence of equality constrained imization problems. We choose to solve a saddle point formulation of each such equality constrained problem for the reason we will explain below. n the second approach, described e.g. in 9, a complementarity problem of the original imization problem is first considered. This problem can be expressed as a single nonlinear equation, and we can use a semismooth Newton method to solve this nonlinear problem. t turns out that the linear problems that are solved in this approach are identical to the saddle point problems of the first approach, and the two approaches differ only in the solution update process. n this paper, we develop the theory only for scalar elliptic problems with inequality constraints and present the following model problem, which is taken from 5, 6, as a motivation: 2 ( ) u i 2 dx fu i dx 2 Ω i Ω i i= 7
8 where u i H (Ω i ),i =,2, Ω = (,) (,),Ω 2 = (,2) (,) u = on u = {} (,) u 2 u on c = {} (,) (3.) The reason we consider only scalar elliptic problems is that the inequality constraints in scalar elliptic problems are much simpler than those in linear elasticity problems and this simplicity allows us to focus on the analysis of the preconditioned operator. We consider a generalization of (3.), a imization problems with multiple bodies Ω i,i =,,N, each of which has many degrees of freedom and is decomposed into subdomains Ω i,j,i =,,N,j =,,N i, constrained by an inequality condition: N ( ) ρ u i 2 dx fu i dx 2 i= Ω i Ω i where u i H (Ω i ), i =,,N, u i = on i u, N B u i (3.2) i= We assume that the boundary of at least one body is clamped, i.e., Ω D := N i= i u. We assume that there are no traction forces. We assume that ρ(x) = ρ i,j ρ >, x Ω i,j, i,j. We also assume the existence of a coefficient C, independent of i, such that ρ i := max j ρ i,j Cρ i,j, i,j. (3.3) We introduce two types of global interfaces: the first one is gl := i j Ω i Ω j, and can be viewed as the potential contact area between the bodies: in the model problem, this is c. The second one, the current contact area, is denoted by k gl, where k gl gl; the superscript k concerns the outer iteration of the active set method and reds us that the current active set/area changes. n each outer iteration of the active set method some of the inequality constraints are adopted as the corresponding equality constraints and the rest are ignored, and k gl,h, the discrete version of k gl, can be viewed as the collection of the nodes at which equality constraints are being imposed. We also introduce the local interfaces loc := j k ( Ω i,j Ω i,k ),i =,,N. We denote the union of the free boundaries by Ω F := ( i Ω i ) \ ( Ω D gl ). We denote the standard finite element space of continuous, piecewise linear functions on Ω i,j by W (i,j). Each W (i,j) is decomposed into a subdomain interior part W (i,j) and a subdomain interface part W (i,j) for functions on Ω i,j gl. We also recall that all functions vanish on the Dirichlet boundary Ω D. We define associated product spaces, W := N i j= W (i,j) and W := N i j= W (i,j). We introduce spaces analogous to W and Ŵ of Section 2. Functions in W are in general discontinuous across the local interface loc, and we define Ŵ as the continuous subspace of W. W := W Ŵ Π is an intermediate space between Ŵ and W. n the two-dimensional case, which is the focus of our analysis in this paper, W consists of functions that are continuous at subdomain vertices. Also, we let W := W W. We introduce the matrices A, which are the direct sums of the stiffness matrices A (i,j),j =,,N i for the individual subdomains: A = A (i,) A (i,ni),i =,,N. (3.4) 8
9 Figure 2: W c for FET-FET, Ŵc for Hybrid H b H s (a) FET-FET H b Hs (b) Hybrid We also introduce the Schur complements S (i,) S =, S (i,j) = A (i,j) S (i,ni) A(i,j)T A (i,j) A (i,j)t,j =,,N i. The finite element formulation of the problem (3.) in the space W c := N W i= and W,c := N i= has been considered in the FET-FET method 2, 8. n the hybrid method, we formulate the problem in the space Ŵ,c := N i= Ŵ. We will briefly describe the FET-FET method here. The problem (3.) can be expressed as u Wc 2 ut à c u f c T u, with Bc u. (3.5) Recalling that we are using an active set method to deal with the inequality conditions, we formulate the imization problem on the current active set: where à c = à () W u Wc 2 ut à c u f c T u, with Z k Bc u =, (3.6) à (N),u = u (). u (N), f c = f (). f (N),u W,i =,,N, and B c = Bloc B gl = B () loc... B (N) loc B () gl B (N) gl, B loc = Z k = Zgl k. B (i,) loc B (i,ni) loc,i =,,N, 9
10 Z k Bc u = in (3.6) indicates the continuity constraint across the local subdomain interface loc,i =,,N, as well as the continuity constraint across the global area of contact k gl. Zk gl is a square matrix obtained by replacing some of the diagonal entries of the identity matrix with zeros; only the entries corresponding to the nodes at which an equality is being imposed are retained. We use the superscript k as an indication that Zgl k and Zk change in each iteration of the active set method. We have B loc u =,u W, exactly when the values associated with more than one subdomain on the body Ω i coincide. Note that B loc has nonzero columns only for the components of W. We also introduce a scaled jump operator, B D,c : B D,c = Bloc,D B gl,d = B () loc,d B (N) loc,d B () gl,d B(N) gl,d and B loc,d = B (i,) loc,d B(i,Ni) loc,d,i =,,N. B loc,d and B gl,d are obtained in the same manner as the B D, of the one-level FET method (see section associated with the Lagrange multipliers for the continuity at the node (x), where N x,loc is the set of indices of the 2.2); the nonzero entry of B (i,j) loc x Ω i,j Ω i,k, multiplied by δ i,k (x) = ργ i,k (x)/ ρ γ s N i,s x,loc subdomains of Ω i with x on their boundary, is the corresponding entry of B loc,d. Similarly, the nonzero entry of B gl associated with the Lagrange multiplier for the continuity at the node x Ω i Ω j, multiplied by δ j (x) = s N (j) x,loc ρ γ j,s (x)/ k N x,gl,t N (k) x,loc ρ γ k,t (x), where N x,gl is the set of indices of the subdomains of any body which share the node x on their boundary, is the corresponding entry of B gl,d. Eliating the interior unknowns in all subdomains of each body, we obtain the following reduced imization problem, u W,c 2 ut S c u g c T u, with Z k B,c u =, (3.7) where S c = S () S (N),u = u ().. u (N), This imization problem has the following KKT system: S c (Z k B,c ) T u Z k B,c λ u = W,i =,,N. g. (3.8) t is natural to reduce this system to an equation for λ as in the one-level FET method and solve it with the PCG method in a proper subspace, using the following preconditioner: M D := Zk BD,c Sc BT D,c Z k. The resulting method, i.e., the FET-FET method, turns out not to be scalable with respect to the number of subdomains; see 8, 2. We now present a scalable alternative, which we name a hybrid method.
11 4 A Hybrid Method The hybrid method introduced in this section is a scalable alternative to the FET-FET method. We use finite element functions in the space Ŵ,c := N i= Ŵ to discretize the contact problem. W denotes a finite element space on Ω i gl. Here, OL stands for one-level; this is because we will use one-level FET type preconditioners, which is obtained by regarding each body as a single subdomain without further decomposition. We introduce the Schur complement Ŝ Ŝ OL on Ŵ, which can be obtained by restricting S to Ŵ : T = R S R,i =,,N, where R : Ŵ W. We also introduce restriction operators R : Ŵ W, R = R (i,).. R (i,ni) R Π, where R (i,j) extracts from a vector in Ŵ defined similarly. We also define the scaled versions (i,j) the part that belongs to W R D, : R D, = R (i,) D,. R (i,ni) D, R Π. and R Π : Ŵ Ŵ Π is Here, R (i,j) D, is obtained as follows: a nonzero entry of R(i,j), which corresponds to a node x Ω i,j,h \ Ω i,h, is multiplied by δ i,j (x), where δ ρ i,j (x) := γ i,j (x) k N x,loc ρ γ i,k (x). The restriction of the imization problem (3.7) in W,c to the subspace Ŵ,c is as follows: u Ŵ,c 2 ut Ŝcu ĝc T u, with Ẑ k B,c u =, (4.) where Ŝ c = Ŝ () Ŝ (N), and where B,c only retains rows and columns of B,c which correspond to Ŵ,c, and Ẑk is obtained by removing irrelevant columns and rows of Z k. n addition to B,c, we define another jump operator B OL which acts on vectors of the space N i= W OL. This operator is needed in the preconditioner for the hybrid method. Recall that B,c acts on vectors of the space Ŵ,c and only has rows corresponding to the Lagrange multipliers enforcing the continuity between different bodies, i.e., continuity across gl. Thus B OL and B,c differ only in that B,c has a number of zero columns which correspond to the nodes on loc,h. We note that B OL can be regarded as the jump operator for the one-level FET method resulting from viewing each body as
12 a subdomain. We can define the scaled jump operator B OL,D in the usual way. ntroducing a vector of Lagrange multipliers λ, we arrive at the following saddle point formulation of (4.): Find (u,λ) Ŵ,c range( B,c ) Ŝ c (Ẑk B,c ) T Ẑ k B,c u λ = ĝc. (4.2) We can solve (4.2) by reducing the system to an equation of λ alone in a proper subspace, but that requires inverting Ŝc, which is expensive. nstead, we keep the saddle point problem (4.2) as is and solve it by a Krylov subspace method which can deal with indefinite systems, such as the preconditioned conjugate residual (PCR) method. Due to the singularity of the matrix Ŝc, the solution of the upper part of the system (4.2) exists if and only if ĝ c (Ẑk B,c ) T λ range(ŝc). Most of the discussion here concerning this issue will be very similar to that of section 2.2 on one-level FET methods. As in the one-level FET method, we introduce a matrix R c such that range(r c ) = ker(ŝc): R c = R () where R consists of the null vectors of Ŝ,i =,,N. n the PCR iterations, we will use an initial vector of Lagrange multipliers λ which satisfies ĝ c (Ẑk B,c ) T λ range(ŝc), and an increment µ with Ẑk BT,c µ range(ŝc). Therefore the space of admissible increments is defined as follows: R (N), V k := {µ range( B,c ) : (Ẑk B,c ) T µ range(ŝc)} = ker(g kt ), where G k := Ẑk B,c R c. We introduce a projection operator P k for the Lagrange multipliers which is an orthogonal projection from U to V k = ker(g kt ): P k := G k (G kt G k ) G kt. We also introduce a subspace Ŵ,R := range(ŝc) of Ŵ,c. We rewrite (4.2) in terms of vectors in the subspace Ŵ,R V. First, noting that any admissible λ has a decomposition of the form λ = λ +µ,µ V, we rewrite the leading equation of (4.2) as Ŝ c u + (Ẑk B,c ) T µ = ĝ c (Ẑk B,c ) T λ. (4.3) Using (4.3) and that P kt µ = P k µ = µ, we can rewrite (4.2): Ŝ c (P k Ẑ k B,c ) T u ĝc Ẑ k B = B,c T λ,c µ. (4.4) The solution of (4.4) satisfies Ŝ c (P k Ẑ k B,c ) T P k Ẑ k B,c u λ ĝc = B,c T λ. (4.5) We use the system (4.5) in order to make sure that our iterates are in the subspace Ŵ,R V. Note that the displacement part of the solution of (4.5), u, does not necessarily satisfy the continuity condition Ẑk B,c u = ; we can recover a solution which satisfies all constraints via the operation u R c (G kt G k ) G kt B,c u, see 5, 2, 8. We now consider a second approach to solve (3.5), following 9; (3.5) can be rewritten as u Ŵ,c 2 ut Ŝcu ĝc T u, with B,c u, (4.6) 2
13 which is equivalent to the following problem { Ŝc u + B T,c λ = ĝ c B,c u, λ, λ T B,c u =. (4.7) The complementarity condition given in the second line is equivalent to C(u,λ) := λ max(,λ + c B,c u ) =, (4.8) for each c >. The system (4.7) can thus be expressed as the following nonlinear system of equations: { Ŝc u + B,c T λ = ĝ c (4.9) C(u,λ) =. t follows that a Newton step for the nonlinear system (4.9) is Ŝ c BT,c δu k ĝc (Ŝcu c B,c,A k k δλ k = k + B,c T λk ) C(u k,λ k ) and where (4.) u k+ = u k + δu k, λ k+ = λ k + δλ k. (4.) k = {i : (λ k + c B,c u k ) i }, A k = {i : (λ k + c B,c u k ) i > }, (4.2) and B,c,A k results from replacing row i of B,c with zeros, for all i A k. k is defined similarly. We can rewrite the second equation of (4.) as follows: and the first equation as (c B,c δu k ) i = (c B,c u k ) i, i A k, and (δλ k ) i = λ k i, i k, (4.3) Ŝ c δu k + B T,c,A k(δλk ) A k + B T,c, k(δλk ) k = ĝ c (Ŝcu k + B T,c,A k(λk ) A k + B T,c, k(λk ) k), (4.4) which is equivalent to Ŝ c δu k + B T,c,A k(δλk ) A k = ĝ c (Ŝcu k + B T,c,A k(λk ) A k), (4.5) due to (4.3). Consequently, we can rewrite the Newton step defined by (4.) and (4.) as Ŝ c BT,c,A δu k k ĝ c (Ŝcu k c B,c,A k δλ k = + B T,c,A (λ k ) k A k) c B,c,A ku k (4.6) and u k+ = u k + δu k, λ k+ = λ k + δλ k, where (δλ k ) i = λ k i, i k. (4.7) Note that with a given pair of (u k,λk ), (4.6) and (4.7) are equivalent to Ŝ c BT,c,A u k+ k c B ĝc,c,a k λ k+ =. (4.8) Notice the similarity between (4.2) and (4.8); in fact, B,c,A k = Ẑk B,c and with c =, (4.2) and (4.8) are the same. n the following we will solve (4.9), with c =, via Newton s method as defined by (4.6) and (4.7). We now discuss our choice of preconditioner. Let A OL denote the stiffness matrix for the entire body 3
14 Ω i : this needs to be distinguished from A, which is a direct sum of stiffness matrices for individual subdomains (see (3.4)). We have A OL = A A A T A, i =,,N, where A corresponds to the nodes on Ω i gl and A to those interior of Ω i, etc. We define the corresponding Schur complement S OL := A A A A T and also a block-diagonal matrix for the entire system as the direct sum of the Schur complements for the individual bodies: S OL := S () OL Since inverting A can be expensive in practice, we need to solve A x = b approximately; we propose a way of doing this in Section 7. We now introduce the following block-diagonal preconditioner for the system: B PR M = BDDC P R P k M D P k (4.9) M BDDC = S (N) OL where P R := R c (Rc T R c ) Rc T is an orthogonal projection operator onto range(ŝc) and R ()T () D, S R () D, where M D We rewrite the KKT system (4.5) as A := Ŝ c (P k Ẑ k B,c ) T P k Ẑ k B,c. R (N) T (N) D, S R (N) D, = Ẑk B OL,DS OL B T OL,DẐkT., Ax = F, (4.2), x := u λ and F := ĝc. (4.2) n our hybrid method, we use a preconditioned conjugate residual (PCR) method with B as the preconditioner. For a description of the preconditioned conjugate residual method, see 26, 8. 5 Convergence Estimates Suppose the following system is solved with the PCR method with the preconditioner M: We define Au = b. (5.) K(M A) = µ max µ = max{ λ : λ σ(m A)} { λ : λ σ(m A)}, (5.2) where σ(m A) is the spectrum of M A. We have the following result, see 26, C.6.2; a proof can be found in 8, Section
15 Lemma 5.. Let A be regular and symmetric and M symmetric and positive definite. Then, after k steps of the PCR algorithm, the norm of the residual is bounded by where ρ = K K+ M /2 r k 2 2ρµ + ρ 2µ M /2 r 2, and µ Z, such that k/2 < µ k/2. According to Lemma 5., we need to study the spectrum of the preconditioned operator B A, which has the same spectrum as B /2 AB /2. General Case We first study a general case, where assug that A B T A = B, B =  Ĉ α u T Âu u T Au α u T Âu, u. (5.3) We assume A, Rn n,ĉ Rm m are real symmetric and positive definite. Then, B /2 AB /2 =  /2A /2  /2 B T Ĉ Ĉ /2. B /2 n the following, we use the notation à :=  /2 A /2 and B := Ĉ /2 B /2. Note that α u T u u T Ãu α u T u, u. (5.4), We study the cases where matrix and  = A and  A separately. When  = A, à is simply the identity B /2 AB /2 BT =. (5.5) B Lemma 5.2. Let B /2 AB /2 be defined as in (5.5). We then have K(B A) = K(B /2 AB /2 ) = /2 + /4 + λ max /2 + /4 + λ, where λ max and λ are the largest and smallest eigenvalues of B T B, respectively. Proof. We consider the following eigenvalue problem: BT u B λ u = t λ, which is equivalent to u + B T λ = tu Bu = tλ (5.6) Notice that t due to the nonsingularity of A and B. Substituting the second equation of (5.6) into the first, we obtain u + t BT Bu = tu. Denoting the eigenvalues of B T B by λi,i =,,n, we obtain ( + λ i /t t)u =, i =,,n. 5
16 Since u = leads to λ =, we need to solve + λ i /t t =,i =,,n, which are equivalent to the quadratic equations t 2 t λ i =. Their solutions are /2 ± /4 + λ i and thus Clearly, σ(b A) = {/2 ± /4 + λ i : i =,,n}. max{ λ : λ σ(b A)} = /2 + /4 + λ max and { λ : λ σ(b A)} = /2 + /4 + λ, where λ max := max n i= λ i and λ := n i=λ i. We now consider the case  A. Then the eigenvalue analysis of A := B /2 AB /2 = à BT B is not as easy, and we left- and right- multiply this symmetrized preconditioned operator with C /2 = à /2 to obtain A 2 := C /2 A C /2 = à /2 BT Bà /2. (5.7) Eigenvalues of A 2 can be analyzed in the same manner as in Lemma 5.2. To relate the spectrum of A to the spectrum of A 2, we use the Courant-Fischer Minimax Theorem. Theorem 5.3 (Courant-Fischer). Let A R n n be a symmetric matrix with real eigenvalues λ i,i =,,n, which are ordered so that λ λ 2 λ n. Then λ k = max dim(v )=k λ k = dim(v )=n k+ x T Ax x V x T x x x max T Ax x V x T x x (5.8) (5.9) Let λ λ 2 λ n denote the eigenvalues of A 2, and λ λ 2 λ n the eigenvalues of A. Suppose λ k > and λ k+ <, where λ k is the smallest positive eigenvalue of A 2 and λ k+ the largest negative eigenvalue of A 2. Also, let q i,i =,,n denote the eigenvectors of A 2 such that A 2 q i = λ i q i and q T i q j = δ ij,i,j =,,n. Using (5.8) and the fact that A = C /2 A 2 C /2 we have λ k = max dim(v )=k x T A x x V x T x = max dim(v )=k x For V := C /2 span{q (),q (2),,q (k) }, we have Noting that due to the definition of C and (5.4), we have x V x x V x (C /2 x) T A 2 (C /2 x) λ x V (C /2 x) T (C /2 k. x) x (C /2 x) T A 2 (C /2 x) (C /2 x) T (C /2 x) (C /2 x) T (C /2 x) x T x α x T x x T Cx α x T x, x (5.) (C /2 x) T A (C /2 x) (C /2 x) T (C /2 x) (C /2 x) T (C /2 x) x T λ k α. x 6
17 Taking the maximum over all k-dimensional subspaces on the left hand side of the previous equation, we obtain λ k λ k α. Similarly, using (5.9), we have λ k+ = dim(v )=n k x max T A x x V x T x = dim(v )=n k x For V := C /2 {q (k+),,q (n) }, we have and max x V x max x V x max (C /2 x) T A 2 (C /2 x) λ x V (C /2 x) T (C /2 k+ x) x (C /2 x) T A 2 (C /2 x) (C /2 x) T (C /2 x) (C /2 x) T (C /2 x) x T x (C /2 x) T A 2 (C /2 x) (C /2 x) T (C /2 x) (C /2 x) T (C /2 x) x T λ k+ α. x Taking the imum on the left hand side of the previous equation, we obtain By a similar argument, λ k+ λ k+ α. λ λ α and λ n λ n α. Letting λ max and λ denote the maximum and the imum eigenvalues of à /2 BT Bà /2, respectively, we obtain K(A ) = max{ λ, λ n } { λ k, λ k+ } α max{λ, λ n } α {λ k, λ k+ } = α K(A 2 ) α /2 + /4 + λ max α α /2 +, (5.) /4 + λ where the second inequality follows from the definition of A 2 in (5.7) and Lemma 5.2. Noticing that and λ max (à /2 BT Bà /2 ) λ max ( B T B)λmax (à ) λ (à /2 BT Bà /2 ) λ ( B T B)λ (à ), α u T u u T à u α u T u, we rewrite (5.) in terms of λ max and λ, the maximum and the imum eigenvalues of B T B and obtain: K(A ) α /2 + /4 + λmax /α α /2 +. (5.2) /4 + λ /α Special Case We now use these results to study the convergence bound of our preconditioned system B A, where B and A are defined in (4.9) and (4.2), respectively. We have A = Ŝ, B = P k Ẑ k B,c,  = P R M BDDC P R, u, Ĉ = P k M D P k. Notice that A,Â, and Ĉ are now singular. However, this does not pose any problem, since in the application of the PCR method our iterates will be in a proper subspace in which those matrices will be nonsingular. From (5.2), we can see that the extreme eigenvalues of BT B and α,α in (5.3) are important parameters, where B T B =  /2 B T Ĉ B /2, which has the same spectrum as B B T Ĉ. n our case, B B T Ĉ = P k Ẑ k B,c P R M BDDC P T R B,cẐkT P k P k M D P k. (5.3) The following lemma indicates the spectral equivalence between the matrices A and Â; for a proof, see 2. 7
18 Lemma 5.4. ( x T Âx x T Ax C + log ( Hs h )) 2 x T Âx, (5.4) for all x range(ŝ). Thus, we can study the spectrum of BA B T Ĉ = P k Ẑ k B,c P R Ŝ P R = P k Ẑ k B,c Ŝ instead of that of the matrix (5.3), where Lemma 5.5. Proof. Let B,c = B T v, where A (i,) B T R (i,)t A (i,) B,c Ŝ Ŝ = T B,cẐkT P k P k M D P k B T,cẐkT P k P k Ẑ k B OL,D S OL B T OL,D Ẑ kt P k (5.5) Ŝ () Ŝ (N). B T,cẐkT P k = B OL S OL BT OL Ẑ kt P k. () (N) B,, B,B OL = B () OL,,B (N) OL. Note that the solution of Ŝ u = v range(ŝ ), can be obtained from the following equation: A (i,2) R (i,2)t A (i,2) A (i,)t R (i,) A (i,2)t R (i,ni)t = A (i,ni) A (i,ni). B T v A (i,ni)t Ni j= R(i,j)T R (i,2) R (i,ni) A (i,j) R(i,j) u (i,) u (i,2)... u (i,ni) û, (5.6) where R (i,j) : Ŵ (i,j) T W is a restriction operator. Noting that all entries of B v, corresponding to the nodes on loc, are zero and eliating those entries results in BT OL v, we can rearrange the system (5.6): A OL u = A A A T A u u = B T OL v where u is the displacement on gl Ω i. The equivalence of (5.6) and (5.7) shows that B OL S OL BT OL Ẑ kt P k. B Ŝ (5.7) B T Ẑ kt P k = Due to Lemma 5.5, the operator (5.5) can be written as P k Ẑ k B OL S OL BT OL Ẑ kt P k P k Ẑ k B OL,D S OL B T OL,D Ẑ kt P k. (5.8) The proof of the following lemma proceeds, line by line, as the proof of 26, Theorem
19 Lemma 5.6. For λ range(p k ), λ,λ P k ẐB k OL S OL BT OL Ẑ kt P k P k ẐB k OL,D S OL B T OL,D Ẑ kt P k λ,λ C( + log(h b /h)) 2 λ,λ. We now can derive a concrete bound for (5.2), using Lemmas 5.4 and 5.6. n our case, λ max = C( + log(h b /h)) 2 ( + log(h s /h)) 2,λ =,α = C( + log(h s /h)) 2, and α =. Assug that H b /h and H s /h are large enough, we have λ max λmax, (5.9) α α and λ = α 2 + ( ( ( + 4 λ = 2 α λ + O 4 λ ) )) 2. (5.2) 2 α α Combining (5.2), (5.9), and (5.2), we have We obtain K(A ) C( + log(h b /h))( + log(h s /h)) 5. (5.2) Theorem 5.7. Let B, A, and K(B A) be defined as in (4.9), (4.2), and (5.2), respectively. Then we have the following bound: 6 Numerical Experiments K(B A) C( + log(h b /h))( + log(h s /h)) 5. Recall that an active set method consists of outer iterations, in which the active set is updated, and inner iterations, in which auxiliary equality constrained problems are solved on the current active set. n this section, we solve the latter by using the hybrid method. We note that such problems were solved using the FET-FET method in 2, 8. We solve the following imization problem: N b N b i= ( ) u i 2 dx fu i dx, (6.) 2 Ω i Ω i where Ω i R 2,i =,,N b N b, are square bodies with side length H b := /N b which form the N b N b system Ω = Ω i =,,. We require u i H (Ω i ),u i Ωi Ω =. Each Ω i is decomposed i= into N s N s square subdomains, each of which is discretized by square bilinear elements of side length h. Also, := i j Ω i Ω j denotes the interface between the bodies. We consider two linearized problems, each with a different contact area between the bodies. n the first problem, the entire is considered as the contact area, i.e., we require the continuity of the displacement vector across the entire. n the second problem, continuity is imposed only on the middle third of the faces between the bodies. We use the preconditioned conjugate residual method. All our experiments have been performed in MATLAB, and the stopping criterion is r n 2 / r 2 < 5, where r n and r are the n th and initial residuals, respectively. n Table, the results obtained with the hybrid method are presented. We have three parameters; the number of bodies across Ω (N b = /H b ), the number of subdomains across each body (N s = H b /H s ), and the number of elements across each subdomain (H s /h). We vary one parameter while keeping the 9
20 Table : Results for the hybrid method. iter denotes the iteration counts. Area on which continuity is imposed between bodies:, i.e., the entire interface for (), and only a proper subset of, for () () () /H b H b /H s H s /h iter iter other two fixed. The results for the first set of experiments, with the entire as the contact surface, are shown in Column (); those for the second set of experiments with a reduced contact area shown in Column (). We observe that the iteration counts are independent of /H b and logarithmically dependent on H s /h. The iteration counts from Table are also plotted in Figure 3. Very similar numerical results have been obtained independently by Klawonn and Rheinbach; see 3 and 4. 7 Active set method combined with the hybrid method t is well known that an active set method can often be slow due to a poor initial guess. However, it has been shown in 9 that a certain primal-dual active set strategy, viewed as a Newton method, has a superlinear convergence provided that the initial point is close to the solution. A very efficient strategy of finding a good initial active set and λ was discussed briefly in 2, Chapter 5 and in full detail in 9. n our experiments this turns out to be a good estimate of the optimal active set. We set u =. The following is the complete algorithm:. Set u = and choose λ as described in 2, Chapter 5, 9. Set k =. 2. Set k = {i : (λ k + B,c u k ) i }, A k = {i : (λ k + B,c u k ) i > }. 3. Solve Ŝ c BT,c,A k B,c,A k u k+ λ k+ = ĝc (7.) and set λ k+ = on k. 2
21 Figure 3: teration counts for the hybrid method for two different contact areas between the bodies, namely, i.e., the entire interface for (), and only a proper subset of, for () 2 teration counts Experiments () Experiments () teration Counts Experiments () Experiments () teration Counts Experiments () Experiments () /H b (a) /H b varies H s /H b (b) H b /H s varies log 2 (H s /h) (c) H s/h varies 4. Stop if A k+ = A k and k+ = k. Otherwise return to 2. We have solved the nonlinear model problem (3.) with the method described above; the results are reported in Table 2. Recall that the application of our preconditioner requires the solution of linear systems with the matrices A,i =,,N, i.e., solving Dirichlet problems for the bodies; this could potentially be expensive and in such a case these Dirichlet problems need to be solved inexactly, e.g., by an iteration. We describe efficient preconditioners Ã, for A, for the preconditioned conjugate gradient method. Let Ŵ and W (i,j) denote the space of continuous finite element functions on loc and Ω i,j loc, which are similar to Ŵ (i,j) and W, respectively. Also, we define a restriction operator R : Ŵ W (i,j). After a symmetric permutation, we can obtain A = A (i,) R (i,)t A (i,) A (i,2) R (i,2)t A (i,2) A (i,)t R (i,) A (i,2)t R (i,ni)t A (i,ni) A (i,ni) A (i,ni)t R (i,j) T R (i,2) R (i,ni) A (i,j) R (i,j) The solution of A x = b can be found by a block factorization. More precisely, with we have where Ŝ := x (j) N i j= = A (i,j) (b (j) Ŝ R (i,j)t (A (i,j) A (i,j) A(i,j) A (i,j)t )R (i,j),. (7.2) A (i,j)t R (i,j) x ), j =,,N i, (7.3) N i x = b R (i,j)t A (i,j) j= A Solving (7.4) can be expensive; the solution of à x = b can be defined as b (j). (7.4) where x (j) x = = A (i,j) (b (j) T R D, S R D, A (i,j)t R (i,j) x ), j =,,N i, (7.5) N i b j= R (i,j)t A (i,j) A b (j), (7.6) 2
22 Table 2: Results: primal-dual active set method + hybrid method. outer it. denotes the number of outer iterations of the active set method; inner it. denotes the number of iterations needed to solve the inner imization problems by the PCR method, until the norm of the residual has been reduced by 5, on the active faces identified in the outer iterations. total it. denotes the total number of inner iterations. N sub (/H) H/h N dof (λ) N dof (total) outer it. inner it. total it. 6(4) (4) (4) (4) (8) (8) (8) (8) (2) (2) (2) (2) (6) (6) (6) (6) R D, and S R D, with defined similarly as and S, respectively. n Table 2, notice that the iteration counts for the inner imizations does not increase rapidly as we increase the number of elements per subdomain or the number of subdomains per body, which is an indication of the scalability of the hybrid algorithm. Also, notice that it takes at most two outer iterations to reach the optimal solution, which is an indication of the effectiveness of the strategy to find an initial active set 9. Figure 4: Solution of the model problem, from different angles. N sub = 6,H/h = References Philip Avery and Charbel Farhat. The FET family of domain decomposition methods for inequality-constrained quadratic programg: Application to contact problems with conforg and nonconforg interfaces. Computer Methods in Applied Mechanics and Engineering, 98(2-26): , 29. Advances in Simulation-Based Engineering Sciences - Honoring J. Tinsley Oden. 22
23 2 Philip Avery, Gert Rebel, Michel Lesoinne, and Charbel Farhat. A numerically scalable dual-primal substructuring method for the solution of contact problems part : the frictionless case. Comput. Methods Appl. Mech. Engrg., 93(23-26): , M. Bergounioux, M. Haddou, M. Hintermüller, and K. Kunisch. A comparison of a Moreau Yosidabased active set strategy and interior point methods for constrained optimal control problems. SAM J. on Optimization, (2):495 52, 2. 4 Maïtine Bergounioux, Kazufumi to, and Karl Kunisch. Primal-dual strategy for constrained optimal control problems. SAM J. Control Optim., 37(4):76 94, Zdeněk Dostál. Optimal quadratic programg algorithms. With applications to variational inequalities., volume 23 of Springer Optimization and ts Applications. Springer, New York, Zdeněk Dostál, David Horák, and Dan Stefanica. A scalable FET-DP algorithm for a semi-coercive variational inequality. Comput. Methods Appl. Mech. Engrg., 96(8): , Charbel Farhat, Po-Shu Chen, and Jan Mandel. A scalable Lagrange multiplier based domain decomposition method for time-dependent problems. nternat. J. Numer. Methods Engrg., 38: , Wolfgang Hackbusch. terative solution of large sparse systems of equations, volume 95 of Applied Mathematical Sciences. Springer-Verlag, New York, 994. Translated and revised from the 99 German original. 9 M. Hintermüller, K. to, and K. Kunisch. The primal-dual active set strategy as a semismooth Newton method. SAM J. on Optimization, 3(3): , 22. Kazufumi to and Karl Kunisch. Augmented Lagrangian methods for nonsmooth, convex optimization in Hilbert spaces. Nonlinear Anal., 4(5-6):59 66, 2. Axel Klawonn and Oliver Rheinbach. A parallel implementation of dual-primal FET methods for three-dimensional linear elasticity using a transformation of basis. SAM J. Sci. Comput., 28(5):886 96, Axel Klawonn and Oliver Rheinbach. Robust FET-DP methods for heterogeneous three dimensional elasticity problems. Comput. Methods Appl. Mech. Engrg., 96(8):4 44, Axel Klawonn and Oliver Rheinbach. A hybrid approach to 3-level FET. PAMM Proc. Appl. Math. Mech., 8():84 843, Axel Klawonn and Oliver Rheinbach. Highly scalable parallel domain decomposition methods with an application to biomechanics. ZAMM Z. Angew. Math. Mech., 9():5 32, 2. 5 Axel Klawonn and Olof B. Widlund. A domain decomposition method with Lagrange multipliers and inexact solvers for linear elasticity. SAM J. Sci. Comput., 22(4):99 29, 2. 6 Axel Klawonn and Olof B. Widlund. Dual-primal FET methods for linear elasticity. Comm. Pure Appl. Math., 59(): , Axel Klawonn, Olof B. Widlund, and Maksymilian Dryja. Dual-primal FET methods for threedimensional elliptic problems with heterogeneous coefficients. SAM J. Numer. Anal., 4():59 79, Jungho Lee. Convergence analysis of a new domain decomposition method for a linearized contact problem. 2. Submitted to SAM J. Sci. Comput. 9 Jungho Lee. A strategy of finding an initial active set for inequality constrained quadratic programg problems. 2. Submitted to Optimization Methods and Software. 23
24 2 Jungho Lee. A Hybrid Domain Decomposition Method and its Applications to Contact Problems. PhD thesis, Courant nstitute of Mathematical Sciences, September Jing Li and Olof B. Widlund. FET-DP, BDDC, and block Cholesky methods. nternat. J. Numer. Methods Engrg., 66(2):25 27, Jan Mandel and Clark R. Dohrmann. Convergence of a balancing domain decomposition by constraints and energy imization. Numer. Linear Algebra Appl., (7): , 23. Dedicated to the 7th birthday of vo Marek. 23 Jan Mandel, Clark R. Dohrmann, and Radek Tezaur. An algebraic theory for primal and dual substructuring methods by constraints. Appl. Numer. Math., 54(2):67 93, Jan Mandel and Radek Tezaur. Convergence of a substructuring method with Lagrange multipliers. Numer. Math., 73(4): , Jan Mandel and Radek Tezaur. On the convergence of a dual-primal substructuring method. Numer. Math., 88(3): , Andrea Toselli and Olof Widlund. Domain decomposition methods algorithms and theory, volume 34 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin, Xue Tu. Three-level BDDC in two dimensions. nternat. J. Numer. Methods Engrg., 69():33 59,
RESEARCH 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions Juan Gabriel Calvo Preprint No. 25-2017 PRAHA 2017 A VIRTUAL
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 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 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 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 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 informationThe antitriangular factorisation of saddle point matrices
The antitriangular factorisation of saddle point matrices J. Pestana and A. J. Wathen August 29, 2013 Abstract Mastronardi and Van Dooren [this journal, 34 (2013) pp. 173 196] recently introduced the block
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 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 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 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 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 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 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
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 informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
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 informationA new algorithm for solving 3D contact problems with Coulomb friction
A new algorithm for solving 3D contact problems with Coulomb friction Radek Kučera VŠB-TU Ostrava, Department of Mathematics and Descriptive Geometry, name@email.address Summary. The paper deals with solving
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 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 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 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 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 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 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 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 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 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 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 information1 Computing with constraints
Notes for 2017-04-26 1 Computing with constraints Recall that our basic problem is minimize φ(x) s.t. x Ω where the feasible set Ω is defined by equality and inequality conditions Ω = {x R n : c i (x)
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 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 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 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 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 informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES On the approximation of a virtual coarse space for domain decomposition methods in two dimensions Juan Gabriel Calvo Preprint No. 31-2017 PRAHA 2017
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 informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More 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 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 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 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 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 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 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 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 informationStructural and Multidisciplinary Optimization. P. Duysinx and P. Tossings
Structural and Multidisciplinary Optimization P. Duysinx and P. Tossings 2018-2019 CONTACTS Pierre Duysinx Institut de Mécanique et du Génie Civil (B52/3) Phone number: 04/366.91.94 Email: P.Duysinx@uliege.be
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 informationMarch 8, 2010 MATH 408 FINAL EXAM SAMPLE
March 8, 200 MATH 408 FINAL EXAM SAMPLE EXAM OUTLINE The final exam for this course takes place in the regular course classroom (MEB 238) on Monday, March 2, 8:30-0:20 am. You may bring two-sided 8 page
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 informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationA DECOMPOSITION PROCEDURE BASED ON APPROXIMATE NEWTON DIRECTIONS
Working Paper 01 09 Departamento de Estadística y Econometría Statistics and Econometrics Series 06 Universidad Carlos III de Madrid January 2001 Calle Madrid, 126 28903 Getafe (Spain) Fax (34) 91 624
More informationElectronic Transactions on Numerical Analysis Volume 49, 2018
Electronic Transactions on Numerical Analysis Volume 49, 2018 Contents 1 Adaptive FETI-DP and BDDC methods with a generalized transformation of basis for heterogeneous problems. Axel Klawonn, Martin Kühn,
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More 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 information1 Constrained Optimization
1 Constrained Optimization Let B be an M N matrix, a linear operator from the space IR N to IR M with adjoint B T. For (column vectors) x IR N and y IR M we have x B T y = Bx y. This vanishes for all y
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationChapter 7. Optimization and Minimum Principles. 7.1 Two Fundamental Examples. Least Squares
Chapter 7 Optimization and Minimum Principles 7 Two Fundamental Examples Within the universe of applied mathematics, optimization is often a world of its own There are occasional expeditions to other worlds
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
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 informationMarch 5, 2012 MATH 408 FINAL EXAM SAMPLE
March 5, 202 MATH 408 FINAL EXAM SAMPLE Partial Solutions to Sample Questions (in progress) See the sample questions for the midterm exam, but also consider the following questions. Obviously, a final
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 informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationTheoretically supported scalable FETI for numerical solution of variational inequalities
Theoretically supported scalable FETI for numerical solution of variational inequalities Zdeněk Dostál and David Horák Abstract The FETI method with a natural coarse grid is combined with recently proposed
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationDomain decomposition on different levels of the Jacobi-Davidson method
hapter 5 Domain decomposition on different levels of the Jacobi-Davidson method Abstract Most computational work of Jacobi-Davidson [46], an iterative method suitable for computing solutions of large dimensional
More informationPart 4: Active-set methods for linearly constrained optimization. Nick Gould (RAL)
Part 4: Active-set methods for linearly constrained optimization Nick Gould RAL fx subject to Ax b Part C course on continuoue optimization LINEARLY CONSTRAINED MINIMIZATION fx subject to Ax { } b where
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationIr O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )
Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
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 informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
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 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 informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More 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 informationAffine covariant Semi-smooth Newton in function space
Affine covariant Semi-smooth Newton in function space Anton Schiela March 14, 2018 These are lecture notes of my talks given for the Winter School Modern Methods in Nonsmooth Optimization that was held
More information