c*m
2 楼
matlab就可以
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
l*i
3 楼
10k variables用Petsc是不是有点奢侈，呵呵。
x*u
4 楼
看看Petsc。
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
l*i
5 楼
ma27 if the matrix is symmetric. look for it in hsl.
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
【在 r*****d 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
m*e
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 的大作中提到】
: 是稀疏的,请教有什么比较稳定的软件可以干这个?
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
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是不错，不过功能太强大了，说明书看着头就昏，而且不太好安装。
【在 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
8 楼
我觉得如果对于sparse matrix 不是很熟悉的话, 有个 叫作sparse 的C程序库很好用,
它也提供了Fortran 的界面. 如果想在自己的程序中用的话, 这个库还是比较容易
上手的. 比较容易找到的是1.3 (netlib上就是这个版本), 大概因为是比较老的C程
, 所以编译常常会有些问题. 幸运的是, 最近作者推出了新版本, 所以编译也不是问
题了. 如果感兴的话, google一下就可以找到了.
它也提供了Fortran 的界面. 如果想在自己的程序中用的话, 这个库还是比较容易
上手的. 比较容易找到的是1.3 (netlib上就是这个版本), 大概因为是比较老的C程
, 所以编译常常会有些问题. 幸运的是, 最近作者推出了新版本, 所以编译也不是问
题了. 如果感兴的话, google一下就可以找到了.
b*e
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.
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
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.
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
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
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?
【在 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
13 楼
Preconditioned Bi-CGSTAB
Really good.
【在 l*****i 的大作中提到】
: There are other Krylov space methods for indefinite systems than CG...
Really good.
【在 l*****i 的大作中提到】
: There are other Krylov space methods for indefinite systems than CG...
相关阅读
Anybody know Lehigh U?CG versus BFGS for iterative linear solversTwo PhD positions are available immediately – computation and simulation of fuel cells求帮忙下载一篇paper 有没有现成的p-median sovler请问怎么增大matlab的可用内存?求助，关于monte carlo计算 (转载)Can we get the data from a figure in the paper?如何在双轴的Maltab图里添加曲线啊。lapack安装问题有人做CFD吗？这是什么算法?哪位高手 熟 通过Matlab 从 MySQL获取数据怎么拟合这样的曲线请问怎么用matlab对N个离散概率密度函数卷积如何实现N层循环结构derivative of the cubic spline interpolation results我来问个问题请教一个网络的问题请教matlab安装