Mercurial > repos > public > sbplib

view spdiag.m @ 1343:09a5783a3d37 feature/D2_boundary_opt

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve on default choice of number of grid points for each block
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 14 Oct 2022 14:42:32 +0200
parents e76321b89c1e
children
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