d********o 发帖数: 59 | 1 A 是sparse, discretized from a PDE. matlab 里的 "\" 自动并行对内存要求比较大
,3D的问题解起来就不快了,bicgstab 没有preconditioner的话,也是很满,gmres好
像就根本不起作用。大家一般都用什么方法呢?有没有好的review的文章呢? |
c*******h 发帖数: 1096 | 2 很多学校、实验室、公司都是自己写code的
在bicgstab, gmres, shur complement, AMG, 以及各种reordering和preconditioner
之间根据自己的喜好选择……
【在 d********o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : A 是sparse, discretized from a PDE. matlab 里的 "\" 自动并行对内存要求比较大 : ,3D的问题解起来就不快了,bicgstab 没有preconditioner的话,也是很满,gmres好 : 像就根本不起作用。大家一般都用什么方法呢?有没有好的review的文章呢?
|
l*****i 发帖数: 3929 | 3 ma27
【在 d********o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : A 是sparse, discretized from a PDE. matlab 里的 "\" 自动并行对内存要求比较大 : ,3D的问题解起来就不快了,bicgstab 没有preconditioner的话,也是很满,gmres好 : 像就根本不起作用。大家一般都用什么方法呢?有没有好的review的文章呢?
|
x*****u 发帖数: 3419 | 4 try PETSc
【在 d********o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : A 是sparse, discretized from a PDE. matlab 里的 "\" 自动并行对内存要求比较大 : ,3D的问题解起来就不快了,bicgstab 没有preconditioner的话,也是很满,gmres好 : 像就根本不起作用。大家一般都用什么方法呢?有没有好的review的文章呢?
|
x*****u 发帖数: 3419 | 5 http://www.greatlakesconsortium.org/events/scaling/files/Drummond.pdf
【在 d********o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : A 是sparse, discretized from a PDE. matlab 里的 "\" 自动并行对内存要求比较大 : ,3D的问题解起来就不快了,bicgstab 没有preconditioner的话,也是很满,gmres好 : 像就根本不起作用。大家一般都用什么方法呢?有没有好的review的文章呢?
|
i**g 发帖数: 134 | |
k*********g 发帖数: 791 | 7 gauss seidel,是最简单、最牢靠的方法;;; |