Mercurial > repos > public > sbplib

view spdiag.m @ 549:ae905a11e32c feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change away from matrix formulation in Rk4 second order form The new formulation can be up to 8 times faster
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 28 Aug 2017 11:06:00 +0200
parents 9c0192cf099f
children 8e0067aef5b0
line wrap: on
line source

function A = spdiag(a,i)
    default_arg('i',0);
    n = length(a)-abs(i);
    A = spdiags(a,i,n,n);
end