I am working on finite element discretization. I want to explore the possibility of using multigrid solvers and preconditioners. Are there any free source codes available which uses multigrid solvers for any general sparse matrix.
Re: multigrid solver
I'm not sure that multigrid actually works for "any sparse" matrix (this is definitely true for SIP type algorithms for example). The reason for this is that prolongation and smoothing operators rely on the smoothness of the solution which is OK for a matrix obtained via FEM and FDM but not in general.
The source code and a description of multigrid can be found in the Numerical Recipes books (FORTRAN 77 & 90, C). The books can be found by doing a web search. The source code from the examples in Peric and Ferziger can also be found on the web.
|All times are GMT -4. The time now is 10:23.|