Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- Pardiso result problem

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

danielsue

Beginner

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

05-24-2013
05:22 PM

47 Views

Pardiso result problem

Hi All,

I have implemented pardiso solver in our model for flow simulation. Pardiso can solve most of our problem well and produce correct results. But recently I have a case that pardiso generates quite different results with our origional solver (ws209), which has been used for over ten years for our model.I guess if there is something wrong with the setting.

Please find the test codes (pardiso_unsym_f.f90), sparse matrices exported from our model (a_i.txt, b_i.txt, ia_i.txt, ja_i.txt), result generated by pardiso (x_out_i.txt) and result generated by ws209 solver (x_ws209_i.txt). The first value in each file is the number of values.

Can anybody help to have a check?

Thanks and regards,

Daniel

Link Copied

2 Replies

Gennady_F_Intel

Moderator

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

05-25-2013
10:13 AM

47 Views

Daniel,

I don't see big differences between x_ws209_1.txt x_out1.txt. Do you estimate the condition number of these matrixes?

danielsue

Beginner

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

05-29-2013
05:33 PM

47 Views

Gennady Fedorov (Intel) wrote:

Daniel,

I don't see big differences between x_ws209_1.txt x_out1.txt. Do you estimate the condition number of these matrixes?

Hi Gennady,

Thanks. I have estimated the condition number, these matrix are ill-conditioned. The condition number is around 1.0E15. The origional solver of our codes can solve but Pardiso can not produce correct result.

I know Pardiso cannot deal with singular matrix, but wonder if it can deal with ill-conditioned matrix?

Thanks,

Daniel

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

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