Mercurial > repos > public > sbplib
changeset 639:021caf7194b8 feature/d1_staggered
Bug fix grid spacing in d1 staggered implementations.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Fri, 10 Nov 2017 15:57:47 -0800 |
| parents | 17dffd56f2fe |
| children | 4046b046c1f7 |
| files | +sbp/+implementations/d1_staggered_2.m +sbp/+implementations/d1_staggered_4.m +sbp/+implementations/d1_staggered_6.m |
| diffstat | 3 files changed, 6 insertions(+), 6 deletions(-) [+] |
line wrap: on
line diff
--- a/+sbp/+implementations/d1_staggered_2.m Fri Nov 10 15:57:21 2017 -0800 +++ b/+sbp/+implementations/d1_staggered_2.m Fri Nov 10 15:57:47 2017 -0800 @@ -4,11 +4,11 @@ % n : number of grid points: n on xp grid and n+1 on xm grid % L : interval length +h = L/(n-1); + n = n - 1; assert(n >= 3, 'Not enough grid points'); -h = L/(n-1); - qm20 = -1/4; qm01 = 1/2; qp11 = -1/4;
--- a/+sbp/+implementations/d1_staggered_4.m Fri Nov 10 15:57:21 2017 -0800 +++ b/+sbp/+implementations/d1_staggered_4.m Fri Nov 10 15:57:47 2017 -0800 @@ -4,11 +4,11 @@ % n : number of grid points: n on xp grid and n+1 on xm grid % L : interval length +h = L/(n-1); + n = n - 1; assert(n >= 7, 'Not enough grid points'); -h = L/(n-1); - % Q+ and Q-, top-left corner QpL = [... -323/378, 2783/3072, -1085/27648, -17/9216, -667/64512;
--- a/+sbp/+implementations/d1_staggered_6.m Fri Nov 10 15:57:21 2017 -0800 +++ b/+sbp/+implementations/d1_staggered_6.m Fri Nov 10 15:57:47 2017 -0800 @@ -4,11 +4,11 @@ % n : number of grid points: n on xp grid and n+1 on xm grid % L : interval length +h = L/(n-1); + n = n - 1; assert(n >= 9, 'Not enough grid points'); -h = L/(n-1); - % Coefficients determined such that the SBP property is satisfied qp44 = -144662434693018337/131093333881896960; pp3 = 793/720;