1

私は現在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_process1秒もかかりません

assemble約5秒かかります

sparse約820秒かかります

時間のかかる主な行は 1 行sparseです。

早く集める方法はありますか？

IA、JAおよびのプロパティAは

  1. のペアは(IA[k],JA[k])、 all に対してすべて異なる値を持ちますk。

  2. 結果の行列Aは下三角であり、ブロック構造を持っています

4

0 に答える 0