Mercurial > repos > public > sbplib
comparison +multiblock/DiffOp.m @ 796:aa1ed37a1b56 feature/poroelastic
Remove hack in multiblock/DiffOp that was used with cell penalties.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Wed, 25 Jul 2018 18:53:07 -0700 |
| parents | 624120275bab |
| children | 5cf9fdf4c98f |
comparison
equal
deleted
inserted
replaced
| 795:1f6b2fb69225 | 796:aa1ed37a1b56 |
|---|---|
| 162 closure = sparse(n,m); | 162 closure = sparse(n,m); |
| 163 penalty = sparse(n,0); | 163 penalty = sparse(n,0); |
| 164 for i = 1:length(boundary) | 164 for i = 1:length(boundary) |
| 165 [closurePart, penaltyPart] = obj.boundary_condition(boundary{i}, type); | 165 [closurePart, penaltyPart] = obj.boundary_condition(boundary{i}, type); |
| 166 closure = closure + closurePart; | 166 closure = closure + closurePart; |
| 167 if iscell(penaltyPart) | |
| 168 penaltyPart = blockmatrix.toMatrix(penaltyPart); | |
| 169 end | |
| 170 penalty = [penalty, penaltyPart]; | 167 penalty = [penalty, penaltyPart]; |
| 171 end | 168 end |
| 172 otherwise | 169 otherwise |
| 173 error('Unknown boundary indentifier') | 170 error('Unknown boundary indentifier') |
| 174 end | 171 end |