Mercurial > repos > public > sbplib

diff +multiblock/DiffOp.m @ 584:b0386d2c180d feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
In case of empty grids, make grid functions with the correct dimensions
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 08 Sep 2017 11:08:50 +0200
parents 55a7777dfcd0
children c360bbecf260
line wrap: on
line diff
--- a/+multiblock/DiffOp.m	Thu Sep 07 14:08:31 2017 +0200
+++ b/+multiblock/DiffOp.m	Fri Sep 08 11:08:50 2017 +0200
@@ -134,7 +134,7 @@
                     op = blockmatrix.toMatrix(blockOp);
                     return
                 case 'multiblock.BoundaryGroup'
-                    op = [];
+                    op = sparse(size(obj.D,1),0);
                     for i = 1:length(boundary)
                         op = [op, obj.getBoundaryOperator(opName, boundary{i})];
                     end
@@ -155,7 +155,7 @@
                 case 'multiblock.BoundaryGroup'
                     [n,m] = size(obj.D);
                     closure = sparse(n,m);
-                    penalty = [];
+                    penalty = sparse(n,0);
                     for i = 1:length(boundary)
                         [closurePart, penaltyPart] = obj.boundary_condition(boundary{i}, type);
                         closure = closure + closurePart;