Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements

被引:39
作者
Anjam, I. [1 ]
Valdman, J. [2 ,3 ]
机构
[1] Univ Jyvaskyla, Dept Math Informat Technol, FI-40014 Jyvaskyla, Finland
[2] Univ South Bohemia, Inst Math & Biomath, Ceske Budejovice, Czech Republic
[3] ASCR, Inst Informat Theory & Automat, Prague, Czech Republic
关键词
MATLAB code vectorization; Finite element method; Edge element; Raviart-Thomas element; Nedelec element;
D O I
10.1016/j.amc.2015.03.105
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We propose an effective and flexible way to assemble finite element stiffness and mass matrices in MATLAB. We apply this for problems discretized by edge finite elements. Typical edge finite elements are Raviart-Thomas elements used in discretizations of H(div) spaces and Nedelec elements in discretizations of H(curl) spaces. We explain vectorization ideas and comment on a freely available MATLAB code which is fast and scalable with respect to time. (C) 2015 Elsevier Inc. All rights reserved.
引用
收藏
页码:252 / 263
页数:12
相关论文
共 18 条
[1]  
Anjam I., 2014, PREPRINT
[2]  
[Anonymous], ADV NUMER ANAL
[3]  
Bahriawati C., 2005, COMPUT METHODS APPL, V5, P333, DOI DOI 10.2478/CMAM-2005-0016)
[4]  
Brezzi Franco, 1991, SPRINGER SERIES COMP, V15
[5]  
Chen L., 2008, IFEM INTEGRATED FINI
[6]  
Cuvelier F., 2013, PREPRINT
[8]  
Funken Stefan, 2011, Computational Methods in Applied Mathematics, V11, P460
[9]  
Hannukainen A., 2012, PREPRINT
[10]   MIXED FINITE-ELEMENTS IN IR3 [J].
NEDELEC, JC .
NUMERISCHE MATHEMATIK, 1980, 35 (03) :315-341