Mercurial > repos > public > sbplib
comparison stencilEquation.m @ 832:5573913a0949 feature/burgers1d
Merged with default, and updated +scheme/Burgers1D accordingly
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Tue, 11 Sep 2018 15:58:35 +0200 |
| parents | 184833fe4c0e |
| children |
comparison
equal
deleted
inserted
replaced
| 831:d0934d1143b7 | 832:5573913a0949 |
|---|---|
| 1 % Find the equation for the stencil for d^k/dx^k | |
| 2 function [A,b] = stencilEquation(k, offsets, order) | |
| 3 q = sym('q', [1, length(offsets)]); | |
| 4 | |
| 5 p = 0:(order-1+k); | |
| 6 | |
| 7 v = vandermonde(offsets, p); | |
| 8 vdiff = vandermonde( 0, p-k); | |
| 9 | |
| 10 eq = q*v == vdiff; | |
| 11 | |
| 12 [A,b] = equationsToMatrix(eq, q); | |
| 13 end |