24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lixy1217

木虫 (著名写手)

应助: 4 (幼儿园)
金币: 4861.1
帖子: 1379
在线: 233.7小时
虫号: 1125891

[交流] 稀疏矩阵的处理

现在感觉数值解PDEs难免要跟稀疏矩阵打交道。所以想跟大家讨论一下是如何在计算中对稀疏矩阵进行操作

偶一直用的是C++，之前曾经尝试用过十字链表，虽然这样可以方便地调用每行每列的数据，但是感觉十字链表的计算效率太低（不知是不是编程技巧的问题）。

由于稀疏矩阵的运算主要是用于矩阵的加减法或矩阵乘向量，而矩阵之间的乘法很少，因此我个人觉得在操作中只需能够方便调用每行就行了，这样便可以用单链表的方式来储存矩阵，那计算效率会高了很多。

不知各位对此的看法如何？

[ 来自科研家族皇家数理科学协会 ]

回复此楼

» 猜你喜欢

北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有4人回复
自荐读博已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有5人回复
论文终于录用啦！满足毕业条件了已经有22人回复
不自信的我已经有5人回复
磺酰氟产物，毕不了业了！已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab稀疏矩阵问题求教已经有3人回复
求大神帮忙解病态线性方程组已经有5人回复
fortran到底怎么弄稀疏矩阵啊？已经有8人回复
matlab中建立稀疏矩阵组并操作已经有4人回复
用fortran存储稀疏矩阵已经有4人回复
【分享】Direct MEthods for Sparse Linear Systems （稀疏矩阵，基本算法）已经有76人回复
图像恢复或去噪中构造方程组的问题（matlab）已经有3人回复
大型非对称稀疏矩阵求解已经有9人回复
matlab 如何处理大规模矩阵？？？大家帮帮忙吧。急已经有5人回复
【求助】请问如何解决insufficient virtual memory 已经有10人回复
【求助】解大型稀疏矩阵的问题已经有6人回复
【求助】急问,如何将一个稀疏的9000阶矩阵,导入 mathlab 进行数值计算?万谢! 已经有3人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2013-04-04 09:40:55

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

hfj1988

新虫 (小有名气)

应助: 4 (幼儿园)
金币: 1057.2
帖子: 185
在线: 115.2小时
虫号: 2138982

★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
lixy1217: 金币+1 2013-04-06 11:03:38

如果要求调用矩阵中（i,j）元素的下面（i+1,j）元素怎么办？
用稀疏存储就是用时间换取空间，没有完全节省的，必然会以另外的东西作为代价。

赞一下

回复此楼

3楼2013-04-05 16:33:52

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

hfj1988

新虫 (小有名气)

应助: 4 (幼儿园)
金币: 1057.2
帖子: 185
在线: 115.2小时
虫号: 2138982

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

4楼: Originally posted by lixy1217 at 2013-04-06 11:03:19
可能还是要具体问题具体方法吧。一般用差分或者有限元解PDE时，每一行的非零元素都很少，从 i+1 行开始寻找到（i+1,j）需要的时间可忽略，相比之下，要在十字链表中把所有元素串起来确实很费时。

不过你说的很有 ...

差分或者有限元解PDE，都是稀疏矩阵，而且非零元素都集中在对角线的一定的带宽里面（如果编码得当的话），可用带宽存储（有变带宽和定带宽），这样存储量少，而且计算速度我认为比链表存储快。
如果所求的稀疏矩阵没有这样的性质（非零元素集中在对角线附近），那么为了节省存储量就不得已用链表了。
我是这样认为的，不知道是否正确。

赞一下

回复此楼

5楼2013-04-07 14:08:35

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