Mercurial > repos > public > sbplib

annotate +time/ExpInt.m @ 503:feebfca90080 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added the new abstract classes to hypsyst3D and Upwind1D to make them work with Jonatans changes
author Ylva Rydin <ylva.rydin@telia.com>
date Tue, 16 May 2017 09:18:10 +0200
parents 2d5dea05456d
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
502
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
1 classdef ExpInt < time.Timestepper
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
2 properties
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
3 F
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
4 k
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
5 t
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
6 v
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
7 m
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
8 n
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
9 end
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
10
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
11
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
12 methods
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
13 % Timesteps v_t = F(v,t), using RK4 fromt t = t0 with timestep k and initial conditions v = v0
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
14 function obj = ExpInt(F, k, t0, v0)
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
15 obj.F = F;
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
16 obj.k = k;
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
17 obj.t = t0;
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
18 obj.v = v0;
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
19 obj.m = length(v0);
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
20 obj.n = 0;
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
21 end
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
22
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
23 function [v,t] = getV(obj)
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
24 v = obj.v;
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
25 t = obj.t;
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
26 end
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
27
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
28 function obj = step(obj)
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
29 [t, res] = exp4(obj.F, [obj.t obj.t+obj.k], obj.v );
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
30 obj.t = t(end);
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
31 obj.v = res(end,:)';
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
32 obj.n = obj.n + 1;
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
33 end
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
34 end
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
35
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
36
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
37
2d5dea05456d An attempt to add the exponential integrators
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
38 end