Hello,
Because of the throuble of PETSC, i want to write a solver by myself, I
try to add a gassian elimination to my
RKDG code and it works very fast, the parallel is only in the
system.update().
Because the PETSC in libmesh solve a big matrix and make the ILU(0)
precondition in the local processor, I think
If I write a GMRES+ILU(0) solver only make the local matrix(mesh in each
processor) and exchange information
by system.update(), Are these two methods equal? I think because the
precondition in the local processor
also cast off some offset information so the local GMRES solver can get
the same result. Can you give some
suggestion about this? Can libmesh run well without PETSC but only mpi?
Thank you very much!
Luyi
2008-5-13
-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2008.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
