Mercurial > repos > public > sbplib
annotate +time/Cdiff.m @ 13:b18d3d201a71
Fixed initialization of step counter in timesteppers.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 22 Sep 2015 08:41:28 +0200 |
| parents | 5ae4f23d9130 |
| children | 484b48e95c83 |
| rev | line source |
|---|---|
| 0 | 1 classdef Cdiff < time.Timestepper |
| 2 properties | |
| 3 D | |
| 4 E | |
| 5 S | |
| 6 k | |
| 7 t | |
| 8 v | |
| 9 v_prev | |
| 10 n | |
| 11 end | |
| 12 | |
| 13 | |
| 14 methods | |
|
13
b18d3d201a71
Fixed initialization of step counter in timesteppers.
Jonatan Werpers <jonatan@werpers.com>
parents:
1
diff
changeset
|
15 function obj = Cdiff(D, E, S, k, t0, n0, v, v_prev) |
| 0 | 16 m = size(D,1); |
| 17 default_arg('E',sparse(m,m)); | |
| 18 default_arg('S',sparse(m,1)); | |
| 19 | |
| 20 if ~(issparse(D) && issparse(E) && issparse(S)) | |
| 21 warning('One of the matrices D, E, S is not sparse!') | |
| 22 print_issparse(D) | |
| 23 print_issparse(E) | |
| 24 print_issparse(S) | |
| 25 end | |
| 26 | |
| 27 obj.D = D; | |
| 28 obj.E = E; | |
| 29 obj.S = S; | |
| 30 obj.k = k; | |
|
1
5ae4f23d9130
Added CdiffNonlin timestepper. Probably fixed a bug with Cdiff. Added default arguments to Rk4SecondOrderNonlin.
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
31 obj.t = t0; |
|
13
b18d3d201a71
Fixed initialization of step counter in timesteppers.
Jonatan Werpers <jonatan@werpers.com>
parents:
1
diff
changeset
|
32 obj.n = n0; |
| 0 | 33 obj.v = v; |
| 34 obj.v_prev = v_prev; | |
| 35 end | |
| 36 | |
| 37 function [v,t] = getV(obj) | |
| 38 v = obj.v; | |
| 39 t = obj.t; | |
| 40 end | |
| 41 | |
| 42 function [vt,t] = getVt(obj) | |
| 43 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) | |
| 44 t = obj.t; | |
| 45 end | |
| 46 | |
| 47 function obj = step(obj) | |
| 48 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E, obj.S); | |
| 49 obj.t = obj.t + obj.k; | |
| 50 obj.n = obj.n + 1; | |
| 51 end | |
| 52 end | |
| 53 end |