Mercurial > repos > public > sbplib
annotate spdiag.m @ 577:e45c9b56d50d feature/grids
Add an Empty grid class
The need turned up for the flexural code when we may or may not have a grid for the open water and want to plot that solution.
In case there is no open water we need an empty grid to plot the empty gridfunction against to avoid errors.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Thu, 07 Sep 2017 09:16:12 +0200 |
parents | 8e0067aef5b0 |
children | e76321b89c1e |
rev | line source |
---|---|
45
b6cfc04e576c
Changed spdiag to be more like matlabs diag.
Jonatan Werpers <jonatan@werpers.com>
parents:
27
diff
changeset
|
1 function A = spdiag(a,i) |
81
9c0192cf099f
Added default argument to spdiag.
Jonatan Werpers <jonatan@werpers.com>
parents:
51
diff
changeset
|
2 default_arg('i',0); |
555
8e0067aef5b0
Make spdiag accept both row and column vectors
Jonatan Werpers <jonatan@werpers.com>
parents:
81
diff
changeset
|
3 |
8e0067aef5b0
Make spdiag accept both row and column vectors
Jonatan Werpers <jonatan@werpers.com>
parents:
81
diff
changeset
|
4 if isrow(a) |
8e0067aef5b0
Make spdiag accept both row and column vectors
Jonatan Werpers <jonatan@werpers.com>
parents:
81
diff
changeset
|
5 a = a'; |
8e0067aef5b0
Make spdiag accept both row and column vectors
Jonatan Werpers <jonatan@werpers.com>
parents:
81
diff
changeset
|
6 end |
8e0067aef5b0
Make spdiag accept both row and column vectors
Jonatan Werpers <jonatan@werpers.com>
parents:
81
diff
changeset
|
7 |
51
0be702829bb9
Generalized spdiag. Added function for createing a banded sparse matrix.
Jonatan Werpers <jonatan@werpers.com>
parents:
45
diff
changeset
|
8 n = length(a)-abs(i); |
45
b6cfc04e576c
Changed spdiag to be more like matlabs diag.
Jonatan Werpers <jonatan@werpers.com>
parents:
27
diff
changeset
|
9 A = spdiags(a,i,n,n); |
27
97a638f91fb8
Added function spdiag(). Fixed a bunch of bugs in the Wave2dCurve scheme.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
10 end |