void hessianAt(double const * p, AddressableMatrix<double> & result) const { checkHessianResultDims(result); result.makeZero(); }
void hessianAt(float const * p, AddressableMatrix<float> & result) const { checkHessianResultDims(result); result.makeZero(); }