Abstract
We conduct simulations for the 3D unsteady state anisotropic diffusion process with DT-MRI data in the human brain by discretizing the governing diffusion equation on Cartesian grid and adopting a high performance differential-algebraic equation (DAE) solver, the parallel version of implicit differential-algebraic (IDA) solver, to tackle the resulting large scale system of DAEs. Parallel preconditioning techniques including sparse approximate inverse and banded-block-diagonal preconditioners are used with the GMRES method to accelerate the convergence rate of the iterative solution. We then investigate and compare the efficiency and effectiveness of the two parallel preconditioners. The experimental results of the diffusion simulations on a parallel supercomputer show that the sparse approximate inverse preconditioning strategy, which is robust and efficient with good scalability, gives a much better overall performance than the banded-block-diagonal preconditioner.
Original language | English |
---|---|
Pages (from-to) | 2389-2399 |
Number of pages | 11 |
Journal | Computers and Structures |
Volume | 82 |
Issue number | 28 SPEC. ISS. |
DOIs | |
State | Published - Nov 2004 |
Bibliographical note
Funding Information:This study was funded by the US Department of Energy Office of Science under the project “Development of a High Performance Anisotropic Diffusion Equation Solver Using the ACTS Toolkit” (DE-FG02-02ER45961). The first two authors would also like to acknowledge DOE’s support of their attending the ACTS Workshop, organized by Drs. Tony Drummond and Osni Marques, at the Lawrence Berkeley National Laboratory, in September 2002. We would also like to thank Dr. Daniel Gembris, at Institute for Medicine, Jülich Research Center, Jülich, Germany, for providing the diffusion tensor data set.
Keywords
- Anisotropic diffusion simulation
- Banded-block-diagonal preconditioner
- Diffusion tensor magnetic resonance imaging (DT-MRI)
- Implicit differential-algebraic solver
- Parallel preconditioning techniques
- Sparse approximate inverse preconditioner
ASJC Scopus subject areas
- Civil and Structural Engineering
- Modeling and Simulation
- General Materials Science
- Mechanical Engineering
- Computer Science Applications