Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dCurvilinearAnisotropic.m @ 1302:a0d615bde7f8 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
comparison
equal deleted inserted replaced
1295:cb053fabbedc 1302:a0d615bde7f8
64 64
65 methods 65 methods
66 66
67 % The coefficients can either be function handles or grid functions 67 % The coefficients can either be function handles or grid functions
68 % optFlag -- if true, extra computations are performed, which may be helpful for optimization. 68 % optFlag -- if true, extra computations are performed, which may be helpful for optimization.
69 function obj = Elastic2dCurvilinearAnisotropic(g, order, rho, C, opSet, optFlag) 69 function obj = Elastic2dCurvilinearAnisotropic(g, order, rho, C, opSet, optFlag, hollow)
70 default_arg('hollow', false);
70 default_arg('rho', @(x,y) 0*x+1); 71 default_arg('rho', @(x,y) 0*x+1);
71 default_arg('opSet',{@sbp.D2VariableCompatible, @sbp.D2VariableCompatible}); 72 default_arg('opSet',{@sbp.D2VariableCompatible, @sbp.D2VariableCompatible});
72 default_arg('optFlag', false); 73 default_arg('optFlag', false);
73 dim = 2; 74 dim = 2;
74 75
186 end 187 end
187 end 188 end
188 end 189 end
189 190
190 gRef = grid.equidistant([m_u, m_v], {0,1}, {0,1}); 191 gRef = grid.equidistant([m_u, m_v], {0,1}, {0,1});
191 refObj = scheme.Elastic2dVariableAnisotropic(gRef, order, rho_tilde, PHI, opSet); 192 refObj = scheme.Elastic2dVariableAnisotropic(gRef, order, rho_tilde, PHI, opSet, [], hollow);
192 193
193 %---- Set object properties ------ 194 %---- Set object properties ------
194 obj.RHO = spdiag(rho); 195 obj.RHO = spdiag(rho);
195 196
196 % Volume quadrature 197 % Volume quadrature