我目前正在使用matlab进行FEM。
我的代码就像这样
function [A,rhs] = Assemble()
pre_process();
[IA,JA,A,rhs] = assemble(pre_processed_parameters); // mex-func
// IA : row index
// JA : col index
// A : value
A = sparse(IA,JA,A);
当A
的大小为ten million x ten million
时,
pre_process
不到1秒
assemble
大约需要5秒
sparse
大约需要820秒
主要耗时行是一行sparse
。
有没有更快的方法来收集它?
IA
,JA
和A
的属性
一对(IA[k],JA[k])
具有所有k
的所有不同值。
结果矩阵A
是下三角形并具有块结构