v0.2.4

Efficiency improvements:

• add f_and_g and f_and_g_prod functions so that drift and diffusion can be computed together

Bug fixes:

• calling solver.integrate twice internally
• adaptive error estimation not recognizing singleton tensors