Mercurial > repos > public > sbplib

view stencilEquation.m @ 975:47b48a97c675 feature/getBoundaryOperator

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Stupid merge makes this branch dangerous. Merging with other branches will result in unexpected changes. Never merge this branch with anything. Closing branch.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 03 Jan 2019 13:24:39 +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