PARDISO produces all-zero result randomly

[crispybits]

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

[crispybits]

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