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.
|
|