Efficient implementation of adaptive P1-FEM in Matlab
Efficient implementation of adaptive P1-FEM in Matlab
复制标题
DOI：
10.2478/cmam-2011-0026
复制
发表时间：
2011
期刊：
影响因子：
--
通讯作者：
S. Funken;D. Praetorius;P. Wissgott
中科院分区：
文献类型：
--
作者：
S. Funken;D. Praetorius;P. Wissgott
摘要我们提供了一个自适应P1有限元法的MATLAB软件包p1afem。这包括用于组装数据的功能、不同的误差估计器和基于指标的自适应网格细化算法。在整个过程中，重点是通过使用MATLAB内置函数和矢量化的有效实现。数值实验强调的效率的代码被观察到几乎是线性的复杂度相对于运行时。虽然本文的范围是AFEM，但其一般思想可以理解为编写高效MATLAB代码的指导方针。
Abstract We provide a MATLAB package p1afem for an adaptive P1-finite element method (AFEM). This includes functions for the assembly of the data, different error estimators, and an indicator-based adaptive meshrefining algorithm. Throughout, the focus is on an efficient realization by use of MATLAB built-in functions and vectorization. Numerical experiments underline the efficiency of the code which is observed to be of almost linear complexity with respect to the runtime. Although the scope of this paper is on AFEM, the general ideas can be understood as a guideline for writing efficient MATLAB code.