Hi!
It is impossible to get LU factorization of the input matrix. It is only possible to work with it through the handle.
Your code at first sight looks fine: each matrix is processed with its own handle, so LU factors are stored separately. The possible reason of the crash could be in allocation of input matrices, in DSS interface or PARDISO. Could you provide us with the test-case including linking line, input matrix generator, number of threads and so on? We will profile it in order to find the root of the problem.
Best regards,
Sergey