Parallel simulation of anisotropic diffusion with human brain DT-MRI data

Ning Kang, Jun Zhang, Eric S. Carlson

Research output: Contribution to journalArticlepeer-review

6 Scopus citations

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 languageEnglish
Pages (from-to)2389-2399
Number of pages11
JournalComputers and Structures
Volume82
Issue number28 SPEC. ISS.
DOIs
StatePublished - 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

Fingerprint

Dive into the research topics of 'Parallel simulation of anisotropic diffusion with human brain DT-MRI data'. Together they form a unique fingerprint.

Cite this