News and Events – NLAFET https://www.nlafet.eu Website for the NLAFET project Thu, 13 Aug 2020 14:18:16 +0000 en-US hourly 1 https://wordpress.org/?v=5.0.19 StarNEig published in Special Issue of CCPE https://www.nlafet.eu/2020/08/13/starneig-published-in-special-issue-of-ccpe/ https://www.nlafet.eu/2020/08/13/starneig-published-in-special-issue-of-ccpe/#respond Thu, 13 Aug 2020 14:16:19 +0000 https://www.nlafet.eu/?p=948 An extended paper describing the task-based StarNEig library for solving nonsymmetric standard and generalized eigenvalue problems has been published in a special issue of the international journal “Concurrency and Computation: Practice and Experience”. The publishers version can be found here. A post print can be found here.

]]>
https://www.nlafet.eu/2020/08/13/starneig-published-in-special-issue-of-ccpe/feed/ 0
New release of the StarNEig library https://www.nlafet.eu/2020/03/26/new-release-of-the-starneig-library/ https://www.nlafet.eu/2020/03/26/new-release-of-the-starneig-library/#respond Thu, 26 Mar 2020 11:51:21 +0000 http://www.nlafet.eu/?p=931 Read more »]]> StarNEig is a new task-based parallel library for solving nonsymmetric standard and generalized eigenvalue problems. A new version (StarNEig 0.1 Beta 4) is now freely available for download from the dedicated website

https://nlafet.github.io/StarNEig/

StarNEig is either comparable to LAPACK and ScaLAPACK or significantly faster depending on the computational step. Moreover, StarNEig realizes new parallel and blocked algorithms for computing eigenvectors without suffering from floating point overflow. In LAPACK the corresponding solvers are sequential scalar codes which compute eigenvectors one by one. In ScaLAPACK the corresponding solvers are vulnerable to overflow.

Eigenvalue problems can be found in every field of natural science. Clear examples are supplied by the analysis of systems of ordinary differential equations. The stability analysis of first order systems produces standard eigenvalue problems which are not necessarily symmetric and the analysis of second order systems produce quadratic eigenvalue problems which are equivalent to nonsymmetric generalized eigenvalue problems. Without the ability to solve eigenvalue problems rapidly and accurately, we would be unable to complete the calculations needed to maintain and advance our civilization. Therefore, it is important that we continue to develop new algorithms and software which maximizes both the performance and the accuracy using existing and emerging hardware

StarNEig is one of very few libraries to offer support for nonsymmetric eigenvalue problems. It is built on top of the runtime system StarPU which is used to schedule the tasks. Currently, StarNEig applies to real problems which have real or complex eigenvalues and eigenvectors. By design, StarNEig applies to both shared and distributed memory machines and it has experimental support for GPU accelerators.

]]>
https://www.nlafet.eu/2020/03/26/new-release-of-the-starneig-library/feed/ 0
Three new papers published https://www.nlafet.eu/2020/03/26/three-new-papers-published/ https://www.nlafet.eu/2020/03/26/three-new-papers-published/#respond Thu, 26 Mar 2020 11:23:44 +0000 http://www.nlafet.eu/?p=924 Read more »]]> Three papers based on research completed at Umeå University during the NLAFET project were presented at the conference PPAM-2019 (Parallel Processing and Applied Mathematics) in Bialystok, Poland in September 2019. They now appear as part of the peer-reviewed conference proceedings published by Springer as LNCS volume 12043. Free preprints are also available through arXiv.

1: Carl Christian Kjelgaard Mikkelsen and Mirko Myllykoski: Parallel Robust Computation of Generalized Eigenvectors of Matrix Pencils
https://link.springer.com/chapter/10.1007/978-3-030-43229-4_6
https://arxiv.org/abs/2003.04776

In this paper, we consider the problem of computing generalized eigenvectors of a matrix pencil in real Schur form. In exact arithmetic, this problem can be solved using substitution. In practice, substitution is vulnerable to floating-point overflow. The robust solvers xtgevc in LAPACK prevent overflow by dynamically scaling the eigenvectors. These subroutines are scalar and sequential codes which compute the eigenvectors one by one. In this paper, we discuss how to derive robust algorithms which are blocked and parallel. The new StarNEig library contains a robust task-parallel solver Zazamoukh which runs on top of StarPU. Our numerical experiments show that Zazamoukh achieves a super-linear speedup compared with dtgevc for sufficiently large matrices.

2: Mirko Myllykoski and Carl Christian Kjelgaard Mikkelsen: Introduction to StarNEig – A Task-based Library for Solving Nonsymmetric Eigenvalue Problems.
https://link.springer.com/chapter/10.1007/978-3-030-43229-4_7
https://arxiv.org/abs/2002.05024

In this paper, we present the StarNEig library for solving dense nonsymmetric (generalized) eigenvalue problems. The library is built on top of the StarPU runtime system and targets both shared and distributed memory machines. Some components of the library support GPUs. The library is currently in an early beta state and only real arithmetic is supported. Support for complex data types is planned for a future release. This paper is aimed at potential users of the library. We describe the design choices and capabilities of the library, and contrast them to existing software such as ScaLAPACK. StarNEig implements a ScaLAPACK compatibility layer that should make it easy for new users to transition to StarNEig. We demonstrate the performance of the library with a small set of computational experiments.

3: Angelika Beatrix Schwarz and Carl Christian Kjelgaard Mikkelsen: Robust Task-Parallel Solution of the Triangular Sylvester Equation.
https://link.springer.com/chapter/10.1007/978-3-030-43229-4_8
https://arxiv.org/abs/1905.10574

The Bartels-Stewart algorithm is a standard approach to solving the dense Sylvester equation. It reduces the problem to the solution of the triangular Sylvester equation. The triangular Sylvester equation is solved with a variant of backward substitution. Backward substitution is prone to overflow. Overflow can be avoided by dynamic scaling of the solution matrix. An algorithm which prevents overflow is said to be robust. The standard library LAPACK contains the robust scalar sequential solver dtrsyl. This paper derives a robust, level-3 BLAS-based task-parallel solver. By adding overflow protection, our robust solver closes the gap between problems solvable by LAPACK and problems solvable by existing non-robust task-parallel solvers. We demonstrate that our robust solver achieves a performance similar to non-robust solvers.

]]>
https://www.nlafet.eu/2020/03/26/three-new-papers-published/feed/ 0
Public deliverables (2.4, 3.7, 4.5) uploaded. https://www.nlafet.eu/2019/04/08/public-deliverables-2-4-3-7-4-5-uploaded/ https://www.nlafet.eu/2019/04/08/public-deliverables-2-4-3-7-4-5-uploaded/#respond Mon, 08 Apr 2019 15:53:17 +0000 http://www.nlafet.eu/?p=843 Public deliverables uploaded:

For the full list of released public deliverables, see this page: http://www.nlafet.eu/public-deliverables/

]]>
https://www.nlafet.eu/2019/04/08/public-deliverables-2-4-3-7-4-5-uploaded/feed/ 0
CCPE – Paper Accepted for special issue https://www.nlafet.eu/2018/10/15/ccpe-paper-accepted-for-special-issue/ https://www.nlafet.eu/2018/10/15/ccpe-paper-accepted-for-special-issue/#respond Mon, 15 Oct 2018 11:38:19 +0000 http://www.nlafet.eu/?p=794 Read more »]]> Paper accepted to a special issue of Concurrency and Computation: Practice and Experience (CCPE). This special issue will consist of extended versions of the best papers accepted by PPAM 2017. The title of the paper is “Parallel robust solution of triangular linear systems”

Robust algorithms do not suffer from overflow and always return a valid result. In LAPACK eigenvectors (standard and generalized) are computed using robust algorithms. The existing algorithms are scalar and sequential. This new work presents algorithms which are blocked and parallel. The analysis is supported by parallel software running on top of StarPU. Further improvements are possible, but the new software is already orders of magnitude faster than the existing software.

The authors are Carl Christian Kjelgaard Mikkelsen, Angelika Beatrix Schwarz and Lars Karlsson

]]>
https://www.nlafet.eu/2018/10/15/ccpe-paper-accepted-for-special-issue/feed/ 0
Working Note 20 have been uploaded https://www.nlafet.eu/2018/10/15/working-note-20-have-been-uploaded/ https://www.nlafet.eu/2018/10/15/working-note-20-have-been-uploaded/#respond Mon, 15 Oct 2018 11:33:55 +0000 http://www.nlafet.eu/?p=792

Working Note 19 have been uploaded:

The full list of published Working Notes can be found on this page: http://www.nlafet.eu/working-notes/

]]>
https://www.nlafet.eu/2018/10/15/working-note-20-have-been-uploaded/feed/ 0
Research poster award for Mawussi Zounon https://www.nlafet.eu/2018/06/28/research-poster-award-for-mawussi-zounon/ https://www.nlafet.eu/2018/06/28/research-poster-award-for-mawussi-zounon/#respond Thu, 28 Jun 2018 13:13:54 +0000 http://www.nlafet.eu/?p=780 Research poster award for Mawussi Zounon on exploiting the low precision floating-point arithmetic available on NVIDIA V100 GPU in solving linear system of equations.

The title of the poster is “Using GPU’s FP16 Tensor Cores Arithmetic to Accelerate Mixed-Precision Iterative Refinement Solvers and Reduce Energy Consumption” with the authors Azzam Haidar, Stan Tomov, Ahmad Abdelfattah, Mawussi Zounon, Jack Dongarra

]]>
https://www.nlafet.eu/2018/06/28/research-poster-award-for-mawussi-zounon/feed/ 0
Public Deliverables uploaded https://www.nlafet.eu/2018/04/18/public-deliverables-uploaded-2/ https://www.nlafet.eu/2018/04/18/public-deliverables-uploaded-2/#respond Wed, 18 Apr 2018 15:28:57 +0000 http://www.nlafet.eu/?p=644

Several deliverables have been uploaded: D3.2, D3.6, D4.3, D6.2, and D6.6.

A full list, including links to all the public deliverables that have been posted until now can be found here: http://www.nlafet.eu/public-deliverables/

]]>
https://www.nlafet.eu/2018/04/18/public-deliverables-uploaded-2/feed/ 0
SIAM PP18 – SIAM Conference on Parallel Processing for Scientific Computing https://www.nlafet.eu/2018/04/13/pp18-siam-conference-on-parallel-processing-for-scientific-computing/ https://www.nlafet.eu/2018/04/13/pp18-siam-conference-on-parallel-processing-for-scientific-computing/#respond Fri, 13 Apr 2018 12:31:04 +0000 http://www.nlafet.eu/?p=516 Read more »]]> NLAFET sent several participants to SIAM PP18 (March 7-10, 2018) at Waseda University in Tokyo, Japan, where they presented some of their ongoing research.  Persons in italic made the presentation at the conference.

Talks

  • Dense Linear Systems for Extreme Scale  Jack J. Dongarra, University of Tennessee, ORNL, and University of Manchester, USA
  • The Batched Blas abstract Sven J. Hammarling and Mawussi Zounon, The University of Manchester, UK
  • Task Based Robust Computation of Eigenvectors abstract Carl Christian Kjelgaard Mikkelsen, Lars Karlsson, Mirko Myllykoski, and Angelika Schwarz, Umeå University, Sweden
  • Sparse Solvers on Extreme-Scale Systems abstract Iain Duff, Science & Technology Facilities Council, United Kingdom and CERFACS, Toulouse, France
  • Task-Based Sparse Direct Solver for Symmetric Indefinite Systems abstract Florent Lopez, Rutherford Appleton Laboratory, United Kingdom; Iain Duff, Science & Technology Facilities Council, United Kingdom and CERFACS, Toulouse, France
  • Robust Algebraic Preconditioners for Large Scale Applications abstract Laura Grigori, Simplice Donfack, and Olivier Tissot, Inria, France
  • Enlarged Conjugate Gradient Method for Reducing Communication abstract Olivier Tissot and Laura Grigori, Inria, France

Posters

  • Efficient Robust Multi-Shift Triangular Solves abstract Angelika Schwarz, Mirko Myllykoski, Carl Christian Kjelgaard Mikkelsen, and Lars Karlsson, Umeå University, Sweden
]]>
https://www.nlafet.eu/2018/04/13/pp18-siam-conference-on-parallel-processing-for-scientific-computing/feed/ 0
Working Note 19 has been uploaded https://www.nlafet.eu/2018/04/13/working-note-19-have-been-uploaded/ https://www.nlafet.eu/2018/04/13/working-note-19-have-been-uploaded/#respond Fri, 13 Apr 2018 11:50:25 +0000 http://www.nlafet.eu/?p=514 Working Note 19 have been uploaded:

The full list of published Working Notes can be found on this page: http://www.nlafet.eu/working-notes/

]]>
https://www.nlafet.eu/2018/04/13/working-note-19-have-been-uploaded/feed/ 0