Mercurial > repos > public > sbplib

comparison +scheme/bcSetup.m @ 759:2645188489f6 feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix? in +scheme/bcSetup.m
author Martin Almquist <malmquist@stanford.edu>
date Mon, 18 Jun 2018 16:32:23 -0700
parents 1de60c4d462d
children 675e571e6f4a
comparison
equal deleted inserted replaced
758:0ef96fcdc028 759:2645188489f6
16 assertType(bc, 'cell'); 16 assertType(bc, 'cell');
17 assert(S_sign == 1 || S_sign == -1, 'S_sign must be either 1 or -1'); 17 assert(S_sign == 1 || S_sign == -1, 'S_sign must be either 1 or -1');
18 18
19 19
20 % Setup storage arrays 20 % Setup storage arrays
21 closure = spzeros(size(diffOp)); 21 closure = spzeros(diffOp.size());
22 gridDataPenalties = {}; 22 gridDataPenalties = {};
23 gridDataFunctions = {}; 23 gridDataFunctions = {};
24 symbolicDataPenalties = {}; 24 symbolicDataPenalties = {};
25 symbolicDataFunctions = {}; 25 symbolicDataFunctions = {};
26 symbolicDataCoords = {}; 26 symbolicDataCoords = {};
56 error('sbplib:scheme:bcSetup:DataWrongNumberOfArguments', 'bc{%d}.data has the wrong number of input arguments. Must be either only time or time and space.', i); 56 error('sbplib:scheme:bcSetup:DataWrongNumberOfArguments', 'bc{%d}.data has the wrong number of input arguments. Must be either only time or time and space.', i);
57 end 57 end
58 end 58 end
59 59
60 % Setup penalty function 60 % Setup penalty function
61 O = spzeros(size(diffOp),1); 61 O = spzeros(diffOp.size(), 1);
62 function v = S_fun(t) 62 function v = S_fun(t)
63 v = O; 63 v = O;
64 for i = 1:length(gridDataFunctions) 64 for i = 1:length(gridDataFunctions)
65 v = v + gridDataPenalties{i}*gridDataFunctions{i}(t); 65 v = v + gridDataPenalties{i}*gridDataFunctions{i}(t);
66 end 66 end