Mercurial > repos > public > sbplib

annotate stripeMatrix.m @ 549:ae905a11e32c feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change away from matrix formulation in Rk4 second order form The new formulation can be up to 8 times faster
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 28 Aug 2017 11:06:00 +0200
parents 5cb99f56f55e
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
51
0be702829bb9 Generalized spdiag. Added function for createing a banded sparse matrix.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 % Creates a matrix of size n,m with the values of val on the diagonals diag.
358
5cb99f56f55e Added some utility functions. Improved some documentation.
Jonatan Werpers <jonatan@werpers.com>
parents: 51
diff changeset
2 % A = stripeMatrix(val,diags,n,m)
51
0be702829bb9 Generalized spdiag. Added function for createing a banded sparse matrix.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 function A = stripeMatrix(val,diags,n,m)
0be702829bb9 Generalized spdiag. Added function for createing a banded sparse matrix.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 default_arg('m',n);
0be702829bb9 Generalized spdiag. Added function for createing a banded sparse matrix.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5
0be702829bb9 Generalized spdiag. Added function for createing a banded sparse matrix.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 D = ones(n,1)*val;
0be702829bb9 Generalized spdiag. Added function for createing a banded sparse matrix.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 A = spdiags(D,diags,n,m);
0be702829bb9 Generalized spdiag. Added function for createing a banded sparse matrix.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8 end