Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d2_variable_periodic_4.m @ 943:21394c78c72e feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default
author Martin Almquist <malmquist@stanford.edu>
date Tue, 04 Dec 2018 15:24:36 -0800
parents bbf303c1f0cf
children
comparison
equal deleted inserted replaced
942:35701c85e356 943:21394c78c72e
28 D1 = HI*(Q - 1/2*e_l*e_l' + 1/2*e_r*e_r'); 28 D1 = HI*(Q - 1/2*e_l*e_l' + 1/2*e_r*e_r');
29 29
30 30
31 scheme_width = 5; 31 scheme_width = 5;
32 scheme_radius = (scheme_width-1)/2; 32 scheme_radius = (scheme_width-1)/2;
33 33
34 r = 1:m; 34 r = 1:m;
35 offset = scheme_width; 35 offset = scheme_width;
36 r = r + offset; 36 r = r + offset;
37 37
38 function D2 = D2_fun(c) 38 function D2 = D2_fun(c)
45 Mp1 = 0 * c(r-2) + 1/6*c(r-1) + 1/2*c(r) + 1/2*c(r+1) + 1/6 *c(r+2); 45 Mp1 = 0 * c(r-2) + 1/6*c(r-1) + 1/2*c(r) + 1/2*c(r+1) + 1/6 *c(r+2);
46 Mp2 = 0 * c(r-2) + 0 * c(r-1) - 1/8*c(r) + 1/6*c(r+1) - 1/8 *c(r+2); 46 Mp2 = 0 * c(r-2) + 0 * c(r-1) - 1/8*c(r) + 1/6*c(r+1) - 1/8 *c(r+2);
47 47
48 vals = -[Mm2,Mm1,M0,Mp1,Mp2]; 48 vals = -[Mm2,Mm1,M0,Mp1,Mp2];
49 diags = -scheme_radius : scheme_radius; 49 diags = -scheme_radius : scheme_radius;
50 M = spdiagsVariablePeriodic(vals,diags); 50 M = spdiagsPeriodic(vals,diags);
51 51
52 M=M/h; 52 M=M/h;
53 D2=HI*(-M ); 53 D2=HI*(-M );
54 54
55 end 55 end