Mercurial > repos > public > sbplib

comparison +rv/+time/getRvTimestepper.m @ 1192:b3c47a716d57 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix indentation
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 29 Jul 2019 16:43:39 +0200
parents 5aa3049a4212
children
comparison
equal deleted inserted replaced
1191:881afc40a3d2 1192:b3c47a716d57
29 fhData = @(v, t) data; 29 fhData = @(v, t) data;
30 end 30 end
31 end 31 end
32 32
33 function F = stabilizedRhs(D, data) 33 function F = stabilizedRhs(D, data)
34 fhData = dataFunctionHandle(data); 34 fhData = dataFunctionHandle(data);
35 F = @(v, t, viscosity) D(v, viscosity) + fhData(v,t); 35 F = @(v, t, viscosity) D(v, viscosity) + fhData(v,t);
36 end 36 end
37 37
38 function F = unstabilizedRhs(D, data) 38 function F = unstabilizedRhs(D, data)
39 fhData = dataFunctionHandle(data); 39 fhData = dataFunctionHandle(data);
40 F = @(v, t) D(v) + fhData(v,t); 40 F = @(v, t) D(v) + fhData(v,t);
66 66
67 function ts = rkRvInstageTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0) 67 function ts = rkRvInstageTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
68 F = stabilizedRhs(diffOpStruct.D_scheme, data); 68 F = stabilizedRhs(diffOpStruct.D_scheme, data);
69 ts = rv.time.RungekuttaRvInstage(F, opt.k, t0, v0, residualViscosity, diffOpStruct.D_t, opt.rkOrder); 69 ts = rv.time.RungekuttaRvInstage(F, opt.k, t0, v0, residualViscosity, diffOpStruct.D_t, opt.rkOrder);
70 end 70 end
71
72
73