Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dVariable.m @ 959:c226fb8c2b8a feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix in Elastic2dVariable.getBoundaryOps, alpha. Adjoint FD conv works now!
author Martin Almquist <malmquist@stanford.edu>
date Wed, 05 Dec 2018 19:30:09 -0800
parents 72cd29107a9a
children ac566f3dc9b3
comparison
equal deleted inserted replaced
958:72cd29107a9a 959:c226fb8c2b8a
547 547
548 d = @kroneckerDelta; % Kronecker delta 548 d = @kroneckerDelta; % Kronecker delta
549 db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta 549 db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta
550 alpha = cell(obj.dim, obj.dim); 550 alpha = cell(obj.dim, obj.dim);
551 551
552 for i = 1:obj.dim 552 alpha_func = @(i,j) tuning*( d(i,j)* a_lambda*LAMBDA ...
553 for j = 1:obj.dim
554 alpha{i,j} = tuning*( d(i,j)* a_lambda*LAMBDA ...
555 + d(i,j)* a_mu_i*MU ... 553 + d(i,j)* a_mu_i*MU ...
556 + db(i,j)*a_mu_ij*MU ); 554 + db(i,j)*a_mu_ij*MU );
555 for i = 1:obj.dim
556 for l = 1:obj.dim
557 alpha{i,l} = d(i,l)*alpha_func(i,j);
557 end 558 end
558 end 559 end
559 560
560 varargout{k} = alpha; 561 varargout{k} = alpha;
561 otherwise 562 otherwise