H*******k
发帖数: 894 | 1 matlab 里面如果矩阵 Ax=b 需要对不同的b解很多次，但是不能把所有b放在一个矩阵
里面用\同时解的话，有没有什么办法对A预处理一下，让它解得更快一些？
A就是一个二维的finite difference PDE矩阵，有五个diagonals,但是没有连在一起。 | c*******h
发帖数: 1096 | 2 do a red-black ordering. the schur complement is pentadiagonal.
【在 H*******k 的大作中提到】
: matlab 里面如果矩阵 Ax=b 需要对不同的b解很多次，但是不能把所有b放在一个矩阵
: 里面用\同时解的话，有没有什么办法对A预处理一下，让它解得更快一些？
: A就是一个二维的finite difference PDE矩阵，有五个diagonals,但是没有连在一起。
| H*******k
发帖数: 894 | 3 thanks! that should reduce the number of unknowns by half, but how can you
make it pentadiagonal? I still see five separated diagonal..
【在 c*******h 的大作中提到】
: do a red-black ordering. the schur complement is pentadiagonal.
|
|
