Mercurial > repos > public > sbplib

comparison +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
comparison
equal deleted inserted replaced
1208:7f427275bc9c 1209:67eee83fd9c9
164 for k = 1:dim 164 for k = 1:dim
165 for l = 1:dim 165 for l = 1:dim
166 PHI{i,j,k,l} = 0*C{i,j,k,l}; 166 PHI{i,j,k,l} = 0*C{i,j,k,l};
167 for m = 1:dim 167 for m = 1:dim
168 for n = 1:dim 168 for n = 1:dim
169 PHI{i,j,k,l} = PHI{i,j,k,l} + J.*K{m,j}.*C{i,m,n,l}.*K{n,k}; 169 PHI{i,j,k,l} = PHI{i,j,k,l} + J.*K{m,i}.*C{m,j,n,l}.*K{n,k};
170 end 170 end
171 end 171 end
172 end 172 end
173 end 173 end
174 end 174 end