问题描述
我有一个带有对称块的通用矩阵,比如说,
M = [ x,v^T
u,A ]
其中 A 是对称块,而 u,v 是一些列向量。我想对 A 执行对称 rank-1 更新,并且我知道 dsyr。然而,似乎 dsyr 只提供更新下三角或上三角的选项,而不是两者。我知道我可以写一些循环来复制矩阵元素,或者使用 dger 代替。但我想知道,有没有更优雅的方法来做到这一点?
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)