Mercurial > repos > public > sbplib

annotate +rv/+time/getRvTimestepper.m @ 1183:27f3b173dcfa feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add RungekuttaRv to getRvTimesteppers
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 05 Jul 2019 17:48:54 +0200
parents cc2fd58989d1
children 5aa3049a4212
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1179
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
1 function ts = getRvTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
2 switch opt.method
1183
27f3b173dcfa Add RungekuttaRv to getRvTimesteppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1179
diff changeset
3 case 'rkRv'
27f3b173dcfa Add RungekuttaRv to getRvTimesteppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1179
diff changeset
4 ts = rkRvTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0);
1179
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
5 case 'rkRvBdf'
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
6 ts = rkRvBdfTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
7 case 'rkRvMs'
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
8 ts = rkRvMsTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
9 case 'rkRvMg'
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
10 ts = rkRvMgTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
11 case 'rkRvInstage'
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
12 ts = rkRvInstageTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
13 otherwise
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
14 error('Timestepping method ''%s'' not supported',method);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
15 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
16 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
17
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
18 function F = stabilizedRhs(D, data)
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
19 if isa(data, 'function_handle')
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
20 F = @(v, t, viscosity) D(v, viscosity) + data(v);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
21 else
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
22 F = @(v, t, viscosity) D(v, viscosity) + data;
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
23 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
24 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
25
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
26 function F = unstabilizedRhs(D, data)
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
27 if isa(data, 'function_handle')
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
28 F = @(v, t) D(v) + data(v);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
29 else
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
30 F = @(v, t) D(v) + data;
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
31 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
32 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
33
1183
27f3b173dcfa Add RungekuttaRv to getRvTimesteppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1179
diff changeset
34 function ts = rkRvTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
27f3b173dcfa Add RungekuttaRv to getRvTimesteppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1179
diff changeset
35 F = stabilizedRhs(diffOpStruct.D_scheme, data);
27f3b173dcfa Add RungekuttaRv to getRvTimesteppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1179
diff changeset
36 ts = rv.time.RungekuttaRv(F, opt.k, t0, v0, residualViscosity, diffOpStruct.D_t, opt.rkOrder);
27f3b173dcfa Add RungekuttaRv to getRvTimesteppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1179
diff changeset
37 end
27f3b173dcfa Add RungekuttaRv to getRvTimesteppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1179
diff changeset
38
1179
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
39 function ts = rkRvBdfTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
40 F = stabilizedRhs(diffOpStruct.D_scheme, data);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
41 ts = rv.time.RungekuttaRvBdf(F, opt.k, t0, v0, residualViscosity, opt.rkOrder, opt.bdfOrders);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
42 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
43
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
44 function ts = rkRvMsTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
45 F = stabilizedRhs(diffOpStruct.D_scheme, data);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
46 F_unstab = unstabilizedRhs(diffOpStruct.D_unstable, data);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
47 ts = rv.time.RungekuttaRvMultiStage(F, F_unstab, opt.k, t0, v0,...
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
48 residualViscosity, diffOpStruct.D_t, opt.rkOrder);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
49 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
50
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
51 function ts = rkRvMgTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
52 F = stabilizedRhs(diffOpStruct.D_scheme, data);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
53 F_coarse = unstabilizedRhs(diffOpStruct.D_coarse, data);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
54 ts = rv.time.RungekuttaRvMultiGrid(F, F_coarse, opt.k, t0, v0,...
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
55 residualViscosity, diffOpStruct.D_t, opt.rkOrder);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
56 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
57
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
58 function ts = rkRvInstageTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
59 F = stabilizedRhs(diffOpStruct.D_scheme, data);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
60 ts = rv.time.RungekuttaRvInstage(F, opt.k, t0, v0, residualViscosity, diffOpStruct.D_t, opt.rkOrder);
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
61 end
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
62
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
63
cc2fd58989d1 Add getTimestepper for the RV time steppers
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
64