Mercurial > repos > public > sbplib
diff +sbp/D2Variable.m @ 670:52a9dedb9171 feature/poroelastic
Add 6th order to D2variable.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Fri, 22 Dec 2017 16:38:08 +0100 |
| parents | 43ea848b6aa1 |
| children | 738de3a4058b |
line wrap: on
line diff
--- a/+sbp/D2Variable.m Fri Dec 22 14:16:17 2017 +0100 +++ b/+sbp/D2Variable.m Fri Dec 22 16:38:08 2017 +0100 @@ -26,6 +26,17 @@ obj.x = linspace(x_l,x_r,m)'; switch order + + case 6 + + [obj.H, obj.HI, obj.D1, obj.D2, ... + ~, obj.e_l, obj.e_r, ~, ~, ~, ~, ~,... + obj.d1_l, obj.d1_r] = ... + sbp.implementations.d4_variable_6(m, obj.h); + obj.borrowing.M.d1 = 0.1878; + obj.borrowing.R.delta_D = 0.3696; + % Borrowing e^T*D1 - d1 from R + case 4 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,... obj.e_r, obj.d1_l, obj.d1_r] = ...