假设一个人在 3D 空间中定义了一个标量场，他想为其梯度找到 Morse-Smale 复合体，以便稍后在吸引子域上执行几个六维函数的积分。要集成的功能大多集中在吸引子周围，在边界处贡献相对较小。同样，边界可以有一个尖尖的形状。什么方法（或算法甚至软件）适合解决这个问题？

目前这个问题分为两个步骤：首先，通过在空间中的许多点测试每个吸引子的边界是否遵循梯度驱动搜索到其他吸引子。沿从吸引子开始的射线选择初始点。这些光线指向由列别捷夫网格在单位半径球体上生成的点的方向。所以后来积分在径向和角度变量中分离。对于径向正交中的每个点（不超过距离吸引子最远的边界点），对相应的球体进行积分，检查每个角点是在积分区域的内部还是外部。

我想知道是否有更有效的方法围绕它具有足够的分辨率来处理尖峰边界并保持角度和径向积分分离的好处。