Mercurial > repos > public > sbplib

diff +scheme/Elastic2dCurvilinearAnisotropic.m @ 1306:633757e582e5 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add transformation gradient as property in CurvilinearAnisotropic
author Martin Almquist <malmquist@stanford.edu>
date Sun, 19 Jul 2020 20:29:48 -0700
parents a0d615bde7f8
children 6c308da9dcbc
line wrap: on
line diff
--- a/+scheme/Elastic2dCurvilinearAnisotropic.m	Mon Jul 13 14:21:41 2020 -0700
+++ b/+scheme/Elastic2dCurvilinearAnisotropic.m	Sun Jul 19 20:29:48 2020 -0700
@@ -23,6 +23,7 @@
 
         D  % Total operator
 
+        K % Transformation gradient
         Dx, Dy % Physical derivatives
         sigma % Cell matrix of physical stress operators
         n_w, n_e, n_s, n_n % Physical normals
@@ -164,6 +165,7 @@
             K{1,2} = -y_u./J;
             K{2,1} = -x_v./J;
             K{2,2} = x_u./J;
+            obj.K = K;
 
             % Physical derivatives
             obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv;