标签: julia julia-jump
如何处理JuMP中的稀疏矩阵?
例如,假设我想强制形式约束:
A * x == 0
其中A是稀疏矩阵,x是变量向量。我假设可以利用A的稀疏性来加快优化。我怎样才能在JuMP中利用这一点?
A
x
答案 0 :(得分:3)
JuMP已经以不同的方式从稀疏矩阵中受益,我没有检查来源,而是从cited paper引用JuMP.jl:
JuMP
在LP的情况下,输入数据结构是向量c和b 和稀疏格式的矩阵A,以及生成这些的例程 数据结构称为矩阵生成器
需要注意的一点是,像JuMP这样的代数建模语言(AML)的主要任务是为求解器生成输入数据结构。像JuMP这样的AML本身并不能解决生成的问题,但是它们会调用标准的适当求解器来完成任务。