Mercurial > repos > public > sbplib

comparison +multiblock/DiffOp.m @ 237:0fc1de5cd85b feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
multigrid.DiffOp: Bug fix.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 15 Jul 2016 13:46:05 +0200
parents b8ee5212f651
children 2d48db7151cf
comparison
equal deleted inserted replaced
236:8592358a3935 237:0fc1de5cd85b
55 for i = 1:nBlocks 55 for i = 1:nBlocks
56 D{i,i} = obj.diffOps{i}.D; 56 D{i,i} = obj.diffOps{i}.D;
57 end 57 end
58 58
59 for i = 1:nBlocks 59 for i = 1:nBlocks
60 for j = i:nBlocks 60 for j = 1:nBlocks
61 intf = grid.connections{i,j}; 61 intf = grid.connections{i,j};
62 if isempty(intf) 62 if isempty(intf)
63 continue 63 continue
64 end 64 end
65
65 66
66 [ii, ij] = obj.diffOps{i}.interface(intf{1}, obj.diffOps{j}, intf{2}); 67 [ii, ij] = obj.diffOps{i}.interface(intf{1}, obj.diffOps{j}, intf{2});
67 D{i,i} = D{i,i} + ii; 68 D{i,i} = D{i,i} + ii;
68 D{i,j} = D{i,j} + ij; 69 D{i,j} = D{i,j} + ij;
69 70