Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d2_variable_periodic_4.m @ 801:bbf303c1f0cf feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Rename spdaigsVariablePeriodic spdiagsPeriodic
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 26 Jul 2018 18:04:29 -0700
parents 5ccf6aaf6d6b
children
comparison
equal deleted inserted replaced
800:87ea9cac3287 801:bbf303c1f0cf
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