Mercurial > repos > public > sbplib

annotate spdiag.m @ 1345:14f44e81e1e3 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add scheme for Elastic2dVariable, not fully compatible.
author Martin Almquist <martin.almquist@it.uu.se>
date Tue, 30 Apr 2024 13:33:39 +0200
parents e76321b89c1e
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
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
798
e76321b89c1e Fix bug in spdiag
Jonatan Werpers <jonatan@werpers.com>
parents: 555
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