[petsc-users] tridiagonal matrix in kspsolve
Shriram Srinivasan
shriram at ualberta.ca
Tue Jun 10 17:54:40 CDT 2014
Apologies if my question is answered elsewhere already. I could not
find this discussed in the manual or mailing list.
I have been using matsolvermumps through ksp to solve my linear system.
I wanted to simply get it working so I was not concerned about the
structure of the matrix.
1) My matrix is tridiagonal, but the elements of the diagonals are not
all equal. Will a dedicated tridiagonal solver be faster than
matsolvermumps in this case ?
2) My matrix is tridiagonal but not in the conventional sense, the minor
diagonals are farther apart. But it is possible to renumber the nodes
and make it a conventional tridiagonal system. In this case, again, am I
losing efficiency using matsolvermumps ?
Please advise if you recommend a completely different course of action.
Thanks,
Shriram
More information about the petsc-users
mailing list