Mercurial > repos > public > sbplib

comparison +rv/+diffops/addClosuresToDiffOp.m @ 1190:79618b58b0a0 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Refactor constructDiffops and remove closures from some of the diffops used to compute the residual.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 08 Jul 2019 15:12:19 +0200
parents 65a577db5ca0
children
comparison
equal deleted inserted replaced
1189:6cb03209f0a7 1190:79618b58b0a0
1 function [D, penalties] = addClosuresToDiffOp(diffOp, D, BCs) 1 function D = addClosuresToDiffOp(D, closures)
2 penalties = cell(size(BCs)); 2 for i = 1:size(closures,1)
3 for i = 1:size(BCs,1) 3 for j = 1:size(closures,2)
4 for j = 1:size(BCs,2) 4 closure = closures{i,j};
5 [closure, penalties{i,j}] = diffOp.boundary_condition(BCs{i,j}.boundary, BCs{i,j}.type);
6 if ~isa(closure, 'function_handle') 5 if ~isa(closure, 'function_handle')
7 closure = @(v) closure*v; 6 closure = @(v) closure*v;
8 end 7 end
9 D = @(v) D(v) + closure(v); 8 D = @(v) D(v) + closure(v);
10 end 9 end