声明稀疏矩阵太慢了

时间:2016-07-12 14:03:19

标签: matlab sparse-matrix finite-element-analysis

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

有没有更快的方法来收集它?

IAJAA的属性

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

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

0 个答案:

没有答案