Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d2_variable_periodic_2.m @ 890:c70131daaa6e feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with feature/poroelastic.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 21 Nov 2018 18:29:29 -0800
parents bbf303c1f0cf
children
comparison
equal deleted inserted replaced
885:18e10217dca9 890:c70131daaa6e
25 D1 = stripeMatrixPeriodic(stencil, diags, m); 25 D1 = stripeMatrixPeriodic(stencil, diags, m);
26 D1 = D1/h; 26 D1 = D1/h;
27 27
28 scheme_width = 3; 28 scheme_width = 3;
29 scheme_radius = (scheme_width-1)/2; 29 scheme_radius = (scheme_width-1)/2;
30 30
31 r = 1:m; 31 r = 1:m;
32 offset = scheme_width; 32 offset = scheme_width;
33 r = r + offset; 33 r = r + offset;
34 34
35 function D2 = D2_fun(c) 35 function D2 = D2_fun(c)
39 M0 = c(r-1)/2 + c(r) + c(r+1)/2; 39 M0 = c(r-1)/2 + c(r) + c(r+1)/2;
40 Mp1 = -c(r)/2 - c(r+1)/2; 40 Mp1 = -c(r)/2 - c(r+1)/2;
41 41
42 vals = [Mm1,M0,Mp1]; 42 vals = [Mm1,M0,Mp1];
43 diags = -scheme_radius : scheme_radius; 43 diags = -scheme_radius : scheme_radius;
44 M = spdiagsVariablePeriodic(vals,diags); 44 M = spdiagsPeriodic(vals,diags);
45 45
46 M=M/h; 46 M=M/h;
47 D2=HI*(-M ); 47 D2=HI*(-M );
48 end 48 end
49 D2 = @D2_fun; 49 D2 = @D2_fun;