Community
cancel
Showing results for 
Search instead for 
Did you mean: 
crispybits
Beginner
31 Views

PARDISO produces all-zero result randomly

Hi. About once in every 1000 solves, I get all zeros as the solution. Any possible reason? It doesn't seem to depend on the input. If I simply run my program again it's fine. The error code is 0 when it happens.

I just noticed that it isn't touching the solution vector. If I put a value in it before calling Pardiso, when it fails that value is still there.
MKL 10.2.2.025
Real symmetric matrix
Compiled with VC++ to a dll
Called by VB.NET app
In core
Reordering type 2 (METIS)
1 RHS vector
2-cores
0 Kudos
1 Reply
crispybits
Beginner
31 Views

I found a workaround. After solving I check to see if the solution was modified, and if not, go and solve again.