[論文レビュー] JuMP 1.0: Recent improvements to a modeling language for mathematical optimization

Miles Lubin, Oscar Dowson+4 more|arXiv (Cornell University)|Jun 1, 2022

Advanced Control Systems Optimization被引用数 8

ひとこと要約

JuMP 1.0 は、Julia における数学的最適化のモデリング言語に大幅な改善をもたらし、パフォーマンス、メモリ効率、拡張性を向上させるために MathOptInterface 抽象レイヤーを活用している。主な進歩には、低レベルのソルバへの直接アクセスを可能にする直接モード、プリコンパイルとキャッシュによるコンパイル遅延の低減、インジケータ制約、補完性制約、一般化不等式といった新しい制約タイプのサポート、および拡張されたソルバ互換性が含まれ、最適化ワークロードにおける使いやすさとスケーラビリティが顕著に向上している。

ABSTRACT

JuMP is an algebraic modeling language embedded in the Julia programming language. JuMP allows users to model optimization problems of a variety of kinds, including linear programming, integer programming, conic optimization, semidefinite programming, and nonlinear programming, and handles the low-level details of communicating with solvers. After nearly 10 years in development, JuMP 1.0 was released in March, 2022. In this short communication, we highlight the improvements to JuMP from recent releases up to and including 1.0.

研究の動機と目的

• JuMP の以前のバージョンにおけるモデル構築およびソルバとの通信におけるパフォーマンスボトルネックを解消すること。
• Julia におけるコンパイル遅延を低減すること。これは、JuMP のインタラクティブおよびコマンドライン利用を妨げる要因である。
• direct_model の導入により、状態なしで低コストのソルバとのやり取りを可能にすることで、メモリ効率を向上させること。
• バイレベル最適化、無限次元最適化、確率的最適化を含む新しい問題タイプへの拡張性とサポートを強化すること。
• MathOptInterface 抽象レイヤーを通じて、ステータス報告とソルバ相互運用性を統一・改善すること。

提案手法

• モデリングとソルバ固有の実装詳細を分離するために、MathOptInterface をコアの抽象レイヤーとして採用すること。
• direct_model を導入し、JuMP が中間のキャッシュを経由せずに問題データをソルバに直接渡す状態なしのラッパーとして機能できるようにすること。
• Julia の JIT コンパイル遅延を緩和するため、プリコンパイルとキャッシュ戦略を実装すること。
• 属性システムと制約 API を拡張し、インジケータ、補完性、一般化不等式制約などの新しいタイプをサポートすること。
• MathOptInterface からの標準化された列挙型を用いて、ソルバに依存しないモデル照会とステータス報告を可能にすること。
• Julia のパッケージマネージャーと統合し、オープンソースのソルバのためのバイナリをプラットフォーム間で自動的にビルドおよび配布すること。

実験結果

リサーチクエスチョン

• RQ1 JuMP は、モデル構築中のメモリオーバーヘッドを低減し、パフォーマンスを向上させるためにどのような手段をとれるか？
• RQ2 科学計算ワークロードにおける最適化モデリングを伴う Julia のコンパイル遅延を緩和する戦略は何か？
• RQ3 JuMP は、ソルバに依存しない状態を保ちながら、より広範な制約タイプと問題クラスをサポートするにはどうすればよいか？
• RQ4 外部ソルバおよびドメイン固有の拡張機能との拡張性と相互運用性を向上させるために、どのような改善が必要か？
• RQ5 さまざまな最適化ソルバ間で、ステータス報告とソルバフィードバックを一貫性があり、情報豊かにすることはどのように実現できるか？

主な発見

• 直接モードにより、メモリ使用量が削減され、Highs_scaleCol などのソルバ C API への直接アクセスが可能になり、最小限の抽象化オーバーヘッドで実現できる。
• プリコンパイルとキャッシュの導入により、単純な線形計画問題においてコンパイル遅延が約 16 秒から 1 秒未満にまで短縮された。
• サポートされるソルバ数は 13 個から 41 個に増加し、Julia のパッケージマネージャーを介した自動バイナリ配布が可能になった。
• JuMP 1.0 は、インジケータ、補完性、一般化不等式制約といった新しい制約タイプを導入し、ネイティブ構文とソルバ統合を実現した。
• MathOptInterface への移行により、プライマル、デュアル、終了ステータスの標準化された列挙型を用いた統一的で洗練されたステータス報告システムが実現された。
• ドキュメントはチュートリアル、ハウツー、リファレンスマニュアルを大幅に強化し、新規ユーザーのアクセス性が向上した。