Mercurial > repos > public > sbplib

view stripeMatrix.m @ 407:e1db62d14835 feature/SBPInTimeGauss

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed bug with initial data and made gauss the default type
author Martin Almquist <martin.almquist@it.uu.se>
date Thu, 02 Feb 2017 19:30:16 +0100
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