Improved parallel solution techniques for the integral transport matrix method
Conference
·
OSTI ID:1041203
- Los Alamos National Laboratory
- NORTH CAROLINA STATE UNIV.
Alternative solution strategies to the parallel block Jacobi (PBJ) method for the solution of the global problem with the integral transport matrix method operators have been designed and tested. The most straightforward improvement to the Jacobi iterative method is the Gauss-Seidel alternative. The parallel red-black Gauss-Seidel (PGS) algorithm can improve on the number of iterations and reduce work per iteration by applying an alternating red-black color-set to the subdomains and assigning multiple sub-domains per processor. A parallel GMRES(m) method was implemented as an alternative to stationary iterations. Computational results show that the PGS method can improve on the PBJ method execution by up to {approx}50% when eight sub-domains per processor are used. However, compared to traditional source iterations with diffusion synthetic acceleration, it is still approximately an order of magnitude slower. The best-performing case are opticaUy thick because sub-domains decouple, yielding faster convergence. Further tests revealed that 64 sub-domains per processor was the best performing level of sub-domain division. An acceleration technique that improves the convergence rate would greatly improve the ITMM. The GMRES(m) method with a diagonal block preconditioner consumes approximately the same time as the PBJ solver but could be improved by an as yet undeveloped, more efficient preconditioner.
- Research Organization:
- Los Alamos National Laboratory (LANL)
- Sponsoring Organization:
- DOE
- DOE Contract Number:
- AC52-06NA25396
- OSTI ID:
- 1041203
- Report Number(s):
- LA-UR-10-07886; LA-UR-10-7886
- Country of Publication:
- United States
- Language:
- English
Similar Records
Iterative stability analysis of spatial domain decomposition based on block Jacobi algorithm for the diamond-difference scheme
AZTEC: A parallel iterative package for the solving linear systems
Implementation of Generalized Coarse-Mesh Rebalance of NEWTRNX for Acceleration of Parallel Block-Jacobi Transport
Journal Article
·
Wed May 27 20:00:00 EDT 2015
· Journal of Computational Physics
·
OSTI ID:1437429
AZTEC: A parallel iterative package for the solving linear systems
Conference
·
Mon Dec 30 23:00:00 EST 1996
·
OSTI ID:433374
Implementation of Generalized Coarse-Mesh Rebalance of NEWTRNX for Acceleration of Parallel Block-Jacobi Transport
Conference
·
Sun Dec 31 23:00:00 EST 2006
·
OSTI ID:958802