Mercurial > repos > public > sbplib

annotate alignedTimestep.m @ 214:f0f4ca946068 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
blockmatrix: added a test to make sure a zero matrix is a blockmatrix. Also added input parameter checking to zero()
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 22 Jun 2016 14:42:01 +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