- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi,
For a project, we woule like to solve multiple independent linear system concurrently. The structure of the code is like this:
omp_set_max_active_levels( 2 );
#pragma omp parallel num_threads( 2 )
{
mkl_set_num_threads_local( 10 );
op A = assemble();
auto pardiso_solver = createSolver(A);
pardiso_solver->solve(rhs, x);
}
When the outer num_threads is set to 1, everything works ok. However, when I increase it to 2. I got error code -2 in solve stage. In my test, the operator A is nothing but a 3x3 matrix:
/*
A = 1 2 3
0 4 5
0 0 6
*/
my solver parameters are set as:
_iparm[1] = 3;
_iparm[7] = 1;
_iparm[34] = 0;
Could I know what could be the cause of the problem?
Link Copied
1 Reply
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hello miaodi,
Thank you for posting here.
Could you tell us your oneMKL version, OS version and hardware platform? As well as providing us a simple reproducer for the issue.
Regards,
Ruqiu
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