Mercurial > repos > public > sbplib

changeset 964:99c2ef883dc6 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix warning in getBoundaryWrapper methods
author Martin Almquist <malmquist@stanford.edu>
date Wed, 19 Dec 2018 07:09:24 +0100
parents c75ddd568fcc
children db3411264b96
files +multiblock/DiffOp.m
diffstat 1 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/+multiblock/DiffOp.m	Wed Dec 19 06:58:10 2018 +0100
+++ b/+multiblock/DiffOp.m	Wed Dec 19 07:09:24 2018 +0100
@@ -150,7 +150,7 @@
 
         % Get a boundary operator specified by opName for the given boundary/BoundaryGroup
         function op = getBoundaryOperatorWrapper(obj, opName, boundary, blockmatrixDiv)
-            default_arg('blockmatrixDiv', obj.blockmatrixDiv);
+            default_arg('blockmatrixDiv', obj.blockmatrixDiv{1});
 
             switch class(boundary)
                 case 'cell'
@@ -163,7 +163,7 @@
                     op = blockmatrix.toMatrix(blockOp);
                     return
                 case 'multiblock.BoundaryGroup'
-                    op = [];
+                    op = sparse(sum(blockmatrixDiv),0);
                     for i = 1:length(boundary)
                         op = [op, obj.getBoundaryOperatorWrapper(opName, boundary{i}, blockmatrixDiv)];
                     end
@@ -174,7 +174,7 @@
 
         % Get a boundary cell of operators, specified by opName for the given boundary/BoundaryGroup
         function opCell = getBoundaryCellOperator(obj, opName, boundary, blockmatrixDiv)
-            default_arg('blockmatrixDiv', obj.blockmatrixDiv);
+            default_arg('blockmatrixDiv', obj.blockmatrixDiv{1});
 
             % Get size of cell
             switch class(boundary)
@@ -204,7 +204,7 @@
                         opCell{i} = op;
 
                     case 'multiblock.BoundaryGroup'
-                        op = sparse(size(obj.D,1),0);
+                        op = sparse(sum(blockmatrixDiv),0);
                         for j = 1:length(boundary)
                             localCell = obj.getBoundaryCellOperator(opName, boundary{j}, blockmatrixDiv);
                             op = [op, localCell{i}];