Mercurial > repos > public > sbplib
diff +time/SBPInTimeImplicitFormulation.m @ 601:5df7f99206b2 feature/grids
SBPInTimeImplicitFormulation: Fix bug in calculation of forcing function
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Thu, 28 Sep 2017 11:14:34 +0200 |
parents | 7dbdf7390265 |
children | 47e86b5270ad |
line wrap: on
line diff
--- a/+time/SBPInTimeImplicitFormulation.m Tue Sep 26 15:38:57 2017 +0200 +++ b/+time/SBPInTimeImplicitFormulation.m Thu Sep 28 11:14:34 2017 +0200 @@ -46,7 +46,6 @@ obj.f = @(t)sparse(length(v0),1); end - obj.k = k; obj.blockSize = blockSize; obj.N = length(v0); @@ -99,7 +98,7 @@ function obj = step(obj) RHS = zeros(obj.blockSize*obj.N,1); - for i = 1:length(obj.blockSize) + for i = 1:obj.blockSize RHS((1 + (i-1)*obj.N):(i*obj.N)) = obj.f(obj.t + obj.nodes(i)); end