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
- solve ill-conditioned linear system

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

Antoine__A_

Beginner

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

01-14-2013
07:16 AM

84 Views

solve ill-conditioned linear system

Hello,

I would like to solve a linear system Ax = b with :

int ia[10] = {1, 5, 9, 11, 15, 19, 21, 23, 25, 26};

int ja[25] = {1, 2, 4, 5, 2, 3, 5, 6, 3, 6, 4, 5, 7, 8, 5, 6, 8, 9, 6, 9, 7, 8, 8, 9, 9};

double a[25] = {

40.2991, 2.5049, 0.929679, 2.44749,

40.2991, 2.5049, 0.929679, 2.44749,

40.2991, 0.929679,

40.2991, 2.5049, 0.929679, 2.44749,

40.2991, 2.5049, 0.929679, 2.44749,

40.2991, 0.929679,

40.2991, 2.5049,

40.2991, 2.5049,

40.2991

};

double rhs[9] = {3.40521, 3.55851, 30.0158, 39300.2, 41805.2, 2506.73, 930.681, 3378.17, 2448.49};

but it returns me a wrong result.

Is there an iparm for detecting ill-conditionned linear systems (to make me sure this is the problem), and in this case how can I solve this problem?

Thanks for any help you can give me.

Link Copied

5 Replies

mecej4

Black Belt

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

01-14-2013
09:04 AM

84 Views

By adapting the provided example dss_sym_c.c in the MKL examples/solverc directory to your example problem above, I found that the residuals never exceed 7e-12 in magnitude. Therefore, I cannot agree that the results are "wrong", nor do I think that ill-conditioning is much of a problem.

Perhaps, you have some errors in the arguments that you passed to the MKL solver routines.

You did not state how you solved the problem, which MKL routines you used, and what the results were. Please provide more details, and state the platform, compiler and MKL versions used.

Antoine__A_

Beginner

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

01-15-2013
01:08 AM

84 Views

I use MKL's Pardiso. The Package ID of mkl is : l_mkl_p_10.0.011.

My compiler is icc and my platform is Redhat.I use the solver's default parameters.

Pardiso returns me this vector :

sol =

-81.5585

-27.2802

-0.657835

906.03

1028.06

60.8015

18.1186

80.0508

60.7579

Alexander_K_Intel2

Employee

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

01-15-2013
01:34 AM

84 Views

Hi Antonie,

Nice to see you here again :) As we can see your matrix has diagonal dominant, so, if you use symmetric type you matrix is positive define. Could provide your input parameters?

With best regards,

Alexander Kalinkin

Antoine__A_

Beginner

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

01-15-2013
06:06 AM

84 Views

Sorry for loosing your time, I made a mistake in checking my result

Sorry

with best regards,

Antoine

Alexander_K_Intel2

Employee

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

01-15-2013
10:01 AM

84 Views

Hi Antoine,

No problem, feel free to ask any questions!

With best regards,

Alexander Kalinkin

- 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.