Mercurial > repos > public > sbplib

annotate alignedTimestep.m @ 834:f1f0bf087e1c feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add support for artificial viscosity - Add option to use upwind operators for discretizing the scheme - Add option to use dissipation operators constructed from undivided differences when discretizing with narrow stencil operators
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 13 Sep 2018 14:21:47 +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