Bisections-Weighted-by-Element-Size-and-Order Algorithm to Optimize Direct Solver Performance on 3D hp-adaptive Grids

Hassan M. AbouEisha, V. M. Calo, K. Jopek, Mikhail Moshkov, A. Paszyńska, M. Paszyński

Research output: Chapter in Book/Report/Conference proceedingConference contribution

1 Scopus citations


The hp-adaptive Finite Element Method (hp-FEM) generates a sequence of adaptive grids with different polynomial orders of approximation and element sizes. The hp-FEM delivers exponential convergence of the numerical error with respect to the mesh size. In this paper, we propose a heuristic algorithm to construct element partition trees. The trees can be transformed directly into the orderings, which control the execution of the multi-frontal direct solvers during the hp refined finite element method. In particular, the orderings determine the number of floating point operations performed by the solver. Thus, the quality of the orderings obtained from the element partition trees is important for good performance of the solver. Our heuristic algorithm has been implemented in 3D and tested on a sequence of hp-refined meshes. We compare the quality of the orderings found by the heuristic algorithm to those generated by alternative state-of-the-art algorithms. We show 50% reduction in flops number and execution time.
Original languageEnglish (US)
Title of host publicationComputational Science – ICCS 2018
PublisherSpringer Nature
Number of pages13
ISBN (Print)9783319937007
StatePublished - Jun 12 2018


Dive into the research topics of 'Bisections-Weighted-by-Element-Size-and-Order Algorithm to Optimize Direct Solver Performance on 3D hp-adaptive Grids'. Together they form a unique fingerprint.

Cite this