Mercurial > repos > public > sbplib

comparison +multiblock/DiffOp.m @ 618:c360bbecf260 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add some todos
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 20 Oct 2017 15:57:49 +0200
parents b0386d2c180d
children 111fcbcff2e9 c29eca357028 aa8cf3851de8 703183ed8c8b
comparison
equal deleted inserted replaced
617:4ced7d47bd1f 618:c360bbecf260
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);