Mercurial > repos > public > sbplib
annotate spyh.m @ 774:66eb4a2bbb72 feature/grids
Remove default scaling of the system.
The scaling doens't seem to help actual solutions. One example that fails in the flexural code.
With large timesteps the solutions seems to blow up. One particular example is profilePresentation
on the tdb_presentation_figures branch with k = 0.0005
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Wed, 18 Jul 2018 15:42:52 -0700 |
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 |