Mercurial > repos > public > sbplib

diff +multiblock/DiffOp.m @ 729:aa8cf3851de8 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update multiblock.DiffOp to work for systems.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 20 Apr 2018 16:56:49 -0700
parents c360bbecf260
children c3b460a7d7ac
line wrap: on
line diff
--- a/+multiblock/DiffOp.m	Fri Apr 20 12:06:38 2018 -0700
+++ b/+multiblock/DiffOp.m	Fri Apr 20 16:56:49 2018 -0700
@@ -53,7 +53,11 @@
 
 
             % Build the differentiation matrix
-            obj.blockmatrixDiv = {grid.Ns, grid.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;