我正在寻找一种对 PDE 进行局部敏感性分析的方法，最好是在 Python 中。

我的印象是，将方程离散化然后将其视为 ODE 是可行的；但是，这在计算上是否易于处理？

保持 PDE 结构而不是恢复到（硬）ODE 不会提高效率吗？特别是，允许​​自适应空间网格（如本文中的自适应网格伴随敏感性）似乎是一个巨大的提升。这种直觉正确吗？如果是，是否有任何现有的 PDE 特定敏感性分析实现，或者我应该调整一个 ODE 还是从头开始编写？