Mercurial > repos > public > sbplib
annotate +time/Cdiff.m @ 289:354e40a8e1a5
Merged in feature/textTable (pull request #3)
Feature/texttable
| author | Martin Almquist <martin.almquist@it.uu.se> |
|---|---|
| date | Mon, 12 Sep 2016 15:18:25 +0200 |
| parents | 484b48e95c83 |
| children | 151ab2b5a686 b5e5b195da1e |
| 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) |
|
142
484b48e95c83
Removed ylim from setup1dPlot added some comments and fixed timestepper paramters.
Jonatan Werpers <jonatan@werpers.com>
parents:
13
diff
changeset
|
16 % Cdiff(D, E, S, k, t0, n0, v, v_prev) |
| 0 | 17 m = size(D,1); |
| 18 default_arg('E',sparse(m,m)); | |
| 19 default_arg('S',sparse(m,1)); | |
| 20 | |
| 21 if ~(issparse(D) && issparse(E) && issparse(S)) | |
| 22 warning('One of the matrices D, E, S is not sparse!') | |
| 23 print_issparse(D) | |
| 24 print_issparse(E) | |
| 25 print_issparse(S) | |
| 26 end | |
| 27 | |
| 28 obj.D = D; | |
| 29 obj.E = E; | |
| 30 obj.S = S; | |
| 31 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
|
32 obj.t = t0; |
|
13
b18d3d201a71
Fixed initialization of step counter in timesteppers.
Jonatan Werpers <jonatan@werpers.com>
parents:
1
diff
changeset
|
33 obj.n = n0; |
| 0 | 34 obj.v = v; |
| 35 obj.v_prev = v_prev; | |
| 36 end | |
| 37 | |
| 38 function [v,t] = getV(obj) | |
| 39 v = obj.v; | |
| 40 t = obj.t; | |
| 41 end | |
| 42 | |
| 43 function [vt,t] = getVt(obj) | |
| 44 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) | |
| 45 t = obj.t; | |
| 46 end | |
| 47 | |
| 48 function obj = step(obj) | |
| 49 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E, obj.S); | |
| 50 obj.t = obj.t + obj.k; | |
| 51 obj.n = obj.n + 1; | |
| 52 end | |
| 53 end | |
| 54 end |