Evaluating Hessians while specifying coefficients in PDE Toolbox

I use the PDE toolbox to solve a partial differential equation of the following form:

-div(c.grad(u))=f,

with c, u and f as in the documentation.

I use the function specifyCoefficients to specify the coefficients c and f. When specifying coefficients, I have access to a structure 'state' with these fields: state.u, state.ux, state.uy. So, I have access to the value and the gradient of u at the given location.

Now the problem is the following: I also need access to the Hessian of u in order to construct the coefficients c and f. Would this somehow be possible?