我已经在 Python 中实现了 Remez 算法，其中所有计算都是使用 Python mpmath 库完成的。我注意到有时和不再单调增加，而是反复迭代相同的值。主要是我的停止条件其中是工作精度的舍入误差，永远不会发生。提高 mpmath 的工作精度并不能解决此问题。$|E_{max}|$$|E_{min}|$$|E_{max}| - |E_{min}| < u$$u$

对于 Remez 算法的所有步骤，我使用默认方法来获得每个步骤的解：用于求解线性方程组的 mpmath 的 LU 分解、用于求根的割线方法和用于极值的黄金分割搜索。对不同的步骤使用不同的方法，可能更精确，会改善这种行为，还是我想念其他东西？我的停止标准是实用的还是我必须放宽它以供一般使用？