blob: 7e3ab4ae0fe691aca269c48dd78f813b96c6c24d (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
#include <limits>
#include <cctk.h>
#include <cctk_Arguments.h>
#include <cctk_Parameters.h>
#include <carpet.hh>
namespace Carpet {
using namespace std;
void CarpetRefineTimeStep (CCTK_ARGUMENTS)
{
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
// Find the smallest CFL factor for all refinement levels
CCTK_REAL min_cfl = numeric_limits<CCTK_REAL>::max();
for (int rl = 0; rl < reflevels; ++ rl) {
CCTK_REAL const level_cfl =
timereffacts.AT(rl) / maxval (spacereffacts.AT(rl));
min_cfl = min (min_cfl, level_cfl);
}
// Reduce the time step correspondingly
cctkGH->cctk_delta_time *= min_cfl;
CCTK_VInfo (CCTK_THORNSTRING,
"Timestep reduced by a factor of %g to %g",
double (1 / min_cfl),
double (cctkGH->cctk_delta_time/cctkGH->cctk_timefac));
}
} // namespace Carpet
|