This paper focuses on the application level improvements in a sparse direct solver specifically used for large-scale unsymmetrical linear equations resulting from unstructured mesh discretization of coupled elliptic/h...
详细信息
This paper focuses on the application level improvements in a sparse direct solver specifically used for large-scale unsymmetrical linear equations resulting from unstructured mesh discretization of coupled elliptic/hyperbolic PDEs. Existing sparse direct solvers are designed for distributed server systems taking advantage of both distributed memory and processing units. We conducted extensive numerical experiments with three state-of-the-art direct linear solvers that can work on distributed-memory parallel architectures;namely, MUMPS(MUMPS solver website, http://***/MUMPS), WSMP (Technical Report TR RC-21886, IBM, Watson Research Center, Yorktown Heights, 2000), and SUPERLU_DIST (ACM Trans Math Softw 29(2): 110-140, 2003). The performance of these solvers was analyzed in detail, using advanced analysis tools such as Tuning and Analysis Utilities (TAU) and Performance Application Programming Interface (PAPI). The performance is evaluated with respect to robustness, speed, scalability, and efficiency in CPU and memory usage. We have determined application level issues that we believe they can improve the performance of a distributed-shared memory hybrid variant of this solver, which is proposed as an alternative solver [SuperLU_MCDT (Many-Core Distributed)] in this paper. The new solver utilizing the MPI/OpenMP hybrid programming is specifically tuned to handle large unsymmetrical systems arising in reservoir simulations so that higher performance and better scalability can be achieved for a large distributed computing system with many nodes of multicore processors. Two main tasks are accomplished during this study: (i) comparisons of public domain solver algorithms;existing state-of-the-art direct sparse linear system solvers are investigated and their performance and weaknesses based on test cases are analyzed, (ii) improvement of direct sparse solver algorithm (SuperLU_MCDT) for many-core distributed systems is achieved. We provided results of numer
We present a method to extract factors of multivariate polynomials with complex coefficients in floating point arithmetic. We establish the connection between the reciprocal of a multivariate polynomial and its Taylor...
详细信息
ISBN:
(纸本)9781605586649
We present a method to extract factors of multivariate polynomials with complex coefficients in floating point arithmetic. We establish the connection between the reciprocal of a multivariate polynomial and its Taylor expansion. Since the multivariate Taylor coefficients are determined by the irreducible factors of the given polynomial, we reconstruct the factors from the Taylor expansion. As each irreducible factor, regardless of its multiplicity, can be separately extracted, our method can lead toward the complete numerical factorization of multivariate polynomials.
We derive a class of iterative formulae to find numerically a factor of arbitrary degree of a polynomial f(x) based on the rational Hermite interpolation. The iterative formula generates the sequence of polynomials wh...
详细信息
暂无评论