Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d2_variable_periodic_6.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
45 Mp2 = c(r-1)/0.20e2 - 0.3e1/0.10e2 * c(r+1) + c(r+3)/0.20e2 + 0.7e1/0.40e2 * c(r) + 0.7e1/0.40e2 * c(r+2); 45 Mp2 = c(r-1)/0.20e2 - 0.3e1/0.10e2 * c(r+1) + c(r+3)/0.20e2 + 0.7e1/0.40e2 * c(r) + 0.7e1/0.40e2 * c(r+2);
46 Mp3 = c(r+1)/0.40e2 + c(r+2)/0.40e2 - 0.11e2/0.360e3 * c(r) - 0.11e2/0.360e3 * c(r+3); 46 Mp3 = c(r+1)/0.40e2 + c(r+2)/0.40e2 - 0.11e2/0.360e3 * c(r) - 0.11e2/0.360e3 * c(r+3);
47 47
48 vals = [Mm3,Mm2,Mm1,M0,Mp1,Mp2,Mp3]; 48 vals = [Mm3,Mm2,Mm1,M0,Mp1,Mp2,Mp3];
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 end 54 end
55 D2 = @D2_fun; 55 D2 = @D2_fun;
56 56
57 57
58 end 58 end