Mercurial > repos > public > sbplib

view 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
line wrap: on
line source

% Creates a matrix of size n,m with the values of val on the diagonals diag.
%   A = stripeMatrix(val,diags,n,m)
function A = stripeMatrix(val,diags,n,m)
    default_arg('m',n);

    D = ones(n,1)*val;
    A = spdiags(D,diags,n,m);
end