Mercurial > repos > public > sbplib
annotate +time/Cdiff.m @ 1:5ae4f23d9130
Added CdiffNonlin timestepper. Probably fixed a bug with Cdiff. Added default arguments to Rk4SecondOrderNonlin.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Thu, 17 Sep 2015 19:19:58 +0200 |
| parents | 48b6fb693025 |
| children | b18d3d201a71 |
| 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 | |
| 15 function obj = Cdiff(D, E, S, k, t0, v, v_prev) | |
| 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; |
| 0 | 32 obj.v = v; |
| 33 obj.v_prev = v_prev; | |
| 34 end | |
| 35 | |
| 36 function [v,t] = getV(obj) | |
| 37 v = obj.v; | |
| 38 t = obj.t; | |
| 39 end | |
| 40 | |
| 41 function [vt,t] = getVt(obj) | |
| 42 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) | |
| 43 t = obj.t; | |
| 44 end | |
| 45 | |
| 46 function obj = step(obj) | |
| 47 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E, obj.S); | |
| 48 obj.t = obj.t + obj.k; | |
| 49 obj.n = obj.n + 1; | |
| 50 end | |
| 51 end | |
| 52 end |