Abstract
We present a multi-threaded solver for symmetric positive definite linear systems where the coefficient matrix of the problem features a bordered-band non-zero pattern. The algorithms that implement this approach heavily rely on a compact storage format, tailored for this type of matrices, that reduces the memory requirements, produces a regular data access pattern, and allows to cast the bulk of the computations in terms of efficient kernels from the Level-3 and Level-2 BLAS. The efficiency of our approach is illustrated by numerical experiments.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
References
Anderson, E., Bai, Z., Bischof, C., Blackford, L.S., Demmel, J., Dongarra, J.J., Du Croz, J., Hammarling, S., Greenbaum, A., McKenney, A., Sorensen, D.: LAPACK Users’ Guide, 3rd edn. Society for Industrial and Applied Mathematics, Philadelphia (1999)
Bollhöfer, M., Kressner, D., Mehl C., Stykel, T., Benner, P. (eds.) Numerical Algebra, Matrix Theory, Differential-Algebraic Equations and Control Theory. Festschrift in Honor of Volker Mehrmann. Springer, Berlin (2015)
Bientinesi, P., Gunnels, J.A., Myers, M.E., Quintana-Ortí, E.S., van de Geijn, R.A.: The science of deriving dense linear algebra algorithms. ACM Trans. Math. Softw. 31(1), 1–26 (2005)
Blobel, V., Kleinwort, C., Meier, F.: Fast alignment of a complex tracking detector using advanced track models. Comput. Phys. Commun. 182(9), 1760–1763 (2011)
de Boor, C., Golub, G.: The numerically stable reconstruction of a Jacobi matrix from spectral data. Linear Algebra Appl. 21(3), 245–260 (1978)
Demmel, J.W.: Applied Numerical Linear Algebra. Society for Industrial and Applied Mathematics, Philadelphia (1997)
Duff, I.S., Scott, J.A.: Stabilized bordered block diagonal forms for parallel sparse solvers. Parallel Comput. 31(3–4), 275–289 (2005)
Golub, G., Welsch, J.H.: Calculation of Gauss quadrature rules. Math. Comput. 23(106), 221–230 (1969)
Intel Corporation. https://software.intel.com/en-us/node/470282
Intel Corporation. http://www.intel.com/
Kleinwort, C.: General broken lines as advanced track fitting method. Nucl. Instrum. Methods Phys. Res. Sect. A Accelerators Spectrometers Detectors Assoc. Equip. 673, 107–110 (2012)
Saak, J.: Effiziente numerische Lösung eines Optimalsteuerungsproblems für die Abkühlung von stahlprofilen. Diplomarbeit, Fachbereich 3/Mathematik und Informatik, Universität Bremen, D-28334 Bremen, September 2003
Van Huffel, S., Park, H.: Efficient reduction algorithms for bordered band matrices. Numer. Linear Algebra Appl. 2(2), 95–113 (1995)
Acknowledgments
The researcher from Universidad Jaime I was supported by the CICYT project TIN2011-23283 of the Ministerio de Economía y Competitividad and FEDER. We would like to thank Matthias Bollhöfer, from TU Braunschweig, for his help during the evaluation of these codes.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2016 Springer International Publishing Switzerland
About this paper
Cite this paper
Benner, P., Ezzatti, P., Quintana-Ortí, E.S., Remón, A. (2016). A Parallel Multi-threaded Solver for Symmetric Positive Definite Bordered-Band Linear Systems. In: Wyrzykowski, R., Deelman, E., Dongarra, J., Karczewski, K., Kitowski, J., Wiatr, K. (eds) Parallel Processing and Applied Mathematics. PPAM 2015. Lecture Notes in Computer Science(), vol 9573. Springer, Cham. https://doi.org/10.1007/978-3-319-32149-3_10
Download citation
DOI: https://doi.org/10.1007/978-3-319-32149-3_10
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-32148-6
Online ISBN: 978-3-319-32149-3
eBook Packages: Computer ScienceComputer Science (R0)