Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Improved parallel solution techniques for the integral transport matrix method

Conference ·
OSTI ID:1041203
 [1];  [2]
  1. Los Alamos National Laboratory
  2. 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
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