Mercurial > repos > public > sbplib

changeset 668:43ea848b6aa1 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add borrowing constants to D2variable
author Martin Almquist <malmquist@stanford.edu>
date Fri, 22 Dec 2017 14:14:24 +0100
parents ed853945ee99
children 17e62551cdc2
files +sbp/D2Variable.m
diffstat 1 files changed, 7 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
diff -r ed853945ee99 -r 43ea848b6aa1 +sbp/D2Variable.m
--- a/+sbp/D2Variable.m	Wed Dec 20 06:52:17 2017 +0100
+++ b/+sbp/D2Variable.m	Fri Dec 22 14:14:24 2017 +0100
@@ -31,17 +31,23 @@
                         obj.e_r, obj.d1_l, obj.d1_r] = ...
                         sbp.implementations.d2_variable_4(m,obj.h);
                     obj.borrowing.M.d1 = 0.2505765857;
+
+                    obj.borrowing.R.delta_D = 0.577587500088313;
+                    % Borrowing e^T*D1 - d1 from R
                 case 2
                     [obj.H, obj.HI, obj.D1, obj.D2, obj.e_l,...
                         obj.e_r, obj.d1_l, obj.d1_r] = ...
                         sbp.implementations.d2_variable_2(m,obj.h);
                     obj.borrowing.M.d1 = 0.3636363636; 
                     % Borrowing const taken from Virta 2014
+
+                    obj.borrowing.R.delta_D = 1.000000538455350;
+                    % Borrowing e^T*D1 - d1 from R
                     
                 otherwise
                     error('Invalid operator order %d.',order);
             end
-
+            obj.borrowing.H11 = obj.H(1,1)/obj.h; % First element in H/h,
             obj.m = m;
             obj.M = [];
         end