CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Question on ILU(0) Preconditioned GMRES

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 5, 2012, 15:21
Default Question on ILU(0) Preconditioned GMRES
  #1
Member
 
Shenren Xu
Join Date: Jan 2011
Location: London, U.K.
Posts: 67
Rep Power: 15
Shenren_CN is on a distinguished road
Dear colleagues,

I have a few questions concerning using ILU(0) for preconditioning the GMRES solver.

As far as I know, ILU(0) decomposes my original system matrix A into

A=L*U+R

where L and R are lower and upper triangular matrices with the same non-zero pattern as my sparse matrix A.

The original linear system

A *x = b

is then preconditioned as

A_new * x = b_new

where A_new=inv(U)*inv(L)*A and b_new = inv(U)*inv(L)*b. The new linear system is solved using GMRES.

My questions are

1) Do I need to explicitly invert L and U? If yes, how do I invert these sparse matrix? Wouldn't the non-zero pattern be destroyed then?

2) In matlab gmres code, one calls it with gmres(... A,..., L,U,...), i.e., only L and U are needed. Then how does matlab invert them, if ever?

3) Anyone has experience on using LASPACK's ILU preconditioned GMRES? Looks like it's ILU is only implemented for symmetric matrix.

Regards,
Shenren
Shenren_CN is offline   Reply With Quote

Old   August 6, 2012, 12:20
Default
  #2
Super Moderator
 
Praveen. C
Join Date: Mar 2009
Location: Bangalore
Posts: 343
Blog Entries: 6
Rep Power: 18
praveen is on a distinguished road
No, you never invert them. Since they are lower/upper matrices, you can solve using forward/backward substitution.
praveen is offline   Reply With Quote

Old   August 6, 2012, 20:19
Default
  #3
Member
 
Shenren Xu
Join Date: Jan 2011
Location: London, U.K.
Posts: 67
Rep Power: 15
Shenren_CN is on a distinguished road
Quote:
Originally Posted by praveen View Post
No, you never invert them. Since they are lower/upper matrices, you can solve using forward/backward substitution.

Thank you Praveen!
Shenren_CN is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Question Re Engineering Data Source imnull ANSYS 0 March 5, 2012 14:51
internal field question - PitzDaily Case atareen64 OpenFOAM Running, Solving & CFD 2 January 26, 2011 16:26
ILU(0) Praveen. C Main CFD Forum 3 April 16, 2007 07:33
Preconditioned Gmres carlos Main CFD Forum 2 April 14, 2003 06:24
Nonlinear/Linear GMRES is... Saebyuk Main CFD Forum 0 May 10, 1999 09:31


All times are GMT -4. The time now is 00:53.