Mercurial > repos > public > sbplib

diff +multiblock/DiffOp.m @ 1071:92cb03e64ca4 feature/grids/LaplaceSquared

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 20 Sep 2018 12:05:20 +0200
parents 5cf9fdf4c98f
children 57760d7088ad 7d4f57725192
line wrap: on
line diff
--- a/+multiblock/DiffOp.m	Tue Jul 24 21:08:08 2018 -0700
+++ b/+multiblock/DiffOp.m	Thu Sep 20 12:05:20 2018 +0200
@@ -53,7 +53,11 @@
 
 
             % Build the differentiation matrix
-            obj.blockmatrixDiv = {g.Ns, g.Ns};
+            Ns = zeros(nBlocks,1);
+            for i = 1:nBlocks
+                Ns(i) = length(obj.diffOps{i}.D);
+            end
+            obj.blockmatrixDiv = {Ns, Ns};
             D = blockmatrix.zero(obj.blockmatrixDiv);
             for i = 1:nBlocks
                 D{i,i} = obj.diffOps{i}.D;
@@ -117,7 +121,7 @@
 
         function ops = splitOp(obj, op)
             % Splits a matrix operator into a cell-matrix of matrix operators for
-            % each g.
+            % each grid.
             ops = sparse2cell(op, obj.NNN);
         end