Bruno Ramon Batista Fernande
2015-10-07 19:13:29 UTC
Hello guys,
I am using petsc is a solver for a finite difference PDE solver.
Basically the coefficient matrix (say A in Ax = b) is well conditioned,
symmetric, and up to 7 diagonals. Also, no pivoting is necessary for such
problem and it is always diagonally dominant.
I am using all default options of petsc but I still feel like I could get
better results by choosing different answers. By default I mean: GMRES
solver preconditioned on left with ILU(0) and a classical Gram-Schmidt to
orthogonalize against krylov vector.
I am using MAT_STRUCTURALLY_SYMMETRIC option and
atol = 1.0e-30
dtol = 1.0e5
You guys have any idea of how to further improve performance of the solver?
Regards,
Bruno
I am using petsc is a solver for a finite difference PDE solver.
Basically the coefficient matrix (say A in Ax = b) is well conditioned,
symmetric, and up to 7 diagonals. Also, no pivoting is necessary for such
problem and it is always diagonally dominant.
I am using all default options of petsc but I still feel like I could get
better results by choosing different answers. By default I mean: GMRES
solver preconditioned on left with ILU(0) and a classical Gram-Schmidt to
orthogonalize against krylov vector.
I am using MAT_STRUCTURALLY_SYMMETRIC option and
atol = 1.0e-30
dtol = 1.0e5
You guys have any idea of how to further improve performance of the solver?
Regards,
Bruno