Mercurial > repos > public > sbplib

diff +multiblock/DiffOp.m @ 868:57760d7088ad bcSetupExperiment

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with feature/grids
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 27 Jul 2018 10:39:12 -0700
parents 1cc5a0d26453 5cf9fdf4c98f
children 386ef449df51 21394c78c72e
line wrap: on
line diff
--- a/+multiblock/DiffOp.m	Thu Jul 26 17:28:48 2018 -0700
+++ b/+multiblock/DiffOp.m	Fri Jul 27 10:39:12 2018 -0700
@@ -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