Mercurial > repos > public > sbplib

view spdiag.m @ 725:e9e15d64f9f9 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Implement symmetric Dirichlet BC with variable tuning parameter in Heat2D scheme.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 30 Mar 2018 15:35:46 -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