Mercurial > repos > public > sbplib

comparison +scheme/bcSetup.m @ 762:675e571e6f4a feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Revert unnecessary change in bcSetup.m
author Martin Almquist <malmquist@stanford.edu>
date Tue, 19 Jun 2018 14:25:43 -0700
parents 2645188489f6
children c70131daaa6e
comparison
equal deleted inserted replaced
761:8ed102db8e9c 762:675e571e6f4a
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(diffOp.size()); 21 closure = spzeros(size(diffOp));
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(diffOp.size(), 1); 61 O = spzeros(size(diffOp),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