Comparison of Parallel Preconditioners in Anisotropic Diffusion Simulation with Human Brain DT-MRI Data

Ning Kang, Jun Zhang
Laboratory for High Performance Scientific Computing and Computer Simulation
Department of Computer Science
University of Kentucky
773 Anderson Hall
Lexington, KY 40506-0046, USA

and
Eric S. Carlson
Department of Chemical Engineering
University of Alabama
P. O. Box 870203
Tuscaloosa, AL 35487-0203, USA

Abstract

We conduct simulations for the 3D unsteady state anisotropic diffusion process 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 precond itioners 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 computational 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.


Key words: anisotropic diffusion simulation, diffusion tensor imaging, differential-algebraic equation, parallel preconditioning techniques

Mathematics Subject Classification:


Download the compressed postscript file brainp.ps.gz, or the PDF file brainp.pdf.
Technical Report No. 376-03, Department of Computer Science, University of Kentucky, Lexington, KY, 2003.
This paper has been accepted for publication in International Journal of Computers and Structures under the title "Parallel simulation of anisotropic diffusion with human brain DT-MRI data".

This research was supported in part by the U.S. National Science Foundation under the grant CCR-9988165, CCR-0092532, and ACR-0202934, in part by the U.S. Department of Energy Office of Science under grant DE-FG02-02ER45961, in part by the Kentucky Science and Engineering Foundation under grant KSEF-02-264-RED-002, in part by the Japanese Research Organization for Information Science & Technology, and in part by the University of Kentucky Research Committee.