plot_dofmap.m

function p = plot_dofmap(df_info,params)
%function p = plot_dofmap(df_info,params)
%
% function plotting the numbers of the dofs for optical consistency
% check across edges

% B. Haasdonk 12.1.2011

p = plot(df_info.grid);
hold on
lagrange_nodes = Fem.Lagrange.nodes_lcoord(df_info.pdeg);
for li = 1:size(lagrange_nodes,1);
  clocal = repmat([1/3,1/3],size(lagrange_nodes,1),1);
  % shrink location a bit towards center for readability;
  local = (lagrange_nodes - clocal)*0.8 + clocal;
  glob = local2global(df_info.grid,1:df_info.grid.nelements,local(li,:));
  nums = num2str(df_info.global_dof_index(:,li));
  p2 = text(glob(:,1),glob(:,2),nums,'Fontsize',5);
end;
p = [p;p2];

end