p1_interpolate.m

来自「解决麦克斯韦的matlab源文件」· M 代码 · 共 20 行

function vals = p1_interpolate(fnct,mesh)
% interpolation of vectorfield with p.w. linear continuous vectorial
% finite elements.
%
% fnct -> string specifying a function of the form y = f(x), where
% both x and y are 2D column vectors
%
% mesh -> 2D triangulation
%

vals = zeros(2*mesh.Nv,1);
fnstr = [fnct,'(x)'];

for i=1:mesh.Nv
  x = mesh.vt(i,:)';
  y = eval(fnstr);
  vals([2*i-1;2*i],1) = y;
end