ed_interpolate.m

解决麦克斯韦的matlab源文件

function vals = ed_interpolate(fnct,mesh)
% Local edge element interpolation of a vectorfield 
%
% 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(mesh.Ne,1);
fnstr = [fnct,'(x)'];

for i=1:mesh.Ne
  vidx = mesh.ep(i,:);
  edgeo = mesh.vt(vidx,:)';
  ed = edgeo(:,2)-edgeo(:,1);
  x = 0.5*sum(edgeo')';
  field = eval(fnstr);
  vals(i) = dot(ed,field);
end