Mercurial > repos > public > sbplib

changeset 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 f35ff0861d5a
children ecc605453733
files +rv/+time/getRvTimestepper.m
diffstat 1 files changed, 7 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- a/+rv/+time/getRvTimestepper.m	Fri Jul 05 17:47:13 2019 +0200
+++ b/+rv/+time/getRvTimestepper.m	Fri Jul 05 17:48:54 2019 +0200
@@ -1,5 +1,7 @@
 function ts = getRvTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
     switch opt.method
+        case 'rkRv'
+            ts = rkRvTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0);
         case 'rkRvBdf'
             ts = rkRvBdfTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0);
         case 'rkRvMs'
@@ -29,6 +31,11 @@
     end
 end
 
+function ts = rkRvTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
+    F = stabilizedRhs(diffOpStruct.D_scheme, data);
+    ts = rv.time.RungekuttaRv(F, opt.k, t0, v0, residualViscosity, diffOpStruct.D_t, opt.rkOrder);
+end
+
 function ts = rkRvBdfTimestepper(opt, diffOpStruct, residualViscosity, data, t0, v0)
     F = stabilizedRhs(diffOpStruct.D_scheme, data);
     ts = rv.time.RungekuttaRvBdf(F, opt.k, t0, v0, residualViscosity, opt.rkOrder, opt.bdfOrders);