Mercurial > repos > public > sbplib

diff +scheme/Elastic2dCurvilinearAnisotropic.m @ 1209:67eee83fd9c9 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Swap indices in anisotropic curvilinear
author Martin Almquist <malmquist@stanford.edu>
date Fri, 20 Sep 2019 15:38:39 -0700
parents 3258dca12af8
children 3d7faa2ca312
line wrap: on
line diff
--- a/+scheme/Elastic2dCurvilinearAnisotropic.m	Fri Sep 20 15:00:19 2019 -0700
+++ b/+scheme/Elastic2dCurvilinearAnisotropic.m	Fri Sep 20 15:38:39 2019 -0700
@@ -166,7 +166,7 @@
                             PHI{i,j,k,l} = 0*C{i,j,k,l};
                             for m = 1:dim
                                 for n = 1:dim
-                                    PHI{i,j,k,l} = PHI{i,j,k,l} + J.*K{m,j}.*C{i,m,n,l}.*K{n,k};
+                                    PHI{i,j,k,l} = PHI{i,j,k,l} + J.*K{m,i}.*C{m,j,n,l}.*K{n,k};
                                 end
                             end
                         end