- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
I have a system of equations I am trying to solve using the FGMRES sparse solver (i.e.- Ax=b). From the system I am solving, I have found in the past that a diagonal preconditioner is optimal (where by a diagonal preconditoner, I just mean the diag. elements of A, and solve P^-1Ax=P^-1b).
I could easily construct this preconditioning matrix for use with the FGMRES solver, but am wondering if there is a built in way to pull this diagonal matrix out? I have noticed that there is a ILUT and ILU0 preconditioner available, and seeing as diagonal preconditioning is quite common, I thought there might be a built in routine to do this as well. (otherwise, I can just write my own).
Thanks!
I could easily construct this preconditioning matrix for use with the FGMRES solver, but am wondering if there is a built in way to pull this diagonal matrix out? I have noticed that there is a ILUT and ILU0 preconditioner available, and seeing as diagonal preconditioning is quite common, I thought there might be a built in routine to do this as well. (otherwise, I can just write my own).
Thanks!
Link Copied
2 Replies
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hello,
MKL Sparse BLAS contain a set of triangular solvers for6 most popular sparse formats (compressed sparse row, compressed sparse columns, block sparse row, coordinate, diagonal and skyline formats). You canperform diagonal preconditioning if you call one of routines mkl_dcsrsv, mkl_dcscsv, mkl_dcoosv and etc depending on sparse representation of your coefficient matrix and set the first parameter of descriptor array matdescra to 'D' or 'd'. The routine automatically extracts the main diagonal of coefficient matrix and performs diagonal preconditioning. Please see the MKL Reference Manual for more details.
All the best
Sergey
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Thanks for the info Sergey! It all looks pretty strait-forward now.
-Mark
-Mark
Reply
Topic Options
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page