Mercurial > repos > public > sbplib

view stripeMatrix.m @ 735:eac5fb4b63db feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make bcSetup work for cell penalties and data from several BC per boundary. Should probably be changed before merge.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 25 Apr 2018 14:47:54 -0700
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