FATODE: A Library for Forward, Adjoint, and Tangent Linear Integration of ODEs
Algorithm 832: UMFPACK V4.3---an unsymmetric-pattern multifrontal method