Mercurial > repos > public > sbplib

view stripeMatrix.m @ 678:06676c40e77f feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add scheme for complete elastic operator. Traction BC working with MMS variable coefficient.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 05 Feb 2018 11:06:15 -0800
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