julia - 带有稀疏矩阵的 JuMP？

Question

我如何处理 JuMP 中的稀疏矩阵?

例如，假设我想施加以下形式的约束:

A * x == 0

哪里 A是一个稀疏矩阵和 x变量向量。我假设 A 的稀疏性可以被利用来使优化更快。如何在 JuMP 中利用这一点？

Answer 1

JuMP已经以不同的方式从稀疏矩阵中受益，我没有检查源，但引用了 cited paper来自 JuMP.jl :

In the case of LP, the input data structures are the vectors c and b and the matrix A in sparse format, and the routines to generate these data structures are called matrix generators

需要注意的一点是，像 JuMP 这样的代数建模语言 (AML) 的主要任务是为求解器生成输入数据结构。像 JuMP 这样的 AML 本身并不解决生成的问题，但它们会调用标准的适当求解器来完成任务。