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

Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- using pardiso with multiple right hand sides

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

gitsnedbutzi

Beginner

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

06-10-2010
06:57 AM

60 Views

I need to set up pardiso for multiple right hand sides. However I failed to do so, because I'm not really a great programmer. For me, it's not quiet clear, in which form x and b have to be.

I tried a two-dimensional dynamic array but it didn't work.

double **x = new double *

for( int j = 0 ; j < n ; j++ )

{

x

}

Is it supposed to bea vector ina row major format?

Thanks for your help!!

Link Copied

Accepted Solutions

Alexander_K_Intel2

Employee

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

06-10-2010
09:11 AM

60 Views

Similar with one right hand side x and b are one dimension double * arrays but therelength become in nrhs (number of rhs) time larger (length of x and b nrhs*n, where n - dimension of problem). For more details please read

IntelMKLibrary ReferenceManual Chapter 8Sparse Solver Routines.

With best regards,

Alexander Kalinkin

5 Replies

Alexander_K_Intel2

Employee

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

06-10-2010
09:11 AM

61 Views

Similar with one right hand side x and b are one dimension double * arrays but therelength become in nrhs (number of rhs) time larger (length of x and b nrhs*n, where n - dimension of problem). For more details please read

IntelMKLibrary ReferenceManual Chapter 8Sparse Solver Routines.

With best regards,

Alexander Kalinkin

gitsnedbutzi

Beginner

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

06-11-2010
07:02 AM

60 Views

Now I have another problem:

My solver works fine for all test examples. However when I try a larger matrix from my FE code, it can get through everything (analysis, factorization)up to the solution phase. At the solution phase (phase 33)pardiso crashes with an access violation when reading.

Any ideas on what might be wrong, or how to check what the problem is?

Other solves work fine on the matrix...

Thanks

Alexander_K_Intel2

Employee

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

06-11-2010
09:58 AM

60 Views

Such situation appears only with many right hand side or with 1 rhs too? Could you describe problem deeply or attach small testcase? :)

With best regards,

Alexander Kalinkin

gitsnedbutzi

Beginner

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

06-13-2010
11:31 PM

60 Views

best regards

Alexander_K_Intel2

Employee

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

06-14-2010
01:15 AM

60 Views

With best regards,

Alexander Kalinkin

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