Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dCurvilinearAnisotropic.m @ 1322:412b8ceafbc6 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add Zt to output args for Elastic2dCurvilinearAnisotropic.interfaceNormalTangential
author Martin Almquist <malmquist@stanford.edu>
date Sat, 24 Oct 2020 20:58:26 -0700
parents b40359c9faed
children df8c71b80c33
comparison
equal deleted inserted replaced
1321:b40359c9faed 1322:412b8ceafbc6
601 end 601 end
602 602
603 % Same interface conditions as in interfaceStandard, but imposed in the normal-tangential 603 % Same interface conditions as in interfaceStandard, but imposed in the normal-tangential
604 % coordinate system. For the isotropic case, the components decouple nicely. 604 % coordinate system. For the isotropic case, the components decouple nicely.
605 % The resulting scheme is not identical to that of interfaceStandard. This appears to be better. 605 % The resulting scheme is not identical to that of interfaceStandard. This appears to be better.
606 function [closure, penalty, forcingPenalties] = interfaceNormalTangential(obj,boundary,neighbour_scheme,neighbour_boundary,type) 606 function [closure, penalty, forcingPenalties, Zt] = interfaceNormalTangential(obj,boundary,neighbour_scheme,neighbour_boundary,type)
607 tuning = type.tuning; 607 tuning = type.tuning;
608 608
609 % u denotes the solution in the own domain 609 % u denotes the solution in the own domain
610 % v denotes the solution in the neighbour domain 610 % v denotes the solution in the neighbour domain
611 611