Linear System solver in Newton method
Hi there,
I am solving an equation with a Newton method. I have a very illconditioned Jacobian matrix with a huge matrix dimension (>10000). The matrix is sparse but not symmetric. At the moment I am using a GMRES (with no restarts because restarts take more time than a larger Krylov space memorization) and an ILU(0) preconditioning. What do you suggest me to try for reduce the computing cost? Thanks in advance and sorry for the beginner's question. Bye, Ale 
Quote:
Really sorry for my bad posed question. It was not so much advanced and your suggest in use direct solvers is enough at the moment for my aims. Thanks a lot for your suggests on how to get a useful answer too. 
No problem, that part of my post really wasn't intended for you. Hopefully the next person will include this sort of information so that they can get a useful answer.

All times are GMT 4. The time now is 16:54. 