## Abstract

The linear solve problems arising in chemical physics and many other fields involve large sparse matrices with a certain block structure, for which special block Jacobi preconditioners are found to be very efficient. In two previous papers [W. Chen, B. Poirier, Parallel implementation of efficient preconditioned linear solver for grid-based applications in chemical physics. I. Block Jacobi diagonalization, J. Comput. Phys. 219 (1) (2006) 185197; W. Chen, B. Poirier, Parallel implementation of efficient preconditioned linear solver for grid-based applications in chemical physics. II. QMR linear solver, J. Comput. Phys. 219 (1) (2006) 198209], a parallel implementation was presented. Excellent parallel scalability was observed for preconditioner construction, but not for the matrixvector product itself. In this paper, we introduce a new algorithm with (1) greatly improved parallel scalability and (2) generalization for arbitrary number of nodes and data sizes.

Original language | English |
---|---|

Pages (from-to) | 779-782 |

Number of pages | 4 |

Journal | Journal of Parallel and Distributed Computing |

Volume | 70 |

Issue number | 7 |

DOIs | |

State | Published - Jul 2010 |

## Keywords

- Block Jacobi
- Chemical physics
- Eigensolver
- Linear solver
- Matrixvector product
- Parallel computing
- Preconditioning
- Sparse matrix