Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d2_variable_periodic_6.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
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