Mercurial > repos > public > sbplib

annotate +noname/testCfl.m @ 819:fdf0ef9150f4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added tag v0.2 for changeset b723495cdb2f
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 07 Sep 2018 14:39:38 +0200
parents 5046ff7d13b8
children de8bcef865b0
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
1 % noname.testCfl(discr, timestepper_method, T, alpha0, tol,threshold, silentFlag)
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
2 % Example:
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
3 % noname.testCfl(Discr(100,4), 'rk4', 1, [0, 1])
39
7249f105e67b Added silent option to noname.testCfl.
Jonatan Werpers <jonatan@werpers.com>
parents: 2
diff changeset
4 function testCfl(discr, timestepper_method, T, alpha0, tol,threshold, silentFlag)
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 default_arg('tol',0.00005);
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 default_arg('threshold',1e2);
39
7249f105e67b Added silent option to noname.testCfl.
Jonatan Werpers <jonatan@werpers.com>
parents: 2
diff changeset
7 default_arg('silentFlag', false);
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8
118
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
9 % TODO:
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
10 % Set threshold from the initial conditions of the pde?
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
11 % Take a set number of steps instead of evolving to a certain time?
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
12 % Stop evolving when it has blown up?
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
13
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
14 testAlpha = getAlphaTester(discr, T, threshold, silentFlag, timestepper_method);
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
16 % Make sure that the upper bound is not working
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
17 ok = testAlpha(alpha0(2));
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
18 if ok % Upper bound too large!
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
19 error('The upper bound on alpha is stable!')
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
20 end
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
21
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
22 % Make sure that the lower bound is ok
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
23 if alpha0(1) ~= 0
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
24 ok = testAlpha(alpha0(1));
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
25 if ~ok
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
26 error('The lower bound on alpha is unstable!');
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
27 end
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
28 end
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
29
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
30 % Use bisection to find sharp estimate
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
31 while( (alpha0(2)-alpha0(1))/alpha0(1) > tol)
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
32 alpha = mean(alpha0);
118
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
33
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
34 if ~silentFlag
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
35 fprintf('[%.3e,%.3e]: ', alpha0(1), alpha0(2));
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
36 end
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
37
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
38 [ok, n_step, maxVal] = testAlpha(alpha);
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
39
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
40 if ok
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
41 alpha0(1) = alpha;
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
42 else
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
43 alpha0(2) = alpha;
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
44 end
118
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
45
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
46 if ~silentFlag
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
47 fprintf('a = %.3e, n_step=%d max = %.2e\n', alpha, n_step, maxVal);
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
48 end
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
49 end
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
50
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
51 if silentFlag
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
52 fprintf('Last calculated: a = %.3e, n_step=%d max = %.2e\n', alpha, n_step, maxVal);
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
53 end
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
54
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
55 fprintf('T = %-3d dof = %-4d order = %d: clf = %.4e\n',T, discr.size(), discr.order, alpha0(1));
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
56
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
57 end
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
58
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
59 function f = getAlphaTester(discr, T, threshold, silentFlag, timestepper_method)
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
60
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
61 % Returns true if cfl was ok
118
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
62 function [ok, n_step, maxVal] = testAlpha(alpha)
101
9933169d2651 Adapted testCfl to new opt style parameters in discrs.
Jonatan Werpers <jonatan@werpers.com>
parents: 39
diff changeset
63 ts = discr.getTimestepper(struct('method', timestepper_method, 'cfl', alpha));
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
64
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
65 warning('off','all')
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
66 ts.evolve(T,true);
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
67 warning('on','all')
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
68
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
69 [v,t] = ts.getV();
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
70 maxVal = max(v);
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
71
118
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
72 if isnan(maxVal) || maxVal == Inf || abs(maxVal) > threshold
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
73 ok = false;
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
74 else
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
75 ok = true;
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
76 end
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
77
118
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
78 n_step = ts.n;
2
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
79 end
bce9e28c1e26 Added a cfl paramater to getTimestepper in Discretization. Added function to get cfl value for a Discr.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
80
116
ee9c03dc7f42 noname.testCfl: Improved relyability and clarity.
Jonatan Werpers <jonatan@werpers.com>
parents: 101
diff changeset
81 f = @testAlpha;
118
5046ff7d13b8 Fixed bug and made more robust.
Jonatan Werpers <jonatan@werpers.com>
parents: 116
diff changeset
82 end