Mercurial > repos > public > sbplib

view spdiag.m @ 761:8ed102db8e9c feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add discont interface in 1D acoustics staggered, using the hat variable interface coupling.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 18 Jun 2018 16:36:24 -0700
parents 8e0067aef5b0
children e76321b89c1e
line wrap: on
line source

function A = spdiag(a,i)
    default_arg('i',0);

    if isrow(a)
        a = a';
    end

    n = length(a)-abs(i);
    A = spdiags(a,i,n,n);
end