VectorizedCodeEllip​ticParabolicPHFEM

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