how to use pardiso to solve a singular matrix problem?

[smth]

Hi, I want to use pardiso to solve a poisson equation with neumann boundary condition.
for instance, Uxx=0 0<1, where partial(U)/parital(x)|x=0,1 equal to zero.
the matrix is just like this :

|1 -1 0 0 |
|-1 2 -1 0 |
|0 -1 2 -1 |
|0 0 -1 1 |

I know in MKL some possion pack can solve this problem. But I want to use paradiso because I intend to use it solve more complex problems.

Does anyone to point me out how to use paradiso to solve this singular problem?

[Sergey_K_Intel1]

PARDISO can only solve non-singular systems of linear equations. So PARDISO can not directly be used for solving rank-deficient problems.

All the best

Sergey