Mercurial > repos > public > sbplib

diff +multiblock/DiffOp.m @ 731:c3b460a7d7ac feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
line wrap: on
line diff
--- a/+multiblock/DiffOp.m	Sat Apr 21 16:04:46 2018 -0700
+++ b/+multiblock/DiffOp.m	Sat Apr 21 16:08:48 2018 -0700
@@ -163,6 +163,9 @@
                     for i = 1:length(boundary)
                         [closurePart, penaltyPart] = obj.boundary_condition(boundary{i}, type);
                         closure = closure + closurePart;
+                        if iscell(penaltyPart)
+                            penaltyPart = blockmatrix.toMatrix(penaltyPart);
+                        end
                         penalty = [penalty, penaltyPart];
                     end
                 otherwise