|
|
[求助]
问一个解线性方程的问题
现在我遇到一个问题,有两个稀疏的大型实对称矩阵A 和 B,其中B是正定的, 对一系列实数x, 只想得到(x-A+i*B)^(-1)的虚部。
我现在用的是最低级最简单的办法,对每个x直接用lapack求逆,然后取虚部。可是效率实在是太低了。
我想应该可以把这个问题转化成求一个线性方程组的解
\right)\left(\begin{matrix}Y\\Z\end{matrix}\right)=\left(\begin{matrix}I\\0\end{matrix}\right))
而且只需要得到Z。
请问有没有什么比较高效而且精确的算法能达到这个目的呢?比如用Krylov subspace? |
» 猜你喜欢
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有4人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
|