Now Trilinos supports Direct solvers with `initialize` and `solve` which are separated. This should allow the use of parallel sparse direct solvers.