Sparse direct methods for model simulation