Publication:
BCYCLIC: A parallel block tridiagonal matrix cyclic solver

Simple item page
dc.affiliation.dptoUC3M. Departamento de Físicaes
dc.affiliation.grupoinvUC3M. Grupo de Investigación: Física de Plasmases
dc.contributor.authorHirshman, S. P.
dc.contributor.authorPerumalla, K. S.
dc.contributor.authorLynch, V. E.
dc.contributor.authorSánchez Fernández, Luis Raúl
dc.date.accessioned2010-06-17T09:37:17Z
dc.date.available2010-06-17T09:37:17Z
dc.date.issued2010-09
dc.description13 pages, 6 figures.
dc.description.abstractA block tridiagonal matrix is factored with minimal fill-in using a cyclic reduction algorithm that is easily parallelized. Storage of the factored blocks allows the application of the inverse to multiple right-hand sides which may not be known at factorization time. Scalability with the number of block rows is achieved with cyclic reduction, while scalability with the block size is achieved using multithreaded routines (OpenMP, GotoBLAS) for block matrix manipulation. This dual scalability is a noteworthy feature of this new solver, as well as its ability to efficiently handle arbitrary (non-powers-of-2) block row and processor numbers. Comparison with a state-of-the art parallel sparse solver is presented. It is expected that this new solver will allow many physical applications to optimally use the parallel resources on current supercomputers. Example usage of the solver in magneto-hydrodynamic (MHD), three-dimensional equilibrium solvers for high-temperature fusion plasmas is cited.
dc.description.sponsorshipThis research has been sponsored by the US Department of Energy under Contract DE-AC05-00OR22725 with UT-Battelle, LLC. This research used resources of the National Center for Computational Sciences at Oak Ridge National Laboratory, which is supported by the Office of Science of the Department of Energy under Contract DE-AC05-00OR22725.
dc.description.statusPublicado
dc.format.mimetypeapplication/pdf
dc.identifier.bibliographicCitationJournal of Computational Physics, 2010, v. 229, n. 18, p. 6392-6404
dc.identifier.doihttps://www.doi.org/10.1016/j.jcp.2010.04.049
dc.identifier.issn0021-9991
dc.identifier.urihttps://hdl.handle.net/10016/8910
dc.language.isoeng
dc.publisherElsevier
dc.relation.publisherversionhttp://dx.doi.org/10.1016/j.jcp.2010.04.049
dc.rights© Elsevier
dc.rights.accessRightsopen access
dc.subject.ecienciaFísica
dc.subject.ecienciaFusión
dc.subject.otherCyclic reduction
dc.subject.otherBlock matrix
dc.subject.otherDense blocks
dc.subject.otherTridiagonal matrix
dc.subject.otherThomas algorithm
dc.subject.otherParallel computing
dc.titleBCYCLIC: A parallel block tridiagonal matrix cyclic solver
dc.typeresearch article*
dc.type.reviewPeerReviewed
dspace.entity.typePublication
Files
Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
Name:
bcyclic_sanchez_jcp_2010_ps.pdf
Size:
1.07 MB
Format:
Adobe Portable Document Format
Description:
post-print version
Download
Collections
DF - GFP - Artículos de revistas