In this paper we analyze the gauss-huard algorithm. From a description of the algorithm in terms of matrix-vector operations we reveal a close relation between the gauss-huard algorithm and an LU factorization as cons...
详细信息
In this paper we analyze the gauss-huard algorithm. From a description of the algorithm in terms of matrix-vector operations we reveal a close relation between the gauss-huard algorithm and an LU factorization as constructed in an ikj variant. (C) 1998 Elsevier Science Inc. All rights reserved.
In 1979, P. huard presented an efficient variant of the gauss-Jordan elimination for the solution of linear systems. In particular, this alternative algorithm exhibits the same computational cost as the traditional LU...
详细信息
ISBN:
(纸本)9783319321493;9783319321486
In 1979, P. huard presented an efficient variant of the gauss-Jordan elimination for the solution of linear systems. In particular, this alternative algorithm exhibits the same computational cost as the traditional LU-based solver, and is considerably cheaper than the gauss-Jordan algorithm, but there exist no recent high performance implementations of the gauss-huard (GH) variant that allow a comparison of these approaches. In this paper we present a reliable GH solver for hybrid platforms equipped with conventional multi-core technology and a graphics processing unit (GPU). The experimental results show that the GH algorithm can beat high performance versions of the LU solver, from tuned libraries for CPU-GPU servers such as MAGMA, for problems of small to moderate scale.
The gauss-huard algorithm (the GHA) is a specialized version of gauss-Jordan elimination for the solution of linear systems that, enhanced with column pivoting, exhibits numerical stability and computational cost clos...
详细信息
ISBN:
(纸本)9783319495835;9783319495828
The gauss-huard algorithm (the GHA) is a specialized version of gauss-Jordan elimination for the solution of linear systems that, enhanced with column pivoting, exhibits numerical stability and computational cost close to those of the conventional solver based on the LU factorization with row pivoting. Furthermore, the GHA can be formulated as a procedure rich in matrix multiplications, so that high performance can be expected on current architectures with multi-layered memories. Unfortunately, in principle the GHA does not admit the introduction of look-ahead, a technique that has been demonstrated to be rather useful to improve the performance of the LU factorization on multi-threaded platforms with high levels of hardware concurrency. In this paper we analyze the effect of this drawback on the implementation of the GHA on systems accelerated with graphics processing units (GPUs), exposing the roles of the CPU-to-GPU and single precision-to-double precision performance ratios, as well as the contribution from the operations in the algorithm's critical path.
The solution of linear systems is a key operation in many scientific and engineering applications. Traditional solvers are based on the LU factorization of the coefficient matrix, and optimized implementations of this...
详细信息
ISBN:
(纸本)9783319269283;9783319269276
The solution of linear systems is a key operation in many scientific and engineering applications. Traditional solvers are based on the LU factorization of the coefficient matrix, and optimized implementations of this method are available in well-known dense linear algebra libraries for most hardware architectures. The gauss-huard algorithm (GHA) is a reliable and alternative method that presents a computational effort close to that of the LU-based approach. In this work we present several implementations of GHA on the Intel Xeon Phi coprocessor. The experimental results show that our solvers based in GHA represent a competitive alternative to LU-based solvers, being an appealing method for the solution of small to medium linear systems, with remarkable reductions in the time-to-solution for systems of dimension n <= 4,000.
暂无评论