Package to compute Hessian matrix numerically

Given numerical zᵢⱼ values defined over a regular 2D grid (xᵢⱼ, yᵢⱼ), is there a package to compute the Hessian matrix at each grid point?

ImageFiltering.jl · ImageFiltering perhaps? Or GitHub - omlins/ParallelStencil.jl: Package for writing high-level code for parallel high-performance stencil computations that can be deployed on both GPUs and CPUs?

1 Like

Thanks Fredrik, I will need to figure out how to use those packages for this task.