# Package For Constructing Finite Difference Operators on Non-Uniform Grid?

**URL:** <https://discourse.julialang.org/t/package-for-constructing-finite-difference-operators-on-non-uniform-grid/138333>\
**Category:** Numerics\
**Tags:** question, package\
**Created:** [July 17, 2026, 8:49pm UTC](https://discourse.julialang.org/t/package-for-constructing-finite-difference-operators-on-non-uniform-grid/138333 "2026-07-17T20:49:23Z")\
**Posts on this page:** 1\
**Showing post:** 9

<div class="post-metadata">

**Author:** ![ducksoverip](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/ducksoverip/32/31967_2.png) [@ducksoverip](https://discourse.julialang.org/u/ducksoverip)\
**Post date:** [July 20, 2026, 4:01pm UTC](https://discourse.julialang.org/t/package-for-constructing-finite-difference-operators-on-non-uniform-grid/138333/9 "2026-07-20T16:01:07Z")

</div>

Right, but I still need to construct the matrix first, even if I wrap it as an operator later. This is because the full matrix has some funky structure due to the underlying equation (see the second bullet point in the [linked post](https://discourse.julialang.org/t/solving-the-helium-eigenvalue-problem-where-to-start/137783/7)). That summed V\_{F0} term, combined with the fact that the operator spans all instances of the index i, means I need to explicitly construct the matrix operator for the PDE. I can do all the steps of the construction as long as I can get the 1D second-derivative operator on one axis of grid points, but doing that with FD methods for a non-uniform grid is where I’m getting hung up.

---

_[View the full topic](https://discourse.julialang.org/t/package-for-constructing-finite-difference-operators-on-non-uniform-grid/138333)._