Autodiff support for higher dimensional Jacobians

[Ebanflo42]

As it is currently written, autodiff expects a scalar output node and the tests which present a vector output fail. We should extend the logic to be able to cope with the case of higher-dimensional output.