I adopted PARDISO to a multiple matrix system, there are hundreds to thousands of matrices to solve. In this situation, each matrix has its own sparsity structure. Thus, each matrix has it own pt values. I tried some benchmarks and got some strange result. In one case, there are thousands of matrix structures. I set the msglvl to 1 and dumped the pt values to see what's going wrong. It turns out the first 14 matrice perform well and their pt values are well-set and different from each other. However, the 15'th matrix (the matrix size is 1x1) has it own pt[i] to be 0. And after the 15'th matrix, all the matrices have pt[i] set to 0 and the PARDISO does not show any messages (even msglvl is 1). I've tried to remove the 1x1 matrix and it turns out to be fine for other matrices until next 1x1 matrix appears. Does 1x1 matrix not set the pt value and affect all other succend matrices? How can I resolve or bypass this issue if I still want the 1x1 matrix to use PARDISO solver (even it's meaningless)? Also, does PARDISO restrict to the number of pt values? I mean does PARDISO has the constraint of the number of different pts. Thanks
ps : In the benchmarking of PARDISO, we force all matrices to use PARDISO solver without regarding the matrix size.