24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1266  |  回复: 18
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 coldwind042 的 19 个金币 ,回帖就立即获得 1 个金币,每人有 1 次机会
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

coldwind042

金虫 (正式写手)

[交流] 【讨论】如何求解大型系数矩阵(10万以上自由度)参与有奖

对于大型的工程问题,(譬如静力学问题),如果计算量过大(如自由度超过10万,甚至达到100万),有什么好的方法来求解此类线性方程组(系数矩阵为带状稀疏对称矩阵)?
前提:我是采用自己编的fortran程序或c程序。
希望大家帮我提点建议和思路,把自己的心得和体会,或者比较好的求解子程序推荐下!感激不尽!
群里面应该有很多航天或者化学或生物模拟方面的高人,应该有这方面类似的问题吧?请不吝赐教!

我目前已经考虑过的技术有:
   1,节点编号优化
   2,采用一维变带宽存储
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redcrag

木虫 (著名写手)

★ ★ ★
csfn(金币+1,VIP+0):thx very much 欢迎常来交流 8-26 10:38
coldwind042(金币+2,VIP+0):请问你的平面波程序是fortran的吗?你求解的是什么性质的稀疏矩阵?是否也是非带状的? 8-27 18:15
只是有收集书的爱好,自己没太大的研究。自己的稀疏矩阵是参照别人的程序用Matlab处理的(有限差分法),也用别人现成的迭代程序(平面波方法)。
引用回帖:
Originally posted by coldwind042 at 2009-8-25 16:19:
追加了金币,希望各位继续参与讨论!

http://redcrag.spaces.live.com/
12楼2009-08-26 08:34:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

xuqingyuan1972

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wangen994(金币+1,VIP+0):欢迎新虫,常来程序软件版 8-22 21:16
采用迭代算法,如pcg等。
2楼2009-07-18 09:31:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wade_thunder

新虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
我也遇到了同样的问题,希望可以共同讨论
3楼2009-07-18 17:31:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wade_thunder

新虫 (小有名气)

关键不是存储量大,而是方程求解太满.谁有好的迭代方法,请赐教.可以考虑收购.
4楼2009-07-18 17:33:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复(可上传附件)
信息提示
请填处理意见