Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dCurvilinearAnisotropic.m @ 1272:15865fbda16e feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Set default stiffness tensors to 0 in CurvilinearAnisotropic schemes
author Martin Almquist <malmquist@stanford.edu>
date Sun, 31 May 2020 21:59:09 -0700
parents e1d4cb8b5309
children a8e730db76e9
comparison
equal deleted inserted replaced
1271:b5025bd67be1 1272:15865fbda16e
71 C_default = cell(dim,dim,dim,dim); 71 C_default = cell(dim,dim,dim,dim);
72 for i = 1:dim 72 for i = 1:dim
73 for j = 1:dim 73 for j = 1:dim
74 for k = 1:dim 74 for k = 1:dim
75 for l = 1:dim 75 for l = 1:dim
76 C_default{i,j,k,l} = @(x,y) 0*x + 1; 76 C_default{i,j,k,l} = @(x,y) 0*x ;
77 end 77 end
78 end 78 end
79 end 79 end
80 end 80 end
81 default_arg('C', C_default); 81 default_arg('C', C_default);