Parallel Direct Solution of Large Sparse Systems in Finite Element Computations