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
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