|May 12, 2003, 02:54||
Implementing an EBE Preconditioner
Having read Saad's "Iterative Methods" and Netlibs "Templates", I'm still unclear on how to implement an Incomplete Lower Upper factorization on an element-by-element basis to precondition (M) a sparse linear system (A).
I'd found coded examples on how to preform an ILU(0) factorization for a general nxn matix. However, I'm not sure how to implement this for iterative methods requiring:
solve w from M*w=A, or w=M^(-1)*v, for example.
How can this be done without assembling a global matrix for the preconditioner M, or then finding the global matrix inverse?
Any help would be much appreciated. Thanks in advance!
|Thread||Thread Starter||Forum||Replies||Last Post|
|thobois class engineTopoChangerMesh error||Peter_600||OpenFOAM||4||August 2, 2014 09:52|
|Help about the preconditioner in gmres method||Dan Gao||Main CFD Forum||3||July 5, 2008 01:18|
|preconditioner of a matrix||khaiching||Main CFD Forum||4||June 18, 2006 20:07|
|Advice on implementing Kass paper||AndyN||Main CFD Forum||0||March 14, 2005 08:41|
|Preconditioner for BiCSTAB||Márcio Ricardo||Main CFD Forum||1||October 12, 2004 07:10|