24小时热门版块排行榜    

查看: 1108  |  回复: 3

sasori

木虫 (小有名气)

[求助] Intel MKL的 BLAS 的SpMV 2求解稀疏矩阵线性方程结果不对 求原因

使用Intel MKL的 BLAS 的SpMV 2中的算子 解线性方程组A*x=y , A为稀疏矩阵 ,y为已知向量,解到的x向量结果不对,   不知道是什么原因。。
A=【1 2 * *
         2 4 * *
         * * 9 12
         * * 12 21】
y=【3 6 21 33】'
x 应该为【1 1 1 1】'才对  但计算结果并非如此

另外用算子计算乘法A*x=y,  A为稀疏矩阵 ,x为已知向量,得到的y是对的,说明对算子含义描述的理解应该没有错。。

求大神指点~~~
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sspa9999

至尊木虫 (著名写手)

2楼2016-10-13 21:17:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sspa9999

至尊木虫 (著名写手)

a不是满秩的,第2行是第一行的2倍,当然结果会有问题了

发自小木虫Android客户端
3楼2016-10-13 21:18:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sasori

木虫 (小有名气)

引用回帖:
3楼: Originally posted by sspa9999 at 2016-10-13 21:18:54
a不是满秩的,第2行是第一行的2倍,当然结果会有问题了

是的 结果发现是用的算子只针对三角矩阵求解的。。
一般稀疏矩阵的线性求解 看网上是说要用pardiso 帅哥会不会用啊?

发自小木虫Android客户端
4楼2016-10-14 19:54:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sasori 的主题更新
信息提示
请填处理意见