Mesh-based multi-frontal solver with reuse of partial LU factorizations for antenna array
MetadataShow full item record
© 2016 Elsevier B.V. There are a number of relevant physical problems in which their problem domains consist of the repetition of a given subdomain. The traditional multi-frontal solver implementations, like MUMPS or SuperLU, get on the input the global sparse linear system of equations. They are not aware of the structure of the computational mesh. They do not know that some parts of the mesh, i.e., some sub-domains are identical. In such a case, some sub-matrices of the global matrix are identical. However, when we assemble the matrices corresponding to identical sub-domains into a global sparse system, they overlap, and we ignore that they corresponded to identical sub-domains. In this paper we advocate another approach to this computational problem, based on the additional knowledge of the structure of the computational mesh. We propose a wrapper over a multi-frontal solver that partitions the computational problem into a cascade of sub-problems, for which a traditional multi-frontal solver is called and asked for the Schur complements. Such solver wrapper can massively reuse computations performed over identical sub-domains, as well as it can propagate this reuse technique towards further elimination steps. We test our reuse solver on a problem consisting of an arrays of antennas and compare it against the execution time of a traditional sparse matrix-based multi-frontal solver called for the entire domain.
Showing items related by title, author, creator and subject.
Parallel simulations of 3d DC borehole resistivity measurements with goaloriented self-adaptive hp finite element methodPaszynski, M.; Pardo, D.; Calo, Victor (2012)In this paper we present a parallel algorithm of the goal-oriented self-adaptive hp Finite Element Method (hp-FEM) with shared data structures and with parallel multi-frontal direct solver. The algorithm generates in a ...
Calo, Victor; Pardo, D.; Paszynski, M. (2011)In this paper we present a goal-oriented self-adaptive hp Finite Element Method (hp-FEM) with shared data structures and a parallel multi-frontal direct solver. The algorithm automatically generates (without any user ...
Computational complexity and memory usage for multi-frontal direct solvers in structured mesh finite elementsCollier, N.; Pardo, D.; Paszynski, M.; Calo, Victor (2012)The multi-frontal direct solver is the state-of-the-art algorithm for the direct solution of sparse linear systems. This paper provides computational complexity and memory usage estimates for the application of the ...