Mercurial > repos > public > sbplib

view spdiag.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 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