Mercurial > repos > public > sbplib
annotate spyh.m @ 1031:2ef20d00b386 feature/advectionRV
For easier comparison, return both the first order and residual viscosity when evaluating the residual. Add the first order and residual viscosity to the state of the RungekuttaRV time steppers
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Thu, 17 Jan 2019 10:25:06 +0100 |
parents | 12db86a8ec7b |
children |
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 |