I noticed this post, where spalloc and sparse are recommanded for efficient assembly in Matlab. I personally use sparse assembling for simple cases.However, when it comes to the case of coupled PDE, say , 3-PDE coupled, then the scalar unknow becomes a 3x3 tensor, in this case I can't figure out a way to exploit sparse(), and have to use for loops to assemble.when the final assembled sparse matrix is as large as 30k*30k(on PC), the assembling process becomes really slow(~10min), while the final matrix solving step is still fast(~less than 3 seconds).is there any suggestion? any generic solution is appreciated(not necessarily Matlab)
Efficient assembly of finite element matrix(coupled equations case)
finite element;matlab;finite difference;linear programming;assembly