声明稀疏矩阵太慢了

Question

我目前正在使用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的属性

1. 一对(IA[k],JA[k])具有所有k的所有不同值。

2. 结果矩阵A是下三角形并具有块结构