This branch interfaces MBDyn with Sandia's Trilinos library , including the NOX nonlinear solver with features like limited memory Broyden updates and a matrix free Newton Krylov solver which uses automatic differentiation to evaluate Jacobian vector products.