Switch to PETSc linear solver library. This offers more flexibility in
the solution method, i.e., users can now select between direct and
iterative solvers, and also accommodates MPI for parallel computing.