Mercurial > repos > public > sbplib

view stencilEquation.m @ 1036:8a9393084b30 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change argument order to the "correct" order, i.e providing diffOp specific parameters before the opSet.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 18 Jan 2019 08:58:26 +0100
parents 184833fe4c0e
children
line wrap: on
line source

% Find the equation for the stencil for d^k/dx^k
function [A,b] = stencilEquation(k, offsets, order)
    q = sym('q', [1, length(offsets)]);

    p = 0:(order-1+k);

    v     = vandermonde(offsets, p);
    vdiff = vandermonde(      0, p-k);

    eq = q*v == vdiff;

    [A,b] = equationsToMatrix(eq, q);
end