add a missing axpyImpl method for hessians.

Code generation is still missing for this function so only the slower fallback implementation is used

Fixes #164 (closed)