Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- Constant residual in GMRES

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

Highlighted
##

Hi all

I am trying to solve a sparse system of linear equations using GMRES (7600 x 7600). Ihave not done any preconditioning yet. The problem I am facing is that after a certain number of iterations, the residual between two consecutive steps does not change. Well, it does change but very very slowly. I am talking about 1d-12 or smaller per iteration.

The system of equations solved nicely with PARDISO. As I need to solve a large system later on, I really need to get GMRES to work. Does anyone know any remedy to this problem? Another thing I should mention is that the iterative solver does not stop even when the maximum number of iterations specified has been exceeded. It just kept on iterating. Why is this so?

One final thing, does the arrangement of the system of equations affect the convergence of the solver? I have tried arranging it both ways and both of them gave the same outcome.

Thank you.

Best wishes,

EH

eh4

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-13-2011
11:58 PM

9 Views

Constant residual in GMRES

I am trying to solve a sparse system of linear equations using GMRES (7600 x 7600). Ihave not done any preconditioning yet. The problem I am facing is that after a certain number of iterations, the residual between two consecutive steps does not change. Well, it does change but very very slowly. I am talking about 1d-12 or smaller per iteration.

The system of equations solved nicely with PARDISO. As I need to solve a large system later on, I really need to get GMRES to work. Does anyone know any remedy to this problem? Another thing I should mention is that the iterative solver does not stop even when the maximum number of iterations specified has been exceeded. It just kept on iterating. Why is this so?

One final thing, does the arrangement of the system of equations affect the convergence of the solver? I have tried arranging it both ways and both of them gave the same outcome.

Thank you.

Best wishes,

EH

1 Reply

Highlighted
##

This post should be moved to the MKL forum.

If you are able to post a small example code that displays the problematic behavior that you described, that code would facilitate identifying the sources of the problems.

The first suspect in such cases, where a complex library is used, is errors in the subroutine argument types/values and proper sequencing of calls.

mecej4

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-14-2011
04:56 AM

9 Views

If you are able to post a small example code that displays the problematic behavior that you described, that code would facilitate identifying the sources of the problems.

The first suspect in such cases, where a complex library is used, is errors in the subroutine argument types/values and proper sequencing of calls.

For more complete information about compiler optimizations, see our Optimization Notice.