哪些 ODE 求解器针对使用稀疏雅可比行列式求解刚性系统进行了优化？例如，当使用典型的有限差分或有限元方法在空间中离散抛物线 PDE 时，就会出现此类系统。