VectorizedCodeEllipticParabolicPHFEM
These programs are supplements to the paper " Vectorized implementation of primal hybrid FEM in MATLAB" by N. Harish et al.
Main_PH.m solves the second order elliptic equation with A=I, p=(1,1) and delta=1:
-nabla.(A nabla u + up)+ delta u = f in (0,1)^2
u=uD on the Dirichlet boundary
(A nabla u + up).n = g on the Neumann boundary.
ParabolicMain.m solves the second order parabolic problem with A=I, p=(1,1) and delta=1:
d/dt u - div(A nabla u+up)+ delta u = f in (0,1)^2,
u = u_D on the Dirichlet boundary
(A nabla u+up).n= g on the Neumann boundary
u0=0 Initial condition
引用格式
Sanjib Acharya (2026). VectorizedCodeEllipticParabolicPHFEM (https://ww2.mathworks.cn/matlabcentral/fileexchange/136359-vectorizedcodeellipticparabolicphfem), MATLAB Central File Exchange. 检索时间: .
