Roc Wang <[email protected]> writes: >> What happens with -pc_type lu? Does it solve the linear system and give >> back answers you expect? > > No, it doesn't work if I use -pc_type lu. The error information shows: > > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [0]PETSC ERROR: No support for this operation for this object type! > [0]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc LU!
If you must run in parallel, either use a third-party solver such as -pc_type lu -pc_factor_mat_solver_package superlu_dist (configure with --download-superlu_dist) or use -pc_type redundant (for testing). I recommend getting your solver working in serial first.
pgpOGqXFp3QUW.pgp
Description: PGP signature
