Matrix solver for Impicit scheme
Hello, I am developing an implicit reimann solver for Euler equations and i am using GMRES algorithm to solve the matrix, but i find that Implicit is taking more CPU time than explicit solver and may be the problem is with the matrix solver. I need some suggestions regarding it as which is the best solver for the implicit schemes and sources from which i can get it?

Take advantage of sparse system. Implement GMRES for sparse system.
GMRES can even solve indefinite system. This much advanced solver in not required for Euler equations. Lot of simple, robust solvers are available and they are fast enough. I recommend sparse PCG solver. 
@duri: thanx a lot for the reply...My guide actually suggested me to try out Gaussseidel point relaxation method..but d problem is i am not able to find the readymade solver..can you suggest me some source from which i can get these readymade PCG or gaussseidel point relaxation solvers n possibly in C language.

Dear Harish,
A good opensource package with a large number of builtin linear solvers is LiS. You can find it with documentation at www.ssisc.org/lis. It supports both F90 and C programming languages, but there is not much usersupport available. The manual however is selfexplanatory and should do the job. Regards, Ganesh 
Thank u very much ganesh. I found it and seems to be helpful..i will try it out.

@ganesh: I an not able to understand how to integrate the linear solvers you mentioned into my code. Can u please just brief me how to do so..

All times are GMT 4. The time now is 22:36. 