SUMMARYA new variant of the FETI method for numerical solution of elliptic PDE is presented. The basic idea is to simplify inversion of the stiffness matrices of subdomains by using Lagrange multipliers not only for gluing the subdomains along the auxiliary interfaces, but also for implementation of the Dirichlet boundary conditions. Results of numerical experiments are presented which indicate that the new method may be even more efficient then the original FETI.
In this paper we discuss a specialization of the augmented Lagrangian-type algorithm of Conn, Gould, and Toint to the solution of strictly convex quadratic programming problems with simple bounds and equality constraints. The new feature of the presented algorithm is the adaptive precision control of the solution of auxiliary problems in the inner loop of the basic algorithm which yields a rate of convergence that does not have any term that accounts for inexact solution of auxiliary problems. Moreover, boundedness of the penalty parameter is achieved for the precision control used. Numerical experiments illustrate the efficiency of the presented algorithm and encourage its usage.
SUMMARYThe direct methods for the solution of systems of linear equations with a symmetric positive-semidefinite (SPS) matrix A usually comprise the Cholesky decomposition of a nonsingular diagonal block A JJ of A and effective evaluation of the action of a generalized inverse of the corresponding Schur complement. In this note we deal with both problems, paying special attention to the stiffness matrices of floating structures without mechanisms. We present a procedure which first identifies a well-conditioned positive-definite diagonal block A JJ of A, then decomposes A JJ by the Cholesky decomposition, and finally evaluates a generalized inverse of the Schur complement S of A JJ . The Schur complement S is typically very small, so the generalized inverse can be effectively evaluated by the singular value decomposition (SVD). If the rank of A or a lower bound on the nonzero eigenvalues of A are known, then the SVD can be implemented without any 'epsilon'. Moreover, if the kernel of A is known, then the SVD can be replaced by effective regularization. The results of numerical experiments show that the proposed method is useful for effective implementation of the FETI-based domain decomposition methods.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.