有人可以指点我从二维矢量场中提取螺线管（又名无散度）场的软件/代码吗？

有大量的论文和文档描述了如何做到这一点，但由于某种原因，没有一个作者（或其他任何人）在网上发布了一段简单的源代码来实现该功能。我发现的只是相当大的包或非常神秘的代码。

更新。更具体地说，对于我的情况：我有两个大小为 64x64 的二维数组，带有双精度浮点数。这两个阵列代表流体流动的 x 和 y 速度。使用周期性边界。我想将速度的螺线管部分提取到两个新的二维阵列中。