由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Computation版 - lapack如何求解XA=B
相关主题
直接用NR的源程序需要注意什么吗?Gfortran library
求矩阵逆的算法C++里用Blas/Lapack的问题
Fast solver of linear equation system我们碰到的大麻烦——急寻计算机工作和博后机会 (转载)
有没有用Lapack的?paralllel eigensolver
mklman is by far the most grandious reference manual i've seen稀疏矩阵压缩
[合集] 这里有lapack用的比较熟悉的人吗?请问如何在C++中使用lapack?
do you use blas/lapack?problem about Fortran 77 in unix
请教一个计算速度的问题LAPACK++ Re: 求解PDE时,用到大型的矩阵,C++里怎么处理?
相关话题的讨论汇总
话题: lapack话题: 求解话题: dtrsm话题: lu话题: dgetrf
进入Computation版参与讨论
1 (共1页)
p*****e
发帖数: 310
1
我用lapack的dgetrf把A进行LU分解,也发现用blas的dtrsm可以进行右边求解,但是发
现dgetrf里的iPiv,即pivoting得信息没法在dtrsm里体现出来。而lapack中的dgetrs没
用右边求解的功能。
我现在只好先用A的LU分解,用dgetri得到逆然后B*A^(-1),这样效率太低了。
另外的方法是用dtrsm求解完以后在对解进行列变换,即X*(PLU)=B => X*P = Y, 可是
发现lapack只提供dlaswp来进行行变换。
另外我还想到因为X*P=Y => P'X'=Y', 这样可以用行变换得到X转置再转置回来,但是这
也太麻烦了,而且还得知道inv(P'),效率如何也不清楚。
不知道有没有好办法?
1 (共1页)
进入Computation版参与讨论
相关主题
LAPACK++ Re: 求解PDE时,用到大型的矩阵,C++里怎么处理?mklman is by far the most grandious reference manual i've seen
复数矩阵对角化有多快?[合集] 这里有lapack用的比较熟悉的人吗?
请教一个fortran 90编程的问题do you use blas/lapack?
lapack++ and blitz++请教一个计算速度的问题
直接用NR的源程序需要注意什么吗?Gfortran library
求矩阵逆的算法C++里用Blas/Lapack的问题
Fast solver of linear equation system我们碰到的大麻烦——急寻计算机工作和博后机会 (转载)
有没有用Lapack的?paralllel eigensolver
相关话题的讨论汇总
话题: lapack话题: 求解话题: dtrsm话题: lu话题: dgetrf