Mercurial > repos > public > sbplib
changeset 959:c226fb8c2b8a feature/poroelastic
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 |
| files | +scheme/Elastic2dVariable.m |
| diffstat | 1 files changed, 4 insertions(+), 3 deletions(-) [+] |
line wrap: on
line diff
--- a/+scheme/Elastic2dVariable.m Wed Dec 05 18:58:10 2018 -0800 +++ b/+scheme/Elastic2dVariable.m Wed Dec 05 19:30:09 2018 -0800 @@ -549,11 +549,12 @@ db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta alpha = cell(obj.dim, obj.dim); - for i = 1:obj.dim - for j = 1:obj.dim - alpha{i,j} = tuning*( d(i,j)* a_lambda*LAMBDA ... + alpha_func = @(i,j) tuning*( d(i,j)* a_lambda*LAMBDA ... + d(i,j)* a_mu_i*MU ... + db(i,j)*a_mu_ij*MU ); + for i = 1:obj.dim + for l = 1:obj.dim + alpha{i,l} = d(i,l)*alpha_func(i,j); end end