Mercurial > repos > public > sbplib

comparison +time/SBPInTimeImplicitFormulation.m @ 462:7dbdf7390265 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in argument to f during time stepping
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 22 Jun 2017 13:06:03 +0200
parents 0b010f8de7cb
children 5df7f99206b2
comparison
equal deleted inserted replaced
461:0b010f8de7cb 462:7dbdf7390265
98 98
99 function obj = step(obj) 99 function obj = step(obj)
100 RHS = zeros(obj.blockSize*obj.N,1); 100 RHS = zeros(obj.blockSize*obj.N,1);
101 101
102 for i = 1:length(obj.blockSize) 102 for i = 1:length(obj.blockSize)
103 RHS((1 + (i-1)*obj.N):(i*obj.N)) = obj.f(obj.nodes(i)); 103 RHS((1 + (i-1)*obj.N):(i*obj.N)) = obj.f(obj.t + obj.nodes(i));
104 end 104 end
105 105
106 RHS = RHS + obj.K*obj.v; 106 RHS = RHS + obj.K*obj.v;
107 107
108 y = obj.L\RHS(obj.p); 108 y = obj.L\RHS(obj.p);