Mercurial > repos > public > sbplib
comparison +scheme/bcSetup.m @ 723:87436a107d8a feature/grids
Fix bug in bcSetup for inhom conditions
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Tue, 20 Mar 2018 14:43:04 +0100 |
parents | facb8ffb5c34 |
children | 1de60c4d462d |
comparison
equal
deleted
inserted
replaced
722:89e14a85a2d0 | 723:87436a107d8a |
---|---|
20 for i = 1:length(bc) | 20 for i = 1:length(bc) |
21 assertType(bc{i}, 'struct'); | 21 assertType(bc{i}, 'struct'); |
22 [localClosure, penalty] = diffOp.boundary_condition(bc{i}.boundary, bc{i}.type); | 22 [localClosure, penalty] = diffOp.boundary_condition(bc{i}.boundary, bc{i}.type); |
23 closure = closure + localClosure; | 23 closure = closure + localClosure; |
24 | 24 |
25 if ~isfield(bc,'data') || isempty(bc{i}.data) | 25 if ~isfield(bc{i},'data') || isempty(bc{i}.data) |
26 continue | 26 continue |
27 end | 27 end |
28 assertType(bc{i}.data, 'function_handle'); | 28 assertType(bc{i}.data, 'function_handle'); |
29 | 29 |
30 coord = diffOp.grid.getBoundary(bc{i}.boundary); | 30 coord = diffOp.grid.getBoundary(bc{i}.boundary); |