Computational Technology Resources - CCP

Keywords: parallel computing, finite element, spectral element, C++, iterative solver, direct solver.

Summary

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.

References

1: 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.
2: 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)

	Computational & Technology Resources an online resource for computational, engineering & technology publications
	not logged in - login
Front Page Browse CCP CSETS CTR IJRT Other Authors Search Purchase Guide FAQ Contact us	Civil-Comp Proceedings ISSN 1759-3433 CCP: 95 PROCEEDINGS OF THE SECOND INTERNATIONAL CONFERENCE ON PARALLEL, DISTRIBUTED, GRID AND CLOUD COMPUTING FOR ENGINEERING Edited by: Paper 68 General Parallel Finite and Spectral-Element Oriented C/C++ Framework A.D. Otero^1,2 and J. Quinteros³ ¹College of Engineering and College of Sciences, University of Buenos Aires, Argentina ²National Council of Scientific and Technical Research, Argentina ³Deutsches GeoForschungsZentrum GFZ, Postdam, Germany doi:10.4203/ccp.95.68 purchase the full-text of this paper 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. Summary 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. References 1 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. 2 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)
Back to top	©Civil-Comp Limited 2023 - terms & conditions