Mercurial > repos > public > sbplib
diff +scheme/Elastic2dCurvilinearAnisotropic.m @ 1302:a0d615bde7f8 feature/poroelastic
Add the hollow option to the anisotropic diffops
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Fri, 10 Jul 2020 20:24:23 -0700 |
| parents | cb053fabbedc |
| children | 633757e582e5 |
line wrap: on
line diff
--- a/+scheme/Elastic2dCurvilinearAnisotropic.m Thu Jul 02 20:05:50 2020 -0700 +++ b/+scheme/Elastic2dCurvilinearAnisotropic.m Fri Jul 10 20:24:23 2020 -0700 @@ -66,7 +66,8 @@ % The coefficients can either be function handles or grid functions % optFlag -- if true, extra computations are performed, which may be helpful for optimization. - function obj = Elastic2dCurvilinearAnisotropic(g, order, rho, C, opSet, optFlag) + function obj = Elastic2dCurvilinearAnisotropic(g, order, rho, C, opSet, optFlag, hollow) + default_arg('hollow', false); default_arg('rho', @(x,y) 0*x+1); default_arg('opSet',{@sbp.D2VariableCompatible, @sbp.D2VariableCompatible}); default_arg('optFlag', false); @@ -188,7 +189,7 @@ end gRef = grid.equidistant([m_u, m_v], {0,1}, {0,1}); - refObj = scheme.Elastic2dVariableAnisotropic(gRef, order, rho_tilde, PHI, opSet); + refObj = scheme.Elastic2dVariableAnisotropic(gRef, order, rho_tilde, PHI, opSet, [], hollow); %---- Set object properties ------ obj.RHO = spdiag(rho);