Dear all, I need to solve Poisson's equation many times with same Matrix A and different RHS b. The most efficient way is do pre-factorization I think. For serial problem, the PCLU seems work as I want to. However, how to do for parallel situation? Should I use an external LU solver such as superlu_dist? Thanks Gong Ding