Support external scalar-valued functions

[Andreas Dutzler]

along with their gradient and hessian.

Note: The math methods inside the function/gradient/hessian must support arbitrary numbers of elementwise-operating trailing axes.