Mercurial > repos > public > sbplib

view spdiag.m @ 695:ba0d31ce4121 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Changed jacobian at right side of the penalty term in 2D, seems like it is working
author Ylva Rydin <ylva.rydin@telia.com>
date Tue, 19 Sep 2017 09:51:08 +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