The extended unsymmetric frontal solution for multiple-point constraints

Purpose – The purpose of this paper is to discuss the linear solution of equality constrained problems by using the Frontal solution method without explicit assembling. Design/methodology/approach – Re-written frontal solution method with a priori pivot and front sequence. OpenMP parallelization, nearly linear (in elimination and substitution) up to 40 threads. Constraints enforced at the local assembling stage. Findings – When compared with both standard sparse solvers and classical frontal implementations, memory requirements and code size are significantly reduced. Research limitations/implications – Large, non-linear problems with constraints typically make use of the Newton method with Lagrange multipliers. In the context of the solution of problems with large number of constraints, the matrix transformation methods (MTM) are often more cost-effective. The paper presents a complete solution, with topological ordering, for this problem. Practical implications – A complete software package in Fortran 2...

[1]  Anshul Gupta,et al.  Recent advances in direct methods for solving unsymmetric sparse systems of linear equations , 2002, TOMS.

[2]  Xiaowei Shi,et al.  An improved GPS method with a new pseudo-peripheral nodes finder in finite element analysis , 2012 .

[3]  W. F. Riley,et al.  Engineering Mechanics: Statics , 1993 .

[4]  S. Valliappan,et al.  A solution algorithm for linear constraint equations in finite element analysis , 1978 .

[5]  John Reid,et al.  Ordering symmetric sparse matrices for small profile and wavefront , 1999 .

[6]  Barbara Chapman,et al.  Using OpenMP - portable shared memory parallel programming , 2007, Scientific and engineering computation.

[7]  P. Areias,et al.  Finite element formulation for modeling nonlinear viscoelastic elastomers , 2008 .

[8]  S. Sloan.  An algorithm for profile and wavefront reduction of sparse matrices , 1986 .

[9]  Iain S. Duff,et al.  The Multifrontal Solution of Unsymmetric Sets of Linear Equations , 1984 .

[10]  Timothy A. Davis,et al.  A combined unifrontal/multifrontal method for unsymmetric sparse matrices , 1999, TOMS.

[11]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[12]  Jennifer A. Scott.  Parallel frontal solvers for large sparse linear systems , 2003, TOMS.

[13]  Jennifer A. Scott.  On ordering elements for a frontal solver , 1999 .

[14]  Iain S. Duff,et al.  Design features of a frontal code for solving sparse unsymmetric linear systems out-of-core , 1984 .

[15]  P. Hood,et al.  Frontal solution program for unsymmetric matrices , 1976 .

[16]  Stuart S. Antman,et al.  Material constraints, lagrange multipliers, and compatibility. Applications to rod and shell theories , 1991 .

[17]  F. Amirouche.  Fundamentals Of Multibody Dynamics. Theory And Applications , 2005 .

[18]  Fred G. Gustavson,et al.  Two Fast Algorithms for Sparse Matrices: Multiplication and Permuted Transposition , 1978, TOMS.

[19]  Timothy A. Davis,et al.  Direct methods for sparse linear systems , 2006, Fundamentals of algorithms.

[20]  Bruce M. Irons,et al.  A frontal solution program for finite element analysis , 1970 .

[21]  M. Ainsworth.  Essential boundary conditions and multi-point constraints in finite element analysis , 2001 .

[22]  M. Shephard.  Linear multipoint constraints applied via transformation as part of a direct stiffness assembly process , 1984 .

[23]  I. Duff,et al.  Direct Methods for Sparse Matrices , 1987 .

[24]  Charles L. Lawson,et al.  Basic Linear Algebra Subprograms for Fortran Usage , 1979, TOMS.

[25]  A. Pothen,et al.  Two improved algorithms for envelope and wavefront reduction , 1997 .

[26]  Nicholas I. M. Gould,et al.  On Modified Factorizations for Large-Scale Linearly Constrained Optimization , 1999, SIAM J. Optim..

[27]  T. Rabczuk,et al.  Implicit solutions with consistent additive and multiplicative components , 2012 .

[28]  Dieter Jungnickel,et al.  Graphs, Networks, and Algorithms , 1980 .

[29]  Jacques Carlier,et al.  Improving the finite element ordering for the frontal solver , 2000 .

[30]  W. Rheinboldt.  Geometric notes on optimization with equality constraints , 1996 .

[31]  Mark S. Shephard,et al.  An algorithm for multipoint constraints in finite element analysis , 1979 .

[32]  S. Antman.  Nonlinear problems of elasticity , 1994 .

[33]  Alan George,et al.  An Implementation of a Pseudoperipheral Node Finder , 1979, TOMS.