Mercurial > repos > public > sbplib
diff +time/SBPInTimeSecondOrderFormImplicit.m @ 1113:47e86b5270ad feature/timesteppers
Change name of property k to dt in time.Timestepper
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Wed, 10 Apr 2019 22:40:55 +0200 |
| parents | f5e14e5986b5 |
| children |
line wrap: on
line diff
--- a/+time/SBPInTimeSecondOrderFormImplicit.m Wed Apr 10 22:22:46 2019 +0200 +++ b/+time/SBPInTimeSecondOrderFormImplicit.m Wed Apr 10 22:40:55 2019 +0200 @@ -5,7 +5,7 @@ n t - k + dt firstOrderTimeStepper end @@ -14,7 +14,7 @@ % Solves A*u_tt + B*u_t + C*u = f(t) % A, B can either both be constants or both be function handles, % They can also be omitted by setting them equal to the empty matrix. - function obj = SBPInTimeSecondOrderFormImplicit(A, B, C, f, k, t0, v0, v0t, do_scaling, TYPE, order, blockSize) + function obj = SBPInTimeSecondOrderFormImplicit(A, B, C, f, dt, t0, v0, v0t, do_scaling, TYPE, order, blockSize) default_arg('f', []); default_arg('TYPE', []); default_arg('order', []); @@ -53,15 +53,15 @@ w0 = [v0; v0t]; - obj.k = k; + obj.dt = dt; obj.t = t0; obj.n = 0; if do_scaling scaling = [ones(m,1); sqrt(diag(C))]; - obj.firstOrderTimeStepper = time.SBPInTimeScaled(obj.AA, obj.BB, obj.ff, obj.k, obj.t, w0, scaling, TYPE, order, blockSize); + obj.firstOrderTimeStepper = time.SBPInTimeScaled(obj.AA, obj.BB, obj.ff, obj.dt, obj.t, w0, scaling, TYPE, order, blockSize); else - obj.firstOrderTimeStepper = time.SBPInTimeImplicitFormulation(obj.AA, obj.BB, obj.ff, obj.k, obj.t, w0, TYPE, order, blockSize); + obj.firstOrderTimeStepper = time.SBPInTimeImplicitFormulation(obj.AA, obj.BB, obj.ff, obj.dt, obj.t, w0, TYPE, order, blockSize); end end