Mercurial > repos > public > sbplib

view spyh.m @ 1295:cb053fabbedc feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add displacement BC on normal tangential components, but it is only provably stable if traction bc are used for the other component.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 02 Jul 2020 20:05:50 -0700
parents 12db86a8ec7b
children
line wrap: on
line source

% spyh mimics the built in spy but returns a handle which allows modifying the plot.
function h = spyh(A)
    [n,m] = size(A);
    [I,J] = find(A);

    if ~isempty(J)
        h = plot(J,I);
        h.LineStyle = 'none';
        h.Marker = '.';
        h.MarkerSize = 14;
    else
        h = [];
    end

    a = gca;
    xlim([0 m+1]);
    ylim([0 n+1]);
    axis square
    a.YDir = 'reverse';
end