Mercurial > repos > public > sbplib

annotate alignedTimestep.m @ 1016:4b42999874c0 feature/advectionRV

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add lower level for boot-strapping to RungeKuttaExteriorRV - Add a lower level to RungeKuttaExteriorRV for which bootstrapping starts, e.g start bootstrapping from time level 3 using a 3rd order BDF - Clean up ResidualViscosity
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 07 Dec 2018 13:11:53 +0100
parents c71095442c57
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
12
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 % Calcualtes the largest timestep smaller than k_max that gives an integer
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 % number of timesteps to time T.
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 % k_max -- largest allowable timestep
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 % T -- time to align with
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 %
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 % Returns:
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 % k -- calculated timestep
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8 % N -- number to of timestep to reach T
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 function [k, N] = alignedTimestep(k_max, T)
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 N = ceil(T/k_max);
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 k = T/N;
c71095442c57 Added function to calculate an aligned timestep.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12 end