首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >是否有计算Cholesky因子更新的BLAS/LAPACK函数?

是否有计算Cholesky因子更新的BLAS/LAPACK函数?
EN

Stack Overflow用户
提问于 2018-05-21 14:07:12
回答 2查看 1.7K关注 0票数 2

A是正定矩阵,A=L*L'是它的cholesky分解,其中L是下三角。设A2 = A + alpha*x*x'是矩阵A的秩-1更新,其中x是适当维数的向量,alpha是标量。

Cholesky因子更新是在不首先计算A2的情况下获得因式分解A2=L2*L2'的过程,在这种低秩矩阵更新的情况下,这对于加快计算非常有用。

我正在使用BLAS/LAPACK库进行初等代数操作。我可以用常规的spptrf计算正定矩阵的Cholesky分解。然而,我一直在寻找,我没有找到一个BLAS/LAPACK功能,执行Cholesky因素更新。是否没有这样做的功能呢?

Additionally:在此旧员额中,讨论了添加此类例程的问题。然而,这是一个非常旧的职位(2013年),我一直未能找到任何更近期的。

EN

回答 2

Stack Overflow用户

发布于 2018-05-27 19:12:22

没有这样的功能。您可以查看我们在SciPy上进行的讨论。我已经编写了一个Python脚本,用相关的文件进行更新。你可以利用这些信息。

https://github.com/scipy/scipy/issues/8188

如果您感到有竞争力,并为此编写了Fortran代码,我将非常感谢您能够将其作为PR https://github.com/Reference-LAPACK/lapack提交给LAPACK回购公司。

票数 2
EN

Stack Overflow用户

发布于 2018-05-24 21:16:07

正如您所指出的,Netlib上的BLAS库,但我怀疑它在网站上。如果你只是在寻找代码,那么就有代码这里。,如果你想要它的话,我会把它转化成朱莉娅。有一本书我从来没有看过,其中可能有这些在它。中也有,请注意,您引用了一篇论文,作者写的代码。你可以直接联系那篇论文的作者。他的网站似乎是这里。,但是这个链接有一个问题。

票数 1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/50450675

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档