Mercurial > repos > public > sbplib

changeset 601:5df7f99206b2 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
SBPInTimeImplicitFormulation: Fix bug in calculation of forcing function
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 28 Sep 2017 11:14:34 +0200
parents 0da1c913a611
children 2dee6c9e827d
files +time/SBPInTimeImplicitFormulation.m
diffstat 1 files changed, 1 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
diff -r 0da1c913a611 -r 5df7f99206b2 +time/SBPInTimeImplicitFormulation.m
--- 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