问一下2d/3d的fast sovler - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Computation - 科学计算

问一下2d/3d的fast sovler

问一下2d/3d的fast sovler# Computation - 科学计算

a*s2003-05-27 07:05

1 楼

求解2d/3d的poisson equation,大家有什么快速的solver?
要求系数不一定是symmetric的,boundary condition可以是
Neumann and Dirichlet boundary condition.
我知道的有CG(conjugate gradient)和multigrid,还有什么
更快,而且3d容易实现的么?

a*s2003-05-27 07:05

2 楼

preconditioned CG

【在 a******s 的大作中提到】

: 求解2d/3d的poisson equation,大家有什么快速的solver?
: 要求系数不一定是symmetric的,boundary condition可以是
: Neumann and Dirichlet boundary condition.
: 我知道的有CG(conjugate gradient)和multigrid,还有什么
: 更快,而且3d容易实现的么?

a*s2003-05-27 07:05

3 楼

我对这个算法很感兴趣,你用这个么?看不少paper都用了这个.
你用这个么? 算3d的case和GS比,要快多少倍?给个大概的概念就行了.
我好看看要不要去试.
thx

【在 a****s 的大作中提到】

: preconditioned CG

a*s2003-05-27 07:05

4 楼

我只用过PCG作练习性的计算, 100*100的二维矩阵.
其速度取决于matrix的condition number k=Lamda(max)/lamda(min), 如果k=1, 只需要
一次iteration就可以了. 你可以想象一下吧. 具体的我也不清楚.
还可以加上residual correction, 只要做一次就可以足够提高速度了.

【在 a******s 的大作中提到】

: 我对这个算法很感兴趣,你用这个么?看不少paper都用了这个.
: 你用这个么? 算3d的case和GS比,要快多少倍?给个大概的概念就行了.
: 我好看看要不要去试.
: thx