r*****d
发帖数: 54 | |
c******m
发帖数: 599 | 2 matlab就可以
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
|
l*****i
发帖数: 3929 | 3 10k variables用Petsc是不是有点奢侈，呵呵。 |
x*****u
发帖数: 3419 | 4 看看Petsc。
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
|
l*****i
发帖数: 3929 | 5 ma27 if the matrix is symmetric. look for it in hsl.
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
|
m********e
发帖数: 5088 | 6 如果是稀疏的：
direct methods:
UMFpack(seriel), Pardiso(OpenMP), SuperLU(MPI), Mumps(MPI)
优点是用起来比较简单，缺点是Factorization耗内存太多，不过才1万阶，用个单线程
我想就行了
iterative methods的软件很多，不过要提供preconditioner，case dependent。所以
我也不大用，我知道的一般的conjugate gradient的方法就很好：http://www.math.ruu.nl/people/sleijpen/
Pets是不错，不过功能太强大了，说明书看着头就昏，而且不太好安装。
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
|
c*******e
发帖数: 8624 | 7 cg记得对spd矩阵才行,还是gmres吧
【在 m********e 的大作中提到】
: 如果是稀疏的：
: direct methods:
: UMFpack(seriel), Pardiso(OpenMP), SuperLU(MPI), Mumps(MPI)
: 优点是用起来比较简单，缺点是Factorization耗内存太多，不过才1万阶，用个单线程
: 我想就行了
: iterative methods的软件很多，不过要提供preconditioner，case dependent。所以
: 我也不大用，我知道的一般的conjugate gradient的方法就很好：http://www.math.ruu.nl/people/sleijpen/
: Pets是不错，不过功能太强大了，说明书看着头就昏，而且不太好安装。
|
w*******e
发帖数: 6 | 8 我觉得如果对于sparse matrix 不是很熟悉的话, 有个 叫作sparse 的C程序库很好用,
它也提供了Fortran 的界面. 如果想在自己的程序中用的话, 这个库还是比较容易
上手的. 比较容易找到的是1.3 (netlib上就是这个版本), 大概因为是比较老的C程
, 所以编译常常会有些问题. 幸运的是, 最近作者推出了新版本, 所以编译也不是问
题了. 如果感兴的话, google一下就可以找到了. |
b***e
发帖数: 38 | 9 Direct method: frontal.
Iterative method: CG for symmetric and GMRES for asymmetric.
If you want to practice, CG is only about 20 lines. GMRES is about 100-200
lines. Frontal method is about a thousand lines if performance is not
strictly required. Once you have Frontal works, pre-conditioner is a piece
of cake. |
l*****i
发帖数: 3929 | 10 Since you mentioned CG, how are you so sure that LZ's problem is positive de
finite?
【在 b***e 的大作中提到】
: Direct method: frontal.
: Iterative method: CG for symmetric and GMRES for asymmetric.
: If you want to practice, CG is only about 20 lines. GMRES is about 100-200
: lines. Frontal method is about a thousand lines if performance is not
: strictly required. Once you have Frontal works, pre-conditioner is a piece
: of cake.
|
b***e
发帖数: 38 | 11
de
I use pre-conditioned CG for near-singular problem. Can you tell me any
iterative methods that are proven more effective for non positive definite
or even singular problems?
【在 l*****i 的大作中提到】
: Since you mentioned CG, how are you so sure that LZ's problem is positive de
: finite?
|
l*****i
发帖数: 3929 | 12 There are other Krylov space methods for indefinite systems than CG...
【在 b***e 的大作中提到】
:
: de
: I use pre-conditioned CG for near-singular problem. Can you tell me any
: iterative methods that are proven more effective for non positive definite
: or even singular problems?
|
l*w
发帖数: 646 | 13 Preconditioned Bi-CGSTAB
Really good.
【在 l*****i 的大作中提到】
: There are other Krylov space methods for indefinite systems than CG...
|
