Commit bf763981 authored by Pierre Cazenave's avatar Pierre Cazenave

Make the debugging plot nicer.

parent 00bdbcbb
...@@ -102,17 +102,20 @@ all_nest_nodes = Mobj.relaxBC_nodes{:}; ...@@ -102,17 +102,20 @@ all_nest_nodes = Mobj.relaxBC_nodes{:};
Mobj.relaxnBC_elems = length(all_nest_elems); Mobj.relaxnBC_elems = length(all_nest_elems);
Mobj.relaxnBC_nodes = length(all_nest_nodes); Mobj.relaxnBC_nodes = length(all_nest_nodes);
% Check it's worked for the first model boundary.
% figure(1) % xc = nodes2elems(Mobj.x, Mobj);
% yc = nodes2elems(Mobj.y, Mobj);
% figure
% clf % clf
% triplot(Mobj.tri,Mobj.x,Mobj.y,'k'); % triplot(Mobj.tri,Mobj.x,Mobj.y,'k');
% hold on % hold on
% %
% plot(Mobj.x( Mobj.relaxBC_nodes{1}), Mobj.y( Mobj.relaxBC_nodes{1}), 'r.') % symbols = {'r.', 'k.', 'rx', 'kx', 'ro', 'ko'};
% plot(Mobj.x( Mobj.relaxBC_nodes{2}), Mobj.y( Mobj.relaxBC_nodes{2}), 'r.') % for nn = 1:length(Mobj.relaxBC_nodes)
% plot(Mobj.x( Mobj.relaxBC_nodes{3}), Mobj.y( Mobj.relaxBC_nodes{3}), 'kx') % plot(Mobj.x(Mobj.relaxBC_nodes{nn}), ...
% plot(Mobj.x( Mobj.relaxBC_nodes{4}), Mobj.y( Mobj.relaxBC_nodes{4}), 'rx') % Mobj.y(Mobj.relaxBC_nodes{nn}), ...
% % symbols{mod(nn, length(Mobj.relaxBC_nodes)) + 1})
% end
% plot(xc(Mobj.relaxBC_elems{3}), yc(Mobj.relaxBC_elems{3}), 'kx') % plot(xc(Mobj.relaxBC_elems{3}), yc(Mobj.relaxBC_elems{3}), 'kx')
% axis('equal', 'tight') % axis('equal', 'tight')
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment