Mercurial > repos > public > sbplib

changeset 1007:9c8ed00732fd feature/getBoundaryOp

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Revert to old code in multiblock.DiffOp.getBoundaryOperator.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 16 Jan 2019 11:31:04 -0800
parents de8c979b9881
children 8d73fcdb07a5
files +multiblock/DiffOp.m
diffstat 1 files changed, 2 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
diff -r de8c979b9881 -r 9c8ed00732fd +multiblock/DiffOp.m
--- a/+multiblock/DiffOp.m	Wed Jan 16 16:50:47 2019 +0100
+++ b/+multiblock/DiffOp.m	Wed Jan 16 11:31:04 2019 -0800
@@ -129,20 +129,19 @@
 
         % Get a boundary operator specified by opName for the given boundary/BoundaryGroup
         function op = getBoundaryOperator(obj, opName, boundary)
-            blockmatrixDiv = obj.blockmatrixDiv{1};
 
             switch class(boundary)
                 case 'cell'
                     blockId = boundary{1};
                     localOp = obj.diffOps{blockId}.getBoundaryOperator(opName, boundary{2});
 
-                    div = {blockmatrixDiv, size(localOp,2)};
+                    div = {obj.blockmatrixDiv{1}, size(localOp,2)};
                     blockOp = blockmatrix.zero(div);
                     blockOp{blockId,1} = localOp;
                     op = blockmatrix.toMatrix(blockOp);
                     return
                 case 'multiblock.BoundaryGroup'
-                    op = sparse(sum(blockmatrixDiv),0);
+                    op = sparse(size(obj.D,1),0);
                     for i = 1:length(boundary)
                         op = [op, obj.getBoundaryOperator(opName, boundary{i})];
                     end