Dear MKL developers,
I am trying to use FEAST eigenvalue solver with CSR formatted matrix, but encountered unexpected error (info = -2).
Both the Table 10 in the FEAST manual (http://arxiv.org/pdf/1203.4031v3.pdf) and MKL manual (https://software.intel.com/en-us/node/521739#GUID-E1DB444D-B362-4DBF-A1D...) said that this happend because "Internal error of the inner system solver. Possible reasons: not enough memory for inner linear system solver or inconsistent input.".
Could you please provide me comments or advices to fix this problem? The icc version is 15.0.3 and MKL version is 11.2, 64-bit system. I attached the sample code with three different CSR matrices. You can choose the matrix (small, medium, large) by commenting in/out the parts at lines 26-46.
Each CSR file contains Hermitian matrix,
1. 'n': dimension of the matrix
2. 'ia': rows, n+1 elements
3. 'n_val': # of non-zero elements
4. 'ja': columns, n_val elements
5. 'val': values, n_val elements
When I choose 'Small' matrix, I could get proper eigenvalues. But, when I choose 'Medium' matrix, it gives info = -2 and unrealistic 'm' (# of eigenvalues obtained). And for the 'Large' matrix, I got a segmentation fault.
When I ran the code with much simpler sample matrices, it worked fine even with very large matrices.
Thank you in advance,
Pilkyung
