24小时热门版块排行榜    

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

coldwind042

金虫 (正式写手)

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

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

我目前已经考虑过的技术有:
   1,节点编号优化
   2,采用一维变带宽存储
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 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分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复(可上传附件)
信息提示
请填处理意见