Mercurial > repos > public > sbplib

view stencilEquation.m @ 770:0090a86d8b72 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve mononomial and vandermonde functions to work in multiple dimension and adapt to sym or double inputs
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 27 Jun 2018 11:10:52 +0200
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