Hi,

just a little question; for a stokes solver,

In almost all the bibliographical FEM articles, I have seen that is mentioned always the GMRES and the BICGSTAB methods as solvers, but in another book I had found the algorithms Orthomin and Orthores, which I did not know them very well, but its algorithms are easy to implement.

Someone can tell me if they are a good alternative to the previous ones, and if they work for not non-symmetric problems?

Thank you very much,

Carlos

A good rule of thumb that I used:

1) Symmetric system with a dominant diagonal: conjugate gradient or orthomin2

2) Non-Symmetric system with a dominant diagonal: orthomin2

3) else: gmres

For the code I wrote when completing me thesis, I always used that rule and it never fails. ( I must be lucky...)