# Package to compute Hessian matrix numerically
**URL:** https://discourse.julialang.org/t/package-to-compute-hessian-matrix-numerically/123151
**Category:** General Usage
**Tags:** finitediff
**Created:** [November 27, 2024, 9:23am UTC](https://discourse.julialang.org/t/package-to-compute-hessian-matrix-numerically/123151 "2024-11-27T09:23:52Z")
**Posts on this page:** 3
**Page:** 1
### Author:  [@rafael.guerra](https://discourse.julialang.org/u/rafael.guerra)
#### Post date: [November 27, 2024, 9:23am UTC](https://discourse.julialang.org/t/package-to-compute-hessian-matrix-numerically/123151/1 "2024-11-27T09:23:52Z")
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?
---
### Author:  [@baggepinnen](https://discourse.julialang.org/u/baggepinnen)
#### Post date: [November 27, 2024, 9:51am UTC](https://discourse.julialang.org/t/package-to-compute-hessian-matrix-numerically/123151/2 "2024-11-27T09:51:46Z")
[ImageFiltering.jl · ImageFiltering](https://juliaimages.org/ImageFiltering.jl/stable/#Linear-filtering:-noteworthy-features) 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](https://github.com/omlins/ParallelStencil.jl)?
---
### Author:  [@rafael.guerra](https://discourse.julialang.org/u/rafael.guerra)
#### Post date: [November 27, 2024, 1:23pm UTC](https://discourse.julialang.org/t/package-to-compute-hessian-matrix-numerically/123151/3 "2024-11-27T13:23:55Z")
Thanks Fredrik, I will need to figure out how to use those packages for this task.