Optional forward mode for jac/jac-vec products

[mlubin]

@amontoison observed that we could use forward mode for the jacobian and jacobian-vector products when the number of variables is less than the number of rows to get better performance than the current approach that uses reverse mode always.

[odow]

This could be best implemented as a new AD backend. It shouldn't be too complicated.