Mercurial > repos > public > sbplib

comparison +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
comparison
equal deleted inserted replaced
728:0aff87f6fb2c 729:aa8cf3851de8
51 end 51 end
52 obj.H = blockmatrix.toMatrix(H); 52 obj.H = blockmatrix.toMatrix(H);
53 53
54 54
55 % Build the differentiation matrix 55 % Build the differentiation matrix
56 obj.blockmatrixDiv = {grid.Ns, grid.Ns}; 56 Ns = zeros(nBlocks,1);
57 for i = 1:nBlocks
58 Ns(i) = length(obj.diffOps{i}.D);
59 end
60 obj.blockmatrixDiv = {Ns, Ns};
57 D = blockmatrix.zero(obj.blockmatrixDiv); 61 D = blockmatrix.zero(obj.blockmatrixDiv);
58 for i = 1:nBlocks 62 for i = 1:nBlocks
59 D{i,i} = obj.diffOps{i}.D; 63 D{i,i} = obj.diffOps{i}.D;
60 end 64 end
61 65