Mercurial > repos > public > sbplib

view stripeMatrix.m @ 1083:a22b23098021 feature/laplace1dUpdate

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Divide interface penalty strength by 4, to make it sharp. Stability verified by eigenvalue computations.
author Martin Almquist <malmquist@stanford.edu>
date Tue, 05 Mar 2019 12:41:50 -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