Mercurial > repos > public > sbplib

changeset 1211:3d7faa2ca312 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add physical derivatives Dx and Dy to ElasticCurvilinearAnisotropic
author Martin Almquist <malmquist@stanford.edu>
date Fri, 27 Sep 2019 13:18:09 -0700
parents 67eee83fd9c9
children 5c5815af4b7a
files +scheme/Elastic2dCurvilinearAnisotropic.m
diffstat 1 files changed, 6 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Elastic2dCurvilinearAnisotropic.m	Fri Sep 20 15:38:39 2019 -0700
+++ b/+scheme/Elastic2dCurvilinearAnisotropic.m	Fri Sep 27 13:18:09 2019 -0700
@@ -23,6 +23,8 @@
 
         D  % Total operator
 
+        Dx, Dy % Physical derivatives
+
         % Boundary operators in cell format, used for BC
         T_w, T_e, T_s, T_n
 
@@ -155,6 +157,10 @@
             K{2,1} = -x_v./J;
             K{2,2} = x_u./J;
 
+            % Physical derivatives
+            obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv;
+            obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;
+
             % Wrap around Aniosotropic Cartesian
             rho_tilde = J.*rho;