annotate spyh.m @ 958:72cd29107a9a feature/poroelastic

Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug!
author Martin Almquist <malmquist@stanford.edu>
date Wed, 05 Dec 2018 18:58:10 -0800
parents 12db86a8ec7b
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
358
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 % spyh mimics the built in spy but returns a handle which allows modifying the plot.
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 function h = spyh(A)
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 [n,m] = size(A);
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 [I,J] = find(A);
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5
443
12db86a8ec7b Fix bug in spyh() when the matix has no non-zero elements.
Jonatan Werpers <jonatan@werpers.com>
parents: 358
diff changeset
6 if ~isempty(J)
12db86a8ec7b Fix bug in spyh() when the matix has no non-zero elements.
Jonatan Werpers <jonatan@werpers.com>
parents: 358
diff changeset
7 h = plot(J,I);
12db86a8ec7b Fix bug in spyh() when the matix has no non-zero elements.
Jonatan Werpers <jonatan@werpers.com>
parents: 358
diff changeset
8 h.LineStyle = 'none';
12db86a8ec7b Fix bug in spyh() when the matix has no non-zero elements.
Jonatan Werpers <jonatan@werpers.com>
parents: 358
diff changeset
9 h.Marker = '.';
12db86a8ec7b Fix bug in spyh() when the matix has no non-zero elements.
Jonatan Werpers <jonatan@werpers.com>
parents: 358
diff changeset
10 h.MarkerSize = 14;
12db86a8ec7b Fix bug in spyh() when the matix has no non-zero elements.
Jonatan Werpers <jonatan@werpers.com>
parents: 358
diff changeset
11 else
12db86a8ec7b Fix bug in spyh() when the matix has no non-zero elements.
Jonatan Werpers <jonatan@werpers.com>
parents: 358
diff changeset
12 h = [];
12db86a8ec7b Fix bug in spyh() when the matix has no non-zero elements.
Jonatan Werpers <jonatan@werpers.com>
parents: 358
diff changeset
13 end
358
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15 a = gca;
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
16 xlim([0 m+1]);
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
17 ylim([0 n+1]);
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
18 axis square
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
19 a.YDir = 'reverse';
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
20 end