comparison +multiblock/DiffOp.m @ 630:31abb4b11133 feature/grids

Merge
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 20 Oct 2017 23:31:20 +0200
parents c360bbecf260
children 111fcbcff2e9 c29eca357028 aa8cf3851de8 703183ed8c8b
comparison
equal deleted inserted replaced
629:2b03ee11e48b 630:31abb4b11133
192 div{2} = size(blockPenalty, 2); % Penalty is a column vector 192 div{2} = size(blockPenalty, 2); % Penalty is a column vector
193 p = blockmatrix.zero(div); 193 p = blockmatrix.zero(div);
194 p{I} = blockPenalty; 194 p{I} = blockPenalty;
195 penalty = blockmatrix.toMatrix(p); 195 penalty = blockmatrix.toMatrix(p);
196 else 196 else
197 % TODO: used by beam equation, should be eliminated. SHould only set one BC per call
197 for i = 1:length(blockPenalty) 198 for i = 1:length(blockPenalty)
198 div{2} = size(blockPenalty{i}, 2); % Penalty is a column vector 199 div{2} = size(blockPenalty{i}, 2); % Penalty is a column vector
199 p = blockmatrix.zero(div); 200 p = blockmatrix.zero(div);
200 p{I} = blockPenalty{i}; 201 p{I} = blockPenalty{i};
201 penalty{i} = blockmatrix.toMatrix(p); 202 penalty{i} = blockmatrix.toMatrix(p);