24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1269  |  回复: 18
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 coldwind042 的 19 个金币 ,回帖就立即获得 1 个金币,每人有 1 次机会

coldwind042

金虫 (正式写手)

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

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

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

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的回帖

雪连天

新虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by wade_thunder at 2009-7-18 17:33:
关键不是存储量大,而是方程求解太满.谁有好的迭代方法,请赐教.可以考虑收购.

收购?你打算出多少银子?
5楼2009-08-02 20:07:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lerel12345

木虫 (职业作家)


小木虫(金币+0.5):给个红包,谢谢回帖交流
我也有同样的问题
月有阴晴圆缺 星星一直在眨眼
6楼2009-08-22 09:33:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redcrag

木虫 (著名写手)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gwdavid(金币+2,VIP+0):谢谢 8-24 15:43
coldwind042(金币+2,VIP+0):谢谢兄弟,看你这几本书就知道你有多牛气啦!我算的是稀疏不对称矩阵,能否更针对性地推荐其中哪本书哪一块啊?感激不尽! 8-25 16:19
找些相关书看看

Numerical Methods for Large Eigenvalue Problems (Manchester) Sadd
Computer Solution of Large Linear Systems (NH) Meurant
Direct Methods for Sparse Linear Systems (SIAM) Davis
Large Scale Eigenvalue Problems (NH) Cullum
Iterative Methods for Sparse Linear Systems (SIAM,2nd) Saad
http://redcrag.spaces.live.com/
7楼2009-08-24 12:01:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magic7004

金虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
coldwind042(金币+1,VIP+0):没什么价值,不过还是谢谢你的回帖! 8-25 16:19
1,写程序
2,找台超级计算机
3,挂机
流氓不可怕,可怕的是流氓有文化,有文化又BH的流氓无敌~~!
8楼2009-08-24 16:31:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coldwind042

金虫 (正式写手)

追加了金币,希望各位继续参与讨论!
9楼2009-08-25 16:19:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxshow

新虫 (初入文坛)

★ ★ ★ ★
csfn(金币+3,VIP+0):thx very much 欢迎常来交流 8-26 10:35
coldwind042(金币+1,VIP+0):谢谢参与讨论。matlab里面有很多好的东西。但直接应用于大问题是不可行的。可惜要将里面的好技术应用于fortran程序还很麻烦,不知道它的内核 8-27 18:13
我是做电磁场有限元分析的
也会遇到这些问题,我一直用的是matlab
可以借用他的存储方式存储
采用一维变带宽存储 是老方法了

解线性系数真的方法也有很多,matlab有很多相应函数
可以结合自己的优化程序求解 很方便

当然需要很大内存和cpu资源

在深一层的分块求解什么的我就不会了
学习知识 结识朋友 提高水平
10楼2009-08-25 19:40:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 coldwind042 的主题更新
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复(可上传附件)
信息提示
请填处理意见