Hello, I would like to use SuperLU as a direct solver for a large size problem running on multiple processors. I have write the analytical Jacobian via FormJacobianLocal(), however, I do not know how to set the direct solver. Any thoughts? Thanks a lot! Rebecca