Mercurial > repos > public > sbplib
annotate +time/Timestepper.m @ 5:5f6b0b6a012b
First try at interface implementation in WaveCurve2s
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Mon, 21 Sep 2015 09:42:15 +0200 |
| parents | 5205251db8c3 |
| children | b18d3d201a71 |
| rev | line source |
|---|---|
| 0 | 1 classdef Timestepper < handle |
| 2 properties (Abstract) | |
| 3 t | |
| 4 k | |
| 5 n | |
| 6 end | |
| 7 | |
| 8 methods (Abstract) | |
| 9 [v,t] = getV(obj) | |
| 10 obj = step(obj) | |
| 11 end | |
| 12 | |
| 13 | |
| 14 methods | |
| 15 function [v,t] = stepN(obj,n,progress_bar) | |
| 16 | |
| 17 if progress_bar && n > 1500 | |
| 18 n1000 = floor(n/1000); | |
| 19 | |
| 20 s = util.replace_string('',' %d %%',0); | |
| 21 for i=1:n | |
| 22 obj.step(); | |
| 23 if mod(i,n1000) == 0 | |
| 24 s = util.replace_string(s,' %.2f %%',i/n*100); | |
| 25 end | |
| 26 end | |
| 27 s = util.replace_string(s,''); | |
| 28 else | |
| 29 for i=1:n | |
| 30 obj.step(); | |
| 31 end | |
| 32 end | |
| 33 v = obj.getV; | |
| 34 t = obj.t; | |
| 35 end | |
| 36 | |
| 37 function [v,t] = evolve(obj, tend, progress_bar) | |
| 38 default_arg('progress_bar',false) | |
| 39 if progress_bar | |
| 40 obj.evolve_with_progress(tend); | |
| 41 else | |
| 42 obj.evolve_without_progress(tend); | |
| 43 end | |
| 44 v = obj.getV; | |
| 45 t = obj.t; | |
| 46 end | |
| 47 | |
| 48 function evolve_with_progress(obj, tend) | |
|
3
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
49 FRAME_RATE = 20; |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
50 |
| 0 | 51 dt = tend-obj.t; |
| 52 | |
|
3
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
53 steps_to_update = 1; |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
54 steps_since_update = 0; |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
55 last_update = tic(); |
| 0 | 56 s = util.replace_string('',' %d %%',0); |
| 57 while obj.t < tend - obj.k/100 | |
| 58 obj.step(); | |
| 59 | |
|
3
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
60 steps_since_update = steps_since_update + 1; |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
61 |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
62 if steps_since_update >= steps_to_update |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
63 s = util.replace_string(s,' %.2f %%',obj.t/tend*100); |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
64 |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
65 time_since_update = toc(last_update); |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
66 time_error = time_since_update - 1/FRAME_RATE; |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
67 time_per_step = time_since_update/steps_since_update; |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
68 |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
69 steps_to_update = max(steps_to_update - 0.9*time_error/time_per_step ,1); |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
70 |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
71 steps_since_update = 0; |
|
5205251db8c3
Added some smartness to evolve with progress.`
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
72 last_update = tic(); |
| 0 | 73 end |
| 74 end | |
| 75 | |
| 76 % if t < tend | |
| 77 % v = rk4.rungekutta_4(v, t, tend-t,F); | |
| 78 % end | |
| 79 | |
| 80 s = util.replace_string(s,''); | |
| 81 end | |
| 82 | |
| 83 function evolve_without_progress(obj, tend) | |
| 84 while obj.t < tend - obj.k/100 | |
| 85 obj.step(); | |
| 86 end | |
| 87 | |
| 88 % if t < tend | |
| 89 % v = rk4.rungekutta_4(v, t, tend-t,F); | |
| 90 % end | |
| 91 end | |
| 92 end | |
| 93 end |