我有一个由一组不等式和等式描述的多面体，这些不等式和等式以numpy数组的形式给出。这些数组通常会很大，因为它们描述了MIPLIB问题的LP松弛。我想计算这个多面体的解析中心。如果可能的话，我想直接在Python代码中完成。但是，我也对这样的解决方案持开放态度，即（例如）将多面体描述写入MPS文件，并使用另一种工具来计算分析中心。
到目前为止，我已经尝试实现了原始-对偶牛顿算法，如“内点算法：Yinyu Ye的《理论与分析》中的方法。但是这种方法太慢了。我尝试过的另一种方法是使用目标为零的路径跟踪内点方法。为此，我使用了scipy的linprog。但是，对简单示例的评估表明，method="interior-point"给出了一个与分析中心不同的内点。即使我设置了options={"resolve": False, "rr": False, "autoscale": False}来防止修改多面体描述。使用method="highs-ipm"将给予顶点解。
我很乐意听听你的想法。