trlib公司

trlib:GLTR方法的无向量实现,用于迭代求解信赖域问题。我们描述了trlib,一个实现了Gould的广义Lanczos方法的变体的库(Gould et al.in SIAM J.Opt。9(2),504-5251999),用于解决信赖域问题。我们的实现有几个不同的特性,这些特性将它与先前存在的特性区分开来。对于Krylov子空间的装配,我们实现了共轭梯度(CG)和Lanczos迭代。无矢量和矩阵的反向通信接口允许使用大多数通用的数据结构,例如函数空间问题离散化后产生的结构。在求解非线性优化问题的序贯方法中,信赖域问题的硬情况经常出现。在这个实现中,我们努力通过考虑所有不变的Krylov子空间来精确地解决硬情况。研究了trlib在最割基准集的无约束问题全子集上的数值性能。除此之外,使用无矢量反向通信接口将PDE离散化工具箱FEniCS与trlib相连接,以解决一系列PDE约束控制信赖域问题。