Computational & Technology Resources
an online resource for computational,
engineering & technology publications
Civil-Comp Proceedings
ISSN 1759-3433
CCP: 95
Edited by:
Paper 68

General Parallel Finite and Spectral-Element Oriented C/C++ Framework

A.D. Otero1,2 and J. Quinteros3

1College of Engineering and College of Sciences, University of Buenos Aires, Argentina
2National Council of Scientific and Technical Research, Argentina
3Deutsches GeoForschungsZentrum GFZ, Postdam, Germany

Full Bibliographic Reference for this paper
A.D. Otero, J. Quinteros, "General Parallel Finite and Spectral-Element Oriented C/C++ Framework", in , (Editors), "Proceedings of the Second International Conference on Parallel, Distributed, Grid and Cloud Computing for Engineering", Civil-Comp Press, Stirlingshire, UK, Paper 68, 2011. doi:10.4203/ccp.95.68
Keywords: parallel computing, finite element, spectral element, C++, iterative solver, direct solver.


The framework includes interfaces to libraries so that the model is isolated from implementation details and those libraries be exchanged. The main idea behind the implemented interfaces is to abstract the model from the libraries that are being used. One class interfaces with the selected Lapack implementation while another class was designed to interface with sparse matrices libraries. These libraries provide fast and tested procedures to solve linear systems stored in sparse matrices.

In the case of this parallel implementation the problem is divided into local (i.e. subdominial) and global problems to be solved. The local problem is solved by means of a direct sparse solver. We analyse iterative as well as direct solvers for the global problem. As a test case, we consider the solution of a general scalar Poisson problem to conduct experiments in order to evaluate the scaling capabilities of the parallel extensions of our code, both using iterative and direct solvers for the global problem.

Results are analysed and conclusions extracted with focus on future lines of development and improvements. The main conclusion regarding the direct solver is the need to implement a parallel global solver and, even further, to work with distributed matrices. The principal conclusion considering the iterative solver are the savings in computing time achieved with respect to the direct solver and the necessity of implementing good preconditioners to reduce the number of iterations.

J. Quinteros, P.M. Jacovkis, V.A. Ramos, "Dise no flexible y modular de modelos numéricos basados en elementos finitos", in S.A. Elaskar, E.A. Pilotta, G.A. Torres, (Editors), "Mecánica Computacional", Volume XXVI, 1724-1740, 2007.
J. Quinteros, V.A. Ramos, P.M. Jacovkis, "An elasto-visco-plastic model using the finite element method for crustal and lithospheric deformation", Journal of Geodynamics, 48(2), 83-94, 2009. doi:10.1016/j.jog.2009.06.006

purchase the full-text of this paper (price £20)

go to the previous paper
go to the next paper
return to the table of contents
return to the book description
purchase this book (price £85 +P&P)