Mercurial > repos > public > sbplib

annotate alignedTimestep.m @ 816:b5e5b195da1e feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add getState to timesteppers, returning the relevant state of the timestepper - Add getState which returns the 'state' of the specialized timestepper.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 10 Sep 2018 16:19:16 +0200
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