Scope for optimizing runtime speed of Jacobian and Hessian generating functions

[escheffel]

The functions which compute the Jacobian and Hessian can surely be sped up a little bit even when more than one core is used. Need to investigate this as it is an important aspect for power users.