Improved parallel solution techniques for the integral transport matrix method
- 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), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC52-06NA25396
- OSTI ID:
- 1041203
- Report Number(s):
- LA-UR-10-07886; LA-UR-10-7886; TRN: US201211%%483
- Resource Relation:
- Conference: ANS International Conference on Mathematics and Computation (M&C) 2011 ; May 8, 2011 ; Rio de Janeiro, RJ, Brazil
- Country of Publication:
- United States
- Language:
- English
Similar Records
Implementation of Generalized Coarse-Mesh Rebalance of NEWTRNX for Acceleration of Parallel Block-Jacobi Transport
Iterative stability analysis of spatial domain decomposition based on block Jacobi algorithm for the diamond-difference scheme