PDLPは、厳密解MIP向きではない

 KKT_toleranceを1e-7にすると、IPMよりもかなり遅くなりますし、それでいて精度も未だ十分ではありません。1e-4にして、推定値を求めるにはよいのですが、PDLPで厳密解を目指して最後まで通すには、不十分であると結論しました。

また、WarmStartもやってみましたが、全く効果がありませんでした。

そこで、Rootで、PDLPで推定⇒IPMに切り替えというのが、最も経済的なアプローチということになると思います。　残念ながらHIGHS　IPMは、未だマルチスレッド化もされていないし、コルスキー分解のGPU化もなされていないので、相当遅いのを覚悟しなければなりません。