Mercurial > repos > public > sbplib

comparison +time/SBPInTimeImplicitFormulation.m @ 461:0b010f8de7cb feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make f parameter optional.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 22 Jun 2017 13:05:32 +0200
parents e0caae9ef6ed
children 7dbdf7390265
comparison
equal deleted inserted replaced
460:e0caae9ef6ed 461:0b010f8de7cb
25 25
26 methods 26 methods
27 function obj = SBPInTimeImplicitFormulation(A, B, f, k, t0, v0, TYPE, order, blockSize) 27 function obj = SBPInTimeImplicitFormulation(A, B, f, k, t0, v0, TYPE, order, blockSize)
28 28
29 default_arg('TYPE','gauss'); 29 default_arg('TYPE','gauss');
30 default_arg('f',[]);
30 31
31 if(strcmp(TYPE,'gauss')) 32 if(strcmp(TYPE,'gauss'))
32 default_arg('order',4) 33 default_arg('order',4)
33 default_arg('blockSize',4) 34 default_arg('blockSize',4)
34 else 35 else
36 default_arg('blockSize',time.SBPInTimeImplicitFormulation.smallestBlockSize(order,TYPE)); 37 default_arg('blockSize',time.SBPInTimeImplicitFormulation.smallestBlockSize(order,TYPE));
37 end 38 end
38 39
39 obj.A = A; 40 obj.A = A;
40 obj.B = B; 41 obj.B = B;
41 obj.f = f; 42
43 if ~isempty(f)
44 obj.f = f;
45 else
46 obj.f = @(t)sparse(length(v0),1);
47 end
48
42 49
43 obj.k = k; 50 obj.k = k;
44 obj.blockSize = blockSize; 51 obj.blockSize = blockSize;
45 obj.N = length(v0); 52 obj.N = length(v0);
46 53