  We apply a flexible numerical integrator to the simulation of adiabatic
quantum computation with nonlinear paths. We find that a nonlinear path may
significantly improve the performance of adiabatic algorithms versus the
conventional straight-line interpolations. The employed integrator is suitable
for solving the time-dependent Schr\"odinger equation for any qubit
Hamiltonian. Its flexible storage format significantly reduces cost for storage
and matrix-vector multiplication in comparison to common sparse matrix schemes.
