Mercurial > repos > public > sbplib
annotate stripeMatrix.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 | 5cb99f56f55e |
children |
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 |