Mercurial > repos > public > sbplib
changeset 357:19495f126f06 feature/beams
Minor clean up of d4_variable_4.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Mon, 05 Dec 2016 13:10:44 +0100 |
parents | 85c2fe06d551 |
children | 0e198af1184b |
files | +sbp/+implementations/d4_variable_4.m |
diffstat | 1 files changed, 5 insertions(+), 5 deletions(-) [+] |
line wrap: on
line diff
--- a/+sbp/+implementations/d4_variable_4.m Wed Nov 16 16:16:44 2016 -0800 +++ b/+sbp/+implementations/d4_variable_4.m Mon Dec 05 13:10:44 2016 +0100 @@ -66,11 +66,11 @@ r = (1+scheme_radius):(m-scheme_radius); function D2 = D2_fun(c) - Mm2 = c(r-2)/0.8e1 - c(r-1)/0.6e1 + c(r) /0.8e1 ; - Mm1 = -c(r-2)/0.6e1 - c(r-1)/0.2e1 - c(r) /0.2e1 - c(r+1)/0.6e1 ; - M0 = c(r-2)/2.4e1 + c(r-1)/1.2e0 + c(r) * 0.3/0.4 + c(r+1)/1.2e0 + c(r+2)/2.4e1; - Mp1 = -c(r-1)/0.6e1 - c(r) /0.2e1 - c(r+1)/0.2e1 - c(r+2)/0.6e1; - Mp2 = c(r) /0.8e1 - c(r+1)/0.6e1 + c(r+2)/0.8e1; + Mm2 = c(r-2)/0.8e1 - c(r-1)/0.6e1 + c(r) /0.8e1 ; + Mm1 = -c(r-2)/0.6e1 - c(r-1)/0.2e1 - c(r) /0.2e1 - c(r+1)/0.6e1 ; + M0 = c(r-2)/2.4e1 + c(r-1)/1.2e0 + c(r)*0.3/0.4 + c(r+1)/1.2e0 + c(r+2)/2.4e1; + Mp1 = -c(r-1)/0.6e1 - c(r) /0.2e1 - c(r+1)/0.2e1 - c(r+2)/0.6e1; + Mp2 = c(r) /0.8e1 - c(r+1)/0.6e1 + c(r+2)/0.8e1; M(r,:) = spdiags([Mm2 Mm1 M0 Mp1 Mp2],0:2*scheme_radius,length(r),m);