Mercurial > repos > public > sbplib
comparison +multiblock/DiffOp.m @ 731:c3b460a7d7ac feature/poroelastic
Hack in multiblockDiffOp, probably not Werpers-compatible.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Sat, 21 Apr 2018 16:08:48 -0700 |
| parents | aa8cf3851de8 |
| children | 624120275bab |
comparison
equal
deleted
inserted
replaced
| 730:9f28cf266f86 | 731:c3b460a7d7ac |
|---|---|
| 161 closure = sparse(n,m); | 161 closure = sparse(n,m); |
| 162 penalty = sparse(n,0); | 162 penalty = sparse(n,0); |
| 163 for i = 1:length(boundary) | 163 for i = 1:length(boundary) |
| 164 [closurePart, penaltyPart] = obj.boundary_condition(boundary{i}, type); | 164 [closurePart, penaltyPart] = obj.boundary_condition(boundary{i}, type); |
| 165 closure = closure + closurePart; | 165 closure = closure + closurePart; |
| 166 if iscell(penaltyPart) | |
| 167 penaltyPart = blockmatrix.toMatrix(penaltyPart); | |
| 168 end | |
| 166 penalty = [penalty, penaltyPart]; | 169 penalty = [penalty, penaltyPart]; |
| 167 end | 170 end |
| 168 otherwise | 171 otherwise |
| 169 error('Unknown boundary indentifier') | 172 error('Unknown boundary indentifier') |
| 170 end | 173 end |