- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hello,
I'm using Pardiso to compute the Schur complement of a symmetric positive-semidefinite matrix with a kernel of dimension 6. I don't know if this is officially supported, but for similar rank-deficient matrices, I sometimes get the correct results. To see whether the Schur complement is right or not, I'm checking in MATLAB the eigenvalues of the result given by Pardiso after the numerical factorization. The first 6 eigenvalues should be 0. With MUMPS, I get the following (correct) output:
>> eig(triu(MUMPS) + tril(MUMPS') - diag(diag(MUMPS)))
ans =
1.0e+11 *
-0.0000
-0.0000
-0.0000
0.0000
0.0000
0.0000
0.0253
0.0412
0.0796
0.0999
0.1078
0.1258
With brute force MATLAB, I get the following (correct) output:
>> sort(real(eig(MATLAB(10:63,10:63) - MATLAB(10:63,1:9) * MATLAB(1:9,1:9)^-1 * MATLAB(1:9,10:63))))
ans =
1.0e+11 *
-0.0000
-0.0000
-0.0000
-0.0000
0.0000
0.0000
0.0253
0.0412
0.0796
0.0999
0.1078
With Pardiso, I get the following (wrong) output:
>> eig(triu(Pardiso) + tril(Pardiso') - diag(diag(Pardiso)))
ans =
1.0e+11 *
-2.5490
-2.1172
-1.5318
-0.8751
-0.4827
-0.3927
-0.3516
-0.2968
-0.1554
-0.1176
I'm attaching a reproducer that factorizes the matrix and output the Schur complement, as well as the original matrix stored in dense format to be used inside MATLAB. Can you reproduce this problem ? Do I need to adjust some parameters, or is there a bug inside the solver ?
Thank you for your help.
Link Copied
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Any help please ?
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Your test returns -1, which means the input inconsistency. I enabled matrix checker and it reported more details, see below
*** Error in PARDISO ( sequence_ido,parameters) error_num= 10
*** Input check: iparam[5] -858993460 (out of bounds)
*** Input parameters: inconsistent error= 10 max_fac_store_in: 1
matrix_number_in: 1 matrix_type_in: -2
ido_in : 12 neqns_in : 63
ia[neqns_in]-1 : 707 nb_in : 1
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi,
I have a problem with Schur complement. In some cases when I calculate Schur complement I get wrong matrix. Values in matrix are correct but on the wrong place in the matrix. Error is appearing randomly.
Any ideas are welcome, thanks.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Nina, could you try the latest update 2 ( MKL 2017)? We have fixed the issue when MKL Pardiso Schur complement returns incorrect results. see MKL v.2017 Bug Fix list - https://software.intel.com/en-us/articles/intel-math-kernel-library-intel-mkl-2017-bug-fixes-list
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
With latest update it works now for symmetric matrix, but for unsymmetrical there is the same bug as before for symmetric.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Nina Z. wrote:
With latest update it works now for symmetric matrix, but for unsymmetrical there is the same bug as before for symmetric.
Hi Nina,
Can you please provide us with a reproducer for nonsymmetric case?
Best regards,
Maria
![](/skins/images/7FC17B7B85029576C25F1E43CE255B51/responsive_peak/images/icon_anonymous_message.png)
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page