Mercurial > repos > public > sbplib
comparison +scheme/elasticShearVariable.m @ 677:eeaf9a00e304 feature/poroelastic
Correct factor of 2 in Dirichlet penalty parameter.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Fri, 26 Jan 2018 15:17:49 -0800 |
| parents | 9e1d2351f539 |
| children |
comparison
equal
deleted
inserted
replaced
| 676:9926efb39330 | 677:eeaf9a00e304 |
|---|---|
| 235 penalties = cell(obj.dim,obj.dim); | 235 penalties = cell(obj.dim,obj.dim); |
| 236 % Loop over components | 236 % Loop over components |
| 237 for i = 1:obj.dim | 237 for i = 1:obj.dim |
| 238 H_gamma_i = obj.H_boundary{i}; | 238 H_gamma_i = obj.H_boundary{i}; |
| 239 sigma_ij = tuning*delta(i,j)*2/(gamma*h(j)) +... | 239 sigma_ij = tuning*delta(i,j)*2/(gamma*h(j)) +... |
| 240 tuning*delta_b(i,j)*(2/(H11*h(j)) + 2/(H11*h(j)*phi)); | 240 tuning*delta_b(i,j)*(2/(H11*h(j)) + 1/(H11*h(j)*phi)); |
| 241 | 241 |
| 242 ci = E{i}*inv(RHO)*nj*Hi*... | 242 ci = E{i}*inv(RHO)*nj*Hi*... |
| 243 ( (e{j}*H_gamma*e{j}'*A*e{j}*d{j}')'*E{i}' + ... | 243 ( (e{j}*H_gamma*e{j}'*A*e{j}*d{j}')'*E{i}' + ... |
| 244 delta(i,j)*(e{j}*H_gamma*e{j}'*A*e{j}*d{j}')'*E{j}' ... | 244 delta(i,j)*(e{j}*H_gamma*e{j}'*A*e{j}*d{j}')'*E{j}' ... |
| 245 ) ... | 245 ) ... |